/home66/gary/public_html/cloudy/c10_branch/source/continuum.h File Reference

#include "thirdparty.h"

Include dependency graph for continuum.h:

This graph shows which files directly or indirectly include this file:

Data Structures
struct	t_continuum
Functions
void	ContCreatePointers (void)
int	ContSetIntensity (void)
void	IncidentContinuumHere (void)
void	ContCreateMesh (void)
void	ContNegative (void)
double	ffun (double anu, double frac_beam_time, double frac_beam_const, double *frac_isotropic)
double	ffun (double anu)
double	ffun1 (double xnu)
void	outsum (double outtot, double outin, double *outout)
double	DrvContPump (transition *t, realnum DopplerWidth)
double	cont_gaunt_calc (double, double, double)
Variables
EXTERN struct t_continuum	continuum

Function Documentation

double cont_gaunt_calc	(	double	,
		double	,
		double
	)

cont_gaunt_calc do table look up of gaunt factor

Parameters:

	temp
	z
	photon

Definition at line 26 of file cont_gaunt.cpp.

References ASSERT, GauntIntegrand(), qg32(), and TE1RYD.

Referenced by dgaunt(), FillGFF(), SanityCheckBegin(), and SaveGaunts().

Here is the call graph for this function:

void ContCreateMesh ( void )

set up continuum energy mesh if first call, otherwise reset to original mesh

Definition at line 38 of file cont_createmesh.cpp.

References t_rfield::anu, t_rfield::anu2, t_rfield::AnuOrg, t_rfield::anusqr, cdEXIT, ChckFill(), t_rfield::ConEmitLocal, t_rfield::ConSourceFcnLocal, continuum, t_rfield::csigc, t_rfield::csigh, DEBUG_ENTRY, t_rfield::egamry, t_rfield::emm, fill(), t_rfield::fine_anu, t_rfield::fine_ener_hi, t_rfield::fine_ener_lo, geometry, ioQQQ, t_rfield::ipnt_coarse_2_fine, t_iterations::iter_malloc, iterations, t_trace::lgTrace, MALLOC, MAX2, NCELL, t_geometry::nend, t_geometry::nend_max, t_rfield::nfine_malloc, t_rfield::nflux, t_continuum::nrange, t_continuum::nStoredBands, t_rfield::nupper, read_continuum_mesh(), rfield, rfield_opac_malloc(), rfield_opac_zero(), t_continuum::StoredEnergy, t_continuum::StoredResolution, trace, t_rfield::trans_coef_total, t_rfield::trans_coef_zone, and t_rfield::widflx.

Referenced by cloudy(), GridGatherInCloudy(), and ParseCompile().

Here is the call graph for this function:

void ContCreatePointers ( void )

ContCreatePointers create pointers for lines and continua, one time per coreload

Todo:: 2 this is redundant with contents of oxygen line arrays use them instead when removing this, make sure all line intensity predictions also go into oi line arrays

Todo:: 2 what are we trying to print here?

Definition at line 59 of file cont_createpointers.cpp.

