#include "cddefines.h"
#include "physconst.h"
#include "rfield.h"
#include "phycon.h"
#include "dense.h"
#include "hmi.h"
#include "thermal.h"
#include "trace.h"
#include "thirdparty.h"
#include "iterations.h"
#include "grainvar.h"
#include "grains.h"

Include dependency graph for grains_qheat.cpp:

Enumerations
enum	QH_Code { QH_OK , QH_ANALYTIC , QH_ANALYTIC_RELAX , QH_DELTA , QH_NEGRATE_FAIL , QH_LOOP_FAIL , QH_ARRAY_FAIL , QH_THIGH_FAIL , QH_RETRY , QH_CONV_FAIL , QH_BOUND_FAIL , QH_DELTA_FAIL , QH_NO_REBIN , QH_LOW_FAIL , QH_HIGH_FAIL , QH_STEP_FAIL , QH_FATAL , QH_WIDE_FAIL , QH_NBIN_FAIL , QH_REBIN_FAIL }

Functions
double	no_atoms (size_t nd)

STATIC void	qheat_init (size_t, vector< double > &, double *)

STATIC void	GetProbDistr_LowLimit (size_t, double, double, double, const vector< double > &, const vector< double > &, vector< double > &, vector< double > &, vector< double > &, long , double , long , QH_Code )

STATIC double	TryDoubleStep (vector< double > &, vector< double > &, vector< double > &, vector< double > &, vector< double > &, const vector< double > &, const vector< double > &, double, double , long, size_t, bool )

STATIC double	log_integral (double, double, double, double)

STATIC void	ScanProbDistr (const vector< double > &, const vector< double > &, long, double, long, double, long , long , long , long , QH_Code *)

STATIC long	RebinQHeatResults (size_t, long, long, vector< double > &, vector< double > &, vector< double > &, vector< double > &, vector< double > &, vector< double > &, vector< double > &, QH_Code *)

STATIC void	GetProbDistr_HighLimit (long, double, double, double, vector< double > &, vector< double > &, vector< double > &, double , long , double , QH_Code )

STATIC double	uderiv (double, size_t)

STATIC double	ufunct (double, size_t, bool *)

STATIC double	inv_ufunct (double, size_t, bool *)

STATIC double	DebyeDeriv (double, long)

void	GrainMakeDiffuse (void)

void	qheat (vector< double > &qtemp, vector< double > &qprob, long int *qnbin, size_t nd)

STATIC void	GetProbDistr_LowLimit (size_t nd, double rel_tol, double Umax, double fwhm, const vector< double > &Phi, const vector< double > &PhiDrv, vector< double > &qtemp, vector< double > &qprob, vector< double > &dPdlnT, long int qnbin, double new_tmin, long nWideFail, QH_Code ErrorCode)

void	InitEnthalpy (void)

Variables
static const long	NQMIN = 10L

static const double	PROB_CUTOFF_LO = 1.e-15

static const double	PROB_CUTOFF_HI = 1.e-20

static const double	SAFETY = 1.e+8

static const long	NSTARTUP = 5L

static const double	MAX_EVENTS = 150.

static const long	LOOP_MAX = 20L

static const double	DEF_FAC = 3.

static const double	PROB_TOL = 0.02

static const long	NQTEST = 500L

static const double	FWHM_RATIO = 0.1

static const double	FWHM_RATIO2 = 0.007

static const long	MAX_LOOP = 2*NQGRID

static const double	QHEAT_TOL = 5.e-3

static const long	WIDE_FAIL_MAX = 3

static const double	STRICT = 1.

static const double	RELAX = 15.

static const double	QT_RATIO = 1.03

static const double	DEN_SIL = 3.30

static const double	MW_SIL = 24.6051

static const double	tlim [5] ={0.,50.,150.,500.,DBL_MAX}

static const double	ppower [4] ={2.00,1.30,0.68,0.00}

static const double	cval [4]

Enumeration Type Documentation

◆ QH_Code

enum QH_Code

