title generate continuum due to hot ism in high Z,z starburst
c
c  TODO the guess of the thickness of the first zone is badly too small,
c  because this model is collisionally ionized, and it used Stromgren
c  length - better to use collisional balance and dr - as result of this
c  the model takes far too many zones
c 
c commands controlling continuum =========
c this sets up coronal equilibrium for temperature and continuum
c cloud is predominantly collisionally ionized
coronal equilibrium, t=1,000,000K
c
c commands for density & abundances =========
c log of hydrogen density (cm-3), by default a constant density model
hden 1
c chemical composition for a well-evolved star cluster with high Z
abundances old solar 84 
abundances starburst Z=10 Zsun
c
c commands controlling geometry  =========
c the log of the total hydrogen column density, cm-2
stop column density 21
c
c other commands for details     =========
c  continuum is plotted in Hzay, we me must iterate to predict this
iterate
c 
c commands controlling output    =========
c normalize to O 7 Lya
normalize to "O  7" 21.60 
c this shortens the printout somewhat
print line faint 1
c output options
punch continuum last "ism_hot_brems.con" no title, units keV
punch overview last "ism_hot_brems.ovr"
punch asserts last "ism_hot_brems.asr"
c this is only to document the problem described above
punch dr last "ism_hot_brems.dr"
c
c commands giving the asserts    =========
// >>chng 01 aug 17, from -5.036 to -4.996, ? ?
// >>chng 04 may 13, from 4.996 to 5.104, dest albedo chng
assert line luminosity "H  1" 4861 -5.104
//
// >>chng 01 mar 17, all lines to luminosity instead of rel intensity
// >>chng 02 may 27, from -1.274 to -1.504, improved ism_hot_brems cooling
// >>chng 02 jun 22, from -1.504 to -1.088, improved ism_hot_brems
assert line luminosity "HFFc" 0 -1.088
//
assert line luminosity "HeFF" 0 -1.605
//
// >>chng 03 nov 23, from -1.581 to -1.606 Kraemer et al. guess of dr
assert line luminosity "MeFF" 0 -1.606
//
assert line luminosity "ToFF" 0 -1.076
//
// >>chng 01 jul 01, from -0.042 to -0.094, change in diff cont transfer
// >>chng 02 feb 14, from -0.094 to -0.072, chng coll rates for he-like ions
// >>chng 02 mar 16, from -0.072 to -0.129, chng he-like ions
// >>chng 02 mar 25, from -0.129 to -0.091, further corrections to above
// >>chng 02 jul 19, from -0.091 to -0.127, better two photon
// >>chng 02 aug 19, from -0.127 to -0.038, better he-like cs
// >>chng 02 dec 20, from -0.038 to -0.0604, improve he-like seq, incr err
// >>chng 03 mar 13, from -0.060 to -0.0200, improve he-like seq, decr err
// >>chng 03 nov 23, from -0.020 to  0.110 Kraemer et al. guess of dr
assert line luminosity "fe10" 171.8  0.110 error 0.1
//
// >>chng 01 jun 15, from -0.073 to -0.120, full CT reaction network
// >>chng 01 jul 01, from -0.120 to -0.147, change in diff cont transfer
// >>chng 01 jul 21, from -0.147 to -0.111, small change, very sensitive
// >>chng 02 aug 19, from -0.119 to -0.094, better he-like cs
// >>chng 03 may 02, from -0.094 to -0.163, several small changes
// >>chng 03 oct 04, from -0.163 to -0.092, atom2 ots
// >>chng 03 nov 23, from -0.109 to -0.151 Kraemer et al. guess of dr
// >>chng 04 may 19, from -0.151 to -0.125, RP update he-like seq
assert line luminosity "fe11" 178.5 -0.125
//
// >>chng 02 jul 19, from -0.081 to -0.-059, better two photon
// >>chng 02 aug 19, from -0.059 to -0.129, better he-like cs
// >>chng 03 jan 24, from -0.129 to -0.166, better he-like cs
// >>chng 03 may 02, from -0.166 to -0.307, several small changes
// >>chng 03 oct 04, from -0.307 to -0.170, atom2 ots
// >>chng 03 nov 23, from -0.170 to -0.560 Kraemer et al. guess of dr
// >>chng 04 may 19, from -0.560 to -0.490, RP update he-like seq
// >>chng 05 jun 05, from -0.490 to -0.455, RP update he-like iso seq
assert line luminosity "fe12" 193.4 -0.455
//
// >>chng 02 jul 19, from -0.111 to -0.200, better two photon
// >>chng 02 aug 19, from -0.200 to -0.030, better he-like cs
// >>chng 03 nov 23, from -0.030 to  0.190 Kraemer et al. guess of dr
// >>chng 05 jun 05, from  0.190 to  0.137, RP update he-like iso seq
assert line luminosity "si 8" 319.1  0.137 error 0.15
//
// >>chng 01 jun 15, from -0.106 to -0.272, full CT reaction network
// >>chng 01 aug 08, from -0.272 to -0.304, ? ? ?
// >>chng 02 feb 10, from -0.304 to -0.254, don't understand, up error to 0.15
// >>chng 02 jul 19, from -0.254 to -0.-0.333, better two photon
// >>chng 02 aug 19, from -0.333 to -0.182, better he-like cs
// >>chng 03 nov 23, from -0.182 to  0.041 Kraemer et al. guess of dr
assert line luminosity "mg 9" 368.5  0.041 error 0.15
//
// >>chng 01 aug 17, -6.832 to -6.738, fe lev2 cs change => ots =>  h ion chng
// >>chng 03 may 02, from -6.738 to -6.761, several small changes
// >>chng 03 oct 08, from -6.761 to -6.712, several small changes
// >>chng 04 may 13, from -6.712 to -7.049, dest albedo chng
assert hydrogen 1 ionization fraction is -7.049
//
// >>chng 01 jun 15, from -5.243 to -5.210, full CT reaction network
// >>chng 01 aug 17, -5.210 to -5.077, fe lev2 cs change => ots =>  h ion chng
// >>chng 03 may 02, from -5.077 to -4.593, several small changes
// >>chng 03 oct 04, from -4.593 to -5.085, atom2 ots
// >>chng 03 nov 23, from -5.085 to -5.052 Kraemer et al. guess of dr
// >>chng 04 may 13, from -5.052 to -5.334, dest albedo chng
assert helium 2 ionization fraction is -5.334
//
// add this, 02 mar 16, o very optically thick, 
// >>chng 02 mar 25, from -0.586 to -0.428, further chngs to he-like
// >>chng 02 may 07, from -0.428 to -0.455, included des of rec photons
// >>chng 02 jun 08, from -0.455 to -0.479, better zoning
// >>chng 02 jun 13, from -0.479 to -0.775, he-like seq
// >>chng 02 aug 19, from -0.775 to -0.362, better he-like cs
// >>chng 02 sep 18, from -0.362 to -0.338, better he-like cs
// >>chng 02 nov 20, from -0.338 to -0.361, new he ots diffuse, incr error
// >>chng 03 jan 24, from -0.361 to -0.244, better he-like cs
// >>chng 04 may 19, from -0.244 to -0.352, RP update he-like seq
// >>chng 05 jun 05, from -0.352 to -0.508, RP update he-like iso seq
assert nitrogen 6 ionization -0.508 error 0.1
//
// >>chng 02 mar 25, from -0.138 to -0.210, further chngs to he-like
// >>chng 02 jun 08, from -0.210 to -0.183, better zoning
// >>chng 02 jun 13, from -0.183 to -0.088, he-like seq
// >>chng 02 aug 19, from -0.088 to -0.255, better he-like cs
// >>chng 02 oct 28, from -0.255 to -0.287, better norm continuum, increase err
// >>chng 03 jan 24, from -0.287 to -0.374, better he-like cs
// >>chng 04 may 19, from -0.374 to -0.263, RP update he-like seq
// >>chng 05 jun 05, from -0.263 to -0.170, RP update he-like iso seq
assert nitrogen 7 ionization -0.170 error 0.1 
//
// add this, 02 mar 16, o very optically thick, all O 7, he like
// >>chng 02 aug 19, from -2.591 to -2.625, better he-like cs
// >>chng 03 may 02, from -2.625 to -2.483, several small changes
// >>chng 03 oct 04, from -2.483 to -2.651, atom2 ots
// >>chng 03 nov 23, from -2.651 to -2.763 Kraemer et al. guess of dr
// >>chng 03 dec 08, from -2.763 to -2.806, chng ots, zones
// >>chng 04 may 13, from -2.806 to -2.871, dest albedo chng
// >>chng 05 jun 05, from -2.871 to -2.837, RP update he-like iso seq
assert oxygen 6 ionization -2.837
//
// >>chng 02 mar 25, from -2.060 to -2.111, further chngs to he-like
// >>chng 02 may 07, from -2.111 to -2.016, included des of rec photons
// >>chng 02 jun 13, from -2.016 to -1.953, he-like seq
// >>chng 02 aug 19, from -1.953 to -2.052, better he-like cs
// >>chng 03 jan 24, from -2.052 to -2.075, better he-like cs
// >>chng 05 jun 05, from -2.075 to -2.037, RP update he-like iso seq
assert oxygen 8 ionization -2.037
//
// >> this changes by about 11% when zoning changes
// >>chng 99 dec 6 caused by changing dr to go from 50 to 10
// >>chng 01 jun 15, from -3.85 to -3.788 full CT reaction network
// >>chng 03 may 02, from -3.788 to -3.268, several small changes
// >>chng 03 oct 04, from -3.268 to -3.808, atom2 ots
// >>chng 03 nov 23, from -3.808 to -3.739 Kraemer et al. guess of dr
// >>chng 04 may 13, from -3.739 to -3.880, dest albedo chng
assert iron 7 ionization fraction is -3.880  0.15
//
// >> this changes by about 7% when zoning changes
// >>chng 99 dec 6 caused by changing dr to go from 50 to 10
// >>chng 02 jul 19, from -2.156 to -2.203, better two photon
// >>chng 02 sep 18, from -2.203 to -2.161, better He-like cs
// >>chng 03 mar 13, from -2.161 to -2.120, better He-like cs and coll ion...how they are split up.
// >>chng 03 may 02, from -2.120 to -1.724, several small changes
// >>chng 03 oct 04, from -1.724 to -2.131, atom2 ots
// >>chng 03 nov 23, from -2.131 to -1.973 Kraemer et al. guess of dr
// >>chng 04 may 13, from -1.973 to -2.047, dest albedo chng
assert iron 8 ionization fraction is -2.047 error 0.1
//
// >>chng 01 jun 15, from -0.852 to -0.885 full CT reaction network
// >>chng 01 jul 01, from -0.885 to -0.945, change in diff cont transfer
// >>chng 02 feb 14, from -0.954 to -0.917, chng coll rates for he-like ions
// >>chng 02 mar 16, from -0.917 to -0.985, chng he-like ions
// >>chng 02 mar 25, from -0.985 to -0.946, further chngs to he-like
// >>chng 02 jul 19, from -0.946 to -0.998, better two photon
// >>chng 02 aug 19, from -0.998 to -0.908, better he-like cs
// >>chng 02 nov 16, from -0.908 to -0.890, increase error
// >>chng 03 may 02, from -0.890 to -0.655, several small changes
// >>chng 03 oct 04, from -0.655 to -0.862, atom2 ots
// >>chng 03 nov 23, from -0.862 to -0.612 Kraemer et al. guess of dr
// >>chng 05 jun 05, from -0.612 to -0.671, RP update he-like iso seq
assert iron 9 ionization fraction is -0.671 error 0.1
//
// >>chng 01 jun 15, from -0.550 to -0.591 full CT reaction network
// >>chng 01 jul 01, from -0.591 to -0.625, change in diff cont transfer
// >>chng 02 feb 14, from -0.625 to -0.604, chng coll rates for he-like ions
// >>chng 02 mar 16, from -0.604 to -0.660, chng he-like ions
// >>chng 02 mar 25, from -0.660 to -0.623, further chngs to he-like
// >>chng 02 jul 19, from -0.623 to -0.658, better two photon
// >>chng 02 aug 19, from -0.658 to -0.580, better he-like cs
// >>chng 03 jan 24, from -0.580 to -0.555, better he-like cs
// >>chng 03 may 02, from -0.555 to -0.529, several small changes
// >>chng 03 oct 04, from -0.529 to -0.553, atom2 ots
// >>chng 03 nov 23, from -0.553 to -0.423 Kraemer et al. guess of dr
// >>chng 05 jun 05, from -0.423 to -0.443, RP update he-like iso seq
assert iron 10 ionization fraction is  -0.443
//
// >>chng 03 may 02, from -0.512 to -0.565, several small changes
// >>chng 03 oct 04, from -0.565 to -0.499, atom2 ots
// >>chng 03 nov 23, from -0.499 to -0.558 Kraemer et al. guess of dr
// >>chng 04 may 19, from -0.558 to -0.532, RP update he-like seq
assert iron 11 ionization fraction is  -0.532
//
// >> next three changed when zoning changes
// >>chng 01 jun 29, from -0.687 to -0.644, small temp chng ?
// >>chng 03 may 02, from -0.644 to -0.823, several small changes
// >>chng 03 oct 04, from -0.823 to -0.695, atom2 ots
// >>chng 03 nov 23, from -0.695 to -1.083 Kraemer et al. guess of dr
// >>chng 05 jun 05, from -1.083 to -0.980, RP update he-like iso seq
assert iron 12 ionization fraction is  -0.980 error 0.2
//
// >>chng 01 feb 21 from -1.247 to -1.194 improved He-seq coll treatment
// >>chng 01 jun 15, from -1.194 to -1.141 full CT reaction network
// >>chng 01 jul 01, from -1.141 to -1.079, change in diff cont transfer
// >>chng 02 feb 14, from -1.079 to -1.155, chng coll rates for he-like ions
// >>chng 02 mar 16, from -1.155 to -1.023, chng he-like ions
// >>chng 02 mar 25, from -1.023 to -1.112, further chngs to he-like
// >>chng 02 jul 19, from -1.112 to -1.016, better two photon
// >>chng 02 aug 19, from -1.016 to -1.217, better he-like cs
// >>chng 03 jan 24, from -1.217 to -1.281, better he-like cs
// >>chng 03 may 02, from -1.281 to -1.435, several small changes
// >>chng 03 oct 04, from -1.435 to -1.289, atom2 ots
// >>chng 03 nov 23, from -1.289 to -2.126 Kraemer et al. guess of dr
// >>chng 04 may 19, from -2.126 to -2.008, RP update he-like seq
// >>chng 05 jun 05, from -2.008 to -1.937, RP update he-like iso seq
assert iron 13 ionization fraction is  -1.937 error 0.1
//
// >>chng 01 feb 21 from -2.224 to -2.128 improved He-seq coll treatment
// >>chng 01 jun 15, from -2.128 to -2.028 full CT reaction network
// >>chng 01 jul 01, from -2.028 to -1.954, change in diff cont transfer
// >>chng 01 jul 21, from -1.954 to -2.015, very sensitive
// >>chng 01 aug 08, from -2.015 to -1.912, ? ? ?
// >>chng 02 feb 14, from -1.912 to -2.105, chng coll rates for he-like ions
// >>chng 02 mar 16, from -2.105 to -1.864, chng he-like ions
// >>chng 02 mar 25, from -1.864 to -2.030, further chngs to he-like
// >>chng 02 jul 19, from -2.030 to -1.865, better two photon
// >>chng 02 aug 19, from -1.865 to -2.219, better he-like cs
// >>chng 03 jan 28, from -2.219 to -2.316, slow drift
// >>chng 03 may 02, from -2.313 to -2.487, several small changes
// >>chng 03 oct 04, from -2.487 to -2.322, atom2 ots
// >>chng 03 nov 23, from -2.322 to -3.580 Kraemer et al. guess of dr
// >>chng 04 may 19, from -3.580 to -3.400, RP update he-like seq
// >>chng 05 jun 05, from -3.400 to -3.285, RP update he-like iso seq
assert iron 14 ionization fraction is  -3.285 error 0.20
//
// >> changed by 5% when zoning changes
// >>chng 99 dec 6 caused by changing dr to go from 50 to 10
// >>chng 01 feb 21 from -3.359 to -3.159 improved He-seq coll treatment
// >>chng 01 jun 15, from -3.159 to -2.973 full CT reaction network
// >>chng 01 jul 01, from -2.973 to -2.880, change in diff cont transfer
// >>chng 02 feb 14, from -2.880 to -3.149, chng coll rates for he-like ions
// >>chng 02 mar 16, from -3.149 to -2.765, chng he-like ions
// >>chng 02 mar 25, from -2.765 to -3.039, further chngs to he-like
// >>chng 02 jul 19, from -3.039 to -2.792, better two photon
// >>chng 02 aug 19, from -2.792 to -3.365, better he-like cs
// >>chng 02 dec 20, from -3.365 to -3.247, better he-like atom
// >>chng 03 jan 24, from -3.247 to -3.508, better he-like cs
// >>chng 03 may 02, from -3.508 to -3.701, several small changes
// >>chng 03 oct 04, from -3.701 to -3.535, atom2 ots
// >>chng 03 nov 23, from -3.535 to -9.824 Kraemer et al. guess of dr
// assert iron 15 ionization fraction is  -9.824 error 0.2
//
// >>chng 02 jun 08, from 75 to 192, smaller rel dr step
// >>chng 03 may 22, from 200 to 20 (??)
// >>chng 03 dec 10, from 20 to 30 dr logic change
// >>chng 04 mar 04, from <30 to <40, chng conver criteria
// >>chng 04 may 19, from -40 to -50 RP update he-like seq
// >>chng 05 jul 16, from 15 to 150, some versions of gcc resolve
// changes in the ionization, that most compilers do not
assert nzone < 150
//
// >>chng 03 nov 23, from 6 to 8 Kraemer et al. guess of dr
// >>chng 04 feb 15, from 8 to 7, conv ionz loop auger
assert itrzn < 7
//
// check that stop column gets exactly right thickness 
assert thickness 20 error 0.001
c ism_hot_brems.in
c class ism 
c ========================================
c 

This model generates a large column constant density cloud
similar to the hot phase of the interstellar medium. 
The continuum is punched to generate one of the figures in Part 2 of Hazy.

There is a strange feature between 1.7e-3A and 2.2e-3A that is the N emission
(head starting at 1.7e-3A with O absorption at 2.2e-3A.  This model is
strongly enriched in heavies so many metal edges, esp O, are optically thick.