t_dynamics Struct Reference

#include <dynamics.h>

Collaboration diagram for t_dynamics:


Data Fields
bool	lgAdvection
double	Cool
double	Heat
double	dCooldT
double	dHeatdT
double	CoolMax
double	HeatMax
double	Rate
double **	Source
double *	H2_molec
double *	CO_molec
bool	lgISO [NISO]
bool	lgMETALS
bool	lgCoolHeat
bool	lg_coronal_time_init
bool	lgStatic
double	time_elapsed
bool	lgRecom
bool	lgStatic_completed
double	AdvecLengthInit
double	FluxCenter
char	chPresMode [20]
double	ShockDepth
double	ShockMach
double	FluxScale
bool	lgFluxDScale
double	FluxIndex
double	dRad
double	oldFullDepth
double	convergence_error
double	convergence_tolerance
double	discretization_error
double	error_scale1
double	error_scale2
bool	lgSetPresMode
float	dDensityDT
long	iteration_first_dynamics

Detailed Description

all of these are initialized in zero

Definition at line 54 of file dynamics.h.

Field Documentation

bool t_dynamics::lgAdvection

is advection turned on ?, set to false in zero

Definition at line 57 of file dynamics.h.

Referenced by advection_set_detault(), CO_solve(), ConvBase(), ConvIterCheck(), DynaZero(), HeLikeLevel(), hmole_step(), HydroLevel(), HydroLevelPop(), ion_solver(), IterRestart(), PrtZone(), radius_increment(), radius_next(), ZoneEnd(), and ZoneStart().

double t_dynamics::Cool

advective cooling minus heating (erg cm^-3 s^-1)

Definition at line 60 of file dynamics.h.

Referenced by CoolEvaluate(), CoolPunch(), DynaEndIter(), DynaIonize(), DynaPrtZone(), DynaPunch(), DynaZero(), HeatPunch(), lines_general(), and PunchDo().

double t_dynamics::Heat

Definition at line 60 of file dynamics.h.

Referenced by CoolEvaluate(), CoolPunch(), DynaIonize(), DynaPrtZone(), DynaPunch(), DynaZero(), HeatPunch(), lines_general(), and PunchDo().

double t_dynamics::dCooldT

Definition at line 60 of file dynamics.h.

Referenced by CoolEvaluate(), DynaIonize(), DynaPunch(), and DynaZero().

double t_dynamics::dHeatdT

Definition at line 60 of file dynamics.h.

Referenced by CoolEvaluate(), DynaIonize(), and DynaZero().

double t_dynamics::CoolMax

largest fraction of cooling and heating

Definition at line 63 of file dynamics.h.

Referenced by DynaZero(), IterStart(), lines_general(), and PrtComment().

double t_dynamics::HeatMax

Definition at line 63 of file dynamics.h.

Referenced by DynaZero(), IterStart(), lines_general(), and PrtComment().

double t_dynamics::Rate

the advection rate (s^-1)

Definition at line 66 of file dynamics.h.

Referenced by CO_solve(), DynaCreateArrays(), DynaIonize(), DynaPunch(), DynaStartZone(), DynaZero(), HeLikeLevel(), hmole_step(), HydroLevelPop(), ion_solver(), and PunchDo().

double** t_dynamics::Source

the advective recombination rate (s^-1)

Definition at line 69 of file dynamics.h.

Referenced by DynaCreateArrays(), DynaIonize(), DynaPunch(), DynaStartZone(), HeLikeLevel(), HydroLevelPop(), ion_solver(), and PunchDo().

double* t_dynamics::H2_molec

save H2 and CO densities

Definition at line 72 of file dynamics.h.

Referenced by DynaCreateArrays(), DynaIonize(), and hmole_step().

double* t_dynamics::CO_molec

Definition at line 73 of file dynamics.h.

Referenced by CO_solve(), DynaCreateArrays(), DynaIonize(), and hmole_step().

bool t_dynamics::lgISO[NISO]

factor to turn off advection for H-like or He-like iso seq, no advection h-like, he-like

Definition at line 76 of file dynamics.h.

Referenced by DynaZero(), HeLikeLevel(), HydroLevelPop(), ParseDont(), and PresTotCurrent().

bool t_dynamics::lgMETALS

factor to turn off advection for rest of ions, "no advection metals"

Definition at line 79 of file dynamics.h.

Referenced by DynaZero(), ParseDont(), and PresTotCurrent().

bool t_dynamics::lgCoolHeat

