t_opac Struct Reference

#include <opacity.h>

Collaboration diagram for t_opac:

Data Fields
realnum **	TauAbsGeo
realnum **	TauScatGeo
realnum **	TauTotalGeo
realnum *	TauAbsFace
realnum *	TauScatFace
double *	opacity_abs
double *	opacity_sct
double *	OldOpacSave
double *	albedo
double *	opacity_abs_savzon1
double *	opacity_sct_savzon1
double *	OpacStatic
double *	FreeFreeOpacity
double *	ExpZone
realnum *	E2TauAbsFace
realnum *	E2TauAbsTotal
realnum *	E2TauAbsOut
realnum *	TauAbsTotal
realnum *	ExpmTau
realnum *	tmn
bool	lgOpacStatic
bool	lgRedoStatic
double *	OpacStack
realnum	taumin
realnum	tlamin
bool	lgCaseB
bool	lgCaseB_HummerStorey
bool	lgCaseB_no_photo
bool	lgCaseB_no_pdest
realnum	telec
realnum	thmin
bool	lgOpacNeg
bool	lgScatON
bool	lgNegOpacIO
bool	lgCompileOpac
bool	lgUseFileOpac
long int	nOpacTot
long int	ipRayScat
long int	iopcom
long int	ippr
long int	ioppr
long int	ipBrems
long int	iphmra
long int	iphmop
long int	ih2pnt [2]
long int	ih2pof
long int	iophe1 [9]
long int	ioptri
long int	ipElement [LIMELM][LIMELM][7][3]
long int	in1 [3]
long int	ipo3exc [3]
long int	ipo3exc3 [3]
long int	ipo1exc [3]
long int	iopo2d
long int	ipmgex
long int	ipOpMgEx
long int	ica2ex [2]
long int	ica2op
long int	ipH2_photo_thresh
long int	ipH2_photo_opac_offset
long int	ipCKshell
realnum	otsmin
realnum	stimax [2]

Detailed Description

Definition at line 75 of file opacity.h.

Field Documentation

double* t_opac::albedo

albedo is local gas albedo

Definition at line 104 of file opacity.h.

Referenced by DumpLine(), emergent_line(), IterRestart(), OpacityAddTotal(), rfield_opac_malloc(), rfield_opac_zero(), RT_DestProb(), RTesc_lya_1side(), and SaveDo().

realnum* t_opac::E2TauAbsFace

saves value of E2(tau), exponential integral, where tau is optical depth to illuminated face

Definition at line 124 of file opacity.h.

Referenced by emergent_line(), IterRestart(), IterStart(), lines_continuum(), rfield_opac_malloc(), rfield_opac_zero(), RT_continuum(), and RT_recom_effic().

realnum* t_opac::E2TauAbsOut

Definition at line 127 of file opacity.h.

Referenced by emergent_line(), IterRestart(), rfield_opac_malloc(), rfield_opac_zero(), and RT_continuum().

realnum* t_opac::E2TauAbsTotal

same things but for outward direction. Only defined on second and later iteratios

Definition at line 126 of file opacity.h.

Referenced by emergent_line(), rfield_opac_malloc(), rfield_opac_zero(), and RT_tau_reset().

realnum* t_opac::ExpmTau

exp(-tau) to illuminated face

Definition at line 132 of file opacity.h.

Referenced by cdSPEC2(), IterRestart(), IterStart(), rfield_opac_malloc(), rfield_opac_zero(), RT_continuum(), RT_diffuse(), and RT_OTS().

double* t_opac::ExpZone

exp(-dTau) for the current zone, evaluated in radinc

Definition at line 120 of file opacity.h.

Referenced by IterRestart(), IterStart(), outline_base_bin(), rfield_opac_malloc(), rfield_opac_zero(), RT_continuum(), and RT_recom_effic().

double* t_opac::FreeFreeOpacity

density/temp factors needed for free-free opacity

Definition at line 117 of file opacity.h.

Referenced by CoolEvaluate(), OpacityAddTotal(), rfield_opac_malloc(), rfield_opac_zero(), and save_opacity().

long int t_opac::ica2ex[2]

calcium excited states

Definition at line 210 of file opacity.h.

Referenced by ContCreatePointers(), OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::ica2op

Definition at line 210 of file opacity.h.

Referenced by OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::ih2pnt[2]

ih2pnt lower, upper limits to bound, ih2pof, opacity offset

