| 1 | title BLR model, density 1e11 cm-3, flux of H-ion phots 1e20 cm2 s-1, Z=20 |
|---|
| 2 | c |
|---|
| 3 | c commands controlling continuum ========= |
|---|
| 4 | table agn |
|---|
| 5 | phi(h) 20 |
|---|
| 6 | c |
|---|
| 7 | c commands for density & abundances ========= |
|---|
| 8 | hden 11 |
|---|
| 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 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 line faint 1 // get rid of lots of faint lines |
|---|
| 23 | punch cooling "blr_n11_p20_Z20.col" |
|---|
| 24 | punch heating "blr_n11_p20_Z20.het" |
|---|
| 25 | punch element calcium "blr_n11_p20_Z20.cal" |
|---|
| 26 | punch asserts "blr_n11_p20_Z20.asr" last |
|---|
| 27 | punch overview "blr_n11_p20_Z20.ovr" |
|---|
| 28 | punch dr "blr_n11_p20_Z20.dr" |
|---|
| 29 | punch temperature history "blr_n11_p20_Z20.tem" |
|---|
| 30 | c |
|---|
| 31 | c commands giving the asserts ========= |
|---|
| 32 | // strong blr lines |
|---|
| 33 | // >>chng 06 jul 18, from 6.004 to 5.978, RP major upgrade VS coll strengths |
|---|
| 34 | assert line luminosity "totl" 1035 5.978 //total of sum of both O VI comp |
|---|
| 35 | // |
|---|
| 36 | // >>chng 04 jun 01, from 8.534 to 8.489, unified calc of Lya line width |
|---|
| 37 | // >>chng 06 jul 18, from 8.489 to 8.529, RP major upgrade VS coll strengths |
|---|
| 38 | // >>chng 06 nov 29, from 8.529 to 8.561, Badnell DR turned on |
|---|
| 39 | // >>chng 07 jul 25, from 8.561 to 8.618, bugfix pvh FeII tau r1290 |
|---|
| 40 | // >>chng 08 mar 17, from 8.618 to 8.575, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes |
|---|
| 41 | assert line luminosity "totl" 1216 8.575 //total Hydrogen Lyman-alpha 1215.67 |
|---|
| 42 | // |
|---|
| 43 | // >>chng 06 nov 29, from 7.853 to 7.831, Badnell DR turned on |
|---|
| 44 | assert line luminosity "totl" 1240 7.831 //total of sum of both N V comp |
|---|
| 45 | // |
|---|
| 46 | // >>chng 06 nov 29, from 8.126 to 8.147, Badnell DR turned on |
|---|
| 47 | assert line luminosity "totl" 1397 8.147 //total of sum of both Si IV comp |
|---|
| 48 | // |
|---|
| 49 | //total O IV] 1402 |
|---|
| 50 | // >>chng 06 nov 24, from 7.177 to 7.199, turn on Badnell RR |
|---|
| 51 | // >>chng 06 nov 29, from 7.199 to 7.256, Badnell DR turned on |
|---|
| 52 | // >>chng 08 mar 17, from 7.256 to 7.226, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes |
|---|
| 53 | assert line luminosity "totl" 1402 7.226 |
|---|
| 54 | // |
|---|
| 55 | assert line luminosity "totl" 1549 8.177 //total of sum of both C IV comp |
|---|
| 56 | // |
|---|
| 57 | // >>chng 06 jul 18, from 7.450 to 7.403, RP major upgrade VS coll strengths |
|---|
| 58 | // >>chng 06 nov 29, from 7.403 to 7.433, Badnell DR turned on |
|---|
| 59 | // >>chng 08 mar 17, from 7.433 to 7.482, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes |
|---|
| 60 | // >>chng 08 nov 08, from 7.482 to 7.456, fix VS80 elec impact coll. str. (r2444) |
|---|
| 61 | assert line luminosity "he 2" 1640 7.456 //total He II Balmer-alpha 1640 |
|---|
| 62 | // |
|---|
| 63 | // >>chng 03 dec 08, from 7.431 to 7.407, chng zones, ots |
|---|
| 64 | // >>chng 06 nov 29, from 7.407 to 7.429, Badnell DR turned on |
|---|
| 65 | // >>chng 08 mar 17, from 7.429 to 7.398, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes |
|---|
| 66 | // >>chng 08 mar 20, from 7.398 to 7.408, dbl precision gets diff answer, take average, widen |
|---|
| 67 | assert line luminosity "o 3" 1666 7.408 0.10 //total O III] 1666 |
|---|
| 68 | // |
|---|
| 69 | // >>chng 03 oct 04, from 7.044 to 7.066, atom2 ots |
|---|
| 70 | // >>chng 03 dec 08, from 7.310 to 7.277, chng zones, ots |
|---|
| 71 | // >>chng 06 jul 18, from 7.277 to 7.305, RP major upgrade VS coll strengths |
|---|
| 72 | // >>chng 06 nov 29, from 7.305 to 7.346, Badnell DR turned on |
|---|
| 73 | // >>chng 07 jul 25, from 7.346 to 7.323, bugfix pvh FeII tau r1290 |
|---|
| 74 | // >>chng 07 oct 02, increase error, slightly different results for all doubles |
|---|
| 75 | assert line luminosity "totl" 1909 7.323 0.10 //total of sum of CIII] 1908.73+ |
|---|
| 76 | // |
|---|
| 77 | //total of sum of Si III] 1892+ |
|---|
| 78 | // >>chng 04 jul 29, from 8.069 to 8.047, drift |
|---|
| 79 | // >>chng 07 jul 25, from 8.047 to 8.013, bugfix pvh FeII tau r1290 |
|---|
| 80 | // >>chng 07 dec 29, from 8.013 to 7.993, drift down |
|---|
| 81 | assert line luminosity "totl" 1888 7.993 |
|---|
| 82 | // |
|---|
| 83 | // >>chng 04 jul 07, from 8.00 to 7.98, oh ct integ |
|---|
| 84 | // >>chng 06 jul 18, from 7.98 to 8.004, RP major upgrade VS coll strengths |
|---|
| 85 | // >>chng 07 mar 20, from 8.004 to 7.984, due to chng in r987 |
|---|
| 86 | // >>chng 07 jul 25, from 7.984 to 7.903, bugfix pvh FeII tau r1290 |
|---|
| 87 | // >>chng 07 sep 15, from 7.903 to 7.881, drift down |
|---|
| 88 | // >>chng 07 dec 02, from 7.881 to 7.861, update Baddnell DR |
|---|
| 89 | // >>chng 08 mar 17, from 7.861 to 7.829, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes |
|---|
| 90 | // >>chng 08 mar 20, from 7.829 to 7.839, dbl precision gets diff answer, take average, widen |
|---|
| 91 | assert line luminosity "totl" 1860 7.839 0.10 //total sum of both Al III comp |
|---|
| 92 | // |
|---|
| 93 | //total C II] 2324.69 + 2328.12 |
|---|
| 94 | // >>chng 03 dec 08, from 7.564 to 7.538, chng zones, ots |
|---|
| 95 | // >>chng 04 jun 01, from 7.538 to 7.597, unified calc of Lya line width |
|---|
| 96 | // >>chng 04 jul 07, from 7.60 to 7.55, oh ct integ |
|---|
| 97 | // >>chng 04 dec 08, from 7.77 to 7.57, co-ion feedback, check converg |
|---|
| 98 | // >>chng 05 mar 25, from 7.57 to 7.54, renorm H pops in level solver |
|---|
| 99 | // >>chng 06 jul 18, from 7.54 to 7.58, RP major upgrade VS coll strengths |
|---|
| 100 | // >>chng 07 jul 25, from 7.58 to 7.51 bugfix pvh FeII tau r1290 |
|---|
| 101 | // >>chng 07 oct 22, from 7.51 to 7.46, resolve l-levels of h-like sequence |
|---|
| 102 | // >>chng 08 nov 08, from 7.46 to 7.418, fix VS80 elec impact coll. str. (r2444) |
|---|
| 103 | assert line luminosity "totl" 2326 7.418 error 0.1 |
|---|
| 104 | // |
|---|
| 105 | // >>chng 04 jun 12, from 8.295 to 8.321 chng cont shield fcn |
|---|
| 106 | // >>chng 06 jan 24, from 8.321 to 8.300, RP merge He-like, drift down |
|---|
| 107 | // >>chng 06 jul 18, from 8.300 to 8.320, RP major upgrade VS coll strengths |
|---|
| 108 | // >>chng 07 mar 20, from 8.320 to 8.300, due to chng in r987 |
|---|
| 109 | // >>chng 07 jul 25, from 8.300 to 8.262, bugfix pvh FeII tau r1290 |
|---|
| 110 | // >>chng 07 oct 22, from 8.262 to 8.239, resolve l-levels of h-like sequence |
|---|
| 111 | // >>chng 08 mar 20, from 8.239 to 8.219, value drifted but did not throw assert |
|---|
| 112 | assert line luminosity "totl" 2798 8.219 //total of sum of both Mg II comp |
|---|
| 113 | // |
|---|
| 114 | // >>chng 03 dec 08, from 6 to 5 chng zones, ots |
|---|
| 115 | // >>chng 04 feb 09, from 5 to 6 chng zones, ots |
|---|
| 116 | // >>chng 04 sep 30, from 6 to 8, conv_base conv heat cool |
|---|
| 117 | // >>chng 08 dec 08, from 8 to 6.925, upper limit -> equals (r2515) |
|---|
| 118 | // >>chng 08 dec 13, from 6.925 to 18.22, merging newsolvers branch |
|---|
| 119 | // >>chng 09 jan 08, from 18.22 to 18.95, minor changes in iter/zn, iterations, and zone that crept in over last month |
|---|
| 120 | assert itrzn 18.95 error 0.02 |
|---|
| 121 | // |
|---|
| 122 | // >>chng 03 dec 08, from 3 to <5 chng zones, ots |
|---|
| 123 | // >>chng 08 dec 08, from 5 to 3, upper limit -> equals (r2515) |
|---|
| 124 | assert niter 3 error 0.01 |
|---|
| 125 | // |
|---|
| 126 | // >>chng 03 nov 10, from 460 to 510, inprove he0 zoning |
|---|
| 127 | // >>chng 03 dec 08, from 510 to 700 chng zones, ots |
|---|
| 128 | // >>chng 05 dec 22, from 700 to 500 (440), rationalize zone asserts |
|---|
| 129 | // >>chng 08 dec 08, from 500 to 467, upper limit -> equals (r2515) |
|---|
| 130 | // >>chng 08 dec 13, from 467 to 462, merging newsolvers branch |
|---|
| 131 | assert nzone 462 error 0.01 |
|---|
| 132 | // |
|---|
| 133 | c blr_n11_p20_Z20.in |
|---|
| 134 | c class blr |
|---|
| 135 | c ======================================== |
|---|
| 136 | |
|---|
| 137 | This is one of the 5 models that sample the LOC plane. |
|---|