/home66/gary/public_html/cloudy/c13_branch/source/atom_feii.cpp File Reference

#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:

Go to the source code of this file.

Defines
#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

Define Documentation

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

Referenced by FeIILevelPops().

#define NLEVDEP 11

Referenced by FeIIPunDepart().

#define NLEVPOP 11

Referenced by FeIIPunPop().

#define NPRADDAT 159

Definition at line 100 of file atom_feii.cpp.

Referenced by FeIICollRatesBoltzmann(), and FeIICreate().

#define USE_OLD true

Referenced by FeIICreate().

Function Documentation

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().

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().

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().

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:

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:

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:

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:

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:

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:

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:

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:

STATIC void FeIICollRatesBoltzmann ( void )

Definition at line 1072 of file atom_feii.cpp.

References ASSERT, EmissionProxy::Aul(), cdEXIT, col_str, CollisionProxy::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:

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, nFeIIConBins, and pow().

Referenced by FeIICreate().

Here is the call graph for this function:

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:

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:

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:

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:

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:

void FeIIPrint ( void )

FeIIPrint

Definition at line 1340 of file atom_feii.cpp.

References DEBUG_ENTRY.

Referenced by PrtZone().

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().

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:

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:

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:

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:

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:

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:

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)