#include <grainvar.h>

Collaboration diagram for GrainVar:

Public Member Functions
	GrainVar ()

	~GrainVar ()

void	clear ()

bool	lgDustOn () const

Data Fields
bool	lgWD01

bool	lgReevaluate

bool	lgGrainPhysicsOn

bool	lgAnyDustVary

bool	lgBakesPAH_heat

bool	lgNegGrnDrg

bool	lgDHetOn

bool	lgQHeatOn

bool	lgDColOn

bool	lgGrainElectrons

vector< string >	ReadRecord

string	chPAH_abundance

realnum	GrainMetal

realnum	elmSumAbund [LIMELM]

enth_type	which_enth [MAT_TOP]

zmin_type	which_zmin [MAT_TOP]

pot_type	which_pot [MAT_TOP]

ial_type	which_ial [MAT_TOP]

pe_type	which_pe [MAT_TOP]

strg_type	which_strg [MAT_TOP]

H2_type	which_H2distr [MAT_TOP]

long	nzone

vector< double >	dstab

vector< double >	dstsc

double	TotalEden

realnum	GrnElecDonateMax

realnum	GrnElecHoldMax

double	GrnRecomTe

long	nChrgRequested

AEInfo *	AugerData [LIMELM]

realnum	GrainChTrRate [LIMELM][LIMELM+1][LIMELM+1]

double	GasCoolColl

double	GasHeatPhotoEl

double	GasHeatTherm

double	GasHeatNet

double	GrainHeatSum

double	GrainHeatLya

double	GrainHeatDif

double	GrainHeatInc

double	GrainHeatCollSum

double	GrainHeatChem

double	dHeatdT

realnum	GrainHeatScaleFactor

realnum	TotalDustHeat

realnum	dphmax

realnum	dclmax

double	dsttmp [NDEMS]

realnum	dstAbundThresholdNear

realnum	dstAbundThresholdFar

bool	lgQHeatAll

bool	lgQHPunLast

FILE *	QHSaveFile

double	rate_h2_form_grains_used_total

vector< realnum >	GrainEmission

vector< realnum >	GraphiteEmission

vector< realnum >	SilicateEmission

vector< GrainBin * >	bin

Private Member Functions
void	p_clear0 ()

void	p_clear1 ()

Detailed Description

NB NB NB NB NB NB

this is the master data structure for grain physics, it is statically allocated

all entries that depend on grain type should go in the GrainBin structure above

NB NB NB NB NB NB

Definition at line 454 of file grainvar.h.

Constructor & Destructor Documentation

GrainVar::GrainVar ( )

inline

Definition at line 460 of file grainvar.h.

References p_clear1().

Here is the call graph for this function:

GrainVar::~GrainVar ( )

inline

Definition at line 464 of file grainvar.h.

References p_clear0().

Here is the call graph for this function:

Member Function Documentation

void GrainVar::clear ( )

inline

Definition at line 468 of file grainvar.h.

References p_clear0(), and p_clear1().

Referenced by GrainZero(), and grid_do().

Here is the call graph for this function:

bool GrainVar::lgDustOn ( ) const

inline

grain logic < have any grains been switched on ?

Definition at line 475 of file grainvar.h.

References bin.

Referenced by AbundancesPrt(), ContRate(), ConvPresTempEdenIoniz(), CoolEvaluate(), fill_array(), GrainDrive(), GrainRestartIter(), GrainsInit(), GrainStartIter(), ion_photo(), lgNetEdenSrcSmall(), t_mole_global::make_species(), mole_h2_grain_form(), mole_update_species_cache(), OpacityAddTotal(), PrtComment(), PrtFinal(), and RT_diffuse().

void GrainVar::p_clear0 ( )

private

Definition at line 339 of file grains.cpp.

References AugerData, bin, dstab, dstsc, GrainEmission, GraphiteEmission, LIMELM, ReadRecord, and SilicateEmission.

Referenced by clear(), and ~GrainVar().

void GrainVar::p_clear1 ( )

private

Definition at line 359 of file grains.cpp.