Definition at line 210 of file opacity.h.

Referenced by ContCreatePointers(), hmole_step(), OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::ih2pof

Definition at line 210 of file opacity.h.

Referenced by hmole_step(), OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::in1[3]

in1 is [NI] excited state

Definition at line 210 of file opacity.h.

Referenced by ContCreatePointers(), IonNitro(), OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::iopcom

iopcom compton scatterin, total recoil

Definition at line 210 of file opacity.h.

Referenced by highen(), OpacityAddTotal(), OpacityCreateAll(), and PresTotCurrent().

long int t_opac::iophe1[9]

iophe1 points for photo from singlet levels

Definition at line 210 of file opacity.h.

Referenced by OpacityCreateAll(), and save_opacity().

long int t_opac::iopo2d

photo to excited O+ levels

Definition at line 210 of file opacity.h.

Referenced by IonOxyge(), and OpacityCreateAll().

long int t_opac::ioppr

Definition at line 210 of file opacity.h.

Referenced by highen(), OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::ioptri

lowest levels of he triplets

Definition at line 210 of file opacity.h.

long int t_opac::ipBrems

ipBrems opacity offset pointer for brems (free-free)

Definition at line 210 of file opacity.h.

Referenced by OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::ipCKshell

index of carbon k-shell threshold in energy array

Definition at line 294 of file opacity.h.

Referenced by ContCreatePointers(), ContSetIntensity(), OpacityCreateAll(), and PrtContinuum().

long int t_opac::ipElement[LIMELM][LIMELM][7][3]

	 * ipElement[nelem][ioniz][shell][purpose],
	 * array index on fortran scale, to energy range are set in ipShells, called by ContCreatePointers
	 * pointers to opacity stack offset defined in OpacityCreate1Element where
	 * opacities are set
	 *
	 * first dim 
	 * [0] is pointer to low-energy threshold in energy array
	 * [1] is highest energy for shell, set by LimitSh
	 * [2] is opacity offset within large opacity stack
	 *
	 * if shell does not exist, set upper limit
	 * to less than lower limit so this never looped upon
	 * these are used as flags by LimitSh to check whether
	 * this is a real shell - 
	 * following code is in ipShells for this non-existant case 
	 * OpacPoint.ipElement[nelem][ion][nshell][0] = 2;
	 * OpacPoint.ipElement[nelem][ion][nshell][1] = 1;
	 * all routines must protect against this since opacities for these
	 * undefined shells is also undefined
	 *
	 * second dim is shell, 0 for k shell, up to 6,
	 * (Shell.chShell[ns] give the label for each of these shells)
	 * These quantities are only defined for the number of shells that exist. 
	 * The number of shells is given by Heavy.nsShells[nelem][ion]
	 *
	 * third dimension is ion stage, 0 for atom
	 * 
	 * last dim is atomic number of element, 0 for H
	 *

Definition at line 210 of file opacity.h.

Referenced by ContCreatePointers(), GammaPrtRate(), GammaPrtShells(), ion_photo(), IonIron(), ipShells(), LimitSh(), OpacityAdd1Element(), OpacityCreate1Element(), OpacityCreateAll(), OpacityValenceRescale(), RT_diffuse(), SanityCheckBegin(), save_opacity(), and SaveDo().

long int t_opac::ipH2_photo_opac_offset

Definition at line 291 of file opacity.h.

Referenced by hmole_step(), OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::ipH2_photo_thresh

threshold for photoionization of H2 - 15.4 eV

Definition at line 291 of file opacity.h.

Referenced by ContCreatePointers(), hmole_step(), OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::iphmop

iphmop H- bound free opacity

Definition at line 210 of file opacity.h.

Referenced by hmole_reactions(), OpacityAddTotal(), OpacityCreateAll(), and RT_diffuse().

long int t_opac::iphmra

iphmra ratio of h minus to neut h ff brems opacity

Definition at line 210 of file opacity.h.

Referenced by CoolEvaluate(), OpacityAddTotal(), OpacityCreateAll(), and RT_diffuse().

long int t_opac::ipmgex

photoionization from upper lev of Mg II 2798

Definition at line 210 of file opacity.h.

Referenced by ContCreatePointers(), IonMagne(), OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::ipo1exc[3]

Definition at line 210 of file opacity.h.

Referenced by ContCreatePointers(), IonOxyge(), OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::ipo3exc[3]

