Changeset 749

Timestamp:

12/26/06 11:40:59 (2 years ago)

Author:

gary

Message:

tsuite - fix formatting of some sims in slow tsuite suite

programs now compile as c++ - mostly TRUE -> true, etc

htm files in programs now have readme_ as start of name

Location:

trunk/tsuite

Files:

• auto/doc_tsuite.htm (modified) (128 diffs)
• auto/runsome_purify.dat (modified) (3 diffs)
• auto/runsome_purify.pl (modified) (1 diff)
• programs/collion/collion.cpp (modified) (2 diffs)
• programs/collion/collion.txt (modified) (43 diffs)
• programs/collion/readme_collion.htm (moved) (moved from trunk/tsuite/programs/collion/collion.htm) (1 diff)
• programs/comp4/comp4.cpp (modified) (2 diffs)
• programs/comp4/readme_comp4.htm (moved) (moved from trunk/tsuite/programs/comp4/comp4.htm)
• programs/hazy_coolingcurve/hazy_coolingcurve.cpp (modified) (1 diff)
• programs/hazy_coolingcurve/hazy_coolingcurve.txt (modified) (1 diff)
• programs/hazy_kmt/hazy_kmt.cpp (modified) (3 diffs)
• programs/hazy_kmt/hazy_kmt.txt (modified) (1 diff)
• programs/hazy_kmt/readme_hazy_kmt.htm (moved) (moved from trunk/tsuite/programs/hazy_kmt/hazy_kmt.htm)
• programs/hizlte/hizlte.cpp (modified) (1 diff)
• programs/hizlte/hizlte.txt (modified) (1 diff)
• programs/hizlte/readme_hizlte.htm (moved) (moved from trunk/tsuite/programs/hizlte/hizlte.htm)
• programs/mpi/mpi.cpp (modified) (1 diff)
• programs/mpi/readme_mpi.htm (moved) (moved from trunk/tsuite/programs/mpi/mpi.htm)
• programs/varyn/readme_varyn.htm (moved) (moved from trunk/tsuite/programs/varyn/varyn.htm)
• programs/varyn/varyn.cpp (modified) (1 diff)
• programs/varyn/varyn.txt (modified) (1 diff)
• programs/vary_nete/readme_vary_nete.htm (moved) (moved from trunk/tsuite/programs/vary_nete/vary_nete.htm)
• programs/vary_nete/vary_nete.cpp (modified) (6 diffs)
• programs/vary_nete/vary_nete.txt (modified) (1 diff)
• slow/doc_tsuite.htm (added)
• slow/doc_tsuite.pl (added)
• slow/feii_blr_n09_p18.in (modified) (1 diff)
• slow/feii_blr_n09_p18_Z20.in (modified) (1 diff)
• slow/feii_blr_n11_p20.in (modified) (1 diff)
• slow/feii_blr_n11_p20_Z20.in (modified) (1 diff)
• slow/feii_blr_n12_p19.in (modified) (1 diff)
• slow/feii_blr_n12_p19_Z20.in (modified) (1 diff)
• slow/feii_blr_n13_p22.in (modified) (1 diff)
• slow/feii_blr_n13_p22_Z20.in (modified) (1 diff)

trunk/tsuite/auto/doc_tsuite.htm

@@ -56 +56 @@
 check that these continua have the expected brightness.
 <p>