References AugerData, bin, chPAH_abundance, dHeatdT, dstAbundThresholdFar, dstAbundThresholdNear, ENTH_CAR, ENTH_CAR2, ENTH_PAH, ENTH_PAH2, ENTH_SIC, ENTH_SIL, ENTH_SIL2, GrainChTrRate, GrainHeatScaleFactor, GrainMetal, H2_CAR, H2_SIL, IAL_CAR, IAL_SIL, lgAnyDustVary, lgBakesPAH_heat, lgDColOn, lgDHetOn, lgGrainElectrons, lgGrainPhysicsOn, lgNegGrnDrg, lgQHeatAll, lgQHeatOn, lgReevaluate, lgWD01, LIMELM, MAT_CAR, MAT_CAR2, MAT_PAH, MAT_PAH2, MAT_SIC, MAT_SIL, MAT_SIL2, nCalledGrainDrive, NCHRG_DEFAULT, nChrgRequested, PE_CAR, PE_SIL, POT_CAR, POT_SIL, STRG_CAR, STRG_SIL, TotalEden, which_enth, which_H2distr, which_ial, which_pe, which_pot, which_strg, which_zmin, ZMIN_CAR, and ZMIN_SIL.

Referenced by clear(), and GrainVar().

Field Documentation

AEInfo* GrainVar::AugerData[LIMELM]

store data about energy spectrum of Auger electrons

Definition at line 536 of file grainvar.h.

Referenced by GrainsInit(), InitBinAugerData(), p_clear0(), p_clear1(), and ReadAugerData().

vector<GrainBin*> GrainVar::bin

per bin grain data pointers to memory allocated for bins

Definition at line 583 of file grainvar.h.

Referenced by AbundancesPrt(), conorm(), elec_esc_length(), GetFracPop(), GetPotValues(), GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), GrainCharge(), GrainChargeTemp(), GrainChrgTransferRates(), GrainCollHeating(), GrainDrift(), GrainDrive(), GrainElecEmis1(), GrainElecRecomb1(), GrainIonColl(), GrainMakeDiffuse(), GrainRestartIter(), GrainScreen(), GrainsInit(), GrainStartIter(), GrainTemperature(), GrainUpdateRadius1(), GrainUpdateRadius2(), GrnStdDpth(), GrnVryDpth(), InitBinAugerData(), InitEmissivities(), InitEnthalpy(), inv_ufunct(), lgCheckMonitors(), lgDustOn(), lines_grains(), mie_read_opc(), diatomics::mole_H2_form(), mole_h2_grain_form(), mole_h_fixup(), mole_update_species_cache(), NewChargeData(), no_atoms(), one_elec(), p_clear0(), p_clear1(), PE_init(), PlanckIntegral(), PrtComment(), PrtFinal(), PrtZone(), qheat(), qheat_init(), radius_increment(), RebinQHeatResults(), SaveDo(), TryDoubleStep(), uderiv(), ufunct(), UpdatePot(), UpdatePot1(), UpdatePot2(), UpdateRecomZ0(), y0b(), y0b01(), y0psa(), y1psa(), and Yfunc().

string GrainVar::chPAH_abundance

which functions describes the PAH abundances, changed with SET PAH

Definition at line 502 of file grainvar.h.

Referenced by GrnStdDpth(), p_clear1(), and ParseSet().

realnum GrainVar::dclmax

largest local cooling of gas by collisions with grains

Definition at line 559 of file grainvar.h.

Referenced by GrainDrive(), GrainStartIter(), lines_grains(), and PrtComment().

double GrainVar::dHeatdT

gas heating derivative for all grains, in erg/cm^3/s/K

Definition at line 555 of file grainvar.h.

Referenced by GrainChargeTemp(), HeatSum(), and p_clear1().

realnum GrainVar::dphmax

largest fraction of local heating due to grain PE heating

Definition at line 559 of file grainvar.h.

Referenced by GrainDrive(), GrainStartIter(), lines_grains(), and PrtComment().

vector<double> GrainVar::dstab

total absorption cross section, current depl is factored in

Definition at line 525 of file grainvar.h.

Referenced by ContRate(), GrainsInit(), GrainUpdateRadius2(), OpacityAddTotal(), p_clear0(), and SaveDo().

realnum GrainVar::dstAbundThresholdFar

min grain abundance for enabling quantum heating in molecular regions

