#include "cddefines.h"
#include "cddrive.h"
#include "thermal.h"
#include "physconst.h"
#include "doppvel.h"
#include "taulines.h"
#include "dense.h"
#include "rfield.h"
#include "radius.h"
#include "lines_service.h"
#include "ipoint.h"
#include "thirdparty.h"
#include "hydrogenic.h"
#include "lines.h"
#include "rt.h"
#include "trace.h"
#include "save.h"
#include "phycon.h"
#include "atomfeii.h"
#include "iso.h"
#include "pressure.h"
#include "parser.h"

Include dependency graph for atom_feii.cpp:

Macros
#define	NPRADDAT 159

#define	USE_OLD true

#define	AMAT(I_, J_) ((amat+(I_)FeII.nFeIILevel_local+(J_)))

#define	NLEVDEP 11

#define	NLEVPOP 11

Functions
STATIC void	FeIIContCreate (double xLamLow, double xLamHigh, long int ncell)

STATIC int	FeIIBandsCreate (const char chFile[])

STATIC void	FeIICollRatesBoltzmann (void)

STATIC void	FeIILyaPump (void)

void	FeII_Colden (const char *chLabel)

void	FeIICreate (void)

void	FeIILevelPops (void)

void	FeIIPrint (void)

double	FeIISumBand (realnum wl1, realnum wl2, double *SumBandInward)

void	FeII_RT_TauInc (void)

void	FeII_RT_tau_reset (void)

void	FeIIPoint (void)

void	FeIIAccel (double *fe2drive)

void	FeII_RT_Make (void)

void	FeIIAddLines (void)

void	FeIIPunchLevels (FILE *ioPUN)

void	FeIIPunchColden (FILE *ioPUN)

void	FeIIPunchOpticalDepth (FILE *ioPUN)

void	FeIISaveLines (FILE *ioPUN)

void	FeII_LineZero (void)

void	FeIIIntenZero (void)

void	FeIIPunData (FILE *ioPUN, bool lgDoAll)

void	FeIIPunDepart (FILE *ioPUN, bool lgDoAll)

void	FeIIPun1Depart (FILE *ioPUN, long int nPUN)

void	FeIIReset (void)

void	FeIIZero (void)

void	FeIIPunPop (FILE *ioPUN, bool lgPunchRange, long int ipRangeLo, long int ipRangeHi, bool lgPunchDensity)

void	AssertFeIIDep (double pred, double BigError, double *StdDev)

void	FeII_OTS (void)

void	FeII_RT_Out (void)

void	ParseAtomFeII (Parser &p)

void	PunFeII (FILE *io)

void	FeIIPunchLineStuff (FILE *io, realnum xLimit, long index)

double	FeIIRadPress (void)

Variables
realnum	CS2SMALL = (realnum)1e-5

int **	ncs1

realnum **	FeII_Bands

realnum **	FeII_Cont

long int	nFeIIBands

long int	nFeIIConBins

static long int *	nnPradDat

static realnum ***	sPradDat

static double **	Fe2SavN

static double **	Fe2A

static double **	Fe2LPump

static double **	Fe2CPump

realnum *	Fe2Energies

static realnum **	Fe2Coll

static double *	Fe2DepCoef

static double *	Fe2LevelPop

static double *	Fe2ColDen

static double *	FeIIBoltzmann

static double	EnerLyaProf1

static double	EnerLyaProf4

static double	PhotOccNumLyaCenter

static double *	yVector

static double **	xMatrix

static double *	amat

STATIC bool	lgFeIIEverCalled =false

Macro Definition Documentation

◆ AMAT

#define AMAT	(	I_,
		J_
	)	((amat+(I_)FeII.nFeIILevel_local+(J_)))

◆ NLEVDEP

#define NLEVDEP 11

◆ NLEVPOP

#define NLEVPOP 11

◆ NPRADDAT

#define NPRADDAT 159

