#include <opacity.h>

Inheritance diagram for t_opac:

Collaboration diagram for t_opac:

Public Member Functions
const char *	chName () const

void	zero ()

void	comment (t_warnings &)

Public Member Functions inherited from module
	module ()

virtual	~module ()

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 *	eeFreeFreeOpacity

double	eeFreeFreeTemp

vector< double, allocator_avx < 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	ih2pnt_ex [2]

long int	ih2pof_ex

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]

double	o3exc

double	o3exc3

long	iopo2d

long	ipmgex

long	ipOpMgEx

long	ica2ex [2]

long	ica2op

long int	ipCKshell

realnum	otsmin

realnum	stimax [2]

Detailed Description

Definition at line 78 of file opacity.h.

Member Function Documentation

const char* t_opac::chName ( ) const

inlinevirtual

Implements module.

Definition at line 80 of file opacity.h.

void t_opac::comment ( t_warnings & )

inlinevirtual

Implements module.

Definition at line 85 of file opacity.h.

void t_opac::zero ( void )

virtual

Todo:: 2 file opacities are disabled for now - reinstate this when arrays settle down

Implements module.

Definition at line 7 of file opacity.cpp.

References DEBUG_ENTRY, lgCompileOpac, lgNegOpacIO, lgOpacNeg, lgOpacStatic, lgScatON, lgUseFileOpac, otsmin, and stimax.

Field Documentation

double* t_opac::albedo

albedo is local gas albedo

Definition at line 113 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 137 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 140 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 139 of file opacity.h.

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

double* t_opac::eeFreeFreeOpacity

Definition at line 127 of file opacity.h.

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

double t_opac::eeFreeFreeTemp

temperature used for the last ee brems evaluation

Definition at line 130 of file opacity.h.

Referenced by InitDefaultsPreparse(), and OpacityAddTotal().

realnum* t_opac::ExpmTau

exp(-tau) to illuminated face

Definition at line 145 of file opacity.h.

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

vector< double, allocator_avx<double> > t_opac::ExpZone

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

Definition at line 133 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 126 of file opacity.h.

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

long t_opac::ica2ex[2]

calcium excited states

Definition at line 300 of file opacity.h.

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

long t_opac::ica2op

Definition at line 300 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 223 of file opacity.h.

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

long int t_opac::ih2pnt_ex[2]

Definition at line 223 of file opacity.h.

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

long int t_opac::ih2pof

Definition at line 223 of file opacity.h.

Referenced by OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::ih2pof_ex

Definition at line 223 of file opacity.h.

Referenced by OpacityAddTotal(), and OpacityCreateAll().

long int t_opac::in1[3]

in1 is [NI] excited state

Definition at line 223 of file opacity.h.

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

long int t_opac::iopcom

iopcom compton scatterin, total recoil

Definition at line 223 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 223 of file opacity.h.

Referenced by OpacityCreateAll(), and save_opacity().

long t_opac::iopo2d

photo to excited O+ levels

Definition at line 300 of file opacity.h.

Referenced by ion_photo(), and OpacityCreateAll().

long int t_opac::ioppr

Definition at line 223 of file opacity.h.

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

long int t_opac::ioptri

lowest levels of he triplets

Definition at line 223 of file opacity.h.

long int t_opac::ipBrems

ipBrems opacity offset pointer for brems (free-free)

Definition at line 223 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 311 of file opacity.h.

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

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 223 of file opacity.h.

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

long int t_opac::iphmop

iphmop H- bound free opacity

Definition at line 223 of file opacity.h.

Referenced by mole_h_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 223 of file opacity.h.

Referenced by t_gaunt::brems_opac(), OpacityCreateAll(), and t_gaunt::p_setup_brems().

long t_opac::ipmgex

photoionization from upper lev of Mg II 2798

Definition at line 300 of file opacity.h.

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

long int t_opac::ipo1exc[3]

Definition at line 223 of file opacity.h.

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

long int t_opac::ipo3exc[3]

pointers to oxy excited states

Definition at line 223 of file opacity.h.

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

long int t_opac::ipo3exc3[3]

Definition at line 223 of file opacity.h.

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

long t_opac::ipOpMgEx

Definition at line 300 of file opacity.h.

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

long int t_opac::ippr

ippr is threshold for pair production, ioppr is opacity offset

Definition at line 223 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 OpacityCreateAllipRayScat opacity pointer for Rayleigh scattering

Definition at line 223 of file opacity.h.

Referenced by OpacityAddTotal(), and OpacityCreateAll().

bool t_opac::lgCaseB

this flag says we are doing case b

Definition at line 174 of file opacity.h.

Referenced by InitDefaultsPreparse(), iso_cascade(), iso_level(), iso_radiative_recomb(), IterStart(), lgConserveEnergy(), ParseCaseB(), ParseCommands(), PrtComment(), 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 178 of file opacity.h.

Referenced by GetHlikeCollisionStrength(), InitDefaultsPreparse(), iso_collide(), iso_get_collision_strength(), iso_get_collision_strength_collapsed_to_collapsed_fast(), and ParseCaseB().

bool t_opac::lgCaseB_no_pdest

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

Definition at line 185 of file opacity.h.

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

bool t_opac::lgCaseB_no_photo

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

Definition at line 182 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 205 of file opacity.h.

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

bool t_opac::lgNegOpacIO

IO unit to save negative opacities

Definition at line 199 of file opacity.h.

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

bool t_opac::lgOpacNeg

flag set if negative opacities occured

Definition at line 192 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 153 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 160 of file opacity.h.

Referenced by ConvBase(), GammaPrtShells(), ion_photo(), iso_photo(), 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 196 of file opacity.h.

Referenced by t_rfield::getCoarseTransCoef(), 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 209 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 214 of file opacity.h.

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

double t_opac::o3exc

threshold energies corresponding to ipo3exc[0] and ipo3exc3[0]

Definition at line 296 of file opacity.h.

Referenced by ContCreatePointers(), and t_mesh::ValidateEdges().

double t_opac::o3exc3

Definition at line 297 of file opacity.h.

Referenced by ContCreatePointers(), and t_mesh::ValidateEdges().

double* t_opac::OldOpacSave

save previous opacity

Definition at line 110 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 104 of file opacity.h.

Referenced by ContRate(), 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_escape(), RT_line_electron_scatter(), 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 117 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 107 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 119 of file opacity.h.

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

double* t_opac::OpacStack

double* t_opac::OpacStatic

these static opacities are only evaluated one time per zone

Definition at line 123 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 314 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 317 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 100 of file opacity.h.

Referenced by ContCreateMesh(), IterRestart(), IterStart(), mole_h_reactions(), 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 91 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().