Definition at line 568 of file grainvar.h.

Referenced by GrainMakeDiffuse(), and p_clear1().

realnum GrainVar::dstAbundThresholdNear

quantum heating physics min grain abundance for enabling quantum heating near I-front

Definition at line 567 of file grainvar.h.

Referenced by GrainMakeDiffuse(), and p_clear1().

vector<double> GrainVar::dstsc

total scattering cross section, current depl and asymmetry factored in

Definition at line 526 of file grainvar.h.

Referenced by GrainsInit(), GrainUpdateRadius2(), OpacityAddTotal(), p_clear0(), and SaveDo().

double GrainVar::dsttmp[NDEMS]

equilibrium temperature grain temperature grid for dstems (in GrainBin)

Definition at line 564 of file grainvar.h.

Referenced by GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), GrainChargeTemp(), GrainTemperature(), InitEmissivities(), InitEnthalpy(), inv_ufunct(), qheat(), RebinQHeatResults(), TryDoubleStep(), and ufunct().

realnum GrainVar::elmSumAbund[LIMELM]

number density of elements summed over all grain bins, at actual depl

Definition at line 510 of file grainvar.h.

Referenced by AbundancesPrt(), GrainsInit(), GrainUpdateRadius1(), and ion_photo().

double GrainVar::GasCoolColl

heating/cooling balance, all entries are valid for current zone, actual depl, and are in erg/cm^3/s cooling of the gas by collisions with all grains, summed in GrGC 0

Definition at line 544 of file grainvar.h.

Referenced by CoolEvaluate(), CoolSave(), GrainChargeTemp(), GrainDrive(), GrainMakeDiffuse(), GrainsInit(), lines_grains(), and SaveDo().

double GrainVar::GasHeatNet

net heating of the gas: PhotoEl + Therm - CoolColl

Definition at line 544 of file grainvar.h.

Referenced by GrainChargeTemp().

double GrainVar::GasHeatPhotoEl

heating of the gas by photoelectric effect of all grains

Definition at line 544 of file grainvar.h.

Referenced by CoolEvaluate(), GrainChargeTemp(), and GrainDrive().

double GrainVar::GasHeatTherm

heating of the gas by thermionic emissions of all grains

Definition at line 544 of file grainvar.h.

Referenced by CoolEvaluate(), and GrainChargeTemp().

realnum GrainVar::GrainChTrRate[LIMELM][LIMELM+1][LIMELM+1]

grain surface recombination rate for element nelem, ionization stage "ion-from" to "ion-to" already multiplied by grain area, actual depletion, GrainChTrRate[nelem][ion-from][ion-to], units s^-1, this can be both ionization and neutralization process, so ion-to can be > or < ion-from

Definition at line 541 of file grainvar.h.

Referenced by fill_array(), GrainChargeTemp(), GrainChrgTransferRates(), GrainDrive(), ion_recomb(), iso_level(), lgNetEdenSrcSmall(), p_clear1(), and PrintRates().

vector<realnum> GrainVar::GrainEmission

grain emission total emission from this zone, per unit vol, for add outward in metfic

Definition at line 578 of file grainvar.h.

Referenced by ContRate(), flxCell(), GrainMakeDiffuse(), GrainsInit(), p_clear0(), and RT_diffuse().

double GrainVar::GrainHeatChem

net amount of energy donated by recombining ions to all grains

Definition at line 544 of file grainvar.h.

Referenced by GrainChargeTemp(), GrainDrive(), and GrainMakeDiffuse().

double GrainVar::GrainHeatCollSum

collisional heating of all grains, added in GraC 0

Definition at line 544 of file grainvar.h.

Referenced by GrainChargeTemp(), GrainDrive(), GrainsInit(), and lines_grains().

double GrainVar::GrainHeatDif

heating of all grains by all diffuse fields (incl Lya), added in GraD 0

Definition at line 544 of file grainvar.h.

Referenced by GrainChargeTemp(), GrainDrive(), GrainsInit(), and lines_grains().

double GrainVar::GrainHeatInc

heating of all grains by incident continuum, added in GraI 0

Definition at line 544 of file grainvar.h.

Referenced by GrainChargeTemp(), GrainDrive(), GrainsInit(), and lines_grains().