Definition at line 100 of file atom_feii.cpp.

◆ USE_OLD

#define USE_OLD true

Function Documentation

◆ AssertFeIIDep()

void AssertFeIIDep	(	double *	pred,
		double *	BigError,
		double *	StdDev
	)

called by assert feii depart coef command

Parameters

*pred
*BigError
*StdDev

Definition at line 2452 of file atom_feii.cpp.

References ASSERT, DEBUG_ENTRY, Fe2DepCoef, FeII, lgFeIIEverCalled, t_FeII::lgSimulate, MAX2, t_FeII::nFeIILevel_local, and POW2.

Referenced by lgCheckMonitors().

◆ FeII_Colden()

void FeII_Colden ( const char * chLabel )

FeII_Colden maintain H2 column densities within X

Parameters

*chLabel

Definition at line 177 of file atom_feii.cpp.

References cdEXIT, DEBUG_ENTRY, t_radius::drad_x_fillfac, EXIT_FAILURE, Fe2ColDen, Fe2LevelPop, FeII, ioQQQ, t_FeII::nFeIILevel_local, t_FeII::nFeIILevel_malloc, and radius.

Referenced by molcol().

◆ FeII_LineZero()

void FeII_LineZero ( void )

initialize optical depth arrays, called by TauOut

Definition at line 1825 of file atom_feii.cpp.

References DEBUG_ENTRY, Fe2LevN, FeII, ipFe2LevN, and t_FeII::nFeIILevel_malloc.

Referenced by RT_tau_init().

◆ FeII_OTS()

void FeII_OTS ( void )

do OTS and outward parts of FeII lines, if large atom is turned on

Definition at line 2504 of file atom_feii.cpp.

References ASSERT, EmissionProxy::Aul(), DEBUG_ENTRY, TransitionProxy::Emis(), Fe2LevN, FeII, TransitionProxy::Hi(), TransitionProxy::ipCont(), ipFe2LevN, t_FeII::nFeIILevel_local, EmissionProxy::ots(), EmissionProxy::Pdest(), and RT_OTS_AddLine().

Referenced by RT_OTS().

Here is the call graph for this function:

◆ FeII_RT_Make()

void FeII_RT_Make ( void )

do line RT for FeII model

Definition at line 1538 of file atom_feii.cpp.

References t_dense::AtomicWeight, DEBUG_ENTRY, dense, Fe2LevN, FeII, GetDopplerWidth(), ioQQQ, TransitionProxy::ipCont(), ipFe2LevN, ipIRON, t_trace::lgTrace, t_FeII::nFeIILevel_local, t_FeII::nFeIILevel_malloc, RT_line_one(), and trace.

Referenced by Fe2_cooling().

Here is the call graph for this function:

◆ FeII_RT_Out()

void FeII_RT_Out ( void )

do outward rates for FeII, called by RT_diffuse

Definition at line 2542 of file atom_feii.cpp.

References DEBUG_ENTRY, dense, Fe2LevN, FeII, TransitionProxy::ipCont(), ipFe2LevN, ipIRON, t_FeII::nFeIILevel_local, TransitionProxy::outline_resonance(), and t_dense::xIonDense.

Referenced by RT_diffuse().

Here is the call graph for this function:

◆ FeII_RT_tau_reset()

void FeII_RT_tau_reset ( void )

FeII_RT_tau_reset reset optical depths for large FeII atom, called by update after each iteration

Definition at line 1425 of file atom_feii.cpp.

References DEBUG_ENTRY, Fe2LevN, FeII, ipFe2LevN, t_FeII::nFeIILevel_malloc, and RT_line_one_tau_reset().

Referenced by RT_tau_reset().

Here is the call graph for this function:

◆ FeII_RT_TauInc()

void FeII_RT_TauInc ( void )

FeII_RT_TauInc called once per zone in RT_tau_inc to increment large FeII atom line optical depths

Definition at line 1396 of file atom_feii.cpp.