Enumerator
QH_OK
QH_ANALYTIC
QH_ANALYTIC_RELAX
QH_DELTA
QH_NEGRATE_FAIL
QH_LOOP_FAIL
QH_ARRAY_FAIL
QH_THIGH_FAIL
QH_RETRY
QH_CONV_FAIL
QH_BOUND_FAIL
QH_DELTA_FAIL
QH_NO_REBIN
QH_LOW_FAIL
QH_HIGH_FAIL
QH_STEP_FAIL
QH_FATAL
QH_WIDE_FAIL
QH_NBIN_FAIL
QH_REBIN_FAIL

Definition at line 27 of file grains_qheat.cpp.

Function Documentation

◆ DebyeDeriv()

STATIC double DebyeDeriv	(	double	x,
		long	n
	)

Definition at line 2545 of file grains_qheat.cpp.

References ASSERT, DEBUG_ENTRY, gauss_init(), gauss_legendre(), MAX2, POW2, POW3, powi(), and TotalInsanity().

Referenced by uderiv().

Here is the call graph for this function:

◆ GetProbDistr_HighLimit()

STATIC void GetProbDistr_HighLimit	(	long	nd,
		double	TolFac,
		double	Umax,
		double	fwhm,
		vector< double > &	qtemp,
		vector< double > &	qprob,
		vector< double > &	dPdlnT,
		double *	tol,
		long *	qnbin,
		double *	new_tmin,
		QH_Code *	ErrorCode
	)

Definition at line 2138 of file grains_qheat.cpp.

References GrainVar::bin, DEBUG_ENTRY, GrainVar::dsttmp, EN1RYD, gv, inv_ufunct(), ioQQQ, t_trace::lgDustBug, t_trace::lgTrace, LN_TWO, log_integral(), MAX2, MIN2, NDEMS, NQGRID, NQMIN, nzone, PI, POW2, PROB_CUTOFF_HI, PROB_CUTOFF_LO, PROB_TOL, QH_ANALYTIC, QH_ANALYTIC_RELAX, QH_BOUND_FAIL, QH_CONV_FAIL, QH_THIGH_FAIL, QH_WIDE_FAIL, QT_RATIO, SAFETY, splint_safe(), STRICT, trace, uderiv(), and ufunct().

Referenced by qheat().

Here is the call graph for this function:

◆ GetProbDistr_LowLimit() [1/2]

STATIC void GetProbDistr_LowLimit	(	size_t	nd,
		double	rel_tol,
		double	Umax,
		double	fwhm,
		const vector< double > &	Phi,
		const vector< double > &	PhiDrv,
		vector< double > &	qtemp,
		vector< double > &	qprob,
		vector< double > &	dPdlnT,
		long int *	qnbin,
		double *	new_tmin,
		long *	nWideFail,
		QH_Code *	ErrorCode
	)

Definition at line 1204 of file grains_qheat.cpp.

References ASSERT, GrainVar::bin, DEBUG_ENTRY, DEF_FAC, GrainVar::dsttmp, EN1RYD, GRAIN_TMIN, gv, ioQQQ, t_trace::lgDustBug, t_trace::lgTrace, LN_TWO, MAX2, MAX3, MAX_LOOP, MIN2, NDEMS, NQGRID, NQMIN, NQTEST, NSTARTUP, nzone, POW3, PROB_CUTOFF_HI, PROB_CUTOFF_LO, PROB_TOL, QH_ARRAY_FAIL, QH_BOUND_FAIL, QH_CONV_FAIL, QH_DELTA_FAIL, QH_HIGH_FAIL, QH_LOOP_FAIL, QH_LOW_FAIL, QH_NBIN_FAIL, QH_NO_REBIN, QH_STEP_FAIL, QH_THIGH_FAIL, QH_WIDE_FAIL, QHEAT_TOL, RebinQHeatResults(), SAFETY, ScanProbDistr(), spldrv_safe(), splint_safe(), trace, TryDoubleStep(), uderiv(), ufunct(), and WIDE_FAIL_MAX.

Here is the call graph for this function:

◆ GetProbDistr_LowLimit() [2/2]

