| 1 | title test very soft continuum, very hot gas |
|---|
| 2 | c |
|---|
| 3 | c commands controlling continuum ========= |
|---|
| 4 | c put in the cosmic background as the only continuum source |
|---|
| 5 | CMB |
|---|
| 6 | c gas has constant temperature of 1e6 K |
|---|
| 7 | constant temperature, t=6 |
|---|
| 8 | c |
|---|
| 9 | c commands for density & abundances ========= |
|---|
| 10 | hden 10. |
|---|
| 11 | c |
|---|
| 12 | c commands controlling geometry ========= |
|---|
| 13 | stop zone 1 |
|---|
| 14 | c specify a thin cell of gas - 1 cm thick |
|---|
| 15 | set dr 0 |
|---|
| 16 | c |
|---|
| 17 | c other commands for details ========= |
|---|
| 18 | c iterate since gas is optically thin |
|---|
| 19 | iterate |
|---|
| 20 | c |
|---|
| 21 | c commands controlling output ========= |
|---|
| 22 | punch asserts last "func_hotgas_coolstar.asr" |
|---|
| 23 | c |
|---|
| 24 | c commands giving the asserts ========= |
|---|
| 25 | // >>chng 08 mar 17, from -6.633 to -6.662, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes |
|---|
| 26 | assert ionization Hydrogen 1 -6.662 |
|---|
| 27 | // >>chng 02 apr 29, from -9.273 to -8.657, chng coll ion hi hei lev |
|---|
| 28 | // >>chng 02 nov 20, from -8.657 to -8.701, new he diff, ots, incr errror |
|---|
| 29 | // since so far from ionization peak |
|---|
| 30 | // >>chng 02 dec 11, from -8.701 to -9.147, new hei |
|---|
| 31 | // >>chng 03 feb 28, from -9.147 to -8.980, new hei |
|---|
| 32 | // >>chng 03 may 02, from -8.980 to -9.066, series of small changes |
|---|
| 33 | // >>chng 04 may 19, from -9.066 to -9.67 RP update helike sequence, incr err |
|---|
| 34 | // >>chng 05 mar 12, from -9.67 to -10.79 RP merge He-like code |
|---|
| 35 | // >>chng 06 jul 18, from -10.79 to -10.62 RP major upgrade VS coll strengths |
|---|
| 36 | // >>chng 07 apr 12, from -10.62 to -9.798 helium-like state-specific DR added |
|---|
| 37 | // >>chng 07 may 21, from -9.798 to -9.909, fix bug in He-like DR extrapolation |
|---|
| 38 | // >>chng 07 nov 29, from -9.909 to -9.699, due to changes in r1579 |
|---|
| 39 | assert ionization Helium 1 -9.699 error 0.2 |
|---|
| 40 | // |
|---|
| 41 | assert ionization Helium 2 -4.509 |
|---|
| 42 | // |
|---|
| 43 | // >>chng 02 aug 19, from -0.605 to -0.747, better he-like cs |
|---|
| 44 | // >>chng 02 sep 18, from -0.747 to -0.693, better he-like cs |
|---|
| 45 | // >>chng 02 dec 20, from -0.693 to -0.671, better he-like atom |
|---|
| 46 | // >>chng 03 jan 24, from -0.671 to -0.732, improv cs |
|---|
| 47 | // >>chng 03 feb 28, from -0.732 to -0.636, new hei |
|---|
| 48 | // >>chng 04 may 19, from -0.636 to -0.667 RP update helike sequence |
|---|
| 49 | // >>chng 05 mar 12, from -0.667 to -0.724 RP merge He-like code |
|---|
| 50 | // >>chng 05 jun 05, from -0.724 to -0.761, RP update he-like iso seq |
|---|
| 51 | // >>chng 06 apr 22, from -0.761 to -0.65, RP use Badnell DR for He-like inc err |
|---|
| 52 | // >>chng 06 jul 18, from -0.65 to -0.53 RP major upgrade VS coll strengths |
|---|
| 53 | // >>chng 07 apr 12, from -0.53 to -0.64, he-like proton collision bugfix |
|---|
| 54 | // >>chng 07 jul 27, from -0.64 to -0.58, repair common error in several collision routines |
|---|
| 55 | assert ionization Carbon 5 -0.58 error 0.1 |
|---|
| 56 | // |
|---|
| 57 | // >>chng 02 aug 19, from -0.222 to -0.183, better he-like cs |
|---|
| 58 | // >>chng 03 feb 28, from -0.183 to -0.212, new hei |
|---|
| 59 | // >>chng 05 mar 12, from -0.212 to -0.188 RP merge He-like code |
|---|
| 60 | // >>chng 06 apr 22, from -0.188 to -0.21, RP use Badnell DR for He-like inc err |
|---|
| 61 | // >>chng 06 jul 18, from -0.21 to -0.25 RP major upgrade VS coll strengths |
|---|
| 62 | // >>chng 07 apr 22, from -0.25 to -0.210, he-like proton collision bugfix |
|---|
| 63 | assert ionization Carbon 6 -0.210 error 0.1 |
|---|
| 64 | // |
|---|
| 65 | // >>chng 02 aug 19, from -0.820 to -0.782, better he-like cs |
|---|
| 66 | // >>chng 03 feb 28, from -0.782 to -0.811, new hei |
|---|
| 67 | // >>chng 05 mar 12, from -0.811 to -0.788 RP merge He-like code |
|---|
| 68 | // >>chng 06 apr 22, from -0.788 to -0.808, RP use Badnell DR for He-like |
|---|
| 69 | // >>chng 06 jul 18, from -0.808 to -0.850 RP major upgrade VS coll strengths |
|---|
| 70 | // >>chng 07 apr 12, from -0.850 to -0.810, he-like proton collision bugfix |
|---|
| 71 | // >>chng 07 oct 03, from -0.810 to -0.831, drift down, inc err |
|---|
| 72 | assert ionization Carbon 7 -0.831 error 0.07 |
|---|
| 73 | // |
|---|
| 74 | // >>chng 06 nov 29, from -2.528 to -2.496, Badnell DR turned on |
|---|
| 75 | assert ionization Oxygen 6 -2.496 |
|---|
| 76 | // |
|---|
| 77 | assert ionization Oxygen 7 -0.004 |
|---|
| 78 | // |
|---|
| 79 | // >>chng 02 aug 19, from -2.205 to -1.972, better he-like cs |
|---|
| 80 | // >>chng 02 sep 18, from -1.972 to -2.108, better he-like cs |
|---|
| 81 | // >>chng 02 oct 29, from -2.108 to -2.187, debug he iso seq |
|---|
| 82 | // >>chng 05 jun 05, from -2.187 to -2.168, RP update he-like iso seq |
|---|
| 83 | // >>chng 06 jul 18, from -2.168 to -2.219 RP major upgrade VS coll strengths |
|---|
| 84 | assert ionization Oxygen 8 -2.219 |
|---|
| 85 | // |
|---|
| 86 | // >>chng 03 nov 23, from -1.560 to -1.482, Kraemer et al. guess of dr |
|---|
| 87 | // >>chng 06 jul 12, from -1.482 to -1.450, use Savin & Gu DR to make mean DR |
|---|
| 88 | // >>chng 06 nov 29, from -1.450 to -1.201, Badnell DR turned on |
|---|
| 89 | assert ionization Iron 8 -1.201 |
|---|
| 90 | // |
|---|
| 91 | // >>chng 03 nov 23, from -0.526 to -0.448, Kraemer et al. guess of dr |
|---|
| 92 | // >>chng 06 jul 12, from -0.448 to -0.416, use Savin & Gu DR to make mean DR |
|---|
| 93 | // >>chng 06 nov 29, from -0.416 to -0.371, Badnell DR turned on |
|---|
| 94 | assert ionization Iron 9 -0.371 |
|---|
| 95 | // |
|---|
| 96 | // >>chng 03 nov 23, from -0.484 to -0.438, Kraemer et al. guess of dr |
|---|
| 97 | // >>chng 06 feb 03, from -0.438 to -0.472, guess dr to Badnell mean |
|---|
| 98 | assert ionization Iron 10 -0.472 |
|---|
| 99 | // |
|---|
| 100 | // >>chng 03 nov 23, from -0.632 to -0.703, Kraemer et al. guess of dr |
|---|
| 101 | // >>chng 06 jul 12, from -0.703 to -0.743, use Savin & Gu DR to make mean DR |
|---|
| 102 | // >>chng 06 nov 29, from -0.743 to -0.864, Badnell DR turned on |
|---|
| 103 | assert ionization Iron 11 -0.864 |
|---|
| 104 | // |
|---|
| 105 | // >>chng 03 nov 23, from -1.017 to -1.355, Kraemer et al. guess of dr |
|---|
| 106 | // >>chng 06 feb 03, from -1.355 to -1.166, guess dr to Badnell mean |
|---|
| 107 | // >>chng 06 jul 12, from -1.166 to -1.262, use Savin & Gu DR to make mean DR |
|---|
| 108 | // >>chng 06 nov 29, from -1.262 to -1.548, Badnell DR turned on |
|---|
| 109 | assert ionization Iron 12 -1.548 |
|---|
| 110 | // |
|---|
| 111 | c func_hotgas_coolstar.in |
|---|
| 112 | c class limit |
|---|
| 113 | c ======================================== |
|---|
| 114 | c |
|---|
| 115 | |
|---|
| 116 | This is a test where the CMB is the only continuum source. |
|---|
| 117 | It does not extend to energies where the |
|---|
| 118 | code needs to work. There are special cases |
|---|
| 119 | used in this situation, for continuum addressing, |
|---|
| 120 | so this checks whether those still function. |
|---|