References t_dense::AtomicWeight, DEBUG_ENTRY, dense, Fe2LevN, FeII, GetDopplerWidth(), TransitionProxy::ipCont(), ipFe2LevN, ipIRON, t_FeII::nFeIILevel_local, t_FeII::nFeIILevel_malloc, and RT_line_one_tauinc().

Referenced by RT_tau_inc().

Here is the call graph for this function:

◆ FeIIAccel()

void FeIIAccel ( double * fe2drive )

called by rt_line_driving to compute radiative acceleration due to FeII lines

Parameters

*fe2drive

Definition at line 1509 of file atom_feii.cpp.

References DEBUG_ENTRY, TransitionProxy::Emis(), TransitionProxy::EnergyErg(), Fe2LevN, FeII, ipFe2LevN, t_FeII::nFeIILevel_local, EmissionProxy::PopOpc(), and EmissionProxy::pump().

Referenced by RT_line_driving().

Here is the call graph for this function:

◆ FeIIAddLines()

void FeIIAddLines ( void )

called by LineSet4, this adds feii line intensities together

Definition at line 1574 of file atom_feii.cpp.

References DEBUG_ENTRY, t_radius::dVeffAper, TransitionList::Emis(), Fe2LevN, Fe2SavN, FeII, ioQQQ, t_LineSave::ipass, ipFe2LevN, LineSave, t_FeII::nFeIILevel_local, t_FeII::nFeIILevel_malloc, nzone, and radius.

Referenced by lines_lv1_k_zn().

Here is the call graph for this function:

◆ FeIIBandsCreate()

STATIC int FeIIBandsCreate ( const char chFile[] )

Definition at line 2293 of file atom_feii.cpp.

References DEBUG_ENTRY, FeII_Bands, FFmtRead(), FILENAME_PATH_LENGTH_2, ioQQQ, lgCalled, t_trace::lgTrace, MALLOC, nFeIIBands, open_data(), read_whole_line(), and trace.

Referenced by FeIICreate().

Here is the call graph for this function:

◆ FeIICollRatesBoltzmann()

STATIC void FeIICollRatesBoltzmann ( void )

Definition at line 1072 of file atom_feii.cpp.

References ASSERT, EmissionProxy::Aul(), cdEXIT, CollisionProxy::col_str(), col_str, TransitionProxy::Coll(), COLL_CONST, CS2SMALL, DEBUG_ENTRY, dsexp(), TransitionProxy::Emis(), TransitionProxy::EnergyWN(), EXIT_FAILURE, Fe2Coll, Fe2LevN, FeII, FeIIBoltzmann, fp_equal(), g, TransitionProxy::Hi(), ioQQQ, ipFe2LevN, TransitionProxy::Lo(), MAX2, MIN2, ncs1, NFE2LEVN, t_FeII::nFeIILevel_local, nnPradDat, NPRADDAT, phycon, POW2, POW3, sPradDat, t_phycon::sqrte, t_phycon::te, and t_phycon::te_wn.

Referenced by FeIILevelPops().

Here is the call graph for this function:

◆ FeIIContCreate()

STATIC void FeIIContCreate	(	double	xLamLow,
		double	xLamHigh,
		long int	ncell
	)

Definition at line 2574 of file atom_feii.cpp.

References DEBUG_ENTRY, FeII_Cont, lgCalled, MALLOC, and nFeIIConBins.

Referenced by FeIICreate().

◆ FeIICreate()

void FeIICreate ( void )

reads in feii data from disk, creates space for main arrays

Definition at line 219 of file atom_feii.cpp.