factor to turn off advective cooling

Definition at line 82 of file dynamics.h.

Referenced by DynaIonize(), DynaZero(), and ParseDont().

bool t_dynamics::lg_coronal_time_init

var set with coronal time init - says to use constant temperature on first relax iterations then let temp run free

Definition at line 86 of file dynamics.h.

Referenced by DynaEndIter(), and ParseCoronal().

bool t_dynamics::lgStatic

set true if static time dependent

Definition at line 89 of file dynamics.h.

Referenced by DynaEndIter(), DynaIonize(), DynaPresChngFactor(), DynaZero(), lgConvPres(), ParseDynaTime(), PresTotCurrent(), PrtComment(), PunchDo(), radius_first(), and radius_next().

double t_dynamics::time_elapsed

elapsed time in time dependent static model

Definition at line 92 of file dynamics.h.

Referenced by DynaEndIter(), DynaPunchTimeDep(), DynaZero(), and PunchDo().

bool t_dynamics::lgRecom

true if recombination logic in place

Definition at line 95 of file dynamics.h.

Referenced by DynaEndIter(), DynaZero(), PrtComment(), radius_next(), and timestep_next().

bool t_dynamics::lgStatic_completed

true if model ends since time dependent model is finished

Definition at line 98 of file dynamics.h.

Referenced by cloudy(), DynaEndIter(), and DynaZero().

double t_dynamics::AdvecLengthInit

the initial value of the advection length, reset with set advection length

Definition at line 101 of file dynamics.h.

Referenced by DynaEndIter(), DynaZero(), and ParseSet().

double t_dynamics::FluxCenter

the center of the particle flux law

Definition at line 104 of file dynamics.h.

Referenced by DynaFlux(), DynaZero(), and ParseDynaWind().

char t_dynamics::chPresMode[20]

flag set by the "set dynamics pressure mode" command

Definition at line 107 of file dynamics.h.

Referenced by DynaPresChngFactor(), and ParseSet().

double t_dynamics::ShockDepth

the shock depth in cm set with "set dynamics shock depth" command

Definition at line 110 of file dynamics.h.

Referenced by DynaPresChngFactor(), and ParseSet().

double t_dynamics::ShockMach

the isothermal Mach number at which to insert an antishock set with "set dynamics antishock mach" command

Definition at line 114 of file dynamics.h.

Referenced by DynaPresChngFactor(), and ParseSet().

double t_dynamics::FluxScale

the scale of the particle flux law

Definition at line 117 of file dynamics.h.

Referenced by DynaFlux(), DynaZero(), and ParseDynaWind().

bool t_dynamics::lgFluxDScale

whether we also need to scale by the face density

Definition at line 120 of file dynamics.h.

Referenced by DynaFlux(), DynaZero(), and ParseDynaWind().

double t_dynamics::FluxIndex

the power law index of the particle flux law

Definition at line 123 of file dynamics.h.

Referenced by DynaFlux(), DynaZero(), and ParseDynaWind().

double t_dynamics::dRad

the proposed thickness for the next zone when advection is included

Definition at line 126 of file dynamics.h.

Referenced by DynaStartZone(), DynaZero(), and radius_next().

double t_dynamics::oldFullDepth

the depth of the last iteration

Definition at line 129 of file dynamics.h.

Referenced by DynaCreateArrays(), DynaIonize(), DynaSaveLast(), and ion_solver().

double t_dynamics::convergence_error

convergence_error and discretization_error give estimates of convergence: :: convergence_error -- change between the last iterations; :: discretization_error -- error in the upstream interpolation. When (and if) discretization_error >> convergence_error, the interpolation length should be decreased.

They should both be based on the same norm of the models, but what norm may be experimented with -- at present, it's H+/Htot just weighted by cell number, which makes the estimates sensitive to the structure of the primary ionization front.

Definition at line 141 of file dynamics.h.

Referenced by ConvIterCheck(), and DynaNewStep().

double t_dynamics::convergence_tolerance

the allowed rel error, by default 0.1

Definition at line 144 of file dynamics.h.

Referenced by ConvIterCheck(), DynaNewStep(), and DynaZero().

double t_dynamics::discretization_error

the error to be expected based on the coarseness of current advection length

Definition at line 147 of file dynamics.h.

Referenced by ConvIterCheck(), DynaEndIter(), DynaNewStep(), DynaStartZone(), and DynaZero().

double t_dynamics::error_scale1

two ways of scaling the error estimate for convergence