| 1 | title BLR model, density 1e12 cm-3, flux of H-ion phots 1e19 cm2 s-1, Z=20 |
|---|
| 2 | c |
|---|
| 3 | c commands controlling continuum ========= |
|---|
| 4 | table agn |
|---|
| 5 | phi(h) 19 |
|---|
| 6 | c |
|---|
| 7 | c commands for density & abundances ========= |
|---|
| 8 | hden 12 |
|---|
| 9 | abundances starburst 20 |
|---|
| 10 | c |
|---|
| 11 | c commands controlling geometry ========= |
|---|
| 12 | stop column density 23 |
|---|
| 13 | c |
|---|
| 14 | c other commands for details ========= |
|---|
| 15 | failures 5 |
|---|
| 16 | iterate convergence |
|---|
| 17 | c these are to try to speed things up |
|---|
| 18 | init "c84.ini" |
|---|
| 19 | no level2 lines |
|---|
| 20 | c |
|---|
| 21 | c commands controlling output ========= |
|---|
| 22 | normalize to "totl" 1216 100 |
|---|
| 23 | print line faint 1 // get rid of lots of faint lines |
|---|
| 24 | punch asserts "blr_n12_p19_Z20.asr" last |
|---|
| 25 | punch convergence reason "blr_n12_p19_Z20.cvr" |
|---|
| 26 | punch overview "blr_n12_p19_Z20.ovr" last |
|---|
| 27 | punch dr "blr_n12_p19_Z20.dr" last |
|---|
| 28 | c |
|---|
| 29 | c commands giving the asserts ========= |
|---|
| 30 | // strong blr lines |
|---|
| 31 | // >>chng 06 jul 07, from 7.581 to 7.560, Stancil N1 CT |
|---|
| 32 | // >>chng 06 aug 02, from 7.560 to 7.609, update O+2, N+2 CT |
|---|
| 33 | // >>chng 06 aug 09, from 7.609 to 7.633 Meudon 06, Bowen change |
|---|
| 34 | // >>chng 07 jul 27, from 7.633 to 7.608, repair common error in several collision routines |
|---|
| 35 | // >>chng 07 dec 03, from 7.608 to 7.587, RP changes iso sequences |
|---|
| 36 | // >>chng 08 mar 17, from 7.587 to 7.561, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes |
|---|
| 37 | assert line luminosity "totl" 1216 7.561 //total Hydrogen Lyman-alpha 1215.67 |
|---|
| 38 | // |
|---|
| 39 | //total of sum of both Si IV comp |
|---|
| 40 | // >>chng 04 sep 06, from 5.393 to 5.345, bug nej Auger |
|---|
| 41 | // >>chng 06 aug 09, from 5.345 to 5.403 Meudon 06, Bowen change |
|---|
| 42 | // >>chng 07 jan 20, from 5.403 to 5.367, e2 bug |
|---|
| 43 | // >>chng 07 jul 25, from 5.367 to 5.334, bugfix pvh FeII tau r1290 |
|---|
| 44 | // >>chng 08 mar 17, from 5.334 to 5.309, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes |
|---|
| 45 | // >>chng 08 nov 04, from 5.309 to 5.329, line RT bug fixes (r2429-35) |
|---|
| 46 | assert line luminosity "totl" 1397 5.329 |
|---|
| 47 | // |
|---|
| 48 | // >>chng 06 aug 02, from 3.076 to 3.112, update O+2, N+2 CT |
|---|
| 49 | // >>chng 06 aug 09, from 3.112 to 3.220 Meudon 06, Bowen change |
|---|
| 50 | // >>chng 07 jan 20, from 3.220 to 3.179, e2 bug |
|---|
| 51 | // >>chng 07 jul 25, from 3.179 to 3.152, bugfix pvh FeII tau r1290 |
|---|
| 52 | // >>chng 08 mar 17, from 3.152 to 3.115, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes |
|---|
| 53 | // >>chng 08 apr 10, from 3.115 to 3.139, change default h-like levels, other minor changes |
|---|
| 54 | assert line luminosity "totl" 1402 3.139 //total O IV] 1402 |
|---|
| 55 | // |
|---|
| 56 | // >>chng 06 aug 02, from 5.224 to 5.244, update O+2, N+2 CT |
|---|
| 57 | // >>chng 06 aug 09, from 5.244 to 5.308 Meudon 06, Bowen change |
|---|
| 58 | // >>chng 06 nov 29, from 5.308 to 5.243, Badnell DR turned on |
|---|
| 59 | // >>chng 07 jan 20, from 5.243 to 5.206, e2 bug |
|---|
| 60 | // >>chng 07 jul 25, from 5.206 to 5.181, bugfix pvh FeII tau r1290 |
|---|
| 61 | assert line luminosity "totl" 1549 5.181 //total of sum of both C IV comp |
|---|
| 62 | // |
|---|
| 63 | // >>chng 06 aug 02, from 6.527 to 6.498, update O+2, N+2 CT |
|---|
| 64 | // >>chng 08 mar 17, from 6.498 to 6.588, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes |
|---|
| 65 | // >>chng 08 apr 30, from 6.588 to 6.609 H, He topoff change |
|---|
| 66 | assert line luminosity "he 2" 1640 6.609 //total He II Balmer-alpha 1640 |
|---|
| 67 | // |
|---|
| 68 | // >>chng 06 aug 02, from 4.709 to 4.734, update O+2, N+2 CT |
|---|
| 69 | // >>chng 06 aug 09, from 4.734 to 4.795 Meudon 06, Bowen change |
|---|
| 70 | // >>chng 06 nov 29, from 4.795 to 4.761, Badnell DR turned on |
|---|
| 71 | // >>chng 07 jul 25, from 4.761 to 4.719, bugfix pvh FeII tau r1290 |
|---|
| 72 | // >>chng 08 mar 17, from 4.719 to 4.656, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes |
|---|
| 73 | assert line luminosity "o 3" 1666 4.656 //total O III] 1666 |
|---|
| 74 | // |
|---|
| 75 | // >>chng 06 aug 09, from 4.551 to 4.591 Meudon 06, Bowen change |
|---|
| 76 | // >>chng 07 jul 25, from 4.591 to 4.552, bugfix pvh FeII tau r1290 |
|---|
| 77 | // >>chng 08 mar 17, from 4.552 to 4.522, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes |
|---|
| 78 | assert line luminosity "totl" 1909 4.522 //total of sum of CIII] 1908.73+ |
|---|
| 79 | // |
|---|
| 80 | //total of sum of Si III] 1892+ |
|---|
| 81 | // >>chng 04 jun 01, from 6.941 to 6.964, unified calc of Lya line width |
|---|
| 82 | // >>chng 04 jul 02, from 6.964 to 6.941, O H ct |
|---|
| 83 | // >>chng 04 jul 03, from 6.941 to 6.962, NA update CO, H chem net |
|---|
| 84 | // >>chng 04 jul 07, from 6.962 to 6.927, oh ct integ |
|---|
| 85 | // >>chng 04 dec 08, from 6.927 to 6.952 co-ion feedback, check converg |
|---|
| 86 | // >>chng 05 mar 25, from 6.952 to 6.913, renorm H pops in level solver |
|---|
| 87 | // >>chng 06 aug 02, from 6.913 to 6.946, update O+2, N+2 CT |
|---|
| 88 | // >>chng 06 aug 25, from 6.946 to 6.967, better T fit to OIII cs |
|---|
| 89 | // >>chng 06 nov 29, from 6.967 to 6.932, Badnell DR turned on |
|---|
| 90 | // >>chng 07 jul 25, from 6.932 to 6.881, bugfix pvh FeII tau r1290 |
|---|
| 91 | // >>chng 07 nov 29, from 6.881 to 6.855, update DR rates |
|---|
| 92 | // >>chng 08 mar 17, from 6.855 to 6.823, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes |
|---|
| 93 | // >>chng 08 apr 20, incr error to more realistic 10% |
|---|
| 94 | assert line luminosity "totl" 1888 6.823 error 0.1 |
|---|
| 95 | // |
|---|
| 96 | // >>chng 06 aug 02, from 6.371 to 6.400, update O+2, N+2 CT |
|---|
| 97 | // >>chng 06 aug 09, from 6.400 to 6.426 Meudon 06, Bowen change |
|---|
| 98 | // >>chng 06 nov 29, from 6.426 to 6.246, Badnell DR turned on |
|---|
| 99 | // >>chng 07 jul 25, from 6.246 to 6.212, bugfix pvh FeII tau r1290 |
|---|
| 100 | // >>chng 07 dec 19, from 6.212 to 6.183, changes in r1698 |
|---|
| 101 | assert line luminosity "totl" 1860 6.183 //total sum of both Al III comp |
|---|
| 102 | // |
|---|
| 103 | //total C II] 2324.69+2328.12 |
|---|
| 104 | // >>chng 03 sep 07, from 5.686 to 5.709, ? ? |
|---|
| 105 | // >>chng 03 sep 10, from 5.709 to 5.771, changes in H2* population |
|---|
| 106 | // >>chng 03 sep 12, from 5.771 to 5.725, changes in H2* population, incr err |
|---|
| 107 | // >>chng 04 jun 12, from 5.649 to 5.816 chng cont shield fcn |
|---|
| 108 | // >>chng 04 jul 07, from 5.816 to 5.963, oh ct integ |
|---|
| 109 | // >>chng 04 jul 16, from 5.963 to 5.902, chng order eval O ion, incr err |
|---|
| 110 | // >>chng 04 jul 29, from 5.902 to 6.050, imp h chem conv, inc err |
|---|
| 111 | // >>chng 04 dec 08, from 6.050 to 5.805 co-ion feedback, check converg |
|---|
| 112 | // >>chng 05 mar 25, from 5.805 to 5.945, renorm H pops in level solver |
|---|
| 113 | // >>chng 05 aug 03, from 5.945 to 5.803, NA update CO, H chem net |
|---|
| 114 | // >>chng 06 aug 02, from 5.803 to 5.892, update O+2, N+2 CT |
|---|
| 115 | // >>chng 07 jan 20, from 5.892 to 5.795, e2 bug |
|---|
| 116 | // >>chng 07 dec 29, from 5.795 to 5.71, drift down |
|---|
| 117 | // >>chng 08 mar 22, from 5.71 to 5.61 bugfix r1879 |
|---|
| 118 | assert line luminosity "totl" 2326 5.61 err 0.2 |
|---|
| 119 | // |
|---|
| 120 | // >>chng 04 jul 07, from 7.807 to 7.833, oh ct integ |
|---|
| 121 | // >>chng 07 jul 25, from 7.833 to 7.803, bugfix pvh FeII tau r1290 |
|---|
| 122 | // >>chng 08 mar 17, from 7.803 to 7.770, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes |
|---|
| 123 | assert line luminosity "totl" 2798 7.770 //total of sum of both Mg II comp |
|---|
| 124 | // |
|---|
| 125 | // >>chng 04 sep 30, from 8 to 10, conv_base conv heat cool |
|---|
| 126 | // >>chng 08 dec 08, from 10 to 9.236, upper limit -> equals (r2515) |
|---|
| 127 | // >>chng 08 dec 13, from 9.236 to 18.39, merging newsolvers branch |
|---|
| 128 | // >>chng 09 jan 08, from 18.39 to 18.79, minor changes in iter/zn, iterations, and zone that crept in over last month |
|---|
| 129 | assert itrzn 18.79 error 0.02 |
|---|
| 130 | // |
|---|
| 131 | // >>chng 04 apr 20, needs one more iteration, H chem has changed |
|---|
| 132 | // >>chng 04 apr 25, needs one more, from 7 to 8 |
|---|
| 133 | // >>chng 04 apr 29, needs one from, to < 10, co col den changes |
|---|
| 134 | // >>chng 08 dec 08, from 10 to 4, upper limit -> equals (r2515) |
|---|
| 135 | assert niter 4 error 0.01 |
|---|
| 136 | // |
|---|
| 137 | // >>chng 03 dec 07, from 450 to 600, dr increase resol metals |
|---|
| 138 | // >>chng 05 dec 22, from 600 to 450 (393), rationalize zone asserts |
|---|
| 139 | // >>chng 08 dec 08, from 450 to 444, upper limit -> equals (r2515) |
|---|
| 140 | // >>chng 08 dec 13, from 444 to 448, merging newsolvers branch |
|---|
| 141 | assert nzone 448 error 0.01 |
|---|
| 142 | // |
|---|
| 143 | c blr_n12_p19_Z20.in |
|---|
| 144 | c class blr |
|---|
| 145 | c ======================================== |
|---|
| 146 | |
|---|
| 147 | This is one of the 5 models that sample the LOC plane. |
|---|