STATIC void GetProbDistr_LowLimit	(	size_t	,
		double	,
		double	,
		double	,
		const vector< double > &	,
		const vector< double > &	,
		vector< double > &	,
		vector< double > &	,
		vector< double > &	,
		long *	,
		double *	,
		long *	,
		QH_Code *
	)

Referenced by qheat().

◆ GrainMakeDiffuse()

void GrainMakeDiffuse ( void )

main routine for generating the grain diffuse emission

Definition at line 170 of file grains_qheat.cpp.

References t_rfield::anu, t_rfield::anu2, ASSERT, GrainVar::bin, CONSERV_TOL, t_thermal::ConstGrainTemp, DEBUG_ENTRY, dense, GrainVar::dstAbundThresholdFar, GrainVar::dstAbundThresholdNear, EN1RYD, FR1RYD, GrainVar::GasCoolColl, GrainVar::GrainEmission, GrainVar::GrainHeatChem, GrainVar::GrainHeatScaleFactor, GrainVar::GrainHeatSum, GrainVar::GraphiteEmission, gv, t_hmi::H2_total, t_thermal::heating, hmi, ipHYDROGEN, lgAbort, GrainVar::lgBakesPAH_heat, GrainVar::lgDColOn, GrainVar::lgDHetOn, MAX2, t_rfield::nflux, NQGRID, PI4, POW2, qheat(), rfield, GrainVar::SilicateEmission, SPEEDLIGHT, STRG_CAR, STRG_SIL, TE1RYD, thermal, TotalInsanity(), GrainVar::which_strg, t_rfield::widflx, and t_dense::xIonDense.

Referenced by RT_diffuse().

Here is the call graph for this function:

◆ InitEnthalpy()

void InitEnthalpy ( void )

initialize interpolation arrays for grain enthalpy

Definition at line 2482 of file grains_qheat.cpp.

References GrainVar::bin, DEBUG_ENTRY, GrainVar::dsttmp, GRAIN_TMIN, gv, log_integral(), NDEMS, spline(), tlim, and uderiv().

Referenced by GrainsInit().

Here is the call graph for this function:

◆ inv_ufunct()

STATIC double inv_ufunct	(	double	enthalpy,
		size_t	nd,
		bool *	lgBoundErr
	)

Definition at line 2456 of file grains_qheat.cpp.

References ASSERT, GrainVar::bin, cdEXIT, DEBUG_ENTRY, GrainVar::dsttmp, EXIT_FAILURE, gv, ioQQQ, NDEMS, and splint_safe().

Referenced by GetProbDistr_HighLimit(), qheat(), and TryDoubleStep().

Here is the call graph for this function:

◆ log_integral()

STATIC double log_integral	(	double	xx1,
		double	yy1,
		double	xx2,
		double	yy2
	)

Definition at line 1778 of file grains_qheat.cpp.

References ASSERT, and DEBUG_ENTRY.

Referenced by GetProbDistr_HighLimit(), InitEnthalpy(), RebinQHeatResults(), and TryDoubleStep().

◆ no_atoms()

double no_atoms ( size_t nd )

inline

Definition at line 17 of file grains_qheat.cpp.

References ATOMIC_MASS_UNIT, GrainVar::bin, and gv.

Referenced by mie_read_form(), and uderiv().

◆ qheat()

void qheat	(	vector< double > &	qtemp,
		vector< double > &	qprob,
		long int *	qnbin,
		size_t	nd
	)

Definition at line 450 of file grains_qheat.cpp.