References TransitionProxy::AddLine2Stack(), multi_arr< T, d, ALLOC, lgBC >::alloc(), AllTransitions, amat, ASSERT, EmissionProxy::Aul(), cdEXIT, CS2SMALL, DEBUG_ENTRY, TransitionProxy::Emis(), TransitionProxy::EnergyWN(), EXIT_FAILURE, Fe2A, Fe2ColDen, Fe2Coll, Fe2CPump, Fe2DepCoef, Fe2Energies, Fe2LevelPop, Fe2LevN, Fe2LPump, Fe2SavN, t_FeII::feconwlHi, t_FeII::feconwlLo, FeII, t_FeII::FeIIAul, FeIIBandsCreate(), FeIIBoltzmann, t_FeII::FeIIColl, FeIIContCreate(), t_FeII::FeIINRGs, t_FeII::FeIISTWT, FILENAME_PATH_LENGTH_2, fp_equal(), EmissionProxy::FracInwd(), g, EmissionProxy::gf(), TransitionProxy::Hi(), ioQQQ, ipFe2LevN, t_FeII::ipRedisFcnResonance, t_FeII::ipRedisFcnSubordinate, EmissionProxy::iRedisFun(), lgFeIIMalloc, t_trace::lgTrace, TransitionProxy::Lo(), MALLOC, max(), MAX2, MIN2, ncs1, t_FeII::nfe2con, NFE2LEVN, t_FeII::nFeIILevel_local, t_FeII::nFeIILevel_malloc, nLine, nnPradDat, NPRADDAT, open_data(), EmissionProxy::ots(), EmissionProxy::phots(), POW2, read_whole_line(), RefIndex(), multi_arr< T, d, ALLOC, lgBC >::reserve(), qList::resize(), TransitionList::resize(), save, TransitionProxy::setHi(), TransitionProxy::setLo(), multi_arr< T, d, ALLOC, lgBC >::size(), sPradDat, TransitionList::states(), trace, TRANS_PROB_CONST, USE_OLD, TransitionProxy::WLAng(), xMatrix, yVector, and TransitionProxy::Zero().

Referenced by atmdat_readin().

Here is the call graph for this function:

◆ FeIIIntenZero()

void FeIIIntenZero ( void )

FeIIIntenZero zero out intensity of FeII atom

Definition at line 1851 of file atom_feii.cpp.

References TransitionProxy::Coll(), EmissionProxy::ColOvTot(), CollisionProxy::cool(), DEBUG_ENTRY, TransitionList::Emis(), TransitionProxy::Emis(), Fe2LevelPop, Fe2LevN, FeII, t_FeII::for7, CollisionProxy::heat(), TransitionProxy::Hi(), ipFe2LevN, ipT191, ipTFe16, ipTFe26, ipTFe34, ipTFe35, ipTFe46, ipTFe56, ipTr48, ipTuv3, TransitionProxy::Lo(), t_FeII::nFeIILevel_malloc, EmissionProxy::ots(), EmissionProxy::phots(), EmissionProxy::PopOpc(), TauLines, and EmissionProxy::xIntensity().

Referenced by Fe2_cooling().

Here is the call graph for this function:

◆ FeIILevelPops()

void FeIILevelPops ( void )

Definition at line 718 of file atom_feii.cpp.

References amat, AMAT, EmissionProxy::Aul(), cdEXIT, EmissionProxy::ColOvTot(), t_FeII::ddT_Fe2_large_cool, DEBUG_ENTRY, dense, t_hydro::dstfe2lya, t_dense::eden, TransitionProxy::Emis(), TransitionProxy::EnergyErg(), EXIT_FAILURE, t_FeII::Fe2_large_cool, t_FeII::Fe2_large_heat, Fe2A, Fe2Coll, Fe2CPump, Fe2DepCoef, Fe2LevelPop, Fe2LevN, Fe2LPump, FeII, FeIIBoltzmann, FeIICollRatesBoltzmann(), FeIILyaPump(), g, getrf_wrapper(), getrs_wrapper(), t_thermal::halfte, TransitionProxy::Hi(), hydro, ioQQQ, ipFe2LevN, ipH1s, ipH2p, ipH_LIKE, ipHYDROGEN, ipIRON, iso_sp, lgFeIIEverCalled, t_FeII::lgFeIILargeOn, t_FeII::lgSimulate, t_trace::lgTrace, TransitionProxy::Lo(), MAX2, NFE2LEVN, t_FeII::nFeIILevel_local, t_FeII::nFeIILevel_malloc, EmissionProxy::Pdest(), EmissionProxy::Pelec_esc(), EmissionProxy::Pesc(), EmissionProxy::phots(), EmissionProxy::PopOpc(), EmissionProxy::pump(), SDIV(), SMALLFLOAT, t_iso_sp::st, thermal, trace, t_iso_sp::trans(), t_thermal::tsq1, EmissionProxy::xIntensity(), t_dense::xIonDense, xMatrix, and yVector.

