title constant pressure H+ region/pdr
c
c commands controlling continuum =========
c the incident continuum is two parts
c star and flux of photons striking it
// >>chng 03 sep 12, now go deep into omc1 and cannot assume plane par geometry
// >>chng 11 dec 21, increase column density to go well past CO dissociation front
c this is the photosphere of the OVI star, its temperature and Q(H)
c this is the result of this command
c table star kurucz 39600 K
c and makes it possible to run these orion sims without
c installing the stellar atmosphere files 
table read "star_kurucz_39600.dat" 
Q(H) 49
c this adds the observed hot brems, its temperature and the flux of 
c photons striking the cloud
brems 6
phi(h) 10
// >>chng 03 apr 30, add cosmic rays to model
c add cosmic rays, which are important in the pdr
cosmic rays, background
c
c commands controlling geometry  =========
c this sets the separation from the star and the face of the cloud
radius 17.4507
// >>chng 00 oct 13, from 50 to 10
// >>chng 04 nov 09, introduce off option
c this turns off the stop temperature opption
c the sim will not stop due to temperature
stop temperature off
c this is typical of a gmc, larson 1981
// >>chng 11 dec 21, 22 to 22.5 to get well past CO dissociation front
stop total hydrogen column density 23
// >>chng 02 dec 21, change to set thickness, with new hmole
// co hits asymptote at 0.9 and so does not stop due to co
// >>chng 03 sep 12, co network cleaned up, reset thickness to
// 1 parsec, roughly the width of the cloud
c this sets the thickness of the HII region & PDR
//>>chng 05 nov 22, from 1 to 0.5 pc, so stop before ices fully form
stop thickness 0.5 linear parsec 
c constant total pressure, so like orion_hii_open
constant pressure 
c this will result in a milli gauss B-field in molecular region
magnetic field -5 gauss
c the orion environment full covers the star, so turn on sphere
sphere 
c we have a spherical geometry but want to simulate observing
c through a spectrometer's slit.  use the aperture
c command for this
aperture beam 
c
c other commands for details     =========
c mimic existence of unmodeled molecular gas
double 
c iterate since lines optically thick
iterate
c the observed microturbulence, partially a flow, so not included in pressure
turbulence 8 km/sec no pressure
c set the line width so lines appear on the save continuum
set save line width 10 km/s
c there should be no failures at all
failures 3
c
c commands for density & abundances =========
c this is the log of the initial H density, cm-3
hden 4
c this will speed up the calculation a bit
init file="ism.ini"
c this uses HII region abundances, but no grains
abundances hii region no grains
c this turns on orion grains, but leaves quantum heating off to save time
grains orion no qheat
c turn on PAHs, with an abundance that depends on H0 fraction,
c as suggested by long-slit observations of Orion bar,
c with an abundance 3x larger than default built into the code
grains PAH function 3
c 
c commands controlling output    =========
normalize to "O  1" 63.17m
c print lots of faint CO lines
print line faint -4
save performance "orion_hii_pdr.per"
save overview last "orion_hii_pdr.ovr"
save heating "orion_hii_pdr.het"
save cooling "orion_hii_pdr.col"
save dr last "orion_hii_pdr.dr"
save results last "orion_hii_pdr.rlt"
save continuum last "orion_hii_pdr.con" units microns
save hydrogen 21 cm last "orion_hii_pdr.21cm"
save hydrogen lya last "orion_hii_pdr.lya" 
save element silicon last "orion_hii_pdr.sil" 
//
// these are for the grains
save grain extinction last  "orion_hii_pdr.grnext"
save grain charge last "orion_hii_pdr.grnchr" 
save grain potential last "orion_hii_pdr.grnpot" 
save grain temperature last "orion_hii_pdr.grntem" 
c
c commands giving the monitors    =========
c
// >>chng 03 sep 12, from plan par to spherical geometry,
// all lines changed, so removed history
//
// >>chng 04 apr 02, from 0.5978 to 0.6171, probably chng in grn abund
// >>chng 04 may 10, from 0.617 to 0.648, ? ?
// >>chng 06 aug 25, from 0.648 to 0.708, add mag field and stop col 22, 
// >>chng 08 mar 17, from 0.708 to 0.732, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 jul 08, from 0.732 to 0.697 line radiation pressure
// >>chng 08 jul 23, fromo 0.697 to 0.729 chng Pesc line radiation pressure
monitor line luminosity "h  1" 4861 0.729
//
// >>chng 04 jan 23, from -1.272 to -1.229, rescale orion dust to gas
// >>chng 04 apr 23, from -1.229 to -1.247, slow drift
// >>chng 04 jul 30, from -1.247 to -1.357, NA implement N chem
// >>chng 04 aug 01, from -1.357 to -1.251, improv C converg
// >>chng 04 sep 30, from -1.251 to -1.175, add pahs 
// >>chng 06 may 30, NA refactor CO, fix H+CH=>C+H2 rate, from -1.175 to-1.153
// >>chng 06 aug 25, from -1.153 to -1.196 add mag field and stop col 22, 
// >>chng 11 dec 12, from -1.196 to -1.237, bug-fix for updating grain opacity (r5759)
monitor line luminosity "c  2" 157.6m -1.237 
//
// >>chng 06 aug 25, from -0.953 to -1.778, add mag field and stop col 22, 
// >>chng 06 oct 25, from -1.778 to -1.733, TE update H2 heating for sphere geo
// >>chng 06 nov 29, from -1.733 to -1.625, Badnell DR turned on
// >>chng 06 dec 25, from -1.625 to -2.905, mole CT reset bug
// >>chng 08 feb 20, from -2.905 to -2.941, grain charge transfer changes in ion_solver and iso*
// >>chng 09 jul 14, from -2.941 to -3.083 update NI atomic data
// >>chng 09 oct 05, from -3.083 to -1.738, [NI] pumping changes
// >>chng 11 dec 12, from -1.738 to -1.789, bug-fix for updating grain opacity (r5759)
// >>chng 13 jun 17, from -1.789 to -1.756, improve in-zone shielding
monitor line luminosity "TOTL"  5199A  -1.756
//
// >>chng 04 jan 23, from -0.0695 to -0.0398, rescale orion dust to gas
// >>chng 04 mar 04, from -0.0398 to -0.0674, improv conv
// >>chng 04 mar 24, from -0.0674 to 0.0221, add PAHs
// >>chng 04 apr 17, from 0.0221 to -0.0667, grn recom, better mole conv
// >>chng 04 jun 28, from -0.0667 to -0.0408, bug in part count for cons pres
// >>chng 04 jul 29, from -0.0408 to -0.017, chng H conv, incr err
// >>chng 04 sep 30, from -0.017 to 0.131, add pahs 
// >>chng 06 aug 25, from 0.131 to 0.253 add mag field and stop col 22, 
// >>chng 11 dec 12, from 0.253 to 0.1942, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from 0.1942 to 0.1443, H column density 1e22 -> 1e23
monitor line luminosity "O  1" 63.17m 0.1443 erro 0.1 
//
// >>chng 04 feb 09, from -1.715 to -1.692, update H2 net
// >>chng 04 mar 04, from -1.692 to -1.688, improv conv
// >>chng 04 apr 20, from -1.688 to -1.721, Nick Abel improve CO post Leiden
// >>chng 04 jul 03, from -1.721 to -1.681, NA update H, CO chem net
// >>chng 06 jul 06, from -1.681 to -1.269, update Ne cs
// >>chng 06 aug 12, from -1.269 to -1.37, dr from tau IR continuum, inc err
// >>chng 06 aug 17, undo above zoning from grains, chng is error
// >>chng 06 aug 25, from -1.27 to -1.40, add mag field and stop col 22, 
// >>chng 06 oct 02, from -1.40 to -1.30 HM update OI cs to Barklem '06
// >>chng 06 dec 25, from -1.30 to -1.37, mole CT reset bug
monitor line luminosity "O  1" 6300 -1.37 error 0.15 
//
// >>chng 06 jul 06, from -0.0201 to -0.0004, update Ne cs
// >>chng 06 jul 25, from -0.0004 to 0.0686 HN update OIII, OII collision data
// >>chng 06 aug 25, from 0.069 to 0.208, add mag field and stop col 22, 
// >>chng 06 nov 29, from 0.208 to 0.175, Badnell DR turned on
// >>chng 07 dec 14, from 0.175 to 0.153 RP update is seq
// >>chng 08 mar 17, from 0.153 to 0.188, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 jul 08, from 0.188 to 0.129 line radiation pressure
// >>chng 08 jul 23, fromo 0.129 to 0.192 chng Pesc line radiation pressure
// >>chng 10 feb 14, from 0.192 to 0.1213, update O atomic data
// >>chng 11 jul 12, from 0.1213 to 0.0798, enable fit to Badnell DR suppression
// >>chng 13 feb 02, from 0.0798 to 0.1184, BUGFIX, largely undo r6703 
monitor line luminosity "totl" 7325 0.1184
//
// >>chng 06 jul 06, from 0.752 to 0.783, update Ne cs
// >>chng 06 jul 25, from 0.783 to 0.808 HN update OIII, OII collision data
// >>chng 06 nov 29, from -0.808 to -0.763, Badnell DR turned on
// >>chng 08 jul 08, from 0.763 to 0.746 line radiation pressure
// >>chng 08 jul 23, from 0.746 to 0.768 chng Pesc line radiation pressure
// >>chng 11 jul 12, from 0.768 to 0.7368, enable fit to Badnell DR suppression
// >>chng 13 feb 02, from 0.7368 to 0.7758, BUGFIX, largely undo r6703 
monitor line luminosity "totl" 3727 0.7758
//
// >>chng 06 nov 29, from -1.248 to -1.273, Badnell DR turned on
// >>chng 08 jul 08, from 1.273 to 1.250 line radiation pressure
// >>chng 08 jul 23, from 1.250 to 1.275 chng Pesc line radiation pressure
// >>chng 11 jul 12, from 1.275 to 1.311, enable fit to Badnell DR suppression
// >>chng 13 feb 02, from 1.311 to 1.288, BUGFIX, largely undo r6703 
monitor line luminositry "o  3" 5007 1.288
//
// >>chng 04 may 10, from -1.094 to -1.155 ? ?
// >>chng 06 jul 25, from -1.155 to -1.123 HN update OIII, OII collision data
// >>chng 06 aug 25, from -1.123 to -1.101, add mag field and stop col 22, 
// >>chng 06 nov 29, from -1.101 to -1.081, Badnell DR turned on
// >>chng 08 jul 08, from -1.081 to -1.116 line radiation pressure
// >>chng 08 jul 23, from -1.116 to -1.083 chng Pesc line radiation pressure
// >>chng 09 feb 23, from -1.083 to -1.063, remove otslin, otscon dampers (PR #84, r2758)
// >>chng 11 nov 16, from -1.063 to -1.043, adjust newmole monitors
monitor line luminosity "totl" 4363 -1.043
//
// >>chng 04 jun 09, from -0.0096 to 0.0073, NA update H chem net
// >>chng 06 aug 25, from 0.007 to -0.046, add mag field and stop col 22, 
// >>chng 06 nov 29, from -0.046 to -0.016, Badnell DR turned on
monitor line luminosity "o  3" 51.80m -0.016
//
//>>chng 06 jul 06, add this line
// >>chng 06 aug 25, from 0.169 to 0.233, add mag field and stop col 22, 
// >>chng 08 jul 08, from 0.233 to 0.210 line radiation pressure
// >>chng 08 jul 23, from 0.210 to 0.248 chng Pesc line radiation pressure
monitor line luminosity "Ne 2" 12.81m 0.248
//
// >>chng 03 sep 27, from -0.863 to -0.891, better conv
// >>chng 03 oct 30, from -0.891 to -0.868, slow drift up
// >>chng 03 nov 16, from -0.868 to -0.829, H2 solomon th85 -> bd96
// >>chng 04 jan 23, from -0.829 to -0.785, rescale orion dust to gas
// >>chng 04 feb 09, from -0.785 to -0.760, update H2 net
// >>chng 04 feb 21, from -0.760 to -0.828, update chem, H- ct, grn recom
// >>chng 04 mar 24, from -0.828 to 0.744, add PAHs
// >>chng 04 apr 17, from -0.744 to -0.835, grn recom, better mole conv
// >>chng 04 jun 28, from -0.835 to -0.808, bug in part count for cons pres
// >>chng 04 sep 30, from -0.808 to -0.644, add pahs 
// >>chng 05 aug 03, from -0.644 to -0.602, NA update CO, H chem net
// >>chng 08 feb 20, from -0.602 to -0.560, grain charge transfer changes in ion_solver and iso*
// >>chng 08 jul 08, from -0.560 to -0.583 line radiation pressure
// >>chng 11 dec 12, from -0.583 to -0.631, bug-fix for updating grain opacity (r5759)
monitor line luminosity "Si 2" 34.81m -0.631 error 0.1
//
// >>chng 04 jan 23, from 3.858 to 3.790, rescale orion dust to gas
// >>chng 04 mar 04, from -3.858 to -3.821 improv conv
// >>chng 04 mar 24, from 3.821 to 3.758, add PAHs
// >>chng 04 apr 17, from 3.758 to 3.818, grn recom, better mole conv
// >>chng 04 may 13, from 3.818 to 3.789, dest albedo chng
// >>chng 04 may 29, from 3.789 to 3.485, NA update H chem
// >>chng 04 jun 28, from 3.485 to 3.630, bug in part count for cons pres
// >>chng 06 aug 12, from 3.630 to 3.653, dr from tau IR continuum, inc err
// >>chng 06 aug 25, from 3.653 to 3.143 add mag field and stop col 22, 
// >>chng 08 jul 08, from 3.143 to 3.122 line radiation pressure
// >>chng 08 jul 23, from 3.122 to 3.149 chng Pesc line radiation pressure
// >>chng 11 dec 21, from 3.149 to 3.276, H column density 1e22 -> 1e23
monitor line luminosity "GraT" 0 3.276
//
// >>chng 04 mar 24, from 24.759 to 24.596 add PAHs
// >>chng 04 apr 17, from 24.596 to 24.771, grn recom, better mole conv
// >>chng 04 may 29, from 24.771 to 24.976, NA update H chem
// >>chng 04 jun 28, from 24.976 to 25.462, bug in part count for cons pres
// >>chng 05 nov 22, from 25.46 to 24.99 stop half thick to avoid grn freeze out
// >>chng 06 aug 25, from 24.99 to 21.13 add mag field and stop col 22, 
// >>chng 11 nov 16, from 21.13 to 21.06, adjust newmole monitors
// >>chng 11 dec 12, from 21.06 to 21.22, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from 21.22 to 22.67, H column density 1e22 -> 1e23
monitor H2 column density 22.67 error 0.1
//
// >>chng 03 sep 27, from 21.391 to 21.527, better conv
// >>chng 04 mar 24, from 21.527 to 21.373 add PAHs
// >>chng 04 apr 17, from 21.373 to 21.549, grn recom, better mole conv
// >>chng 04 may 29, from 21.549 to 21.750, NA update H chem
// >>chng 04 jun 28, from 21.750 to 22.234, bug in part count for cons pres
// >>chng 05 jun 07, from 22.23 to 21.90, NA add CO grain condensation
// >>chng 05 nov 22, from 21.90 to 21.62 stop half thick to avoid grn freeze out
// >>chng 05 dec 17, from 21.62 to 21.70, NA increase C chemistry
// >>chng 06 aug 25, from 21.70 to 16.30 add mag field and stop col 22, 
// >>chng 06 nov 29, from 16.30 to 16.45 Badnell DR turned on
// >>chng 07 jan 15, from 16.45 to 16.37, small model H2 bugfix
// >>chng 07 dec 06, from 16.37 to 16.46, changes in h-like collisions
// >>chng 08 feb 20, from 16.46 to 16.25, grain charge transfer changes in ion_solver and iso*
// >>chng 08 mar 17, from 16.25 to 16.17, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 jul 08, from 16.17 to 16.27 line radiation pressure
// >>chng 08 jul 23, from 16.27 to 16.17 chng Pesc line radiation pressure
// >>chng 09 oct 05, from 16.17 to 16.28, [NI] pumping changes
// >>chng 10 feb 14, from 16.28 to 16.34, update O atomic data
// >>chng 10 mar 14, from 16.34 to 16.26, update As to NIST recommended values
// >>chng 10 jul 26, from 16.26 to 16.29, icc mac, inc error
// >>chng 11 jul 12, from 16.29 to 16.1, enable fit to Badnell DR suppression
// >>chng 11 aug 05, from 16.10 to 16.17, add helike quadrupole transitions
// >>chng 11 nov 16, from 16.17 to 14.16, adjust newmole monitors
// >>chng 11 dec 12, from 14.16 to 14.49, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from 14.49 to 19.39, H column density 1e22 -> 1e23
monitor CO column density 19.39 error 0.20
//
// >>chng 04 apr 20, add these four
// >>chng 04 may 19, from 18.11 to 18.04, NA add several moles, reactions
// >>chng 04 may 29, from 18.04 to 17.99, NA update H chem
// >>chng 04 jun 28, from 17.99 to 18.14, bug in part count for cons pres
// >>chng 04 jul 03, from 18.14 to 18.31, NA update H, CO chem net
// >>chng 04 jul 17, from 18.31 to 18.51, NA update chem net
// >>chng 04 jul 30, from 18.51 to 18.65, NA implement N chem
// >>chng 04 aug 01, from 18.65 to 18.77, improv C converg
// >>chng 04 aug 09, from 18.77 to 18.81, small chng grn eden conv, inc err
// >>chng 04 aug 15, from 18.81 to 18.52, NA adjust ion conv, inc err
// >>chng 04 sep 08, from 18.52 to 18.42, rewrite auger/grain rec ion solver
// >>chng 04 oct 19, from 18.42 to 18.50, small change in d Mole/ d Zone
// >>chng 05 mar 25, from 18.50 to 18.59 renorm H pops in level solver
// >>chng 05 jun 07, from 18.59 to 19.39, NA add CO grain condensation
// >>chng 05 aug 03, from 19.39 to 18.91, NA update CO, H chem net
// >>chng 05 nov 22, from 18.91 to 17.64 stop half thick to avoid grn freeze out
// >>chng 06 mar 02, from 17.64 to 17.74, NA Fe0, Mg0, <-> Si+, S+, C+ ct,
// also energy barrier for H + (CH and CH2) reactions
// >>chng 06 aug 12, from 17.74 to 17.83 dr from tau IR continuum, inc err
// >>chng 06 aug 25, from 17.83 to 11.24 add mag field and stop col 22, 
// >>chng 11 dec 21, from 11.24 to 16.01, H column density 1e22 -> 1e23
// >>chng 12 apr 09, from 16.01 to 16.31, update CR backbround to Indriolo+07, ~8x larger
monitor H2O column density 16.31 error 0.20
// 
// >>chng 04 may 19, from 20.27 to 19.29, NA add several moles, reactions
// >>chng 04 may 29, from 19.29 to 18.85, NA update H chem
// >>chng 04 jun 28, from 18.85 to 19.00, bug in part count for cons pres
// >>chng 04 jul 03, from 19.00 to 19.27, NA update H, CO chem net
// >>chng 04 jul 17, from 19.27 to 19.22, NA update chem net
// >>chng 04 jul 30, from 19.22 to 19.08, NA implement N chem
// >>chng 04 aug 01, from 19.08 to 19.16, improv C converg
// >>chng 04 aug 15, from 19.16 to 19.20, NA adjust ion conv, inc err
// >>chng 05 jun 07, from 19.20 to 19.87, NA add CO grain condensation
// >>chng 05 jul 04, from 19.87 to 18.81 turn cl on
// >>chng 05 aug 03, from 19.39 to 19.34, NA update CO, H chem net, inc err
// >>chng 05 nov 22, from 19.34 to 18.12 stop half thick to avoid grn freeze out
// >>chng 06 aug 25, from 18.12 to 10.71 add mag field and stop col 22, 
// >>chng 07 jan 07, from 10.71 to 10.77, enable McCall et al. H3+ rate
// >>chng 11 jul 12, from 10.77 to 10.67, enable fit to Badnell DR suppression
// >>chng 11 nov 16, from 10.67 to 10.43, adjust newmole monitors
// >>chng 11 dec 21, from 10.43 to 16.21, H column density 1e22 -> 1e23
// >>chng 12 apr 09, from 16.21 to 16.41, update CR backbround to Indriolo+07, ~8x larger
monitor O2 column density 16.41 error 0.2
// 
// >>chng 04 may 19, from 18.41 to 17.95, NA add several moles, reactions
// >>chng 04 jun 09, from 17.95 to 17.87 NA update H chem net
// >>chng 04 jun 28, from 17.87 to 18.21, bug in part count for cons pres
// >>chng 04 jul 03, from 18.21 to 18.33, NA update H, CO chem net
// >>chng 04 jul 17, from 18.33 to 18.62, NA update chem net
// >>chng 04 jul 30, from 18.62 to 19.27, NA implement N chem
// >>chng 04 aug 01, from 19.27 to 19.03, improv C converg
// >>chng 04 aug 04, from 19.03 to 18.98, finer zones, smaller d CO net, inc err
// >>chng 04 aug 09, from 18.98 to 19.09, small chng grn eden conv, inc err
// >>chng 04 aug 15, from 19.09 to 18.68, NA adjust ion conv, inc err
// >>chng 05 mar 25, from 18.68 to 18.84 renorm H pops in level solver
// >>chng 05 jun 07, from 18.84 to 19.24, NA add CO grain condensation
// >>chng 05 nov 22, from 19.24 to 18.16 stop half thick to avoid grn freeze out
// >>chng 06 apr 15, from 18.16 to 18.27 NA added N2H+ to the chemistry
// >>chng 06 aug 25, from 18.27 to 10.10 add mag field and stop col 22, 
// >>chng 06 nov 29, from 10.10 to 10.25 Badnell DR turned on
// >>chng 08 feb 20, from 10.25 to 10.11, grain charge transfer changes in ion_solver and iso*
// >>chng 08 jul 23, from 10.11 to 9.99 chng Pesc line radiation pressure
// >>chng 09 feb 23, from 9.99 to 10.09, remove otslin, otscon dampers (PR #84, r2758)
// >>chng 10 feb 14, from 10.09 to 10.18, update O atomic data
// >>chng 11 jul 12, from 10.18 to 9.957, enable fit to Badnell DR suppression
// >>chng 11 nov 16, from 9.957 to 9.56, adjust newmole monitors
// >>chng 11 dec 12, from 9.56 to 9.796, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from 9.796 to 15.87, H column density 1e22 -> 1e23
// >>chng 12 apr 09, from 15.87 to 16.32, update CR backbround to Indriolo+07, ~8x larger
monitor SiO column density 16.32 error 0.25
// 
// >>chng 04 may 19, from 15.60 to 15.43, NA add several moles, reactions
// >>chng 04 may 29, from 15.43 to 15.28, NA update H chem
// >>chng 04 jun 28, from 15.28 to 15.37, bug in part count for cons pres
// >>chng 04 jul 03, from 15.37 to 15.61, NA update H, CO chem net
// >>chng 04 jul 17, from 15.61 to 15.66, NA update chem net
// >>chng 04 jul 30, from 15.66 to 15.49, NA implement N chem
// >>chng 04 aug 01, from 15.49 to 15.62, improv C converg
// >>chng 04 oct 19, from 15.62 to 15.68, larger step size in dMole/dZone
// >>chng 05 jun 07, from 15.68 to 16.20, NA add CO grain condensation
// >>chng 05 nov 22, from 16.20 to 14.94 stop half thick to avoid grn freeze out
// >>chng 06 aug 12, from 14.94 to 15.01 dr from tau IR continuum, inc err
// >>chng 06 aug 17, undo above zoning from grains, chng is error
// >>chng 06 aug 25, from 14.94 to 12.71 add mag field and stop col 22, 
// >>chng 08 mar 17, from 12.71 to 12.78, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 jul 08, from 12.78 to 12.70 line radiation pressure
// >>chng 08 jul 23, from 12.70 to 12.77 chng Pesc line radiation pressure
// >>chng 11 dec 12, from 12.77 to 12.55, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from 12.55 to 13.03, H column density 1e22 -> 1e23
// >>chng 12 apr 09, from 13.03 to 13.56, update CR backbround to Indriolo+07, ~8x larger
monitor OH column density 13.56 error 0.15
//
// >>chng 03 oct 30, from -3.206 to -3.224, slow drift up
// >>chng 04 mar 04, from -3.224 to -3.302 improv conv, incr error
// >>chng 04 mar 24, from -3.302 to -3.128 add PAHs
// >>chng 04 apr 02, from -3.138 to -3.094, update co photo dissoc rate
// >>chng 04 apr 17, from -3.094 to -3.272, grn recom, better mole conv
// >>chng 04 jun 09, from -3.272 to -3.242 NA update H chem net
// >>chng 04 jun 28, from -3.24 to -3.13, bug in part count for cons pres
// >>chng 04 jul 03, from -3.13 to -3.20, NA update H, CO chem net
// >>chng 05 feb 18, from -3.20 to -3.17, incld outwrd cont in line pump
// >>chng 06 may 30, NA refactor CO, fix H+CH=>C+H2 rate, from -3.17 to -3.08
// >>chng 06 aug 25, from -3.08 to -4.01 add mag field and stop col 22, 
// >>chng 04 sep 11, from -4.01 to -4.07, turn on CO net when H+/H<0.5, inc err
// >>chng 06 nov 24, from -4.07 to -4.03, split diff betweeen icl, vs05, inc er
// >>chng 06 nov 29, from -4.04 to -3.89 Badnell DR turned on
// >>chng 08 feb 20, from -3.89 to -4.06, grain charge transfer changes in ion_solver and iso*
// >>chng 08 mar 20, from -4.06 to -4.15 value drifted but did not throw assert
// >>chng 08 jul 08, from -4.15 to -4.02 line radiation pressure
// >>chng 08 jul 23, from -4.02 to -4.15 chng Pesc line radiation pressure
// >>chng 09 oct 05, from -4.15 to -4.052, [NI] pumping changes
// >>chng 10 mar 02, from -4.052 to -3.97, fix botches after CO model switch
// >>chng 11 jul 12, from -3.97 to -4.185, enable fit to Badnell DR suppression
// >>chng 11 aug 05, from -4.185 to -4.103, add helike quadrupole transitions
// >>chng 11 nov 16, from -4.103 to -6.125, adjust newmole monitors
// >>chng 11 dec 12, from -6.125 to -5.741, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from -5.741 to -2.921, H column density 1e22 -> 1e23
monitor line luminosity "CO  " 650.1m  -2.921 error 0.20
//
// >>chng 04 jan 23, from -2.935 to -2.960 rescale orion dust to gas
// >>chng 04 mar 04, from -2.960 to -3.029 improv conv, incr error
// >>chng 04 mar 24, from -3.029 to -2.838 add PAHs
// >>chng 04 apr 02, from -2.838 to -2.806, update co photo dissoc rate
// >>chng 04 apr 17, from -2.804 to -3.004, grn recom, better mole conv
// >>chng 04 jun 28, from -3.00 to -2.86, bug in part count for cons pres
// >>chng 04 jul 03, from -2.86 to -2.93, NA update H, CO chem net
// >>chng 05 aug 03, from -2.93 to -2.90, NA update CO, H chem net, inc err
// >>chng 06 mar 02, from -2.90 to -2.87, NA Fe0, Mg0, <-> Si+, S+, C+ ct,
// also energy barrier for H + (CH and CH2) reactions
// >>chng 06 may 30, NA refactor CO, fix H+CH=>C+H2 rate, from -2.87 to -2.80
// >>chng 06 aug 25, from -2.80 to -3.72 add mag field and stop col 22, 
// >>chng 04 sep 11, from -3.72 to -3.77, turn on CO net when H+/H<0.5, inc err
// >>chng 06 nov 24, from -3.77 to -3.73, split diff betweeen icl, vs05, inc er
// >>chng 06 nov 29, from -3.73 to -3.60 Badnell DR turned on
// >>chng 08 feb 20, from -3.60 to -3.76, grain charge transfer changes in ion_solver and iso*
// >>chng 08 apr 03, from -3.76 to -3.90 include 10830 in fine opacity array
// >>chng 08 jul 08, from -3.90 to -3.73 line radiation pressure
// >>chng 08 jul 23, from -3.72 to -3.84 chng Pesc line radiation pressure
// >>chng 09 oct 05, from -3.84 to -3.747, [NI] pumping changes
// >>chng 10 mar 02, from -3.747 to -3.622, fix botches after CO model switch
// >>chng 11 jul 12, from -3.622 to -3.83, enable fit to Badnell DR suppression
// >>chng 11 aug 05, from -3.83 to -3.754, add helike quadrupole transitions
// >>chng 11 nov 16, from -3.754 to -5.764, adjust newmole monitors
// >>chng 11 dec 12, from -5.764 to -5.429, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from -5.429 to -2.715, H column density 1e22 -> 1e23
monitor line luminosity "CO  " 520.1m  -2.715 error 0.20
//
// >>chng 04 mar 04, from -2.723 to -2.810 improv conv, incr error
// >>chng 04 mar 24, from -2.810 to -2.607 add PAHs
// >>chng 04 apr 02, from -2.607 to -2.585, update co photo dissoc rate
// >>chng 04 apr 17, from -2.585 to -2.799, grn recom, better mole conv
// >>chng 04 jun 28, from -2.80 to -2.66, bug in part count for cons pres
// >>chng 04 jul 03, from -2.66 to -2.72, NA update H, CO chem net
// >>chng 06 mar 02, from -2.72 to -2.66, NA Fe0, Mg0, <-> Si+, S+, C+ ct,
// also energy barrier for H + (CH and CH2) reactions
// >>chng 06 may 30, NA refactor CO, fix H+CH=>C+H2 rate, from -2.66 to -2.59
// >>chng 06 aug 25, from -2.59 to -3.52 add mag field and stop col 22, 
// >>chng 04 sep 11, from -3.52 to -3.57, turn on CO net when H+/H<0.5, inc err
// >>chng 06 nov 24, from -3.57 to -3.54, split diff betweeen icl, vs05, inc er
// >>chng 06 nov 29, from -3.54 to -3.43 Badnell DR turned on
// >>chng 08 feb 20, from -3.43 to -3.56, grain charge transfer changes in ion_solver and iso*
// >>chng 08 apr 03, from -3.56 to -3.69 include 10830 in fine opacity array
// >>chng 08 jul 08, from -3.69 to -3.55 line radiation pressure
// >>chng 08 jul 23, from -3.55 to -3.63 chng Pesc line radiation pressure
// >>chng 09 oct 05, from -3.63 to -3.544, [NI] pumping changes
// >>chng 10 mar 02, from -3.544 to -3.385, fix botches after CO model switch
// >>chng 11 jul 12, from -3.385 to -3.587, enable fit to Badnell DR suppression
// >>chng 11 nov 16, from -3.587 to -5.514, adjust newmole monitors
// >>chng 11 dec 12, from -5.514 to -5.236, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from -5.236 to -2.587, H column density 1e22 -> 1e23
monitor line luminosity "CO  " 433.4m  -2.587 error 0.20
//
// >>chng 03 oct 03, from -2.565 to -2.529, integ co ion solvers
// >>chng 04 jan 23, from -2.529 to -2.558 rescale orion dust to gas
// >>chng 04 mar 04, from -2.558 to -2.623 improv conv, incr error
// >>chng 04 mar 24, from -2.623 to -2.415 add PAHs
// >>chng 04 apr 17, from -2.415 to -2.626, grn recom, better mole conv
// >>chng 04 may 13, from -2.626 to -2.605, dest albedo chng, incr err
// >>chng 04 may 29, from -2.605 to -2.646, NA update H chem
// >>chng 04 jun 28, from -2.65 to -2.51, bug in part count for cons pres
// >>chng 04 jul 03, from -2.51 to -2.56, NA update H, CO chem net
// >>chng 05 feb 18, from -2.56 to -2.53, incld outwrd cont in line pump
// >>chng 06 may 30, NA refactor CO, fix H+CH=>C+H2 rate, from -2.53 to -2.44
// >>chng 06 aug 25, from -2.44 to -3.42 add mag field and stop col 22, 
// >>chng 04 sep 11, from -3.42 to -3.46, turn on CO net when H+/H<0.5, inc err
// >>chng 06 nov 24, from -3.46 to -3.44, split diff betweeen icl, vs05, inc er
// >>chng 06 nov 29, from -3.44 to -3.33 Badnell DR turned on
// >>chng 08 feb 20, from -3.33 to -3.44, grain charge transfer changes in ion_solver and iso*
// >>chng 08 apr 03, from -3.44 to -3.55 include 10830 in fine opacity array
// >>chng 08 jul 08, from -3.55 to -3.44 line radiation pressure
// >>chng 08 jul 23, from -3.44 to -3.50 chng Pesc line radiation pressure
// >>chng 09 oct 06, from -3.5 to -3.422, [NI] pumping changes
// >>chng 10 mar 02, from -3.422 to -3.232, fix botches after CO model switch
// >>chng 11 jul 12, from -3.232 to -3.428, enable fit to Badnell DR suppression
// >>chng 11 nov 16, from -3.428 to -5.348, adjust newmole monitors
// >>chng 11 dec 12, from -5.348 to -5.134, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from -5.134 to -2.526, H column density 1e22 -> 1e23
monitor line luminosity "CO  " 371.5m  -2.526 error 0.20
//
// >>chng 03 oct 03, from -2.429 to -2.377, integ co ion solvers
// >>chng 04 jan 23, from -2.377 to -2.409 rescale orion dust to gas
// >>chng 04 mar 04, from -2.409 to -2.473 improv conv, incr error
// >>chng 04 mar 24, from -2.473 to -2.249 add PAHs
// >>chng 04 apr 17, from -2.249 to -2.490, grn recom, better mole conv
// >>chng 04 may 29, from -2.490 to -2.52, NA update H chem
// >>chng 04 jun 28, from -2.52 to -2.40, bug in part count for cons pres
// >>chng 06 may 30, NA refactor CO, fix H+CH=>C+H2 rate, from -2.40 to -2.33
// >>chng 06 aug 25, from -2.33 to -3.38 add mag field and stop col 22, 
// >>chng 06 nov 24, from -3.38 to -3.40, split diff betweeen icl, vs05, inc er
// >>chng 06 nov 29, from -3.40 to -3.31 Badnell DR turned on
// >>chng 07 dec 06, from -3.31 to -3.25, changes in h-like collisions
// >>chng 08 feb 20, from -3.25 to -3.38, grain charge transfer changes in ion_solver and iso*
// >>chng 08 mar 20, from -3.38 to -3.44 value drifted but did not throw assert
// >>chng 08 jul 23, from -3.44 to -3.44 chng Pesc line radiation pressure
// >>chng 09 oct 05, from -3.44 to -3.368, [NI] pumping changes
// >>chng 10 mar 02, from -3.368 to -3.148, fix botches after CO model switch
// >>chng 11 jul 12, from -3.148 to -3.337, enable fit to Badnell DR suppression
// >>chng 11 aug 05, from -3.337 to -3.281, add helike quadrupole transitions
// >>chng 11 nov 16, from -3.281 to -5.248, adjust newmole monitors
// >>chng 11 dec 12, from -5.248 to -5.103, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from -5.103 to -2.546, H column density 1e22 -> 1e23
monitor line luminosity "CO  " 325.1m  -2.546 error 0.15
//
// >>chng 03 oct 03, from -2.316 to -2.263, integ co ion solvers
// >>chng 04 jan 23, from -2.263 to -2.295 rescale orion dust to gas
// >>chng 04 feb 21, from -2.295 to -2.270, update chem, H- ct, grn recom
// >>chng 04 mar 04, from -2.270 to -2.365 improv conv, incr error
// >>chng 04 mar 24, from -2.365 to -2.110 add PAHs
// >>chng 04 apr 17, from -2.110 to -2.380, grn recom, better mole conv
// >>chng 04 may 13, from -2.380 to -2.330, dest albedo chng, incr err
// >>chng 04 may 29, from -2.33 to -2.41, NA update H chem
// >>chng 04 jun 28, from -2.41 to -2.29, bug in part count for cons pres
// >>chng 06 may 30, NA refactor CO, fix H+CH=>C+H2 rate, from -2.29 to -2.24
// >>chng 06 aug 25, from -2.24 to -3.40 add mag field and stop col 22, 
// >>chng 06 nov 24, from -3.40 to -3.42, split diff betweeen icl, vs05, inc er
// >>chng 06 nov 29, from -3.42 to -3.35 Badnell DR turned on
// >>chng 07 dec 06, from -3.35 to -3.28, changes in h-like collisions
// >>chng 08 feb 20, from -3.28 to -3.39, grain charge transfer changes in ion_solver and iso*
// >>chng 08 apr 03, from -3.39 to -3.48 include 10830 in fine opacity array
// >>chng 09 feb 23, from -3.48 to -3.422, remove otslin, otscon dampers (PR #84, r2758)
// >>chng 10 feb 14, from -3.422 to -3.338, update O atomic data
// >>chng 10 mar 02, from -3.338 to -3.124, fix botches after CO model switch
// >>chng 10 aug 12, from -3.124 to -3.188, fix botches related to r4125; star botches from recent continuum mesh changes
// >>chng 11 jul 12, from -3.188 to -3.306, enable fit to Badnell DR suppression
// >>chng 11 nov 16, from -3.306 to -5.201, adjust newmole monitors
// >>chng 11 dec 12, from -5.201 to -5.127, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from -5.127 to -2.757, H column density 1e22 -> 1e23
// >>chng 12 oct 22, from -2.757 to -2.834, zone tau in rt_line_one
monitor line luminosity "CO  " 289.0m  -2.834 error 0.15
//
// >>chng 03 oct 03, from -2.246 to -2.196, integ co ion solvers
// >>chng 04 jan 23, from -2.196 to -2.224 rescale orion dust to gas
// >>chng 04 feb 21, from -2.224 to -2.197, update chem, H- ct, grn recom
// >>chng 04 mar 04, from -2.197 to -2.315 improv conv, incr error
// >>chng 04 mar 24, from -2.315 to -2.002 add PAHs
// >>chng 04 apr 17, from -2.002 to -2.323, grn recom, better mole conv
// >>chng 04 may 13, from -2.323 to -2.251, dest albedo chng, incr err
// >>chng 04 may 29, from -2.25 to -2.35, NA update H chem
// >>chng 04 jun 28, from -2.35 to -2.21, bug in part count for cons pres
// >>chng 06 aug 25, from -2.21 to -3.48 add mag field and stop col 22, 
// >>chng 07 jan 15, from -3.48 to -3.41, small model H2 bugfix
// >>chng 08 mar 17, from -3.41 to -3.50, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 jul 23, from -3.50 to -3.50 chng Pesc line radiation pressure
// >>chng 09 oct 05, from -3.5 to -3.436, [NI] pumping changes
// >>chng 10 mar 02, from -3.436 to -3.153, fix botches after CO model switch
// >>chng 10 jun 17, from -3.153 to -3.22, fix botches due to recent dBase cooling changes
// >>chng 11 jul 12, from -3.22 to -3.328, enable fit to Badnell DR suppression
// >>chng 11 nov 16, from -3.328 to -5.196, adjust newmole monitors
// >>chng 11 dec 21, from -5.196 to -3.278, H column density 1e22 -> 1e23
// >>chng 12 apr 09, from -3.278 to -3.352, update CR backbround to Indriolo+07, ~8x larger
monitor line luminosity "CO  " 260.2m  -3.352 error 0.15
//
// >>chng 03 oct 03, from -2.245 to -2.201, integ co ion solvers
// >>chng 04 mar 04, from -2.201 to -2.349 improv conv, incr error
// >>chng 04 mar 24, from -2.349 to -1.930 add PAHs
// >>chng 04 apr 17, from -1.930 to -2.353, grn recom, better mole conv
// >>chng 04 may 13, from -2.353 to -2.239, dest albedo chng, incr err
// >>chng 04 may 29, from -2.24 to -2.38, NA update H chem
// >>chng 04 jun 28, from -2.38 to -2.19, bug in part count for cons pres
// >>chng 06 aug 25, from -2.19 to -3.61 add mag field and stop col 22, 
// >>chng 07 jan 15, from -3.61 to -3.53, small model H2 bugfix
// >>chng 08 mar 17, from -3.53 to -3.60, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 jul 08, from -3.60 to -3.76 line radiation pressure
// >>chng 08 jul 23, from -3.76 to -3.60 line radiation pressure bug
// >>chng 10 feb 14, from -3.6 to -3.528, update O atomic data
// >>chng 10 mar 02, from -3.528 to -3.231, fix botches after CO model switch
// >>chng 10 jun 17, from -3.231 to -3.307, fix botches due to recent dBase cooling changes
// >>chng 11 jul 12, from -3.307 to -3.397, enable fit to Badnell DR suppression
// >>chng 11 nov 16, from -3.397 to -5.222, adjust newmole monitors
// >>chng 11 dec 16, from -5.222 to -5.29, adjust botches due to recent iso recomb changes
// >>chng 11 dec 21, from -5.29 to -3.984, H column density 1e22 -> 1e23
// >>chng 12 apr 09, from -3.984 to -4.07, update CR backbround to Indriolo+07, ~8x larger
monitor line luminosity "CO  " 236.5m  -4.07 error 0.15
//
// >>chng 03 oct 03, from -2.310 to -2.285, integ co ion solvers
// >>chng 04 mar 04, from -2.285 to -2.508 improv conv, incr error
// >>chng 04 mar 24, from -2.508 to -1.904 add PAHs
// >>chng 04 apr 17, from -1.904 to -2.518, grn recom, better mole conv
// >>chng 04 may 13, from -2.518 to -2.312, dest albedo chng, incr err
// >>chng 04 may 29, from -2.31 to -2.68, NA update H chem
// >>chng 04 jun 28, from -2.68 to -2.27, bug in part count for cons pres
// >>chng 05 feb 18, from -2.27 to -2.20, incld outwrd cont in line pump
// >>chng 06 aug 25, from -2.20 to -3.79 add mag field and stop col 22, 
// >>chng 07 jan 15, from -3.79 to -3.71, small model H2 bugfix
// >>chng 08 apr 03, from -3.71 to -3.78 include 10830 in fine opacity array
// >>chng 08 jul 08, from -3.78 to -3.85 line radiation pressure
// >>chng 08 jul 23, from -3.85 to -3.76 line radiation pressure bug
// >>chng 10 feb 14, from -3.76 to -3.699, update O atomic data
// >>chng 10 mar 02, from -3.699 to -3.354, fix botches after CO model switch
// >>chng 10 jun 17, from -3.354 to -3.437, fix botches due to recent dBase cooling changes
// >>chng 11 jul 12, from -3.437 to -3.509, enable fit to Badnell DR suppression
// >>chng 11 nov 16, from -3.509 to -5.271, adjust newmole monitors
// >>chng 11 dec 12, from -5.271 to -5.4, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from -5.4 to -4.67, H column density 1e22 -> 1e23
// >>chng 12 apr 18, from -4.67 to -4.737, stout merge
monitor line luminosity "CO  " 216.9m  -4.737 error 0.15
// 
// >>chng 05 sep 30, add these 5
// >>chng 05 nov 22, from 11.19 to 10.99 stop half thick to avoid grn freeze out
// >>chng 06 aug 25, from 10.99 to 10.77 add mag field and stop col 22, 
// >>chng 08 feb 20, from 10.77 to 10.89, grain charge transfer changes in ion_solver and iso*
// >>chng 08 jun 25, from 10.89 to 10.77 (12.1) add source/sink to ioniz ladder
// >>chng 08 dec 13, from 10.77 to 10.7, merging newsolvers branch
// >>chng 09 oct 05, from 10.7 to 11.39, [NI] pumping changes
// >>chng 11 dec 12, from 11.39 to 11.34, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from 11.34 to 11.38 H column density 1e22 -> 1e23
monitor column density H2+  11.38 error 0.1
// 
// >>chng 05 nov 22, from 14.63 to 13.78 stop half thick to avoid grn freeze out
// >>chng 06 jan 04, from 13.78 to 13.82, NA C chem restore missed block reacts
// >>chng 06 apr 15, from 13.82 to 13.77 NA added N2H+ to the chemistry
// >>chng 06 aug 25, from 13.77 to 9.04 add mag field and stop col 22, 
// >>chng 07 jan 07, from 9.04 to 9.41, enable McCall et al. H3+ rate
// >>chng 08 jun 25, from 9.41 to 9.36 (12.1) add source/sink to ioniz ladder
// >>chng 08 jul 08, from 9.36 to 9.44 line radiation pressure
// >>chng 08 jul 23, from 9.44 to 9.36 chng Pesc line radiation pressure
// >>chng 09 feb 24, from 9.36 to 9.463, botches due to r2763
// >>chng 09 oct 05, from 9.463 to 9.758, [NI] pumping changes
// >>chng 11 nov 16, from 9.758 to 9.684, adjust newmole monitors
// >>chng 11 dec 12, from 9.684 to 9.755, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from 9.755 to 11.46, H column density 1e22 -> 1e23
// >>chng 12 apr 09, from 11.46 to 12.33, update CR backbround to Indriolo+07, ~8x larger
monitor column density H3+  12.33 error 0.1
// 
// >>chng 05 nov 22, from 25.47 to 24.99 stop half thick to avoid grn freeze out
// >>chng 06 aug 25, from 24.99 to 21.13 add mag field and stop col 22, 
// >>chng 11 nov 16, from 21.13 to 21.06, adjust newmole monitors
// >>chng 11 dec 12, from 21.06 to 21.22, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from 21.22 to 22.67, H column density 1e22 -> 1e23
monitor column density H2g  22.67 error 0.1
// 
// >>chng 05 dec 08, from 17.25 to 16.63, update to Elwert Solomon rate
// >>chng 06 may 30, NA refactor CO, fix H+CH=>C+H2 rate, from 16.63 to 16.83
// >>chng 06 aug 25, from 16.83 to 16.95 add mag field and stop col 22, 
// >>chng 11 nov 16, from 16.95 to 16.28, adjust newmole monitors
// >>chng 11 dec 16, from 16.28 to 16.24, adjust botches due to recent iso recomb changes
// >>chng 11 dec 21, from 16.24 to 16.2, H column density 1e22 -> 1e23
monitor column density H2*  16.2 error 0.1
// 
// >>chng 05 dec 08, add this temperature
monitor temperature hydrogen 2 8995
// 
// >>chng 06 may 30, NA refactor CO, fix H+CH=>C+H2 rate, from 667 to 635
// >>chng 06 aug 25, from 635 to 765 add mag field and stop col 22, 
// >>chng 08 mar 19, from 765 to 790, open error, this value varies
// over about 4% depending on compiler - asserted value is now near middle
// of range
// >>chng 08 dec 13, from 790 to 694.4, merging newsolvers branch
// >>chng 09 feb 08, error 0.1 -> 0.05
// >>chng 09 oct 05, from 694.4 to 1060, [NI] pumping changes
// >>chng 11 dec 12, from 1060 to 997.5, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from 997.5 to 1077, H column density 1e22 -> 1e23
monitor temperature hydrogen 1 1077
// 
// >>chng 06 aug 12, from 24.5 to 25.8 dr from tau IR continuum, inc err
// >>chng 06 aug 17, undo above zoning from grains, chng is error
// >>chng 06 aug 25, from 24.5 to 138 add mag field and stop col 22, 
// >>chng 11 dec 12, from 138 to 113.9, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from 113.9 to 30.24, H column density 1e22 -> 1e23
monitor temperature H2 30.24 error 0.1
// 
// >>chng 06 aug 25, from 11.25 to 11.30 add mag field and stop col 22, 
// >>chng 08 feb 20, from 11.30 to 11.50, grain charge transfer changes in ion_solver and iso*
// >>chng 08 jun 25, from 11.50 to 11.26 (12.1) add source/sink to ioniz ladder
// >>chng 08 jul 08, from 11.26 to 11.21 line radiation pressure
// >>chng 08 jul 23, from 11.21 to 11.26 chng Pesc line radiation pressure
// >>chng 09 feb 24, from 11.26 to 11.21, botches due to r2763
// >>chng 09 oct 05, from 11.21 to 11.28, [NI] pumping changes
// >>chng 11 nov 16, from 11.28 to 11.2, adjust newmole monitors
// >>chng 12 aug 06, from 11.2 to 11.12, bugfix - include Pelec_esc in several places
monitor column density HeH+ 11.12 error 0.1
//
// >>chng 04 mar 04, from < 1200 to <900 improv conv
// >>chng 04 mar 31, from < 900 to <950 alter dr logic-test all chng hv mole
// >>chng 04 apr 17, from 950 to 1050 grn recom, better mole conv
// >>chng 04 jun 28, from 1050 to 1200, bug in part count for cons pres
// >>chng 04 jul 17, from 1200 to 1400 NA update chem net
// >>chng 05 jul 04, from 1400 to 1600 turn cl on
// >>chng 05 nov 22, from 1600 to 1800 stop half thick to avoid grn freeze out
// >>chng 06 aug 12, from 1800 to 2200 (2082) dr from tau IR continuum, inc err
// >>chng 07 oct 29, from 2200 to 719, change from "<" to "="
// >>chng 08 dec 08, from 719 to 705, upper limit -> equals (r2515)
// >>chng 09 jan 08, from 705 to 713, minor changes in iter/zn, iterations, and zone that crept in over last month
// >>chng 09 jan 16, from 713 to 711, average of 10 runs
// >>chng 09 jun 01, from 711 to 719, drift up
// >>chng 11 nov 16, from 719 to 646, adjust newmole monitors
// >>chng 11 dec 12, from 646 to 692, bug-fix for updating grain opacity (r5759)
// >>chng 11 dec 21, from 692 to 1211 H column density 1e22 -> 1e23
monitor nzone 1211 
//
// >>chng 03 oct 30, from 11 to 15, many fewer zone, larger itr/zn, faster
// >>chng 03 nov 16, from 15 to 9, actual 7.44
// >>chng 04 jun 28, from 9 to 14, bug in part count for cons pres
// >>chng 04 aug 04, from 14 to 12, finer zones
// >>chng 04 sep 16, from 12 to 16, incre edentrue convergence
// >>chng 05 feb 06, from 16 to 18, actul 16.1
// >>chng 08 dec 08, from 18 to 12.1, upper limit -> equals (r2515)
// >>chng 08 dec 13, from 12.1 to 47.06, merging newsolvers branch
// >>chng 09 jan 08, from 47.06 to 46.19, minor changes in iter/zn, iterations, and zone that crept in over last month
// >>chng 09 jan 16, from 46.19 to 46.14 (error 0.04), average of 10 runs
// >>chng 11 dec 21, from 46.14 to 30.36, H column density 1e22 -> 1e23
monitor itrzn 30.36
//
// >>chng 08 jul 07, add radiation pressure
// >>chng 08 jul 08, set to 0.163
// >>chng 08 jul 23, from 0.163 to 0.0184 chng Pesc line radiation pressure
// >>chng 12 oct 09, from 0.0184 to 0.0175, fix CRD esc prob
// >>chng 13 jan 14, from 0.0175 to 0.0166, Pelec_esc improvements
monitor Pradmax 0.0166
//
c orion_hii_pdr.in
c class hii pdr  
c ========================================

This extends BFM from the H+ region into the PDR as in Abel et al 2005.
This is the correct way to do a PDR calcualtion.

>>refer	Orion	model	Baldwin, J., Ferland, G. J., 
>>refercon	Martin, P. G., Corbin, M., Cota, S., Peterson, 
>>refercon	B. M., & Slettebak, A. 1991, ApJ, 374, 580

>>refer	physics	HII/PDR	Abel, N.P., Ferland, G.J., Shaw, G. & 
>>refercon	van Hoof, P.A.M. 2005, ApJS, 161, 65