+</p><hr><kbd><h2>agn_lex00_u0.in  <br>stop column density 16
+<br>print line sort wavelength 
+<br>*constant temperature 720,000
+<br>no induced processes
+<br>print lines column linear 
+<br>normalise to "H  1" 1216
+<br>print line faint -1
+<br>print last iteration
+<br>hden 5
+<br>iterate
+<br>phi(h) 15.477 range 7.353 to 735.3
+<br>interpolate (-8 -3)
+<br>continue (-6 1.)
+<br>continue (-2, 14.5)
+<br>continue (0.4771 12.7) 
+<br>continue (0.8663 10.6)
+<br>continue (3.8663 7.6)
+<br>continue (6 1)
+<br>continue (9 -3)
+<br>element abundance helium  -1
+<br>element abundance carbon -3.432
+<br>element abundance nitrogen -3.959
+<br>element abundance oxygen -3.097
+<br>element abundance neon -3.959
+<br>element abundance magnesium -4.4318
+<br>element abundance silicon -4.456
+<br>element abundance sulphur -4.7959
+<br>element abundance argon -5.4318
+<br>element abundance iron -4.398
+<br>element lithium off
+<br>element beryllium off
+<br>element boron off
+<br>element fluorine off
+<br>element sodium off
+<br>element aluminium off
+<br>element phosphrous off
+<br>element chlorine off
+<br>element potasium off
+<br>element calcium off
+<br>element scandium off
+<br>element titanium off
+<br>element vanadium off
+<br>element chromium off
+<br>element manganese off
+<br>element cobalt off
+<br>element nickel off
+<br>element copper off
+<br>element zinc off
+<br>c
+<br>c agn_lex00_u0.in
+</kbd><p>
+<p>
+<p>
+</p><hr><kbd><h2>agn_lex00_u1.in  <br>stop column density 16
+<br>*constant temper 720,000K
+<br>no induced processes
+<br>print lines column linear 
+<br>normalise to "H  1" 1216
+<br>print line faint -1
+<br>print last iteration
+<br>hden 5
+<br>iterate
+<br>phi(h) 16.477 range 7.353 to 735.3
+<br>interpolate (-8 -3)
+<br>continue (-6 1.)
+<br>continue (-2, 14.5)
+<br>continue (0.4771 12.7) 
+<br>continue (0.8663 10.6)
+<br>continue (3.8663 7.6)
+<br>continue (6 1)
+<br>continue (9 -3)
+<br>element abundance helium  -1
+<br>element abundance carbon -3.432
+<br>element abundance nitrogen -3.959
+<br>element abundance oxygen -3.097
+<br>element abundance neon -3.959
+<br>element abundance magnesium -4.4318
+<br>element abundance silicon -4.456
+<br>element abundance sulphur -4.7959
+<br>element abundance argon -5.4318
+<br>element abundance iron -4.398
+<br>element lithium off
+<br>element beryllium off
+<br>element boron off
+<br>element fluorine off
+<br>element sodium off
+<br>element aluminium off
+<br>element phosphrous off
+<br>element chlorine off
+<br>element potasium off
+<br>element calcium off
+<br>element scandium off
+<br>element titanium off
+<br>element vanadium off
+<br>element chromium off
+<br>element manganese off
+<br>element cobalt off
+<br>element nickel off
+<br>element copper off
+<br>element zinc off
+<br>c
+<br>c agn_lex00_u1.in 
+</kbd><p>
+<p>
+</p><hr><kbd><h2>agn_lex00_um1.in  <br>stop column density 16
+<br>no induced processes
+<br>print lines column linear
+<br>print lines sort intensity 
+<br>normalise to "H  1" 1216
+<br>print line faint -2
+<br>print last iteration
+<br>hden 5
+<br>iterate
+<br>phi(h) 14.477 range 7.353 to 735.3
+<br>interpolate (-8 -3)
+<br>continue (-6 1.)
+<br>continue (-2, 14.5)
+<br>continue (0.4771 12.7) 
+<br>continue (0.8663 10.6)
+<br>continue (3.8663 7.6)
+<br>continue (6 1)
+<br>continue (9 -3)
+<br>element abundance helium  -1
+<br>element abundance carbon -3.432
+<br>element abundance nitrogen -3.959
+<br>element abundance oxygen -3.097
+<br>element abundance neon -3.959
+<br>element abundance magnesium -4.4318
+<br>element abundance silicon -4.456
+<br>element abundance sulphur -4.7959
+<br>element abundance argon -5.4318
+<br>element abundance iron -4.398
+<br>element lithium off
+<br>element beryllium off
+<br>element boron off
+<br>element fluorine off
+<br>element sodium off
+<br>element aluminium off
+<br>element phosphrous off
+<br>element chlorine off
+<br>element potasium off
+<br>element calcium off
+<br>element scandium off
+<br>element titanium off
+<br>element vanadium off
+<br>element chromium off
+<br>element manganese off
+<br>element cobalt off
+<br>element nickel off
+<br>element copper off
+<br>element zinc off
+<br>c 
+<br>c
+<br>c agn_lex00_um1.in 
+</kbd><p>
+<p>
 </p><hr><kbd><h2>agn_reflector.in  </kbd><I>model of Compton reflector
 </I></h2><kbd><br>title model of Compton reflector
@@ -237 +393 @@
 <br>stop zone 1
 <br>set dr 0
+<br>c
+<br>c commands controlling output    =========
+<br>punch heating "agn_S_u150.het"
+<br>punch cooling "agn_S_u150.col"
 <br>c
 <br>c
