cloudy
trunk
|
#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"
Go to the source code of this file.
Macros | |
#define | NO_ATOMS(ND) (gv.bin[ND]->AvVol*gv.bin[ND]->dustp[0]/ATOMIC_MASS_UNIT/gv.bin[ND]->atomWeight) |
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 | |
STATIC void | qheat_init (long, double[], double *) |
STATIC void | GetProbDistr_LowLimit (long, double, double, double, double[], double[], double[], double[], double[], long *, double *, long *, QH_Code *) |
STATIC double | TryDoubleStep (double[], double[], double[], double[], double[], double[], double[], double, double *, long, long, bool *) |
STATIC double | log_integral (double, double, double, double) |
STATIC void | ScanProbDistr (double[], double[], long, double, long, long *, long *, long *, long *, QH_Code *) |
STATIC long | RebinQHeatResults (long, long, long, double[], double[], double[], double[], double[], double[], double[], QH_Code *) |
STATIC void | GetProbDistr_HighLimit (long, double, double, double, double[], double[], double[], double *, long *, double *, QH_Code *) |
STATIC double | uderiv (double, long) |
STATIC double | ufunct (double, long, bool *) |
STATIC double | inv_ufunct (double, long, bool *) |
STATIC double | DebyeDeriv (double, long) |
void | GrainMakeDiffuse (void) |
void | qheat (double qtemp[], double qprob[], long int *qnbin, long int nd) |
STATIC void | GetProbDistr_LowLimit (long int nd, double rel_tol, double Umax, double fwhm, double Phi[], double PhiDrv[], double qtemp[], double qprob[], double dPdlnT[], long int *qnbin, double *new_tmin, long *nWideFail, QH_Code *ErrorCode) |
STATIC double | uderiv (double temp, long int nd) |
STATIC double | ufunct (double temp, long int nd, bool *lgBoundErr) |
STATIC double | inv_ufunct (double enthalpy, long int nd, bool *lgBoundErr) |
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] |
#define NO_ATOMS | ( | ND) | (gv.bin[ND]->AvVol*gv.bin[ND]->dustp[0]/ATOMIC_MASS_UNIT/gv.bin[ND]->atomWeight) |
Definition at line 17 of file grains_qheat.cpp.
Referenced by uderiv().
enum QH_Code |
Definition at line 24 of file grains_qheat.cpp.
STATIC double DebyeDeriv | ( | double | x, |
long | n | ||
) |
Definition at line 2556 of file grains_qheat.cpp.
References ASSERT, DEBUG_ENTRY, gauss_init(), gauss_legendre(), MALLOC, MAX2, POW2, POW3, powi(), and TotalInsanity().
Referenced by uderiv().
STATIC void GetProbDistr_HighLimit | ( | long | nd, |
double | TolFac, | ||
double | Umax, | ||
double | fwhm, | ||
double | qtemp[], | ||
double | qprob[], | ||
double | dPdlnT[], | ||
double * | tol, | ||
long * | qnbin, | ||
double * | new_tmin, | ||
QH_Code * | ErrorCode | ||
) |
Definition at line 2146 of file grains_qheat.cpp.
References t_gv::bin, GrainBin::chDstLab, GrainBin::cnv_GR_pCM3, GrainBin::cnv_GR_pH, GrainBin::cnv_H_pGR, DEBUG_ENTRY, GrainBin::dstems, GrainBin::dstslp, GrainBin::dstslp2, t_gv::dsttmp, GrainBin::DustEnth, EN1RYD, GrainBin::GrainHeat, 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().
STATIC void GetProbDistr_LowLimit | ( | long | , |
double | , | ||
double | , | ||
double | , | ||
double | [], | ||
double | [], | ||
double | [], | ||
double | [], | ||
double | [], | ||
long * | , | ||
double * | , | ||
long * | , | ||
QH_Code * | |||
) |
Referenced by qheat().
STATIC void GetProbDistr_LowLimit | ( | long int | nd, |
double | rel_tol, | ||
double | Umax, | ||
double | fwhm, | ||
double | Phi[], | ||
double | PhiDrv[], | ||
double | qtemp[], | ||
double | qprob[], | ||
double | dPdlnT[], | ||
long int * | qnbin, | ||
double * | new_tmin, | ||
long * | nWideFail, | ||
QH_Code * | ErrorCode | ||
) |
Definition at line 1220 of file grains_qheat.cpp.
References ASSERT, t_gv::bin, GrainBin::chDstLab, GrainBin::cnv_GR_pCM3, GrainBin::cnv_H_pGR, DEBUG_ENTRY, DEF_FAC, GrainBin::dstems, GrainBin::dstslp2, t_gv::dsttmp, EN1RYD, GRAIN_TMIN, GrainBin::GrainHeat, gv, ioQQQ, t_trace::lgDustBug, t_trace::lgTrace, LN_TWO, MAX2, MAX3, MAX_LOOP, MIN2, t_gv::nBin, 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, GrainBin::qtmin, RebinQHeatResults(), SAFETY, ScanProbDistr(), spldrv_safe(), splint_safe(), GrainBin::tedust, trace, TryDoubleStep(), uderiv(), ufunct(), and WIDE_FAIL_MAX.
void GrainMakeDiffuse | ( | void | ) |
main routine for generating the grain diffuse emission
Definition at line 165 of file grains_qheat.cpp.
References t_rfield::anu, t_rfield::anu2, ASSERT, t_gv::bin, GrainBin::BolFlux, cdEXIT, GrainBin::cnv_H_pCM3, CONSERV_TOL, t_thermal::ConstGrainTemp, DEBUG_ENTRY, dense, GrainBin::dstab1, GrainBin::dstAbund, t_gv::dstAbundThresholdFar, t_gv::dstAbundThresholdNear, EN1RYD, FR1RYD, t_gv::GasCoolColl, t_gv::GrainEmission, GrainBin::GrainHeat, t_gv::GrainHeatChem, t_gv::GrainHeatScaleFactor, t_gv::GrainHeatSum, t_gv::GraphiteEmission, gv, t_hmi::H2_total, t_thermal::heating, hmi, ioQQQ, ipHYDROGEN, lgAbort, t_gv::lgBakesPAH_heat, GrainBin::lgChrgConverged, t_gv::lgDColOn, t_gv::lgDHetOn, GrainBin::lgQHeat, GrainBin::lgTdustConverged, GrainBin::lgUseQHeat, MALLOC, GrainBin::matType, MAX2, t_gv::nBin, t_rfield::nflux, NQGRID, PI4, POW2, qheat(), GrainBin::RateDn, GrainBin::RateUp, rfield, t_gv::SilicateEmission, SPEEDLIGHT, STRG_CAR, STRG_SIL, TE1RYD, GrainBin::tedust, thermal, GrainBin::Tsublimat, t_gv::which_strg, t_rfield::widflx, and t_dense::xIonDense.
Referenced by RT_diffuse().
void InitEnthalpy | ( | void | ) |
initialize interpolation arrays for grain enthalpy
Definition at line 2490 of file grains_qheat.cpp.
References t_gv::bin, DEBUG_ENTRY, t_gv::dsttmp, GrainBin::DustEnth, GrainBin::EnthSlp, GrainBin::EnthSlp2, GRAIN_TMIN, gv, log_integral(), t_gv::nBin, NDEMS, spline(), tlim, and uderiv().
Referenced by GrainsInit().
STATIC double inv_ufunct | ( | double | , |
long | , | ||
bool * | |||
) |
Referenced by GetProbDistr_HighLimit(), qheat(), and TryDoubleStep().
STATIC double inv_ufunct | ( | double | enthalpy, |
long int | nd, | ||
bool * | lgBoundErr | ||
) |
Definition at line 2464 of file grains_qheat.cpp.
References ASSERT, t_gv::bin, cdEXIT, DEBUG_ENTRY, t_gv::dsttmp, GrainBin::DustEnth, GrainBin::EnthSlp2, gv, ioQQQ, t_gv::nBin, NDEMS, and splint_safe().
STATIC double log_integral | ( | double | xx1, |
double | yy1, | ||
double | xx2, | ||
double | yy2 | ||
) |
Definition at line 1788 of file grains_qheat.cpp.
References ASSERT, and DEBUG_ENTRY.
Referenced by GetProbDistr_HighLimit(), InitEnthalpy(), RebinQHeatResults(), and TryDoubleStep().
void qheat | ( | double | qtemp[], |
double | qprob[], | ||
long int * | qnbin, | ||
long int | nd | ||
) |
Definition at line 459 of file grains_qheat.cpp.
References t_rfield::anu, ASSERT, t_gv::bin, GrainBin::chDstLab, GrainBin::chrg, GrainBin::cnv_H_pCM3, GrainBin::cnv_H_pGR, CONSERV_TOL, DEBUG_ENTRY, DEF_FAC, GrainBin::dstems, GrainBin::dstslp2, t_gv::dsttmp, GrainBin::DustEnth, EN1RYD, ChargeBin::FracPop, FWHM_RATIO, FWHM_RATIO2, GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), GRAIN_TMIN, GrainBin::GrainCoolTherm, GrainBin::GrainHeat, GrainBin::GrainHeatColl, gv, GrainBin::HeatingRate1, ChargeBin::HeatingRate2, inv_ufunct(), ioQQQ, iterations, t_trace::lgDustBug, GrainBin::lgEverQHeat, t_iterations::lgLastIt, GrainBin::lgQHeat, t_gv::lgQHPunLast, GrainBin::lgQHTooWide, t_trace::lgTrace, GrainBin::lgUseQHeat, LN_TWO, LOOP_MAX, MALLOC, MAX2, MAX_EVENTS, MIN2, t_gv::nBin, GrainBin::nChrg, 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(), GrainBin::QHeatFailures, t_gv::QHPunchFile, GrainBin::qnflux, GrainBin::qtmin, GrainBin::qtmin_zone1, RELAX, rfield, spldrv_safe(), STRICT, GrainBin::tedust, TorF(), trace, uderiv(), ufunct(), WIDE_FAIL_MAX, and t_rfield::widflx.
STATIC void qheat_init | ( | long | nd, |
double | phiTilde[], | ||
double * | check | ||
) |
Definition at line 884 of file grains_qheat.cpp.
References t_rfield::anu, t_gv::anumax, t_gv::anumin, ASSERT, t_gv::bin, BOLTZMANN, GrainBin::chrg, GrainBin::cnv_H_pCM3, CONSERV_TOL, DEBUG_ENTRY, dprintf(), GrainBin::dstab1, ChargeBin::DustZ, ChargeBin::ehat, Ehi, EN1RYD, ChargeBin::fac1, ChargeBin::FracPop, GrainBin::GrainHeat, gv, GrainBin::HeatingRate1, ChargeBin::HeatingRate2, hunt_bisect(), ioQQQ, ChargeBin::ipThresInf, ChargeBin::ipThresInfVal, GrainBin::lgQHeat, t_gv::lgWD01, MALLOC, MAX2, MIN2, t_gv::nBin, GrainBin::nChrg, t_rfield::nflux, phycon, ChargeBin::PotSurf, ChargeBin::PotSurfInc, GrainBin::qnflux, GrainBin::qnflux2, rfield, sign(), t_rfield::SummedCon, t_phycon::te, GrainBin::tedust, ChargeBin::ThresInfInc, ChargeBin::ThresSurf, ChargeBin::ThresSurfVal, t_rfield::widflx, ChargeBin::yhat, and ChargeBin::yhat_primary.
Referenced by qheat().
STATIC long RebinQHeatResults | ( | long | nd, |
long | nstart, | ||
long | nend, | ||
double | p[], | ||
double | qtemp[], | ||
double | qprob[], | ||
double | dPdlnT[], | ||
double | u1[], | ||
double | delu[], | ||
double | Lambda[], | ||
QH_Code * | ErrorCode | ||
) |
Definition at line 1918 of file grains_qheat.cpp.
References ASSERT, t_gv::bin, GrainBin::cnv_GR_pCM3, GrainBin::cnv_GR_pH, CONSERV_TOL, DEBUG_ENTRY, GrainBin::dstems, GrainBin::dstslp, t_gv::dsttmp, EN1RYD, GrainBin::GrainHeat, gv, ioQQQ, t_trace::lgDustBug, t_trace::lgTrace, log_integral(), MALLOC, MAX2, MIN2, t_gv::nBin, 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().
STATIC void ScanProbDistr | ( | double | u1[], |
double | dPdlnT[], | ||
long | nbin, | ||
double | maxVal, | ||
long | nmax, | ||
long * | nstart, | ||
long * | nstart2, | ||
long * | nend, | ||
long * | nWideFail, | ||
QH_Code * | ErrorCode | ||
) |
Definition at line 1817 of file grains_qheat.cpp.
References ASSERT, DEBUG_ENTRY, 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().
STATIC double TryDoubleStep | ( | double | u1[], |
double | delu[], | ||
double | p[], | ||
double | qtemp[], | ||
double | Lambda[], | ||
double | Phi[], | ||
double | PhiDrv[], | ||
double | step, | ||
double * | cooling, | ||
long | index, | ||
long | nd, | ||
bool * | lgBoundFail | ||
) |
Definition at line 1615 of file grains_qheat.cpp.
References t_rfield::anu, ASSERT, t_gv::bin, GrainBin::cnv_H_pGR, DEBUG_ENTRY, GrainBin::dstems, GrainBin::dstslp2, t_gv::dsttmp, EN1RYD, gv, inv_ufunct(), log_integral(), MAX2, NDEMS, PROB_CUTOFF_LO, QHEAT_TOL, GrainBin::qnflux, rfield, and splint_safe().
Referenced by GetProbDistr_LowLimit().
STATIC double uderiv | ( | double | , |
long | |||
) |
Referenced by GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), InitEnthalpy(), qheat(), and RebinQHeatResults().
STATIC double uderiv | ( | double | temp, |
long int | nd | ||
) |
Definition at line 2320 of file grains_qheat.cpp.
References ASSERT, t_gv::bin, BOLTZMANN, cdEXIT, cval, DEBUG_ENTRY, DebyeDeriv(), EN1RYD, ENTH_CAR, ENTH_CAR2, ENTH_PAH, ENTH_PAH2, ENTH_SIL, ENTH_SIL2, gv, ioQQQ, GrainBin::matType, MAX2, MIN2, t_gv::nBin, NO_ATOMS, POW2, ppower, tlim, and t_gv::which_enth.
STATIC double ufunct | ( | double | , |
long | , | ||
bool * | |||
) |
Referenced by GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), qheat(), and RebinQHeatResults().
STATIC double ufunct | ( | double | temp, |
long int | nd, | ||
bool * | lgBoundErr | ||
) |
Definition at line 2438 of file grains_qheat.cpp.
References ASSERT, t_gv::bin, cdEXIT, DEBUG_ENTRY, t_gv::dsttmp, GrainBin::DustEnth, GrainBin::EnthSlp, gv, ioQQQ, t_gv::nBin, NDEMS, and splint_safe().
|
static |
|
static |
Definition at line 66 of file grains_qheat.cpp.
Referenced by GetProbDistr_LowLimit(), and qheat().
|
static |
Definition at line 115 of file grains_qheat.cpp.
|
static |
Definition at line 77 of file grains_qheat.cpp.
Referenced by qheat().
|
static |
Definition at line 80 of file grains_qheat.cpp.
Referenced by qheat().
|
static |
Definition at line 63 of file grains_qheat.cpp.
Referenced by cnewton(), lgCoolNetConverge(), and qheat().
|
static |
Definition at line 59 of file grains_qheat.cpp.
Referenced by qheat().
|
static |
Definition at line 83 of file grains_qheat.cpp.
Referenced by GetProbDistr_LowLimit().
|
static |
Definition at line 118 of file grains_qheat.cpp.
|
static |
Definition at line 47 of file grains_qheat.cpp.
Referenced by GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), RebinQHeatResults(), and ScanProbDistr().
|
static |
Definition at line 73 of file grains_qheat.cpp.
Referenced by GetProbDistr_LowLimit().
|
static |
Definition at line 55 of file grains_qheat.cpp.
Referenced by GetProbDistr_LowLimit().
|
static |
Definition at line 122 of file grains_qheat.cpp.
Referenced by uderiv().
|
static |
Definition at line 51 of file grains_qheat.cpp.
Referenced by GetProbDistr_HighLimit(), and GetProbDistr_LowLimit().
|
static |
Definition at line 50 of file grains_qheat.cpp.
Referenced by GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), qheat(), RebinQHeatResults(), and TryDoubleStep().
|
static |
Definition at line 69 of file grains_qheat.cpp.
Referenced by GetProbDistr_HighLimit(), and GetProbDistr_LowLimit().
|
static |
Definition at line 86 of file grains_qheat.cpp.
Referenced by GetProbDistr_LowLimit(), and TryDoubleStep().
|
static |
Definition at line 100 of file grains_qheat.cpp.
Referenced by GetProbDistr_HighLimit(), and RebinQHeatResults().
|
static |
Definition at line 93 of file grains_qheat.cpp.
Referenced by qheat().
|
static |
Definition at line 52 of file grains_qheat.cpp.
Referenced by gauss_legendre(), GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), RebinQHeatResults(), spldrv_safe(), and splint_safe().
|
static |
Definition at line 92 of file grains_qheat.cpp.
Referenced by GetProbDistr_HighLimit(), and qheat().
|
static |
Definition at line 121 of file grains_qheat.cpp.
Referenced by InitEnthalpy(), and uderiv().
|
static |
Definition at line 89 of file grains_qheat.cpp.
Referenced by GetProbDistr_LowLimit(), qheat(), and ScanProbDistr().