Referenced by Fe2_cooling().

Here is the call graph for this function:

◆ FeIILyaPump()

STATIC void FeIILyaPump ( void )

Todo:: 2 change this number to speed of light.

Definition at line 2941 of file atom_feii.cpp.

References t_dense::AtomicWeight, EmissionProxy::Aul(), DEBUG_ENTRY, dense, TransitionProxy::Emis(), TransitionProxy::EnergyWN(), EnerLyaProf1, EnerLyaProf4, Fe2LevN, Fe2LPump, FeII, g, GetDopplerWidth(), TransitionProxy::Hi(), t_hydro::HLineWidth, hydro, ipFe2LevN, ipH1s, ipH2p, ipH_LIKE, ipHYDROGEN, ipIRON, iso_sp, iteration, t_FeII::lgLyaPumpOn, TransitionProxy::Lo(), MAX2, t_FeII::nFeIILevel_local, PhotOccNumLyaCenter, powi(), SMALLFLOAT, t_iso_sp::st, EmissionProxy::TauIn(), and EmissionProxy::TauTot().

Referenced by FeIILevelPops().

Here is the call graph for this function:

◆ FeIIPoint()

void FeIIPoint ( void )

FeIIPoint called by ContCreatePointers to create pointers for lines in large FeII atom

Definition at line 1456 of file atom_feii.cpp.

References abscf(), EmissionProxy::Aul(), t_rfield::chLineLabel, EmissionProxy::dampXvel(), DEBUG_ENTRY, TransitionProxy::Emis(), TransitionProxy::EnergyRyd(), TransitionProxy::EnergyWN(), Fe2LevN, FeII, g, EmissionProxy::gf(), TransitionProxy::ipCont(), ipFe2LevN, EmissionProxy::ipFine(), ipFineCont(), ipoint(), TransitionProxy::Lo(), t_FeII::nFeIILevel_malloc, EmissionProxy::opacity(), PI4, and rfield.

Referenced by ContCreatePointers().

Here is the call graph for this function:

◆ FeIIPrint()

void FeIIPrint ( void )

FeIIPrint

Definition at line 1340 of file atom_feii.cpp.

References DEBUG_ENTRY.

Referenced by PrtZone().

◆ FeIIPun1Depart()

void FeIIPun1Depart	(	FILE *	ioPUN,
		long int	nPUN
	)

send the departure coef for physical level nPUN to unit ioPUN

Parameters

ioPUN	the io unit where the print should be directed
nPUN	the physical (not c) number of the level

Definition at line 2082 of file atom_feii.cpp.

References ASSERT, DEBUG_ENTRY, Fe2DepCoef, FeII, and t_FeII::nFeIILevel_local.

Referenced by FeIIPunDepart().

◆ FeIIPunchColden()

void FeIIPunchColden ( FILE * ioPUN )

FeIIPunchColden save level energies, stat weights, column density

Parameters

*ioPUN file we will save to

Definition at line 1647 of file atom_feii.cpp.

References DEBUG_ENTRY, TransitionProxy::EnergyWN(), Fe2ColDen, Fe2LevN, FeII, g, TransitionProxy::Hi(), ipFe2LevN, and t_FeII::nFeIILevel_malloc.