@@ -854 +1014 @@
 <br>punch dr "blr_n09_p18_Z20.dr"
 <br>punch convergence reason "blr_n09_p18_Z20.cvr"
+<br>punch convergence error "blr_n09_p18_Z20.cve"
 <br>c
 <br>c blr_n09_p18_Z20.in 
@@ -1501 +1662 @@
 <br>c
 <br>c commands controlling output    =========
-<br>normalize to "c  5" 40.27A 
 <br>c output files
 <br>punch overview "coll_heat_only.ovr"
+<br>punch heating "coll_heat_only.het"
+<br>punch cooling "coll_heat_only.col"
 <br>punch continuum "coll_heat_only.con"  units angstroms 
 <br>c
-</kbd><p>
-c coll_heat_only.in
-c class coronal 
-c ========================================
-<p>
+<br>c  
+<br>c coll_heat_only.in
+<br>c class coronal 
+<br>c ========================================
+</kbd><p>
 This test is an optically thin collisionally ionized gas with no
 photoionization at all.  Before c9601 the code would not do the ionization
@@ -1624 +1786 @@
 Test with only collisional ionization at a
 high gas high temperature. 
+</p><hr><kbd><h2>dynamics_orion_flow.in  <br>* stop zone 1000 1 1 
+<br>* no co molecules
+</kbd><I>Orion nebula blister with wind
+</I></h2><kbd><br>title Orion nebula blister with wind
+<br>c
+<br>c commands controlling continuum =========
+<br>blackbody 40,000
+<br>phi(h) 13.0
+<br>brems 6
+<br>phi(h) 10
+<br>c
+<br>c commands for density & abundances =========
+<br>hden 4
+<br>abundances hii region no grains
+<br>grains orion no qheat single 
+<br>c
+<br>c commands controlling geometry  =========
+<br>wind -7 km/s advection
+<br>c this is the account for continued matter beyond i front 
+<br>double optical depths
+<br>sphere 
+<br>stop AV 15 
+<br>stop temperature linear 5  
+<br>c
+<br>c other commands for details     =========
+<br>iterate 3
+<br>magnetic field -4 
+<br>cosmic rays background 
+<br>c this has no effect on dynamics, but is to desaturate the lines as per flow
+<br>turbulence 8 km/sec no pressure
+<br>c try to speed up simulation without hurting dynamics
+<br>init "fast.ini"
+<br>c use Bakes & Tielens heating
+<br>set PAH Bakes 
+<br>* init "c84.ini"
+<br>* failures 5
+<br>c
+<br>c commands controlling output    =========
+<br>c want to print line intensities are surface brightness, per arcsec^2
+<br>print lines surface brightness arcsec
+<br>print diffuse continuum 
+<br>print line faint -1.5
+<br>punch overview "dynamics_orion_flow.ovr" no hash
+<br>punch hydrogen 21 cm "dynamics_orion_flow.21cm" no hash
+<br>punch dr "dynamics_orion_flow.dr" no hash
+<br>punch continuum last "dynamics_orion_flow.con" units microns no hash
+<br>punch pressure "dynamics_orion_flow.pre" no hash
+<br>punch heating "dynamics_orion_flow.het" no hash
+<br>punch cooling "dynamics_orion_flow.col" no hash
+<br>punch molecules "dynamics_orion_flow.mol" no hash
+<br>punch wind "dynamics_orion_flow.wnd" last 
+<br>c
+<br>c dynamics_orion_flow.in
+<br>c class wind hii 
+<br>c ========================================
+<br>c 
+</kbd><p>
+This is a model similar in spirit to the blister geometry HII region
+model computed by Baldwin et al. (1991), but with a D-critical flow.
+<p>
+</p><hr><kbd><h2>dynamics_veryfast.in  </kbd><I>very fast dusty windy model
+</I></h2><kbd><br>title very fast dusty windy model
+<br>c
+<br>c commands controlling continuum =========
+<br>table agn
+<br>ionization parameter -2
+<br>c
+<br>c commands for density & abundances =========
+<br>hden 9
+<br>c this is a set of commands to speed up calc
+<br>init file "fast.ini"
+<br>c
+<br>c commands controlling geometry  =========
+<br>stop thickness 8.4
+<br>c this is a subsonic wind with advection 
+<br>wind -5 advection
+<br>c
+<br>c other commands for details     =========
+<br>c absolutely no convergence problems should occur 
+<br>failures 1
+<br>c this degrades the continuum resolution by factor of 3, to
+<br>c speed up calc 
+<br>set continuum resolution 3
+<br>iterate
+<br>c
+<br>c commands controlling output    =========
+<br>punch pressure last "dynamics_veryfast.pre"
+<br>punch total opacity last "dynamics_veryfast.opc"
+<br>punch ionizing continuum last "dynamics_veryfast.ion"
+<br>punch continuum last "dynamics_veryfast.con"
+<br>c
+<br>c 
+<br>c dynamics_veryfast.in 
+<br>c class dynamics  
+<br>c ========================================
+<br>c 
+</kbd><p>
+This is meant to be a very fast calculation to 
+use when running extensive debug-enabled runtimes.
+<p>
+<p>
+</p><hr><kbd><h2>dynamics_wind.in  </kbd><I>test of equations of motion in a wind
+</I></h2><kbd><br>title test of equations of motion in a wind
+<br>c 
+<br>c test of wind code
+<br>c radiative acceleration (e- only) is 9.54E-7 cm s^-2
+<br>c terminal velocity (e- only) is 7.6 km s^-1 
+<br>c
+<br>c commands controlling continuum =========
+<br>table agn
+<br>luminosiy (total) 45
+<br>c 
+<br>c commands for density & abundances =========
+<br>hden 4
+<br>c this only includes H and He
+<br>init file="hheonly.ini"
+<br>c
+<br>c commands controlling geometry  =========
+<br>radius (parsecs) 1
+<br>stop thickness (parsecs) -1
+<br>wind 0.1
+<br>c
+<br>c other commands for details     =========
+<br>c this will speed things up a bit 
+<br>no level2
+<br>no radiation pressure
+<br>constant temperature 8
+<br>c
+<br>c commands controlling output    =========
+<br>punch assets "dynamics_wind.asr"
+<br>punch pressure "dynamics_wind.pre"
+<br>punch dr "dynamics_wind.dr"
+<br>c
+<br>c dynamics_wind.in
+<br>c class limit dynamics 
+<br>c ========================================
+<br>c 
+</kbd><p>
+This tests the management of the radiative acceleration of an
+electron scattering wind, and the resulting velocity.  The parameters
+were chosen so that electron scattering is the dominant opacity source,
+so that the equations can be solved both numerically (in the example)
+and analytically (the expected solution given above).  In a realistic wind
+the gas would be more neutral and line driving would dominate.  The actual
+acceleration is slightly below that due to Thomson scattering alone because
+the gas opacity at high energies is slightly below Thomson. 
+<p>
+<br>Checks:
+<br>-   The radiative acceleration is correct (e- 9.543910-7 cm s-2).
+<br>-   The terminal velocity should be 7.57 km s-1.
+<br>-   Force multiplier near unity (no line driving since so highly ionized).
+<br>-   Thickness of cloud correct (R-Ro + dr/2 should be 3.086391017 cm).
+<p>
 </p><hr><kbd><h2>feii_hin.in  </kbd><I>test feii in high density limit
 </I></h2><kbd><br>title test feii in high density limit