References t_rfield::anu, ASSERT, GrainVar::bin, CONSERV_TOL, DEBUG_ENTRY, DEF_FAC, GrainVar::dsttmp, EN1RYD, FWHM_RATIO, FWHM_RATIO2, GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), GRAIN_TMIN, gv, inv_ufunct(), ioQQQ, iterations, t_trace::lgDustBug, t_iterations::lgLastIt, GrainVar::lgQHPunLast, t_trace::lgTrace, LN_TWO, LOOP_MAX, MAX2, MAX_EVENTS, MIN2, NDEMS, NQGRID, t_rfield::nupper, nzone, POW2, PROB_CUTOFF_LO, QH_ANALYTIC, QH_CONV_FAIL, QH_DELTA, QH_DELTA_FAIL, QH_FATAL, QH_HIGH_FAIL, QH_LOW_FAIL, QH_NEGRATE_FAIL, QH_NO_REBIN, QH_OK, QH_RETRY, QH_WIDE_FAIL, qheat_init(), GrainVar::QHSaveFile, RELAX, rfield, spldrv_safe(), STRICT, TorF(), trace, uderiv(), ufunct(), WIDE_FAIL_MAX, and t_rfield::widflx.

Referenced by GrainMakeDiffuse().

Here is the call graph for this function:

◆ qheat_init()

STATIC void qheat_init	(	size_t	nd,
		vector< double > &	phiTilde,
		double *	check
	)

Definition at line 866 of file grains_qheat.cpp.

References t_rfield::anu, GrainVar::anumax, GrainVar::anumin, ASSERT, GrainVar::bin, BOLTZMANN, CONSERV_TOL, DEBUG_ENTRY, dprintf(), ChargeBin::DustZ, ChargeBin::ehat, EN1RYD, ChargeBin::fac1, ChargeBin::FracPop, gv, ChargeBin::HeatingRate2, hunt_bisect(), ioQQQ, ChargeBin::ipThresInf, ChargeBin::ipThresInfVal, GrainVar::lgWD01, MAX2, MIN2, t_rfield::nflux, phycon, ChargeBin::PotSurf, ChargeBin::PotSurfInc, rfield, sign(), t_rfield::SummedCon, t_phycon::te, ChargeBin::ThresInfInc, ChargeBin::ThresSurf, ChargeBin::ThresSurfVal, t_rfield::widflx, ChargeBin::yhat, and ChargeBin::yhat_primary.

Referenced by qheat().

Here is the call graph for this function:

◆ RebinQHeatResults()

STATIC long RebinQHeatResults	(	size_t	nd,
		long	nstart,
		long	nend,
		vector< double > &	p,
		vector< double > &	qtemp,
		vector< double > &	qprob,
		vector< double > &	dPdlnT,
		vector< double > &	u1,
		vector< double > &	delu,
		vector< double > &	Lambda,
		QH_Code *	ErrorCode
	)

Definition at line 1913 of file grains_qheat.cpp.

References ASSERT, GrainVar::bin, CONSERV_TOL, DEBUG_ENTRY, GrainVar::dsttmp, EN1RYD, gv, ioQQQ, t_trace::lgDustBug, t_trace::lgTrace, log_integral(), MAX2, MIN2, NDEMS, NQGRID, NQMIN, PROB_CUTOFF_LO, QH_CONV_FAIL, QH_REBIN_FAIL, QT_RATIO, SAFETY, splint_safe(), trace, uderiv(), and ufunct().

Referenced by GetProbDistr_LowLimit().

Here is the call graph for this function:

◆ ScanProbDistr()

STATIC void ScanProbDistr	(	const vector< double > &	u1,
		const vector< double > &	dPdlnT,
		long	nbin,
		double	maxVal,
		long	nmax,
		double	qtmin1,
		long *	nstart,
		long *	nstart2,
		long *	nend,
		long *	nWideFail,
		QH_Code *	ErrorCode
	)

Definition at line 1807 of file grains_qheat.cpp.

References ASSERT, DEBUG_ENTRY, GRAIN_TMIN, ioQQQ, t_trace::lgDustBug, t_trace::lgTrace, MAX2, NQMIN, QH_HIGH_FAIL, QH_LOW_FAIL, QH_NBIN_FAIL, QH_NO_REBIN, QH_WIDE_FAIL, trace, and WIDE_FAIL_MAX.

Referenced by GetProbDistr_LowLimit().

◆ TryDoubleStep()

