title Orion nebula blister with wind
c
c commands controlling continuum =========
blackbody 40000
phi(h) 13.0
brems 6
phi(h) 10
c
c commands for density & abundances =========
hden 4
abundances hii region no grains
grains orion no qheat single 
c
c commands controlling geometry  =========
wind -7 km/s advection
c this accounts for matter beyond end of the model
double optical depths
sphere 
// >>chng 06 jul 05, stop on AV rather than temperature */
// >>chng 08 jul 01, from 15 to 5 
stop AV 5 
stop temperature linear 5  
c
c other commands for details     =========
iterate 4
magnetic field -4 
cosmic rays background 
c this has no effect on dynamics, but is to desaturate the lines as per flow
turbulence 8 km/sec no pressure
c try to speed up simulation without hurting dynamics
init "fast.ini"
c use Bakes & Tielens heating
set PAH Bakes 
* init "c84.ini"
* failures 5
c
c commands controlling output    =========
c want to print line intensities are surface brightness, per arcsec^2
print lines surface brightness arcsec
print line faint -1.5
save performance "dynamics_orion_flow.per"
save overview "dynamics_orion_flow.ovr" no hash
save hydrogen 21 cm "dynamics_orion_flow.21cm" no hash
save dr "dynamics_orion_flow.dr" no hash
save continuum last "dynamics_orion_flow.con" units microns no hash
save monitors last "dynamics_orion_flow.asr" no hash
save pressure "dynamics_orion_flow.pre" no hash
save heating "dynamics_orion_flow.het" no hash
save cooling "dynamics_orion_flow.col" no hash
save molecules "dynamics_orion_flow.mol" no hash
save wind "dynamics_orion_flow.wnd" last 
c
c commands giving the monitors    =========
// >>chng 06 mar 11, from -11.241 to -11.175, chng small H from 4 to 5 lev
// >>chng 07 jun 25, from -11.175 to -11.328, old h-like topoff added indirect to n=4, new only to n=6, huge effect on 4861.
// >>chng 07 oct 26, from -11.328 to -11.350, very close to botch after recent changes
monitor line "totl" 4861 intensity -11.350
//
monitor line "CaBo" 4861 1.0
//
// >>chng 05 oct 01, from 0.103 to 0.082, RP update He levels
// >>chng 07 dec 29, from 0.082 to 0.857, ddrift up
// >>chng 08 apr 30, from 0.0857 to 0.0908 H, He topoff change
// >>chng 12 nov 01, from 0.0908 to 0.095, Chianti experimental levels, revised Badnell DR suppression
monitor line "he 1" 5876 0.095
// 
monitor line "Ca B" 5876 0.127
// 
monitor line "He 1" 3889 0.0801
monitor line "Ca B" 3889 0.105
// 
monitor line "He 1" 7065 0.0273
monitor line "Ca B" 7065 0.0331
//
// >>chng 08 apr 30, from 0.913 to 0.842 H, He topoff change
monitor line "totl" 10830 0.842
//
// >>chng 06 mar 11, from 0.241 to 0.263, chng small H from 4 to 5 lev
// >>chng 06 nov 29, from 0.263 to 0.221, Badnell DR turned on
// >>chng 08 apr 30, from 0.221 to 0.250 H, He topoff change
monitor line "totl" 2326 0.250 error 0.1
//
// >>chng 05 oct 01, from 0.200 to 0.223, RP update He levels
// >>chng 07 oct 22, from 0.223 to 0.246, resolve l-levels of h-like sequence
// >>chng 09 feb 23, from 0.246 to 0.2746, remove otslin, otscon dampers (PR #84, r2758)
monitor line "totl" 1909 0.2746 error 0.1
//
// >>chng 05 oct 01, from 0.505 to 0.480, RP update He levels
// >>chng 05 dec 20, from 0.480 to 0.516, RP change default helium levels
// >>chng 06 nov 29, from 0.516 to 0.449, Badnell DR turned on
// >>chng 08 apr 30, from 0.449 to 0.497 H, He topoff change
// >>chng 10 mar 05, from 0.497 to 0.5436, changes caused by r3848 atomic data update
monitor line "n  2" 6584 0.5436 error 0.1
//
// >>chng 05 dec 20, from 0.621 to 0.655, RP change default helium levels
// >>chng 06 nov 29, from 0.655 to 0.610, Badnell DR turned on
// >>chng 08 apr 30, from 0.610 to 0.668 H, He topoff change
monitor line "totl" 3727 0.668 error 0.1
//
// >>chng 05 oct 01, from 3.707 to 3.871, RP update He levels
// >>chng 07 dec 06, from 3.871 to 4.133, changes in h-like collisions
monitor line "o  3" 5007 4.133
//
// >>chng 06 jul 06, from 0.159 to 0.175, update Ne cs, CO chem now
// includes dynamics but no cross talk between CO and H2 in H2 part
// >>chng 07 dec 14, from 0.175 to 0.167 RP update is seq
// >>chng 08 nov 04, from 0.167 to 0.1746, line RT bug fixes (r2429-35)
// >>chng 11 jul 12, from 0.1746 to 0.1629, enable fit to Badnell DR suppression
// >>chng 13 feb 02, from 0.1629 to 0.1768, BUGFIX, largely undo r6703 
monitor line "ne 2" 12.81m 0.1768
//
// >>chng 05 dec 08, from 0.342 to 0.309 update to Elwert Solomon rate
// >>chng 06 jul 06, from 0.309 to 0.484, sim goes much deeper, CO chem now
// includes dynamics but no cross talk between CO and H2 in H2 part
// >>chng 06 nov 29, from 0.484 to 0.561, Badnell DR turned on, no co in advec
// >>chng 07 may 08, from 0.561 to 0.543, dyn logic changed to set stop 
// thickness after relaxing 2 iter, increased to 4 interation in this
// sim to get more dynamics logic tested
// >>chng 08 jul 01, from 0.561 to 0.0161, Av from 15 to 5
monitor HHeicf 0.0161 error 0.02
//
// some dynamic stuff
// >>chng 12 oct 22, from -0.13 to -0.1362, zone tau in rt_line_one
monitor velocity -0.1362
//
// assert dynamics convergence error < xx 
//
// >>chng 05 mar 28, from 80 to 50, impr mol converg
// >>chng 08 jul 01, from 50 to 65 (58) Av from 15 to 5
// >>chng 08 dec 08, from 65 to 57.33, upper limit -> equals (r2515)
// >>chng 08 dec 13, from 57.33 to 100.8, merging newsolvers branch
// >>chng 09 jan 16, from 100.8 to 100.1, average of 10 runs
// >>chng 09 apr 19, from 100.1 to 101.5 refactor p esc and dest with cooling
// >>chng 10 jun 17, from 101.5 to 91.29, fix botches due to recent dBase cooling changes
// >>chng 11 aug 06, from 91.29 to 100.8, add default min dr rel to depth
// >>chng 11 nov 16, from 100.8 to 75.31, adjust newmole monitors
// >>chng 12 mar 07, from 75.31 to 58.32, renorm to r5961 results
// >>chng 12 aug 16, from 58.32 to 55.3, further fix of Rawlings H2 rates
monitor itrzn 55.3 error 0.05
//
// >>chng 06 jul 06, from 900 to 1000 sim goes much deeper, CO chem now
// includes dynamics but no cross talk between CO and H2 in H2 part
// >>chng 08 dec 08, from 1000 to 578, upper limit -> equals (r2515)
// >>chng 09 jan 16, from 578 to 582, average of 10 runs
// >>chng 10 jun 17, from 582 to 658, fix botches due to recent dBase cooling changes
// >>chng 11 aug 05, from 658 to 688, fix mistake in r5333
// >>chng 11 aug 06, from 688 to 408, add default min dr rel to depth
monitor nzone 408
//
c dynamics_orion_flow.in
c class dynamics hii pdr
c ========================================
c 

This is a model similar in spirit to the blister geometry H+ region
model computed by Baldwin et al. (1991), but with a D-critical flow.
Many physical processes have been disabled to make this simulationn faster.
Grain physics is not done so the gas temperature is incorrect.  The
main purpose is to do a quick test of the dynamical flow with grain
opacities included.  The slow directory contains a full simulation
of a flow like Orion.