@@ -1682 +1997 @@
 <br>c commands controlling output    =========
 <br>print last iteration
+<br>punch feii relative populations range 0 200 "feii_hirad.lv1" last
+<br>punch feii relative populations range 201 371 "feii_hirad.lv2" last
 <br>c
 <br>c
@@ -1715 +2032 @@
 <br>c
 <br>c commands controlling output    =========
+<br>punch feii continuum "feii_pump.con"
 <br>c
 <br>c feii_pump.in
@@ -1802 +2120 @@
 <br>punch feii continuum last "feii_ste.con"
 <br>punch feii column density last "feii_ste.col"
+<br>punch feii populations all last "feii_ste.pop"
 <br>print last iteration
 <br>c
@@ -1929 +2248 @@
 the globule model.
 <p>
+</p><hr><kbd><h2>func_grid.in  </kbd><I>test grid command and fits output
+</I></h2><kbd><br>title test grid command and fits output
+<br>c
+<br>c commands controlling continuum =========
+<br>c following gives range for variation of blackbody
+<br>c initial value ignored, will vary through above range
+<br>c value on blackbody command is ignored but needed to pass parser
+<br>blackbody 100,000K vary
+<br>grid, range from 4 to 6 with increment 1 
+<br>luminosity 39
+<br>c
+<br>c commands for density & abundances =========
+<br>hden 0 vary
+<br>grid, range from 2 to 3 with increment 1 
+<br>init "ism.ini"
+<br>c
+<br>c commands controlling geometry  =========
+<br>stop zone 1
+<br>radius 4 parsecs
+<br>c
+<br>c other commands for details     =========
+<br>constant temperature 4
+<br>c
+<br>c commands controlling output    =========
+<br>punch xspec atable "func_grid.fit"
+<br>punch overview "func_grid.ovr" no clobber 
+<br>punch linelist out to "func_grid.lin" read "LineList_HII.dat" no clobber
+<br>punch grid "func_grid.grd"
+<br>c
+<br>c
+<br>c func_grid.in
+<br>c class functionality 
+<br>c ========================================
+<br>c 
+</kbd><p>
+this model tests the grid command and fits output.
+<p>
+<p>
 </p><hr><kbd><h2>func_hotgas_coolstar.in  </kbd><I>test very soft continuum, very hot gas
 </I></h2><kbd><br>title test very soft continuum, very hot gas