References abscf(), multi_arr< T, d, ALLOC, lgBC >::alloc(), t_rfield::anu, t_rfield::anu2, t_rfield::AnuOrg, t_rfield::anusqr, t_iso::As2nu, ASSERT, atmdat_2phot_shapefunction(), atoms, t_emission::Aul, cdEXIT, t_rfield::chContLabel, t_elementnames::chElementSym, chIonLbl(), t_elementnames::chIonStage, t_rfield::chLineLabel, chLineLbl(), multi_arr< T, d, ALLOC, lgBC >::clone(), ContBandsCreate(), continuum, t_emission::damp, t_emission::dampXvel, dBaseLines, DEBUG_ENTRY, dense, eina(), ELECTRON_MASS, elementnames, transition::Emis, EN1RYD, t_rfield::EnerGammaRay, transition::EnergyErg, transition::EnergyK, t_continuum::EnergyKshell, transition::EnergyWN, ERG1CM, EVRYD, ExtraLymanLines, fe, FeIIPoint(), fiddle(), t_continuum::filbnd, t_rfield::fine_anu, fixit(), t_quantumState::g, GetGF(), t_emission::gf, H2_ContPoint(), H2_Create(), he, Heavy, HFLines, transition::Hi, hmi, t_oxy::i2d, t_oxy::i2p, t_opac::ica2ex, t_opac::ih2pnt, t_hmi::iheh1, t_hmi::iheh2, t_opac::in1, t_fe2ovr_la::init_pointers(), Singleton< t_yield >::Inst(), Singleton< T >::Inst(), Singleton< t_fe2ovr_la >::Inst(), t_yield::ion_emit(), ionbal, ioQQQ, t_rfield::ip1000A, t_he::ip374, t_he::ip660, t_rfield::ipB_filter, t_opac::ipCKshell, t_ionbal::ipCompRecoil, transition::ipCont, ipContEnergy(), t_opac::ipElement, t_rfield::ipEnerGammaRay, t_fe::ipfe10, t_emission::ipFine, ipFineCont(), t_rfield::ipG0_DB96_hi, t_rfield::ipG0_DB96_lo, t_rfield::ipG0_spec_hi, t_rfield::ipG0_spec_lo, t_rfield::ipG0_TH85_hi, t_rfield::ipG0_TH85_lo, ipH1s, t_opac::ipH2_photo_thresh, ipH2p, ipH2s, ipH_LIKE, t_iso::ipHalfTwoPhoE, ipHe1s1S, ipHe2s1S, ipHE_LIKE, t_Heavy::ipHeavy, ipHELIUM, t_hmi::iphmin, ipHYDROGEN, t_iso::ipIsoLevNIonCon, ipLineEnergy(), t_opac::ipmgex, t_opac::ipo1exc, t_opac::ipo3exc, t_opac::ipo3exc3, t_atoms::ipoiex, ipoint(), ipOXYGEN, t_opac::ippr, t_secondaries::ipSecIon, ipShells(), t_iso::ipSym2nu, t_iso::ipTwoPhoE, t_rfield::ipV_filter, t_rt::ipxry, iso, iso_create(), t_continuum::KshellLimit, t_trace::lgConBug, t_iso::lgDielRecom, t_dense::lgElmtOn, t_rt::lgFstOn, lgLinesAdded, t_trace::lgPointBug, lgStatesAdded, t_trace::lgTrace, t_trace::lgTrLine, LIMELM, linesAdded2, transition::Lo, MAX2, NCRS, t_yield::nelem(), nHFLines, NISO, nLevel1, t_yield::nlines(), t_iso::nLyman_malloc, t_Heavy::nsShells, t_iso::numLevels_max, t_rfield::nupper, nUTA, nWindLine, opac, t_emission::opacity, oxy, t_fe::pfe10, t_fe::pfe11a, t_fe::pfe11b, t_fe::pfe14, PI, PI4, pow(), multi_arr< T, d, ALLOC, lgBC >::reserve(), rfield, rt, RYDLAM, SatelliteLines, secondaries, t_yield::set_ipoint(), t_iso::SmallA, SMALLFLOAT, SPEEDLIGHT, T1CM, TauDummy, TauLine2, TauLines, trace, t_emission::tran, Transitions, TwoS, UTALines, t_Heavy::Valence_IP_Ryd, WAVNRYD, t_rfield::widflx, WL_B_FILT, WL_V_FILT, transition::WLAng, and t_iso::xIsoLevNIonRyd.

Referenced by cloudy(), and ParseCompile().

Here is the call graph for this function:

void ContNegative ( void )

ContNegative sanity check for negative continuum intensities

Definition at line 10 of file cont_negative.cpp.

References t_rfield::anu, t_rfield::chContLabel, t_rfield::chLineLabel, t_rfield::ConInterOut, DEBUG_ENTRY, dense, t_dense::eden, t_rfield::flux, ioQQQ, t_rfield::nflux, nzone, t_rfield::otscon, t_rfield::otslin, t_rfield::outlin, phycon, rfield, t_rfield::SummedCon, and t_phycon::te.

Referenced by iso_level(), and PrintRates().

int ContSetIntensity ( void )

ContSetIntensity derive intensity of incident continuum,

Returns:: returns 0 if continuum is ok, 1 if failure of some kind

Definition at line 100 of file cont_setintensity.cpp.