double GrainVar::GrainHeatLya

heating of all grains by Lya, added in GraL 1216

Definition at line 544 of file grainvar.h.

Referenced by GrainChargeTemp(), GrainDrive(), GrainsInit(), and lines_grains().

realnum GrainVar::GrainHeatScaleFactor

scale factor for PE heating as per Allers et al. 2005 (SET GRAINS HEAT)

Definition at line 557 of file grainvar.h.

Referenced by GrainDrive(), GrainMakeDiffuse(), GrainTemperature(), p_clear1(), and ParseSet().

double GrainVar::GrainHeatSum

total heating of all grain types, added in GraT 0

Definition at line 544 of file grainvar.h.

Referenced by GrainChargeTemp(), GrainDrive(), GrainMakeDiffuse(), GrainsInit(), and lines_grains().

realnum GrainVar::GrainMetal

grain depletion chemical abundances of grains are normalized such that in any given zone the total number density of an element nelem locked up in grains is given by elmSumAbund(nelem) = Sum_over_nd ( elmAbund[nd][nelem]*dstAbund[nd]*dense.hden ) grain depletion factor, from METALS xxx GRAINS command, usually 1

Definition at line 510 of file grainvar.h.

Referenced by GrainsInit(), GrainUpdateRadius1(), p_clear1(), and ParseMetal().

vector<realnum> GrainVar::GraphiteEmission

graphite emission from this zone only

Definition at line 579 of file grainvar.h.

Referenced by flxCell(), GrainMakeDiffuse(), GrainsInit(), and p_clear0().

realnum GrainVar::GrnElecDonateMax

largest local fraction of electrons donated by grains

Definition at line 531 of file grainvar.h.

Referenced by GrainDrive(), GrainStartIter(), lines_grains(), and PrtComment().

realnum GrainVar::GrnElecHoldMax

largest local fraction of electrons contained by grains

Definition at line 532 of file grainvar.h.

Referenced by GrainDrive(), GrainStartIter(), lines_grains(), and PrtComment().

double GrainVar::GrnRecomTe

the electron temperature that was used to calculate the recom rates

Definition at line 533 of file grainvar.h.

Referenced by GrainChargeTemp(), GrainsInit(), and NewChargeData().

bool GrainVar::lgAnyDustVary

do any of the grain abundances vary with depth

Definition at line 479 of file grainvar.h.

Referenced by GrainDrive(), GrainsInit(), and p_clear1().

bool GrainVar::lgBakesPAH_heat

turn on simple formula for PAH heating of gas

Definition at line 479 of file grainvar.h.

Referenced by CoolEvaluate(), GrainDrive(), GrainMakeDiffuse(), GrainTemperature(), p_clear1(), and ParseSet().

bool GrainVar::lgDColOn

turn off grain gas collisional energy exchange, set false with NO GRAIN COLLISIONAL ENERGY EXCHANGE default true, turned off with GRAIN NO COOLING

Definition at line 494 of file grainvar.h.

Referenced by GrainCollHeating(), GrainMakeDiffuse(), p_clear1(), ParseCommands(), ParseDont(), and ParseGrain().

bool GrainVar::lgDHetOn

test logic default true, turned off with GRAIN NO HEATING

Definition at line 489 of file grainvar.h.

Referenced by CoolEvaluate(), GrainMakeDiffuse(), GrainTemperature(), p_clear1(), and ParseGrain().

bool GrainVar::lgGrainElectrons

should electrons from/to grains be included in the total electron sum? del true, set false with no grain electrons command

Definition at line 498 of file grainvar.h.

Referenced by ConvBase(), eden_sum(), p_clear1(), and ParseDont().

bool GrainVar::lgGrainPhysicsOn

this is option to turn off all grain physics while leaving the opacity in, set false with no grain physics command

Definition at line 479 of file grainvar.h.

Referenced by ContRate(), ConvPresTempEdenIoniz(), CoolEvaluate(), fill_array(), GrainDrive(), GrainRestartIter(), GrainStartIter(), lines_grains(), mole_h2_grain_form(), p_clear1(), ParseDont(), PrtFinal(), PrtZone(), and RT_diffuse().

bool GrainVar::lgNegGrnDrg