@@ -1934 +2291 @@
 <br>c commands controlling continuum =========
 <br>c put in the cosmic background as the only continuum source 
-<br>blackbody 3, lte  
+<br>CMB
+<br>c gas has constant temperature of 1e6 K
 <br>constant temperature, t=6
 <br>c
@@ -1991 +2349 @@
 <br>c
 <br>c func_ion_increase.in
-<br>c class limit 
+<br>c class functionality 
 <br>c ========================================
 <br>c 
@@ -2059 +2417 @@
 <br>c
 <br>c commands controlling output    =========
+<br>normalize to "O  3" 5007 
 <br>punch dr "func_set_ion.dr" 
 <br>c
@@ -2130 +2489 @@
 This is a test of the highest temperature the code can do.
 <p>
-</p><hr><kbd><h2>func_t30.in  </kbd><I>test very soft continuum, very hot gas
-</I></h2><kbd><br>title test very soft continuum, very hot gas
+</p><hr><kbd><h2>func_t3.in  </kbd><I>test low temperature limit of code, 3K
+</I></h2><kbd><br>title test low temperature limit of code, 3K
 <br>c
 <br>c commands controlling continuum =========
@@ -2155 +2514 @@
 <br>c
 <br>c 
-<br>c func_t30.in
+<br>c func_t3.in
 <br>c class limit 
 <br>c ========================================
@@ -2161 +2520 @@
 </kbd><p>
 This is a test of the lowest temperature the code can do.
+It runs a constant temperature of 3K
 <p>
 </p><hr><kbd><h2>func_test.in  </kbd><I>this runs the standard, one command, test, which contains many asserts
@@ -2176 +2536 @@
 <br>c following will print physical constants used by the code
 <br>print constants
+<br>c print the past to the data directory
+<br>print path 
 <br>c following prints column densities 
 <br>print column density
@@ -2209 +2571 @@
 <br>punch continuum last "func_trans_punch.con"
 <br>c this continuum will be used by transread.in 
-<br>punch transmitted continuum last "func_trans_punch.trn"
+<br>punch last transmitted continuum "func_trans_punch.trn"
 <br>c
 <br>init file "trans.dat"
@@ -2301 +2663 @@
 <br>eden 5
 <br>init file "hheonly.ini"
+<br>element oxygen on 
 <br>grain abund 1 type 1
 <br>grain abund 1 type 2
@@ -2321 +2684 @@
 <br>c
 <br>c commands controlling output    =========
+<br>normalize to "O  1" 63.17m
 <br>c
 <br>c grains_lte.in
@@ -2367 +2731 @@
 <br>print diffuse continuum
 <br>print line faint 1
-<br>punch grains_qheat last "grains_qheat.qht"
+<br>punch grains temperature last "grains_qheat.qht"
+<br>punch continuum  "grains_qheat.con" units microns
 <br>c
 <br>c
@@ -2620 +2985 @@
 <br>c
 <br>c other commands for details     =========
+<br>c turn on the large H2 model
 <br>atom h2 
+<br>c turn down the number of convergence failures - there should be none
 <br>failures 3
 <br>c use leiden initialization file
@@ -2626 +2993 @@
 <br>c This command defines the grain temperature to be a constant 20 Kelvin
 <br>constant grain temperature 20 
-<br>c This sets the temperature to a constant 50 Kelvin
+<br>c This sets the gas kinetic temperature to a constant 50 Kelvin
 <br>constant temperature 50 linear
 <br>c
 <br>c commands controlling output    =========