STATIC double TryDoubleStep	(	vector< double > &	u1,
		vector< double > &	delu,
		vector< double > &	p,
		vector< double > &	qtemp,
		vector< double > &	Lambda,
		const vector< double > &	Phi,
		const vector< double > &	PhiDrv,
		double	step,
		double *	cooling,
		long	index,
		size_t	nd,
		bool *	lgBoundFail
	)

Definition at line 1605 of file grains_qheat.cpp.

References t_rfield::anu, ASSERT, GrainVar::bin, DEBUG_ENTRY, GrainVar::dsttmp, EN1RYD, gv, inv_ufunct(), log_integral(), MAX2, NDEMS, NQGRID, PROB_CUTOFF_LO, QHEAT_TOL, rfield, and splint_safe().

Referenced by GetProbDistr_LowLimit().

Here is the call graph for this function:

◆ uderiv()

STATIC double uderiv	(	double	temp,
		size_t	nd
	)

Definition at line 2312 of file grains_qheat.cpp.

References ASSERT, GrainVar::bin, BOLTZMANN, cdEXIT, cval, DEBUG_ENTRY, DebyeDeriv(), EN1RYD, ENTH_CAR, ENTH_CAR2, ENTH_PAH, ENTH_PAH2, ENTH_SIL, ENTH_SIL2, EXIT_FAILURE, gv, ioQQQ, MAX2, MIN2, no_atoms(), POW2, ppower, tlim, and GrainVar::which_enth.

Referenced by GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), InitEnthalpy(), qheat(), and RebinQHeatResults().

Here is the call graph for this function:

◆ ufunct()

STATIC double ufunct	(	double	temp,
		size_t	nd,
		bool *	lgBoundErr
	)

Definition at line 2430 of file grains_qheat.cpp.

References ASSERT, GrainVar::bin, cdEXIT, DEBUG_ENTRY, GrainVar::dsttmp, EXIT_FAILURE, gv, ioQQQ, NDEMS, and splint_safe().

Referenced by GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), qheat(), and RebinQHeatResults().

Here is the call graph for this function:

Variable Documentation

◆ cval

const double cval[4]

static

Initial value:

={
        1.40e3/DEN_SIL*ATOMIC_MASS_UNIT*MW_SIL/EN1RYD,
        2.20e4/DEN_SIL*ATOMIC_MASS_UNIT*MW_SIL/EN1RYD,
        4.80e5/DEN_SIL*ATOMIC_MASS_UNIT*MW_SIL/EN1RYD,
        3.41e7/DEN_SIL*ATOMIC_MASS_UNIT*MW_SIL/EN1RYD}

Definition at line 126 of file grains_qheat.cpp.

Referenced by uderiv().

◆ DEF_FAC

const double DEF_FAC = 3.

static

Definition at line 69 of file grains_qheat.cpp.

Referenced by GetProbDistr_LowLimit(), and qheat().

◆ DEN_SIL

const double DEN_SIL = 3.30

static

Definition at line 118 of file grains_qheat.cpp.

◆ FWHM_RATIO

const double FWHM_RATIO = 0.1

static

Definition at line 80 of file grains_qheat.cpp.

Referenced by qheat().

◆ FWHM_RATIO2

const double FWHM_RATIO2 = 0.007

static

Definition at line 83 of file grains_qheat.cpp.

Referenced by qheat().

◆ LOOP_MAX

const long LOOP_MAX = 20L

static

Definition at line 66 of file grains_qheat.cpp.

Referenced by cnewton(), lgCoolNetConverge(), and qheat().

◆ MAX_EVENTS

const double MAX_EVENTS = 150.

static

Definition at line 62 of file grains_qheat.cpp.

Referenced by qheat().

◆ MAX_LOOP

const long MAX_LOOP = 2*NQGRID

static

Definition at line 86 of file grains_qheat.cpp.

Referenced by GetProbDistr_LowLimit().

◆ MW_SIL

const double MW_SIL = 24.6051

static

Definition at line 121 of file grains_qheat.cpp.

◆ NQMIN

const long NQMIN = 10L

static

Definition at line 50 of file grains_qheat.cpp.