flag set if negative grin drag force encountered

Definition at line 486 of file grainvar.h.

Referenced by GrainDrift(), GrainStartIter(), p_clear1(), and PrtComment().

bool GrainVar::lgQHeatAll

use quantum heating for all grains, WILL BECOME DEFAULT

Definition at line 569 of file grainvar.h.

Referenced by GrainsInit(), and p_clear1().

bool GrainVar::lgQHeatOn

default true, turned off with GRAIN NO QHEAT

Definition at line 489 of file grainvar.h.

Referenced by GrainsInit(), p_clear1(), and ParseDont().

bool GrainVar::lgQHPunLast

only save quantum heating information on last iteration

Definition at line 570 of file grainvar.h.

Referenced by qheat().

bool GrainVar::lgReevaluate

reevaluate keyword on grains cmd => always reevalutae grain conditions during calculation

Definition at line 479 of file grainvar.h.

Referenced by GrainDrive(), p_clear1(), and ParseGrain().

bool GrainVar::lgWD01

use physics as described in Weingartner & Draine (2001)

Definition at line 479 of file grainvar.h.

Referenced by GrainsInit(), InitBinAugerData(), ion_photo(), p_clear1(), ParseDont(), qheat_init(), y0b(), and y2s().

long GrainVar::nChrgRequested

number of charge states requested by user, default value is 2

Definition at line 534 of file grainvar.h.

Referenced by GrainBin::p_clear1(), p_clear1(), and SetNChrgStates().

long GrainVar::nzone

grain opacities remember in what zone grain quantities were last updated

Definition at line 524 of file grainvar.h.

Referenced by GrainChargeTemp(), GrainsInit(), and NewChargeData().

FILE* GrainVar::QHSaveFile

file pointer for PUNCH QHEAT command

Definition at line 571 of file grainvar.h.

Referenced by CloseSaveFiles(), qheat(), and SaveFilesInit().

double GrainVar::rate_h2_form_grains_used_total

H2 physics rate H2 forms on grains, summed over bins, units s^-1, actual depl when multiplied with hden, this is formation rate in H2-molecules/cm^3/s

Definition at line 574 of file grainvar.h.

Referenced by diatomics::H2_PunchDo(), lgCheckMonitors(), diatomics::mole_H2_form(), mole_h2_grain_form(), and SaveDo().

vector<string> GrainVar::ReadRecord

record of all the files read by mie_read_opc

Definition at line 500 of file grainvar.h.

Referenced by mie_read_opc(), and p_clear0().

vector<realnum> GrainVar::SilicateEmission

silicate emission from this zone only

Definition at line 580 of file grainvar.h.

Referenced by flxCell(), GrainMakeDiffuse(), GrainsInit(), and p_clear0().

realnum GrainVar::TotalDustHeat

total PE heating integrated over model, erg/s or erg/cm^2/s, actual depl

Definition at line 559 of file grainvar.h.

Referenced by GrainDrive(), GrainStartIter(), lines_grains(), and PrtComment().

double GrainVar::TotalEden

grain charging contribution to eden from all grain species, a positive number means that the grains contribute to the free electron pool, in cm^-3

Definition at line 529 of file grainvar.h.

Referenced by ConvBase(), eden_sum(), GrainChargeTemp(), GrainDrive(), lines_grains(), p_clear1(), and SaveDo().

enth_type GrainVar::which_enth[MAT_TOP]

these arrays define the physical properties of material MAT_PAH, MAT_CAR, MAT_SIL, etc... they are initialized in the constructor defines expression for the enthalpy function

Definition at line 515 of file grainvar.h.

Referenced by p_clear1(), and uderiv().

H2_type GrainVar::which_H2distr[MAT_TOP]

defines expression for H2 ro-vib distribution at formation

Definition at line 521 of file grainvar.h.

Referenced by GrainCollHeating(), diatomics::mole_H2_form(), and p_clear1().

ial_type GrainVar::which_ial[MAT_TOP]

defines expression for the inverse attenuation length

Definition at line 518 of file grainvar.h.

Referenced by mie_write_opc(), and p_clear1().

pe_type GrainVar::which_pe[MAT_TOP]

defines expression for the photo-electric effect