| 1 | title BLR model, density 1e09 cm-3, flux of H-ion phots 1e18 cm2 s-1, Z=20 |
|---|
| 2 | c |
|---|
| 3 | c commands controlling continuum ========= |
|---|
| 4 | table agn |
|---|
| 5 | phi(h) 18 |
|---|
| 6 | c |
|---|
| 7 | c commands for density & abundances ========= |
|---|
| 8 | hden 9 |
|---|
| 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 | iterate to convergence |
|---|
| 16 | c these are to try to speed things up |
|---|
| 17 | init "c84.ini" |
|---|
| 18 | no level2 lines |
|---|
| 19 | c |
|---|
| 20 | c commands controlling output ========= |
|---|
| 21 | normalize to "totl" 1216 100 |
|---|
| 22 | print h-like hydrogen departure coefficients |
|---|
| 23 | print line faint 1 // get rid of lots of faint lines |
|---|
| 24 | punch asserts "blr_n09_p18_Z20.asr" last |
|---|
| 25 | punch overview "blr_n09_p18_Z20.ovr" |
|---|
| 26 | punch dr "blr_n09_p18_Z20.dr" |
|---|
| 27 | punch convergence reason "blr_n09_p18_Z20.cvr" |
|---|
| 28 | punch convergence error "blr_n09_p18_Z20.cve" |
|---|
| 29 | c |
|---|
| 30 | c commands giving the asserts ========= |
|---|
| 31 | // strong blr lines |
|---|
| 32 | // >>chng 04 mar 04, from 3.696 to 3.806, chng convergence criteria |
|---|
| 33 | // >>chng 04 sep 06, from 3.806 to 3.792, bug nej Auger |
|---|
| 34 | // >>chng 04 oct 02, from 3.792 to 3.778, drift down |
|---|
| 35 | assert line luminosity "totl" 1035 3.778 //total of sum of both O VI comp |
|---|
| 36 | // |
|---|
| 37 | //total Hydrogen Lyman-alpha 1215.67 |
|---|
| 38 | // >>chng 04 jun 01, from 6.796 to 6.900, unified calc of Lya line width |
|---|
| 39 | // >>chng 06 nov 29, from 6.900 to 6.944, Badnell DR turned on |
|---|
| 40 | // >>chng 07 jul 25, from 6.944 to 6.987, bugfix pvh FeII tau r1290 |
|---|
| 41 | assert line luminosity "totl" 1216 6.987 |
|---|
| 42 | // |
|---|
| 43 | //total of sum of both N V comp |
|---|
| 44 | // >>chng 06 nov 29, from 5.621 to 5.566, Badnell DR turned on |
|---|
| 45 | assert line luminosity "totl" 1240 5.566 |
|---|
| 46 | // |
|---|
| 47 | //total of sum of both Si IV comp |
|---|
| 48 | // >>chng 04 jul 28, from 5.731 to 5.753, drift |
|---|
| 49 | // >>chng 06 aug 09, from 5.753 to 5.777 Meudon 06, Bowen change |
|---|
| 50 | // >>chng 06 nov 29, from 5.777 to 5.817, Badnell DR turned on |
|---|
| 51 | assert line luminosity "totl" 1397 5.817 |
|---|
| 52 | // |
|---|
| 53 | // >>chng 06 nov 29, from 5.163 to 5.214, Badnell DR turned on |
|---|
| 54 | assert line luminosity "totl" 1402 5.214 //total O IV] 1402 |
|---|
| 55 | // |
|---|
| 56 | // >>chng 07 jan 20, from 5.980 to 5.955, e2 bug |
|---|
| 57 | assert line luminosity "totl" 1549 5.955 //total of sum of both C IV comp |
|---|
| 58 | // |
|---|
| 59 | // >>chng 06 nov 29, from 5.394 to 5.415, Badnell DR turned on |
|---|
| 60 | // >>chng 07 oct 22, from 5.415 to 5.388, resolve l-levels of h-like sequence |
|---|
| 61 | // >>chng 08 feb 13, from 5.388 to 5.449, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes |
|---|
| 62 | assert line luminosity "he 2" 1640 5.449 //total He II Balmer-alpha 1640 |
|---|
| 63 | // |
|---|
| 64 | // >>chng 06 jul 25, from 5.464 to 5.382 HN update OIII collision data |
|---|
| 65 | // >>chng 06 nov 29, from 5.382 to 5.418, Badnell DR turned on |
|---|
| 66 | assert line luminosity "o 3" 1666 5.418 //total O III] 1666 |
|---|
| 67 | // |
|---|
| 68 | // >>chng 04 sep 23, from 6.062 to 6.041, slow ddrift down |
|---|
| 69 | // >>chng 06 nov 29, from 6.041 to 6.113, Badnell DR turned on |
|---|
| 70 | assert line luminosity "totl" 1909 6.113 //total of sum of CIII] 1908.73+ |
|---|
| 71 | // |
|---|
| 72 | // >>chng 04 jun 01, from 6.378 to 6.343, unified calc of Lya line width |
|---|
| 73 | // >>chng 07 jul 25, from 6.343 to 6.321, bugfix pvh FeII tau r1290 |
|---|
| 74 | // >>chng 07 oct 22, from 6.321 to 6.278, resolve l-levels of h-like sequence |
|---|
| 75 | // >>chng 07 oct 26, from 6.278 to 6.323, due to changes in r1472 |
|---|
| 76 | // >>chng 07 nov 17, from 6.323 to 6.282, changes in r1512 |
|---|
| 77 | // >>chng 07 nov 23, from 6.282 to 6.316, result with 32 bit vs03 |
|---|
| 78 | assert line luminosity "totl" 1888 6.316 //total of sum of Si III] 1892+ |
|---|
| 79 | // |
|---|
| 80 | // >>chng 04 jun 01, from 5.432 to 5.372, unified calc of Lya line width |
|---|
| 81 | // >>chng 07 jul 25, from 5.372 to 5.352, bugfix pvh FeII tau r1290 |
|---|
| 82 | assert line luminosity "totl" 1860 5.352 //total sum of both Al III comp |
|---|
| 83 | // |
|---|
| 84 | //total C II] 2324.69 + 2328.12 |
|---|
| 85 | // >>chng 04 may 19, from 5.715 to 5.686 RP update helike sequence |
|---|
| 86 | // >>chng 04 jun 01, from 5.686 to 5.634, unified calc of Lya line width |
|---|
| 87 | // >>chng 04 jul 02, from 5.634 to 5.914, O H ct |
|---|
| 88 | // >>chng 04 jun 03, from 5.914 to 5.622, NA update CO, H chem net, incr err |
|---|
| 89 | // >>chng 04 jul 07, from 5.62 to 6.12, oh ct integ |
|---|
| 90 | // >>chng 04 jul 16, from 6.12 to 5.95, chng order of eval of heav elem, inc er |
|---|
| 91 | // >>chng 04 jul 29, from 5.96 to 6.22, impr conv h chem, incr error |
|---|
| 92 | // >>chng 04 sep 25, from 6.22 to 6.21, incr temp toler 0.02 to 0.005 |
|---|
| 93 | // >>chng 05 jul 17, had been missing, reset |
|---|
| 94 | // >>chng 07 nov 17, from 5.926 to 6.006, changes in r1512 |
|---|
| 95 | // >>chng 07 nov 23, from 6.006 to 5.917, result with 32 bit vs03 |
|---|
| 96 | // >>chng 08 apr 30, from 5.92 to 5.77, H, He topoff change |
|---|
| 97 | // >>chng 08 nov 13, from 5.77 to 5.683, fix Boltzmann factor in CT rates (r2455) |
|---|
| 98 | assert line luminosity "totl" 2326 5.683 error 0.2 |
|---|
| 99 | // |
|---|
| 100 | //total of sum of both Mg II comp |
|---|
| 101 | // >>chng 04 jun 01, from 6.629 to 6.578, unified calc of Lya line width |
|---|
| 102 | // >>chng 04 jun 12, from 6.578 to 6.616 chng cont shield fcn |
|---|
| 103 | // >>chng 04 jul 02, from 6.615 to 6.745, O H ct |
|---|
| 104 | // >>chng 04 jun 03, from 6.745 to 6.614, NA update CO, H chem net, incr err |
|---|
| 105 | // >>chng 04 jul 07, from 6.61 to 6.88, oh ct integ |
|---|
| 106 | // >>chng 04 jul 16, from 6.88 to 6.76, chng order of eval of heav elem, inc er |
|---|
| 107 | // >>chng 04 jul 29, from 6.76 to 6.94, impr conv h chem, incr error |
|---|
| 108 | // >>chng 04 sep 25, from 6.94 to 6.79, incr temp toler 0.02 to 0.005 |
|---|
| 109 | // >>chng 04 dec 08, from 6.79 to 6.61 co-ion feedback, check converg |
|---|
| 110 | // >>chng 05 mar 25, from 6.61 to 6.74 renorm H pops in level solver |
|---|
| 111 | // >>chng 08 apr 30, from 6.74 to 6.63, H, He topoff change |
|---|
| 112 | assert line luminosity "totl" 2798 6.63 error 0.2 |
|---|
| 113 | // |
|---|
| 114 | // >>chng 04 sep 25, from 6.6 to 8.5 incr temp toler 0.02 to 0.005 |
|---|
| 115 | // >>chng 08 dec 08, from 8.5 to 6.686, upper limit -> equals (r2515) |
|---|
| 116 | // >>chng 08 dec 13, from 6.686 to 18.5, merging newsolvers branch |
|---|
| 117 | assert itrzn 18.5 error 0.02 |
|---|
| 118 | // >>chng 07 dec 06, from 5 to 7, changes in h-like collisions |
|---|
| 119 | // >>chng 07 dec 19, from 7 to 9, changes in r1698 |
|---|
| 120 | // >>chng 08 dec 08, from 9 to 3, upper limit -> equals (r2515) |
|---|
| 121 | assert niter 3 error 0.01 |
|---|
| 122 | // |
|---|
| 123 | // >>chng 03 nov 10, from 390 to 450, improve he0 zoning |
|---|
| 124 | // >>chng 03 dec 07, from 450 to 550, dr improve metals resolution |
|---|
| 125 | // >>chng 05 dec 22, from 550 to 420 (371), rationalize zone asserts |
|---|
| 126 | // >>chng 08 apr 30, from 420 to 460 (444) H, He topoff change |
|---|
| 127 | // >>chng 08 dec 08, from 460 to 455, upper limit -> equals (r2515) |
|---|
| 128 | // >>chng 08 dec 13, from 455 to 400, merging newsolvers branch |
|---|
| 129 | // >>chng 09 jan 08, from 400 to 407, minor changes in iter/zn, iterations, and zone that crept in over last month |
|---|
| 130 | assert nzone 407 error 0.01 |
|---|
| 131 | // |
|---|
| 132 | c blr_n09_p18_Z20.in |
|---|
| 133 | c class blr |
|---|
| 134 | c ======================================== |
|---|
| 135 | |
|---|
| 136 | This is one of the 5 models that sample the LOC plane. |
|---|