Referenced by GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), RebinQHeatResults(), and ScanProbDistr().

◆ NQTEST

const long NQTEST = 500L

static

Definition at line 76 of file grains_qheat.cpp.

Referenced by GetProbDistr_LowLimit().

◆ NSTARTUP

const long NSTARTUP = 5L

static

Definition at line 58 of file grains_qheat.cpp.

Referenced by GetProbDistr_LowLimit().

◆ ppower

const double ppower[4] ={2.00,1.30,0.68,0.00}

static

Definition at line 125 of file grains_qheat.cpp.

Referenced by uderiv().

◆ PROB_CUTOFF_HI

const double PROB_CUTOFF_HI = 1.e-20

static

Definition at line 54 of file grains_qheat.cpp.

Referenced by GetProbDistr_HighLimit(), and GetProbDistr_LowLimit().

◆ PROB_CUTOFF_LO

const double PROB_CUTOFF_LO = 1.e-15

static

Definition at line 53 of file grains_qheat.cpp.

Referenced by GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), qheat(), RebinQHeatResults(), and TryDoubleStep().

◆ PROB_TOL

const double PROB_TOL = 0.02

static

Definition at line 72 of file grains_qheat.cpp.

Referenced by GetProbDistr_HighLimit(), and GetProbDistr_LowLimit().

◆ QHEAT_TOL

const double QHEAT_TOL = 5.e-3

static

Definition at line 89 of file grains_qheat.cpp.

Referenced by GetProbDistr_LowLimit(), and TryDoubleStep().

◆ QT_RATIO

const double QT_RATIO = 1.03

static

Definition at line 103 of file grains_qheat.cpp.

Referenced by GetProbDistr_HighLimit(), and RebinQHeatResults().

◆ RELAX

const double RELAX = 15.

static

Definition at line 96 of file grains_qheat.cpp.

Referenced by qheat().

◆ SAFETY

const double SAFETY = 1.e+8

static

Definition at line 55 of file grains_qheat.cpp.

Referenced by gauss_legendre(), GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), radius_next(), RebinQHeatResults(), spldrv_safe(), and splint_safe().

◆ STRICT

const double STRICT = 1.

static

Definition at line 95 of file grains_qheat.cpp.

Referenced by GetProbDistr_HighLimit(), and qheat().

◆ tlim

const double tlim[5] ={0.,50.,150.,500.,DBL_MAX}

static

Definition at line 124 of file grains_qheat.cpp.

Referenced by InitEnthalpy(), and uderiv().

◆ WIDE_FAIL_MAX

const long WIDE_FAIL_MAX = 3

static

Definition at line 92 of file grains_qheat.cpp.

Referenced by GetProbDistr_LowLimit(), qheat(), and ScanProbDistr().

Enumerations

Functions

Variables

Enumeration Type Documentation

◆ QH_Code

Function Documentation

◆ DebyeDeriv()

◆ GetProbDistr_HighLimit()

◆ GetProbDistr_LowLimit() [1/2]

◆ GetProbDistr_LowLimit() [2/2]

◆ GrainMakeDiffuse()

◆ InitEnthalpy()

◆ inv_ufunct()

◆ log_integral()

◆ no_atoms()

◆ qheat()

◆ qheat_init()

◆ RebinQHeatResults()

◆ ScanProbDistr()

◆ TryDoubleStep()

◆ uderiv()

◆ ufunct()

Variable Documentation

◆ cval

◆ DEF_FAC

◆ DEN_SIL

◆ FWHM_RATIO

◆ FWHM_RATIO2

◆ LOOP_MAX

◆ MAX_EVENTS

◆ MAX_LOOP

◆ MW_SIL

◆ NQMIN

◆ NQTEST

◆ NSTARTUP

◆ ppower

◆ PROB_CUTOFF_HI

◆ PROB_CUTOFF_LO

◆ PROB_TOL

◆ QHEAT_TOL

◆ QT_RATIO

◆ RELAX

◆ SAFETY

◆ STRICT

◆ tlim

◆ WIDE_FAIL_MAX