pointers to oxy excited states

Definition at line 210 of file opacity.h.

Referenced by ContCreatePointers(), ContSetIntensity(), IonOxyge(), OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::ipo3exc3[3]

Definition at line 210 of file opacity.h.

Referenced by ContCreatePointers(), IonOxyge(), OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::ipOpMgEx

Definition at line 210 of file opacity.h.

Referenced by IonMagne(), OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::ippr

ippr is threshold for pair production, ioppr is opacity offset

Definition at line 210 of file opacity.h.

Referenced by ContCreatePointers(), highen(), OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::ipRayScat

NBNBNBthis must exactly parallel the read/write statements in OpacityCreateAll ipRayScat opacity pointer for Rayleigh scattering

Definition at line 210 of file opacity.h.

Referenced by OpacityAddTotal(), and OpacityCreateAll().

bool t_opac::lgCaseB

this flag says we are doing case b

Definition at line 161 of file opacity.h.

Referenced by InitDefaultsPreparse(), iso_cascade(), iso_radiative_recomb(), IterStart(), lgConserveEnergy(), ParseCaseB(), ParseCommands(), RT_tau_init(), RT_tau_reset(), and SanityCheckBegin().

bool t_opac::lgCaseB_HummerStorey

this flag turns off n=2 collisional excitations for comparison with their paper, keyword HUMMER on caseb command

Definition at line 165 of file opacity.h.

Referenced by HydroCSInterp(), InitDefaultsPreparse(), iso_collide(), and ParseCaseB().

bool t_opac::lgCaseB_no_pdest

another case b option, turn off background opacities, no Pdest

Definition at line 172 of file opacity.h.

Referenced by InitDefaultsPreparse(), ParseCaseB(), and RT_line_all().

bool t_opac::lgCaseB_no_photo

this flag turns off excited state photoionization, keyword NO PHOTO on caseb command

Definition at line 169 of file opacity.h.

Referenced by InitDefaultsPreparse(), iso_photo(), and ParseCaseB().

bool t_opac::lgCompileOpac

variable dealing with the option to generate and use an ancillary file containing a stored opacity set. These are all initialized in zero1 lgCompileOpac flag saying to compile opacity only

Definition at line 192 of file opacity.h.

Referenced by OpacityCreateAll(), ParseCompile(), and zero().

bool t_opac::lgNegOpacIO

IO unit to save negative opacities

Definition at line 186 of file opacity.h.

Referenced by ParseSet(), pnegopc(), and zero().

bool t_opac::lgOpacNeg

flag set if negative opacities occured

Definition at line 179 of file opacity.h.

Referenced by PrtComment(), RT_continuum(), and zero().

bool t_opac::lgOpacStatic

true by default, and set false with the *no static opacities* command. When false always update all opacities

Definition at line 140 of file opacity.h.

Referenced by ConvBase(), ParseDont(), and zero().

bool t_opac::lgRedoStatic

this flag is set true in OpacityZero when the OpacStatic array is zeroed, and is false if the array has been left alone. all later opacities must be reevaluated when this flag is true this is the stack used to hold opacities - entered one time when code is initialized, in routine OpacityCreateAll

Definition at line 149 of file opacity.h.

Referenced by ConvBase(), GammaPrtShells(), ion_photo(), iso_photo(), iso_update_rates(), OpacityAdd1Subshell(), OpacityAdd1SubshellInduc(), OpacityZero(), SanityCheckBegin(), and save_opacity().

bool t_opac::lgScatON

flag saying whether (true) or not (false) scattering opacity is enabled

Definition at line 183 of file opacity.h.

Referenced by ParseDont(), ParseSphere(), RT_continuum(), and zero().

bool t_opac::lgUseFileOpac

set false if no file opacity command entered, option to ignore opacity file

Definition at line 196 of file opacity.h.

Referenced by ParseDont(), and zero().

long int t_opac::nOpacTot

total number of opacity cells used in opacity stack in OpacityCreateAll used as a counter to remember where next opacity goes

Definition at line 201 of file opacity.h.

Referenced by OpacityCreate1Element(), OpacityCreateAll(), OpacityCreatePowerLaw(), and OpacityCreateReilMan().

double* t_opac::OldOpacSave

save previous opacity

Definition at line 101 of file opacity.h.

Referenced by IterRestart(), OpacityZero(), OpacityZeroOld(), rfield_opac_malloc(), and rfield_opac_zero().

