Example 2: Photoionized Grid

In this example a grid of photoionization models is calculated with varying ionization parameter, column density, and element abundances. The ionizing spectrum is a power law with index -1.

cfrac,r,a,0.,0.,1.,"covering fraction soft maximum"
cfractyp,i,h,0,0,2,"covering fraction variation type"
cfracint,i,a,1,0,1,"covering fraction interpolation type"
cfracsof,r,a,0.,0.,1.,"covering fraction soft minimum"
cfracnst,i,a,1,1,20,"covering fraction number of steps"
temperature,r,h,1.,0.,1.E4,"temperature soft maximum (/10**4K)"
temperaturetyp,i,h,0,0,2,"temperature variation type"
temperatureint,i,a,1,0,1,"temperature interpolation type"
temperaturesof,r,a,0.,0.,1.,"temperature soft minimum"
temperaturenst,i,a,1,1,20,"temperature number of steps"
pressure,r,h,0.03,0.,1.,"pressure soft maximum (dyne/cm**2)"
pressuretyp,i,h,0,0,2,"pressure variation type"
pressureint,i,a,1,0,1,"pressure interpolation type"
pressuresof,r,a,0.,0.,1.,"pressure soft minimum"
pressurenst,i,a,1,1,20,"pressure number of steps"
density,r,a,1.E+10,0.,1.E18,"density soft maximum (cm**-3)"
densitytyp,i,h,0,0,2,"density variation type"
densityint,i,a,1,0,1,"density interpolation type"
densitysof,r,a,0.,0.,1.,"density soft minimum"
densitynst,i,a,1,1,20,"density number of steps"
trad,r,a,-1.,,,"radiation temperature or alpha soft maximum?"
tradtyp,i,h,0,0,2,"radiation temperature variation type"
tradint,i,a,1,0,1,"radiation temperature interpolation type"
tradsof,r,a,0.,0.,1.,"radiation temperature soft minimum"
tradnst,i,a,1,1,20,"radiation temperature number of steps"
rlrad38,r,a,1.e+6,0.,1.E10,"luminosity soft maximum (/10**38 erg/s)"
rlrad38typ,i,h,0,0,2,"luminosity variation type"
rlrad38int,i,a,1,0,1,"luminosity interpolation type"
rlrad38sof,r,a,0.,0.,1.,"luminosity soft minimum"
rlrad38nst,i,a,1,1,20,"luminosity number of steps"
column,r,a,1.E23,0.,1.E25,"column density soft maximum (cm**-2)"
columntyp,i,h,2,0,2,"column density variation type"
columnint,i,a,1,0,1,"column density interpolation type"
columnsof,r,a,1.E20,1.,1.E25,"column density soft minimum"
columnnst,i,a,7,1,20,"column density number of steps"
rlogxi,r,a,4.5,-10.,+10.,"log(ionization parameter) soft maximum (erg cm/s)"
rlogxityp,i,h,2,0,2,"log(ionization parameter) variation type"
rlogxiint,i,a,0,0,1,"log(ionization parameter) interpolation type"
rlogxisof,r,a,-1.,-10.0,+10.0,"log(ionization parameter) soft minimum"
rlogxinst,i,a,20,1,20,"log(ionization parameter) number of steps"
habund,r,h,1.,0.,100.,"hydrogen abundance soft maximum"
habundtyp,i,h,0,0,2,"hydrogen abundance variation type"
habundint,i,a,1,0,1,"hydrogen abundance interpolation type"
habundsof,r,a,0.,0.,1.,"hydrogen abundance soft minimum"
habundnst,i,a,1,1,20,"hydrogen abundance number of steps"
heabund,r,h,1.,0.,100.,"helium abundance soft maximum"
heabundtyp,i,h,0,0,2,"helium abundance variation type"
heabundint,i,a,1,0,1,"helium abundance interpolation type"
heabundsof,r,a,0.,0.,1.,"helium abundance soft minimum"
heabundnst,i,a,1,1,20,"helium abundance number of steps"
cabund,r,h,1.,0.,100.,"carbon abundance soft maximum"
cabundtyp,i,h,1,0,2,"carbon abundance variation type"
cabundint,i,a,1,0,1,"carbon abundance interpolation type"
cabundsof,r,a,0.,0.,1.,"carbon abundance soft minimum"
cabundnst,i,a,1,1,20,"carbon abundance number of steps"
nabund,r,h,1.,0.,100.,"nitrogen abundance soft maximum"
nabundtyp,i,h,1,0,2,"nitrogen abundance variation type"
nabundint,i,a,1,0,1,"nitrogen abundance interpolation type"
nabundsof,r,a,0.,0.,1.,"nitrogen abundance soft minimum"
nabundnst,i,a,1,1,20,"nitrogen abundance number of steps"
oabund,r,h,1.,0.,100.,"oxygen abundance soft maximum"
oabundtyp,i,h,1,0,2,"oxygen abundance variation type"
oabundint,i,a,1,0,1,"oxygen abundance interpolation type"
oabundsof,r,a,0.,0.,1.,"oxygen abundance soft minimum"
oabundnst,i,a,1,1,20,"oxygen abundance number of steps"
neabund,r,h,1.,0.,100.,"neon abundance soft maximum"
neabundtyp,i,h,1,0,2,"neon abundance variation type"
neabundint,i,a,1,0,1,"neon abundance interpolation type"
neabundsof,r,a,0.,0.,1.,"neon abundance soft minimum"
neabundnst,i,a,1,1,20,"neon abundance number of steps"
mgabund,r,h,1.,0.,100.,"magnesium abundance soft maximum"
mgabundtyp,i,h,1,0,2,"magnesium abundance variation type"
mgabundint,i,a,1,0,1,"magnesium abundance interpolation type"
mgabundsof,r,a,0.,0.,1.,"magnesium abundance soft minimum"
mgabundnst,i,a,1,1,20,"magnesium abundance number of steps"
siabund,r,h,1.,0.,100.,"silicon abundance soft maximum"
siabundtyp,i,h,1,0,2,"silicon abundance variation type"
siabundint,i,a,1,0,1,"silicon abundance interpolation type"
siabundsof,r,a,0.,0.,1.,"silicon abundance soft minimum"
siabundnst,i,a,1,1,20,"silicon abundance number of steps"
sabund,r,h,1.,0.,100.,"sulfur abundance soft maximum"
sabundtyp,i,h,1,0,2,"sulfur abundance variation type"
sabundint,i,a,1,0,1,"sulfur abundance interpolation type"
sabundsof,r,a,0.,0.,1.,"sulfur abundance soft minimum"
sabundnst,i,a,1,1,20,"sulfur abundance number of steps"
arabund,r,h,1.,0.,100.,"argon abundance soft maximum"
arabundtyp,i,h,1,0,2,"argon abundance variation type"
arabundint,i,a,1,0,1,"argon abundance interpolation type"
arabundsof,r,a,0.,0.,1.,"argon abundance soft minimum"
arabundnst,i,a,1,1,20,"argon abundance number of steps"
caabund,r,h,1.,0.,100.,"calcium abundance soft maximum"
caabundtyp,i,h,1,0,2,"calcium abundance variation type"
caabundint,i,a,1,0,1,"calcium abundance interpolation type"
caabundsof,r,a,0.,0.,1.,"calcium abundance soft minimum"
caabundnst,i,a,1,1,20,"calcium abundance number of steps"
feabund,r,h,1.,0.,100.,"iron abundance soft maximum"
feabundtyp,i,h,1,0,2,"iron abundance variation type"
feabundint,i,a,1,0,1,"iron abundance interpolation type"
feabundsof,r,a,0.,0.,1.,"iron abundance soft minimum"
feabundnst,i,a,1,1,20,"iron abundance number of steps"
niabund,r,h,0.,0.,100.,"nickel abundance soft maximum"
niabundtyp,i,h,0,0,2,"nickel abundance variation type"
niabundint,i,a,1,0,1,"nickel abundance interpolation type"
niabundsof,r,a,0.,0.,1.,"nickel abundance soft minimum"
niabundnst,i,a,1,1,20,"nickel abundance number of steps"
spectrum,s,a,"pow",,,"spectrum type?"
spectrum_file,s,a,"spct.dat",,,"spectrum file?"
spectun,i,a,0,0,1,"spectrum units? (0=energy, 1=photons)"
redshift,i,h,1,0,1,"Is redshift a parameter? (0=no, 1=yes)"
nsteps,i,h,3,1,1000,"number of steps"
niter,i,h,99,,,"number of iterations"
lwrite,i,h,0,0,1,"write switch (1=yes, 0=no)"
lprint,i,h,0,0,1,"print switch (1=yes, 0=no)"
lstep,i,h,0,,,"step size choice switch"
npass,i,h,1,1,10000,"number of passes"
lcpres,i,h,0,0,1,"constant pressure switch (1=yes, 0=no)"
emult,r,h,0.5,1.e-6,1.e+6,"Courant multiplier"
taumax,r,h,5.,1.,10000.,"tau max for courant step"
xeemin,r,h,0.1,1.e-6,0.5,"minimum electron fraction"
critf,r,h,1.e-14,1.e-24,0.1,"critical ion abundance"
vturbi,r,h,1.,0.,30000.,"turbulent velocity (km/s)"
modelname,s,a,"photoionized grid",,,"model name"
loopcontrol,i,h,0,0,30000,"loop control (0=standalone)"
elow,r,h,1.0E+2,0.,5.11E+5,"energy band low end (eV)"
ehigh,r,h,2.0E+4,0.,5.11E+5,"energy band high end (eV)"
mode,s,h,"ql",,,"mode"