Referenced by SaveDo().

Here is the call graph for this function:

◆ FeIIPunchLevels()

void FeIIPunchLevels ( FILE * ioPUN )

called by parse_punch, save level energies and stat weights

Parameters

ioPUN

Definition at line 1625 of file atom_feii.cpp.

References DEBUG_ENTRY, TransitionProxy::EnergyWN(), Fe2LevN, FeII, g, TransitionProxy::Hi(), ipFe2LevN, and t_FeII::nFeIILevel_malloc.

Referenced by SaveDo().

Here is the call graph for this function:

◆ FeIIPunchLineStuff()

void FeIIPunchLineStuff	(	FILE *	io,
		realnum	xLimit,
		long	index
	)

include FeII lines in punched optical depths, etc, called from SaveLineStuff

Parameters

io
xLimit
index

Definition at line 2820 of file atom_feii.cpp.

References t_dense::AtomicWeight, DEBUG_ENTRY, dense, Fe2LevN, FeII, GetDopplerWidth(), ipFe2LevN, ipIRON, t_FeII::nFeIILevel_local, and Save1Line().

Referenced by SaveLineStuff().

Here is the call graph for this function:

◆ FeIIPunchOpticalDepth()

void FeIIPunchOpticalDepth ( FILE * ioPUN )

called by FeIIPunchLevels, this creates the save feii line optical depths

Parameters

ioPUN

Definition at line 1673 of file atom_feii.cpp.

References DEBUG_ENTRY, TransitionProxy::Emis(), Fe2LevN, FeII, ipFe2LevN, t_FeII::nFeIILevel_malloc, EmissionProxy::TauIn(), and TransitionProxy::WLAng().

Referenced by SaveDo().

Here is the call graph for this function:

◆ FeIIPunData()

void FeIIPunData	(	FILE *	ioPUN,
		bool	lgDoAll
	)

save line data for FeII atom

Parameters

ioPUN	io unit for save
lgDoAll	save all levels if true, only subset if false

Definition at line 1954 of file atom_feii.cpp.

References EmissionProxy::Aul(), cdEXIT, DEBUG_ENTRY, TransitionProxy::Emis(), EXIT_FAILURE, Fe2LevN, FeII, ioQQQ, ipFe2LevN, lgFeIIEverCalled, MIN2, ncs1, t_FeII::nFeIILevel_local, and Save1LineData().

Referenced by SaveLineData().

Here is the call graph for this function:

◆ FeIIPunDepart()

void FeIIPunDepart	(	FILE *	ioPUN,
		bool	lgDoAll
	)

save some departure coef for large atom, set with save feii departure command

Parameters

ioPUN
lgDoAll	option to save all dep coef if true

Definition at line 2029 of file atom_feii.cpp.

References DEBUG_ENTRY, FeII, FeIIPun1Depart(), t_FeII::nFeIILevel_local, and NLEVDEP.

Referenced by SaveDo().

Here is the call graph for this function:

◆ FeIIPunPop()

void FeIIPunPop	(	FILE *	ioPUN,
		bool	lgPunchRange,
		long int	ipRangeLo,
		long int	ipRangeHi,
		bool	lgPunchDensity
	)

save some level pops for large atom, set with save feii level populations command

Parameters

ioPUN
lgPunchRange	save range of levels if true, only selected subset if false
ipRangeLo	if ipPunchRange is true, this is lower bound of range on C scale
ipRangeHi	if ipPunchRange is true, this is upper bound of range on C scale
lgPunchDensity	flag saying whether to save density (cm-3, true) or relative population (flase)

Definition at line 2189 of file atom_feii.cpp.

References ASSERT, DEBUG_ENTRY, dense, Fe2LevelPop, FeII, ipIRON, MIN2, t_FeII::nFeIILevel_local, NLEVPOP, SDIV(), and t_dense::xIonDense.