double* t_opac::opacity_abs

local absorption opacity, cm-1

Definition at line 95 of file opacity.h.

Referenced by ContRate(), ConvBase(), CoolCalc(), DumpLine(), GrainRateDr(), IterEnd(), IterRestart(), IterStart(), OpacityAdd1Subshell(), OpacityAdd1SubshellInduc(), OpacityAddTotal(), OpacityCreate1Element(), OpacityCreateAll(), OpacityZero(), OpacityZeroOld(), pltcon(), pltopc(), pnegopc(), PresTotCurrent(), radius_first(), radius_next(), rfield_opac_malloc(), rfield_opac_zero(), RT_continuum(), RT_DestProb(), RT_diffuse(), RT_line_all(), RT_line_pumping(), RT_OTS(), RT_OTS_AddCont(), RT_OTS_AddLine(), RT_OTS_PrtRate(), RT_OTS_Update(), RT_recom_effic(), RTesc_lya(), SanityCheckBegin(), Save_Line_RT(), save_opacity(), SaveDo(), and ZoneStart().

double* t_opac::opacity_abs_savzon1

initial opacities from zone 1, used to reset opacity at restart and the saved value of local opacity

Definition at line 108 of file opacity.h.

Referenced by IterRestart(), IterStart(), pltopc(), and rfield_opac_malloc().

double* t_opac::opacity_sct

local scattering opacity, cm-1

Definition at line 98 of file opacity.h.

Referenced by IterRestart(), IterStart(), OpacityAddTotal(), OpacityZero(), pltopc(), PresTotCurrent(), rfield_opac_malloc(), rfield_opac_zero(), RT_continuum(), Save_Line_RT(), save_opacity(), SaveDo(), and SaveSpecial().

double* t_opac::opacity_sct_savzon1

opacity_sct_savzon1 save local opacity at start of calculations

Definition at line 110 of file opacity.h.

Referenced by IterRestart(), IterStart(), pltopc(), and rfield_opac_malloc().

double* t_opac::OpacStack

[NOPSV]

Definition at line 151 of file opacity.h.

Referenced by CoolEvaluate(), csphot(), GammaBn(), GammaK(), GammaPrt(), highen(), opacity_more_memory(), OpacityAdd1Subshell(), OpacityAdd1SubshellInduc(), OpacityAddTotal(), OpacityCreate1Element(), OpacityCreateAll(), OpacityCreatePowerLaw(), OpacityCreateReilMan(), OpacityValenceRescale(), PresTotCurrent(), RT_diffuse(), SanityCheckBegin(), and save_opacity().

double* t_opac::OpacStatic

these static opacities are only evaluated one time per zone

Definition at line 114 of file opacity.h.

Referenced by OpacityAdd1Subshell(), OpacityAdd1SubshellInduc(), OpacityAddTotal(), OpacityZero(), rfield_opac_malloc(), rfield_opac_zero(), SanityCheckBegin(), and save_opacity().

realnum t_opac::otsmin

smallest ots rate, set to 0 in scalar

Definition at line 297 of file opacity.h.

Referenced by ParseDont(), RT_recom_effic(), and zero().

realnum t_opac::stimax[2]

max correction for stim emission in continuum opacities at Lyman and Balmer edges

Definition at line 300 of file opacity.h.

Referenced by IterRestart(), OpacityAddTotal(), PrtComment(), and zero().

realnum* t_opac::TauAbsFace

these are the integrated absorption and scattering optical depths to the illuminated face of the cloud

Definition at line 91 of file opacity.h.

Referenced by hmole_reactions(), IterRestart(), IterStart(), radius_increment(), rfield_opac_malloc(), RT_continuum(), RT_tau_reset(), and SaveDo().

realnum** t_opac::TauAbsGeo

[0] is optical depth for current position, [1] is total optical depth from previous iteration TauAbsGeo, TauScatGeo, absorption, scattering, optical depths, at current position, includes through far side of slab

Definition at line 82 of file opacity.h.

Referenced by iso_radiative_recomb(), iter_end_check(), PrtAllTau(), PrtComment(), PrtFinal(), PrtZone(), radius_next(), rfield_opac_malloc(), rfield_opac_zero(), RT_continuum(), RT_recom_effic(), RT_tau_init(), RT_tau_reset(), SaveDo(), state_get_put(), and tauff().