Changeset 1993

Show
Ignore:
Timestamp:
04/30/08 09:10:12 (2 weeks ago)
Author:
gary
Message:

continue rationalizing case a and case b

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved

  • trunk/tsuite/auto/h_caseb_hs87.in

    r1953 r1993  
    1414print populaions h-like 
    1515print departure h-like 
    16 assert Case B H-like hydrogen  
     16c  
     17assert Case B H-like hydrogen error 0.15 range 1,000 30,000A 
    1718assert line "TOTL" 4861    1.0000 error 0.02 
    1819assert line "Ca B" 4861    1.0000 error 0.02 

  • trunk/tsuite/auto/h_caseb_lot.in

    r1488 r1993  
    1 title log density case B, T=5000, log n=2 
     1title log density case B, T=500 log n=2 
    22c 
    33c commands controlling continuum ========= 
    …… 
    77c commands for density & abundances ========= 
    88hden = 2 
    9 init file="hheonly.ini" 
     9init file="honly.ini" 
    1010c 
    1111c commands controlling geometry  ========= 
    …… 
    3232punch asserts "h_caseb_lot.asr" last 
    3333c 
    34 c commands giving the asserts    ========= 
    35 // >>chng 99Jun15 next two changed, by .2 and .08 dex respectively, 
    36 // when ots oscil fix put in place 92.20 
    37 // >>chng 99jun18 increased error on H1 and He1, since case b does not 
    38 // really define this, and they are so affected by lyman line transport 
    39 // >>chng 99nov11, changed H and He to HS case B expectations 
    40 // >>chng 01 sep 24, from -3.44 to -2.994, chng def case b tau,  
    41 // >>chng 06 aug 09, from -2.994 to -3.13 Meudon 06, Bowen change 
    42 // >>chng 07 oct 22, from -3.130 to -2.894, resolve l-levels of h-like sequence 
    43 assert ionization hydrogen 1 -2.894 error 0.1 
    44 // 
    45 // >>chng 01 sep 24, from -5.686 to -5.270 chng def case b tau,  
    46 // >>chng 02 jan 19, from -5.270 to -5.306, had not blocked grnd rec in 
    47 //      case b for he 
    48 // >>chng 02 feb 03, from -5.306 to -5.145, turned off excit state photo 
    49 // >>chng 02 dec 11, from -5.145 to -5.240 new hei 
    50 // >>chng 06 aug 09, from -5.24 to -5.39 Meudon 06, Bowen change 
    51 // >>chng 07 oct 22, from -5.39 to -5.12, resolve l-levels of h-like sequence 
    52 assert ionization helium 1 -5.12 error 0.1 
    53 // 
    54 assert ionizaiton helium 2 -1.749 
    55 assert line "H  1" 4861 1.00 error 0.07 
    56 // 
    57 // >>chng 03 jun 08, from xx to 3.69, incr error 
    58 // >>chng unknown  , from 3.69 to 3.97, ?? */ 
    59 // >>chng 07 oct 22, from 3.97 to 3.62, resolve l-levels of h-like sequence 
    60 assert line "H  1" 6563 3.62 error 0.08 
    61 // 
    62 // >>chng 03 jun 02, from 0.752 to 0.665, turn on large H atom 
    63 // >>chng 06 may 08, from 0.665 to 0.700, RP update high-n As H-like 
    64 // >>chng 07 oct 22, from 0.700 to 0.625, resolve l-levels of h-like sequence 
    65 assert line "H  1"  1.875m 0.625 
    66 // 
    67 // >>chng 03 jun 02, from 0.249 to 0.228, turn on large H atom 
    68 // >>chng 06 may 08, from 0.228 to 0.247, RP update high-n As H-like 
    69 // >>chng 07 oct 22, from 0.247 to 0.234, resolve l-levels of h-like sequence 
    70 assert line "H  1"  1.282m 0.234 
    71 // 
    72 // >>chng 03 jun 02, from 0.233 to 0.197, turn on large H atom 
    73 // >>chng 06 may 08, from 0.197 to 0.214, RP update high-n As H-like 
    74 // >>chng 07 oct 22, from 0.214 to 0.181, resolve l-levels of h-like sequence 
    75 assert line "H  1"  4.051m 0.181 
    76 // 
    77 // >>chng 03 jun 02, from 0.095 to 0.0836, turn on large H atom 
    78 // >>chng 06 may 08, from 0.0836 to 0.0933, RP update high-n As H-like 
    79 // >>chng 07 oct 22, from 0.0933 to 0.0845, resolve l-levels of h-like sequence 
    80 assert line "H  1"  2.625m 0.0845 
    81 // 
    82 // >>chng 03 jun 02, from 8.745 to 8.117, turn on large H atom 
    83 // >>chng 06 may 08, from 8.12 to 8.48 RP update high-n As H-like 
    84 // >>chng 07 oct 22, from 8.48 to 7.48, resolve l-levels of h-like sequence 
    85 assert line "He 2" 1640 7.48 
    86 // 
    87 // >>chng 03 jun 02, from 1.933 to 1.760, turn on large H atom 
    88 // >>chng 06 may 08, from 1.76 to 1.93 RP update high-n As H-like 
    89 // >>chng 07 oct 22, from 1.93 to 1.39, resolve l-levels of h-like sequence 
    90 assert line "He 2" 4686 1.39 
    91 // 
    92 assert line "Bac " 3646 144.4 
    93 assert line "Ba C" 0 1.93 
    94 assert line "Ca B" 4861 1.00 error 0.07 
    95 assert line "Ca B" 6563 3.97 error 0.08 
    96 assert line "Ca B"  1.875m 0.752 
    97 assert line "Ca B"  1.282m 0.247 
    98 assert line "Ca B"  4.051m 0.233 
    99 assert line "Ca B"  2.625m 0.0933 
    100 assert line "Ca B" 1640 8.48 
    101 assert line "Ca B" 4686 1.93 
     34assert Case B H-like element hydrogen  0.05 range 1,000 40,000A 
     35 
    10236c h_caseb_lot.in 
    10337c class limit  

  • trunk/tsuite/auto/h_caseb_n8.in

    r1935 r1993  
    77c commands for density & abundances ========= 
    88hden = 8 
    9 init file="hheonly.ini" 
     9init file="honly.ini" 
    1010c 
    1111c commands controlling geometry  ========= 
    …… 
    1414c 
    1515c other commands for details     ========= 
    16 case b hummer no photoionization 
     16case b hummer no photoionization no Pdest  
     17no scattering escape 
    1718no induced processes (OK) 
    1819constant temperature = 20000 
    19 atom h-like levels resolved 15  
    20 atom h-like levels collapsed 35  
    2120iterate  
    2221c 
    …… 
    2726punch asserts "h_caseb_n8.asr" 
    2827c 
    29 c commands giving the asserts    ========= 
    30 // 99 jun 18, increase error on H1 and He1 since small chagnes 
    31 // in lyman line dest probs affect HeII Lya and ionization  
    32 // >>chng 01 sep 24, put back in after lowering lya tay to 1e5 
    33 // >>chng 02 nov 20, from -4.379 to -4.442, new he ots diff 
    34 // >>chng 03 mar 13, from -4.442 to -4.419, H- rates better 
    35 // >>chng 03 dec 08, from -4.419 to -4.444, chng zones, ots 
    36 // >>chng 06 aug 09, from -4.444 to -4.960 Meudon 06, Bowen change 
    37 // >>chng 07 mar 13, from -4.960 to -4.935, change in lines_service to fix cache errors 
    38 // >>chng 07 jul 27, from -4.935 to -4.896, repair common error in several collision routines 
    39 // >>chng 07 oct 22, from -4.896 to -4.961, resolve l-levels of h-like sequence 
    40 // >>chng 07 nov 17, from -4.961 to -4.820, changes in r1512 
    41 // >>chng 07 nov 23, from -4.820 to -4.784, this is result w/32 bit vs03 
    42 // >>chng 07 nov 29, from -4.784 to -4.731, due to changes in r1579 
    43 // >>chng 08 feb 05, from -4.731 to -4.820, the nov 17 value, bi-stable? 
    44 // >>chng 08 apr 07, from -4.820 to -4.627, a number of mostly minor iso changes 
    45 // >>chng 08 apr 10, from -4.627 to -4.459, change default h-like levels, other minor changes 
    46 assert ionization hydrogen 1 -4.459 
    47 // 
    48 assert ionizaiton helium 2 -2.859 
    49 // 
    50 // 99jun18, changed error to default since no longer 
    51 // agreed to better than 3%.  change was when lyman 
    52 // line dest probs cal with damper - lyman lines are escape route? 
    53 // >>chng 01 aug 23, from 11.98 to 10.04, better continuum resolution 
    54 assert line "BaC " 3646 10.04  error 0.02 
    55 assert line "H  1" 6563  2.6632 
    56 assert line "H  1" 4861  1.0000 
    57 assert line "H  1" 4340  0.4893 
    58 assert line "H  1"  1.875m  0.2594 
    59 assert line "H  1"  1.282m  0.1380 
    60 assert line "H  1"  1.094m  0.0813 
    61 assert line "H  1"  4.051m  0.0547 
    62 assert line "H  1"  2.625m  0.0340 
    63 // >>chng 07 jul 27, from 0.0221 to 0.0233, repair common error in several collision routines 
    64 assert line "H  1"  2.166m  0.0233 
    65 assert line "H  1"  7.458m  0.0163 
    66 // >>chng 07 jul 27, from 0.0112 to 0.0118, repair common error in several collision routines 
    67 // >>chng 07 oct 26, from 0.0118 to 0.0113, very close to botch after recent changes 
    68 assert line "H  1"  4.653m  0.0113 
    69 // >>chng 07 jul 27, from 0.0078 to 0.0083, repair common error in several collision routines 
    70 // >>chng 07 oct 26, from 0.0083 to 0.0079, very close to botch after recent changes 
    71 // >>chng 08 mar 17, from 0.0079 to 0.0083, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes 
    72 assert line "H  1"  3.740m  0.0083 
    73 assert line "He 2" 1640  7.8736 
    74 assert line "He 2" 1215  2.7056 
    75 assert line "He 2" 1085  1.2602 
    76 assert line "He 2" 4686  0.9804 
    77 assert line "He 2" 3203  0.4699 
    78 assert line "He 2" 2733  0.2595 
    79 // >>chng 07 oct 22, from 0.2351 to 0.2233, resolve l-levels of h-like sequence 
    80 assert line "He 2"  1.012m  0.2233 
    81 assert line "He 2" 6560  0.1340 
    82 assert line "He 2" 5412  0.0821 error 0.06 
    83 c h_caseb_n8.in 
     28assert case b h-like element hydrogen 0.10 range 1,000A to 40,000A 
    8429c class limit  
    8530c ======================================== 
    …… 
    9540-       Relative line intensities 
    9641 
     42high density causes disagreemeent with HS - collisions 

  • trunk/tsuite/auto/limit_caseb_he_den.in

    r1992 r1993  
    2424c commands controlling output    ========= 
    2525atom H-like print levels  
    26 punch asserts "grid_caseb_he.asr" 
     26punch asserts "limit_caseb_he_den.asr" 
    2727c  
    2828c commands giving the asserts    ========= 
    2929assert Case B H-like Helium   error 0.09 range 1,000A to 10,000A 
    3030c  
    31 c grid_caseb_hhe.in  
     31c limit_caseb_he_den.in  
    3232c class limit  
    3333c ======================================== 
    3434c  
    3535 
    36 this tests the predicted H I, He I, and He II spectra in the 
    37 Case B limit. 
     36this tests the predicted He II spectra in the Case B limit. 
    3837 
    3938this effectively turns off hydrogen to avoid the problem with every 
    …… 
    4140the number of levels for H I. 
    4241 
    43 this asserts the values are within 15%. 
    44 Actually they are all nearly within a few percent except at the lowest 
    45 temperature.  The error greater than 10% occurs at the lowest temperature 
    46 of 5,000K. 
     42this asserts the values are within 9% for the standard T = 1e4K and 
     43a range of densities.  Actually they are all nearly within a few 
     44percent except at the lowest temperature of 5,000K. 

  • trunk/tsuite/auto/limit_caseb_he_den_temp.in

    r1992 r1993  
    1 title Test model H and He atoms in Case B limit 
     1title Test model He atoms in Case B limit 
    22c grid over nebular range of density and temperature to check Case B H & He 
    33c 
    …… 
    2525c commands controlling output    ========= 
    2626atom H-like print levels  
    27 punch asserts "grid_caseb_he.asr" 
     27punch asserts "limit_caseb_he_den_temp.asr" 
    2828c  
    2929c commands giving the asserts    ========= 
    3030assert Case B H-like Helium   error 0.13 range 1,000A to 10,000A 
    3131c  
    32 c grid_caseb_hhe.in  
     32c limit_caseb_he_den_temp.in  
    3333c class limit  
    3434c ======================================== 
    3535c  
    3636 
    37 this tests the predicted H I, He I, and He II spectra in the 
    38 Case B limit. 
     37this tests the predicted He II spectra in the Case B limit. 
    3938 
    4039this effectively turns off hydrogen to avoid the problem with every 
    …… 
    4241the number of levels for H I. 
    4342 
    44 this asserts the values are within 15%. 
     43this asserts the values are within 15% for a range of density and temperature. 
    4544Actually they are all nearly within a few percent except at the lowest 
    4645temperature.  The error greater than 10% occurs at the lowest temperature 

  • trunk/tsuite/auto/limit_caseb_hhe_den.in

    r1992 r1993  
    2121c commands controlling output    ========= 
    2222atom H-like print levels 
    23 punch asserts "grid_caseb_hhe.asr" 
     23punch asserts "limit_caseb_hhe_den.asr" 
    2424c  
    2525c commands giving the asserts    ========= 
    …… 
    2727assert Case B He-like Helium   error 0.05 range 1,000A to 40,000A 
    2828c  
    29 c grid_caseb_hhe.in  
     29c limit_caseb_hhe_den.in  
    3030c class limit  
    3131c ======================================== 
    3232c  
    3333 
    34 this tests the predicted H I and He I spectra in the 
    35 Case B limit. 
     34this tests the predicted H I and He I spectra in the Case B limit. 
    3635 
     36the grid is over density at the standard temperature of 1e4 K. 

  • trunk/tsuite/auto/limit_caseb_hhe_den_temp.in

    r1992 r1993  
    2222c commands controlling output    ========= 
    2323atom H-like print levels 
    24 punch asserts "grid_caseb_hhe.asr" 
     24punch asserts "limit_caseb_hhe_den_temp.asr" 
    2525c  
    2626c commands giving the asserts    ========= 
    …… 
    2828assert Case B He-like Helium   error 0.10 range 1,000A to 40,000A 
    2929c  
    30 c grid_caseb_hhe.in  
     30c limit_caseb_hhe_den_temp.in  
    3131c class limit  
    3232c ======================================== 
    3333c  
    3434 
    35 this tests the predicted H I and He I spectra in the 
    36 Case B limit. 
     35this tests the predicted H I and He I spectra in the Case B limit. 
     36The grid is over both density and temperature. 
    3737