References t_rfield::anu, t_rfield::anu2, t_rfield::anu3, t_rfield::anulog, t_rfield::AnuOrg, ASSERT, aweigh, BIGFLOAT, called, cdEXIT, t_rfield::chContLabel, t_rfield::chLineLabel, t_rfield::chSpType, t_continuum::cn1216, t_continuum::cn4861, t_rfield::ConEmitOut, t_rfield::ConEmitReflec, conorm(), t_rfield::ConRefIncid, t_thermal::ConstTemp, t_rfield::ContBoltz, continuum, t_rfield::convoc, t_ionbal::CosRayIonRate, t_hextra::cryden, t_rfield::csigc, t_rfield::csigh, t_secondaries::csupra, DEBUG_ENTRY, dense, t_dense::density_low_limit, t_dense::eden, t_dense::EdenExtra, t_dense::EdenHCorr, t_dense::EdenSet, t_dense::EdenTrue, t_hextra::effneu, t_rfield::egamry, ELECTRON_MASS, ELEM_CHARGE_ESU, t_rfield::emm, EN1RYD, extin(), t_continuum::fbeta, ffun(), t_rfield::fine_anu, t_rfield::flux, t_rfield::flux_beam_const, t_rfield::flux_beam_const_save, t_rfield::flux_beam_time, t_rfield::flux_isotropic, t_rfield::flux_isotropic_save, t_rfield::flux_time_beam_save, t_rfield::flux_total_incident, t_rfield::FluxFaint, t_continuum::fluxv, FR1RYD, t_hextra::frcneu, fweigh, t_prt::fx1ryd, t_prt::GammaLumin, t_dense::gas_phase, Heavy, hextra, HNU3C2, HPLANCK, hydro, t_oxy::i2d, t_oxy::i2p, t_rfield::Illumination, Singleton< T >::Inst(), ionbal, t_dense::IonHigh, t_dense::IonLow, ioQQQ, ipCARBON, t_opac::ipCKshell, t_ionbal::ipCompRecoil, t_prt::ipeak, t_rfield::ipEnerGammaRay, ipH1s, ipH2p, ipH_LIKE, ipHE_LIKE, t_Heavy::ipHeavy, ipHELIUM, ipHYDROGEN, t_iso::ipIsoLevNIonCon, ipLITHIUM, t_opac::ipo3exc, ipoint(), ipOXYGEN, t_rfield::ipPlasma, t_rfield::ipPlasmax, t_rfield::ipSpec, iso, t_rfield::lgBeamed, t_trace::lgComBug, t_ionbal::lgCompRecoil, t_continuum::lgCon0, t_trace::lgConBug, t_dense::lgElmtOn, t_rfield::lgHabing, t_hydro::lgHInducImp, t_rfield::lgHionRad, t_hextra::lgNeutrnHeatOn, t_noexec::lgNoExec, t_rfield::lgOcc1Hi, t_rfield::lgPlasNu, t_trace::lgPtrace, t_dense::lgSetIoniz, t_called::lgTalk, t_thermal::lgTemperatureConstant, t_rfield::lgTimeVary, t_trace::lgTrace, LIMELM, max(), MAX2, MIN2, t_rfield::nfine, t_rfield::nfine_malloc, t_rfield::nflux, noexec, t_rfield::nShape, t_iso::numLevels_max, t_rfield::nupper, t_rfield::occ1nu, t_rfield::occmax, t_rfield::occmnu, t_rfield::OccNumbIncidCont, opac, t_rfield::OpticalDepthScaleFactor, oxy, t_prt::pbal, phycon, PI, PI4, t_rfield::plsfrq, t_rfield::plsfrqmax, pow(), POW2, pow2(), t_prt::powion, t_prt::pradio, prt, ptrcer(), t_prt::q, t_rfield::qbal, t_prt::qgam, t_rfield::qhe, t_rfield::qheii, t_rfield::qhtot, t_rfield::qrad, t_rfield::qtot, t_prt::qx, rfield, RT_OTS_Zero(), secondaries, t_secondaries::SetCsupra, t_dense::SetIoniz, ShowMe(), t_rfield::slope, SMALLFLOAT, SPEEDLIGHT, StatesElemNEW, sumcon(), t_continuum::sv1216, t_continuum::sv4861, t_rfield::tbr4nu, t_rfield::tbrmax, t_rfield::tbrmnu, t_phycon::te, TE1RYD, t_phycon::telogn, t_phycon::TEnerDen, thermal, t_timesc::TimeErode, timesc, TorF(), TotalInsanity(), t_continuum::TotalLumin, t_continuum::totlsv, t_hextra::totneu, trace, t_rfield::uh, t_rfield::uheii, t_rfield::widflx, t_dense::xIonDense, and t_prt::xpow.

Referenced by cloudy().

Here is the call graph for this function:

double DrvContPump	(	transition *	t,
		realnum	DopplerWidth
	)

DrvContPump local continuum pumping rate radiative transfer for all lines

Parameters:

	*t
	DopplerWidth

Todo:: 2 if used, add damp as arg since calling routine probably evaluated it

Definition at line 32 of file cont_pump.cpp.

References a0, BREAK_, t_emission::damp, my_Integrand_con_pump_op::damp, DEBUG_ENTRY, t_radius::dRNeff, transition::Emis, FITTED, t_emission::iRedisFun, t_rfield::lgInducProcess, t_emission::opacity, t_emission::PopOpc, my_Integrand_con_pump_op::PumpTau, radius, rfield, Integrator< Integrand, Method >::sum(), and t_emission::TauIn.

Referenced by ParseDrive().

Here is the call graph for this function:

double ffun ( double anu )

ffun version without fractions

Definition at line 17 of file cont_ffun.cpp.

References DEBUG_ENTRY, and ffun().

Referenced by ContSetIntensity(), ffun(), and IncidentContinuumHere().

Here is the call graph for this function:

double ffun	(	double	anu,
		double *	frac_beam_time,
		double *	frac_beam_const,
		double *	frac_isotropic
	)

ffun evaluate total flux for sum of all continuum sources

Parameters:

	anu	photon energy (Rydberg) where continuum is evaluated
	frac_beam_time	fraction of beamed continuum that is varies with time
	frac_beam_const	fraction of beamed continuum that is constant
	frac_isotropic	fraction of continuum that is isotropic