+<br>c default normalization line is Hbeta, which is not produced by this gas
 <br>normalize to "C  2" 157.6m
 <br>print line faint -4 
 <br>punch overview "h2_pdr_leiden_f1.ovr"
+<br>punch line emissivity "h2_pdr_leiden_f1.ems" 
+<br>H2   2.121m
+<br>12CO 647.2m
+<br>C  2 157.6m
+<br>O  1 63.17m
+<br>Si 2 34.81m
+<br>end of lines
 <br>punch leiden lines "h2_pdr_leiden_f1.lin"
 <br>punch leiden "h2_pdr_leiden_f1.lei"
@@ -2792 +3167 @@
 <br>iterate
 <br>case b
-<br>cosmic rays
+<br>cosmic rays background
 <br>c
 <br>c commands controlling output    =========
@@ -3654 +4029 @@
 <br>c
 <br>c commands controlling output    =========
+<br>c this increases the number of wavelength digits in the line output
+<br>c to make it easier to identify highly excited Rydberg lines
+<br>set line precision 6
 <br>print he-like departure helium
-<br>normalise to "He 1" 4471
+<br>normalise to "He 1" 4471.47
 <br>print line faint 0.01
 <br>c
@@ -3666 +4044 @@
 This is one of the standard comparisons with the Benjamin et al.
 (1999), ApJ, 514, 307, paper.  The asserts give their answers.  
+<p>
+the line wavelength precision is increased to differentiate between
+various HeI lines
 <p>
 </p><hr><kbd><h2>he1n2t4_best.in  </kbd><I>test hei atom vs Benjamin et al. 99
@@ -3787 +4168 @@
 <br>c
 <br>c commands controlling output    =========
+<br>punch fits "he1n4t4.fit"
 </kbd><I>punch diffuse continuum "he1n4t4.dif" last no units microns
 </I></h2><kbd><br>punch diffuse continuum "he1n4t4.dif" last no title units microns
@@ -4266 +4648 @@
 <br>hden 7
 <br>set dr -12
+<br>init file "ism.ini"
+<br>c must turn on since ism.ini turned it off
+<br>element cobalt on
+<br>c then set abundance and ionization after turning it on
 <br>element cobalt abundance -2
 <br>element cobalt ionization -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 0 -5
-<br>init file "ism.ini"
+<br>c now turn lots of elements off to save time
 <br>element helium off
 <br>element nitrogen off
@@ -4279 +4665 @@
 <br>element argon off
 <br>element iron off
-<br>element cobalt on
 <br>laser 1300
 <br>ionization -1
@@ -4294 +4679 @@
 <br>hden 7
 <br>set dr -12
+<br>init file "ism.ini"
+<br>c must turn on since ism.ini turned it off
+<br>element copper on
+<br>c then set abundance and ionization after turning it on
 <br>element copper abundance -2
 <br>element copper ionization -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 0 -5
-<br>init file "ism.ini"
+<br>c now turn lots of elements off to save time
 <br>element helium off
 <br>element nitrogen off
@@ -4307 +4696 @@
 <br>element argon off
 <br>element iron off
-<br>element copper on
 <br>laser 1300
 <br>ionization -1
@@ -4514 +4902 @@
 <br>hden 7
 <br>set dr -12
+<br>init file "ism.ini"
+<br>c must turn on since ism.ini turned it off
+<br>element nickel on
+<br>c then set abundance and ionization after turning it on
 <br>element nickel abundance -2
 <br>element nickel ionization -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 0 -5
-<br>init file "ism.ini"
+<br>c now turn lots of elements off to save time
 <br>element helium off
 <br>element nitrogen off
@@ -4527 +4919 @@
 <br>element argon off
 <br>element iron off
-<br>element nickel on
 <br>laser 1300
 <br>ionization -1
@@ -4566 +4957 @@
 <br>c commands controlling output    =========
 <br>print he-like departure oxygen
-<br>print line faint -2
+<br>print line sort wavelength 
+<br>print line faint -5
 <br>normalize to "o  7" 21.60A
 <br>c
@@ -4634 +5026 @@
 <br>hden 7
 <br>set dr -12
+<br>init file "ism.ini"
+<br>c must turn on since ism.ini turned it off
+<br>element zinc on
+<br>c then set abundance and ionization after turning it on
 <br>element zinc abundance -2
 <br>element zinc ionization -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 -5 0 -5
-<br>init file "ism.ini"
+<br>c now turn lots of elements off to save time