#include <dynamics.h>
Public Member Functions | |
double | Cool () |
double | Heat () |
double | dCooldT () |
Data Fields | |
bool | lgAdvection |
double | Cool_r |
double | Heat_v |
double | dHeatdT |
double | CoolMax |
double | HeatMax |
double | Rate |
double ** | Source |
double *** | StatesElem |
double * | molecules |
bool | lgISO [NISO] |
bool | lgMETALS |
bool | lgCoolHeat |
bool | lg_coronal_time_init |
bool | lgTimeDependentStatic |
double | time_elapsed |
bool | lgRecom |
bool | lgStatic_completed |
double | AdvecLengthInit |
double | FluxCenter |
char | chPresMode [20] |
double | ShockDepth |
double | ShockMach |
long int | n_initial_relax |
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 |
realnum | Upstream_density |
realnum | DivergePresInteg |
bool | lgEquilibrium |
bool | lgTracePrint |
double | timestep_init |
double | timestep |
double | timestep_stop |
double | timestep_factor |
all of these are initialized in zero
Definition at line 57 of file dynamics.h.
double t_dynamics::Cool | ( | ) |
Definition at line 2182 of file dynamics.cpp.
References Cool_r, t_phycon::EnthalpyDensity, Heat_v, MERGE, phycon, and scalingDensity().
Referenced by CoolEvaluate(), CoolSave(), DynaIterEnd(), DynaPrtZone(), DynaSave(), lines_general(), SaveDo(), and SaveHeat().
double t_dynamics::dCooldT | ( | ) |
Definition at line 2197 of file dynamics.cpp.
References Cool_r, phycon, t_pressure::PresGasCurr, pressure, and t_phycon::te.
Referenced by CoolEvaluate(), and DynaSave().
double t_dynamics::Heat | ( | ) |
Definition at line 2168 of file dynamics.cpp.
References Cool_r, t_phycon::EnthalpyDensity, Heat_v, MERGE, phycon, and scalingDensity().
Referenced by CoolEvaluate(), CoolSave(), DynaPrtZone(), DynaSave(), lines_general(), SaveDo(), and SaveHeat().
double t_dynamics::AdvecLengthInit |
the initial value of the advection length, reset with set advection length
Definition at line 108 of file dynamics.h.
Referenced by DynaIterEnd(), DynaZero(), and ParseSet().
char t_dynamics::chPresMode[20] |
flag set by the "set dynamics pressure mode" command
Definition at line 114 of file dynamics.h.
Referenced by ParseSet(), PresMode::set(), and stepDensity().
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 153 of file dynamics.h.
Referenced by ConvIterCheck(), and DynaNewStep().
the allowed rel error, by default 0.1
Definition at line 156 of file dynamics.h.
Referenced by ConvIterCheck(), DynaNewStep(), and DynaZero().
double t_dynamics::Cool_r |
advective cooling minus heating (erg cm^-3 s^-1)
Definition at line 63 of file dynamics.h.
Referenced by Cool(), dCooldT(), DynaIonize(), DynaZero(), and Heat().
double t_dynamics::CoolMax |
largest fraction of cooling and heating
Definition at line 68 of file dynamics.h.
Referenced by DynaZero(), IterStart(), lines_general(), and PrtComment().
double t_dynamics::dHeatdT |
Definition at line 63 of file dynamics.h.
Referenced by CoolEvaluate(), DynaIonize(), and DynaZero().
the error to be expected based on the coarseness of current advection length
Definition at line 159 of file dynamics.h.
Referenced by ConvIterCheck(), DynaIterEnd(), DynaNewStep(), DynaStartZone(), and DynaZero().
Definition at line 171 of file dynamics.h.
Referenced by DynaEndZone(), DynaIterEnd(), DynaZero(), pressureZone(), and stepDensity().
double t_dynamics::dRad |
the proposed thickness for the next zone when advection is included
Definition at line 138 of file dynamics.h.
Referenced by DynaStartZone(), DynaZero(), and radius_next().
double t_dynamics::error_scale1 |
two ways of scaling the error estimate for convergence
Definition at line 162 of file dynamics.h.
Referenced by DynaNewStep().
double t_dynamics::error_scale2 |
Definition at line 162 of file dynamics.h.
Referenced by ConvIterCheck(), DynaIterEnd(), DynaNewStep(), DynaStartZone(), and DynaZero().
double t_dynamics::FluxCenter |
the center of the particle flux law
Definition at line 111 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 135 of file dynamics.h.
Referenced by DynaFlux(), DynaZero(), and ParseDynaWind().
double t_dynamics::FluxScale |
the scale of the particle flux law
Definition at line 129 of file dynamics.h.
Referenced by DynaFlux(), DynaZero(), and ParseDynaWind().
double t_dynamics::Heat_v |
Definition at line 63 of file dynamics.h.
Referenced by Cool(), DynaIonize(), DynaZero(), and Heat().
double t_dynamics::HeatMax |
Definition at line 68 of file dynamics.h.
Referenced by DynaZero(), IterStart(), lines_general(), and PrtComment().
var set with coronal time init - says to use constant temperature on first relax iterations then let temp run free
Definition at line 93 of file dynamics.h.
Referenced by DynaIterEnd(), InitDefaultsPreparse(), and ParseCoronal().
is advection turned on ?, set to false in zero
Definition at line 60 of file dynamics.h.
Referenced by AbundChange(), advection_set_default(), atom_levelN(), ConvBase(), ConvIterCheck(), DynaZero(), funjac(), HomogeneousSource(), iso_level(), IterRestart(), IterStart(), mole_eval_dynamic_balance(), PrintRates(), PrtZone(), radius_next(), ZoneEnd(), and ZoneStart().
factor to turn off advective cooling
Definition at line 89 of file dynamics.h.
Referenced by DynaIonize(), DynaZero(), and ParseDont().
Enforce equilibrium populations
Definition at line 174 of file dynamics.h.
Referenced by DynaZero(), HomogeneousSource(), iso_level(), and ParseSet().
whether we also need to scale by the face density
Definition at line 132 of file dynamics.h.
Referenced by DynaFlux(), DynaZero(), and ParseDynaWind().
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 83 of file dynamics.h.
Referenced by DynaZero(), iso_level(), ParseDont(), and PresTotCurrent().
bool t_dynamics::lgMETALS |
factor to turn off advection for rest of ions, "no advection metals"
Definition at line 86 of file dynamics.h.
Referenced by DynaZero(), ParseDont(), and PresTotCurrent().
bool t_dynamics::lgRecom |
true if recombination logic in place
Definition at line 102 of file dynamics.h.
Referenced by DynaIterEnd(), DynaIterStart(), DynaZero(), PrtComment(), radius_next(), and timestep_next().
flag set true if set dynamics flow type was set - this means to use the specified option, and not to derive one
Definition at line 166 of file dynamics.h.
Referenced by DynaZero(), ParseSet(), and PresMode::set().
true if model ends since time dependent model is finished
Definition at line 105 of file dynamics.h.
Referenced by cloudy(), DynaIterEnd(), DynaZero(), and SaveDo().
set true if time dependent static simulation
Definition at line 96 of file dynamics.h.
Referenced by check_grid_file(), ConvInitSolution(), CoolEvaluate(), DynaIonize(), DynaIterEnd(), DynaZero(), IterEnd(), lgConserveEnergy(), ParseDynaTime(), PresTotCurrent(), PrtComment(), radius_first(), radius_next(), and SaveDo().
Definition at line 177 of file dynamics.h.
Referenced by DynaEndZone(), DynaIonize(), DynaIterEnd(), DynaStartZone(), InitDefaultsPreparse(), ParseCoronal(), ParseDynaTime(), ParseDynaWind(), and stepDensity().
double* t_dynamics::molecules |
save molecular network densities
Definition at line 80 of file dynamics.h.
Referenced by DynaCreateArrays(), DynaIonize(), and mole_eval_dynamic_balance().
long int t_dynamics::n_initial_relax |
set how many iterations we will start with, before allowing changes. This allows the solution to relax to an equilibrium set with "set dynamics relax" command
Definition at line 126 of file dynamics.h.
Referenced by atom_levelN(), ConvInitSolution(), CoolEvaluate(), DynaIonize(), DynaIterEnd(), DynaZero(), funjac(), HomogeneousSource(), iso_level(), IterEnd(), mole_eval_dynamic_balance(), ParseSet(), PrtFinal(), and radius_next().
double t_dynamics::oldFullDepth |
the depth of the last iteration
Definition at line 141 of file dynamics.h.
Referenced by DynaCreateArrays(), DynaIonize(), and DynaSaveLast().
double t_dynamics::Rate |
the advection rate (s^-1)
Definition at line 71 of file dynamics.h.
Referenced by DynaCreateArrays(), DynaIonize(), DynaSave(), DynaStartZone(), DynaZero(), funjac(), HomogeneousSource(), iso_level(), mole_eval_dynamic_balance(), PrintRates(), and SaveDo().
double t_dynamics::ShockDepth |
the shock depth in cm set with "set dynamics shock depth" command
Definition at line 117 of file dynamics.h.
Referenced by ParseSet(), and PresMode::set().
double t_dynamics::ShockMach |
the isothermal Mach number at which to insert an antishock set with "set dynamics antishock Mach" command
Definition at line 121 of file dynamics.h.
Referenced by ParseSet(), and PresMode::set().
double** t_dynamics::Source |
the advective recombination rate (cm^-3 s^-1)
Definition at line 74 of file dynamics.h.
Referenced by DynaCreateArrays(), DynaIonize(), DynaSave(), DynaStartZone(), HomogeneousSource(), iso_level(), mole_eval_dynamic_balance(), PrintRates(), and SaveDo().
double*** t_dynamics::StatesElem |
the advective isolevel balance terms
Definition at line 77 of file dynamics.h.
Referenced by DynaCreateArrays(), DynaIonize(), and iso_level().
double t_dynamics::time_elapsed |
elapsed time in time dependent static model
Definition at line 99 of file dynamics.h.
Referenced by ConvInitSolution(), DynaIterEnd(), DynaIterStart(), DynaPunchTimeDep(), DynaZero(), and SaveDo().
double t_dynamics::timestep |
Definition at line 181 of file dynamics.h.
Referenced by atom_levelN(), CoolEvaluate(), DynaIonize(), DynaIterEnd(), DynaIterStart(), DynaPunchTimeDep(), IterEnd(), ParseDynaTime(), and timestep_next().
double t_dynamics::timestep_factor |
Definition at line 181 of file dynamics.h.
Referenced by DynaIterEnd(), and DynaZero().
double t_dynamics::timestep_init |
Definition at line 181 of file dynamics.h.
Referenced by DynaZero(), ParseDynaTime(), and timestep_next().
double t_dynamics::timestep_stop |
Definition at line 181 of file dynamics.h.
Referenced by DynaIterEnd(), and ParseDynaTime().
Upstream density
Definition at line 169 of file dynamics.h.
Referenced by CoolEvaluate(), DynaStartZone(), and zero().