#include "module.h"

Include dependency graph for dynamics.h:

This graph shows which files directly or indirectly include this file:

Data Structures
struct	t_dynamics

Functions
void	DynaIterStart (void)

void	DynaIterEnd (void)

void	DynaStartZone (void)

void	DynaEndZone (void)

void	DynaIonize (void)

void	DynaCreateArrays (void)

void	ParseDynaWind (Parser &p)

void	ParseDynaTime (Parser &p)

void	DynaPrtZone (void)

void	DynaSave (FILE *ipPnunit, char chJob)

void	DynaPunchTimeDep (FILE ipPnunit, const char chJob)

realnum	DynaFlux (double depth)

Variables
t_dynamics	dynamics

Function Documentation

void DynaCreateArrays ( void )

DynaCreateArrays allocate some space needed to save the dynamics structure variables, called from atmdat_readin

Definition at line 1398 of file dynamics.cpp.

References DEBUG_ENTRY, dense, EnthalpyDensity, fixit, ipH_LIKE, ipHYDROGEN, iso_sp, t_dense::lgElmtOn, LIMELM, MALLOC, mole_global, t_dynamics::molecules, NISO, t_mole_global::num_calc, t_iso_sp::numLevels_max, t_struc::nzlim, Old_DenMass, Old_density, Old_depth, Old_ednstr, Old_EnthalpyDensity, Old_gas_phase, Old_hiistr, Old_histr, Old_molecules, Old_pressure, Old_StatesElem, Old_xIonDense, Old_xLyman_depth, t_dynamics::oldFullDepth, t_dynamics::Rate, t_dynamics::Source, t_dynamics::StatesElem, struc, Upstream_molecules, UpstreamElem, UpstreamIon, and UpstreamStatesElem.

Referenced by atmdat_readin().

void DynaEndZone ( void )

DynaEndZone called at end of iteration when advection is turned on

Definition at line 833 of file dynamics.cpp.

References DEBUG_ENTRY, t_radius::depth, t_dynamics::DivergePresInteg, t_radius::drad, DynaFlux(), fprintf(), ioQQQ, t_dynamics::lgTracePrint, radius, wind, and Wind::windv.

Referenced by ZoneEnd().

Here is the call graph for this function:

realnum DynaFlux ( double depth )

Definition at line 1290 of file dynamics.cpp.

References DEBUG_ENTRY, dense, t_dynamics::FluxCenter, t_dynamics::FluxIndex, t_dynamics::FluxScale, t_dynamics::lgFluxDScale, pow(), and t_dense::xMassDensity0.

Referenced by DynaEndZone(), PresTotCurrent(), and stepDensity().

Here is the call graph for this function:

void DynaIonize ( void )

DynaIonize, called from ionize to evaluate advective terms for current conditions

Definition at line 160 of file dynamics.cpp.

References AdvecSpecificEnthalpy, ASSERT, t_dynamics::Cool_r, DEBUG_ENTRY, dense, t_radius::depth, t_dynamics::dHeatdT, Dyn_dr, t_phycon::EnthalpyDensity, EnthalpyDensity, fprintf(), t_dense::gas_phase, t_dynamics::Heat_v, t_thermal::htot, t_dense::IonHigh, t_dense::IonLow, ioQQQ, ipCARBON, ipH_LIKE, ipHYDROGEN, ipUpstream, iso_sp, iteration, t_dynamics::lgCoolHeat, t_dense::lgElmtOn, t_dynamics::lgTimeDependentStatic, t_dynamics::lgTracePrint, LIMELM, mole_global, t_dynamics::molecules, t_dynamics::n_initial_relax, NISO, t_mole_global::num_calc, t_iso_sp::numLevels_local, t_struc::nzlim, nzone, Old_depth, Old_xIonDense, t_dynamics::oldFullDepth, phycon, POW2, t_pressure::PresGasCurr, pressure, radius, t_dynamics::Rate, scalingDensity(), t_dynamics::Source, t_dynamics::StatesElem, struc, t_phycon::te, thermal, t_dynamics::timestep, Upstream_molecules, UpstreamElem, UpstreamIon, UpstreamStatesElem, wind, Wind::windv, t_dense::xIonDense, and t_dense::xMassDensity.

Referenced by ConvBase().

Here is the call graph for this function:

void DynaIterEnd ( void )

routine called at end of iteration when advection is turned on

Definition at line 854 of file dynamics.cpp.

References t_dynamics::AdvecLengthInit, ASSERT, t_thermal::ConstTemp, t_dynamics::Cool(), cosmology, t_thermal::ctot, DEBUG_ENTRY, t_radius::depth, t_dynamics::discretization_error, t_dynamics::DivergePresInteg, t_radius::dr_max_last_iter, Dyn_dr, DynaNewStep(), DynaSaveLast(), t_dynamics::error_scale2, fprintf(), GetHubbleFactor(), InitDynaTimestep(), ioQQQ, iphUpstream, ipUpstream, ipyUpstream, t_iterations::iter_malloc, iteration, iterations, t_dynamics::lg_coronal_time_init, t_cosmology::lgDo, lgNeedTimestep(), lgPrintDynamics, t_dynamics::lgRecom, t_radius::lgSdrmaxRel, Wind::lgStatic(), t_dynamics::lgStatic_completed, t_thermal::lgTemperatureConstant, lgtime_dt_specified, lgtime_Recom, t_dynamics::lgTimeDependentStatic, t_dynamics::lgTracePrint, t_dynamics::n_initial_relax, nTime_flux, phycon, t_pressure::PresGasCurr, pressure, radius, t_cosmology::redshift_current, t_cosmology::redshift_step, rfield, scalingDensity(), t_radius::sdrmax, StopCalc, t_iterations::StopThickness, t_phycon::te, t_StopCalc::TempHiStopIteration, t_StopCalc::TempLoStopIteration, thermal, t_rfield::time_continuum_scale, time_dt, time_dt_scale_factor, t_dynamics::time_elapsed, time_elapsed_time, t_dynamics::timestep, t_dynamics::timestep_factor, timestep_next(), t_dynamics::timestep_stop, wind, and Wind::windv0.

Referenced by IterRestart().

Here is the call graph for this function:

void DynaIterStart ( void )

routine called at start of iteration when advection is turned on

Definition at line 2225 of file dynamics.cpp.

References cdEXIT, DEBUG_ENTRY, EXIT_FAILURE, fprintf(), hextra, ioQQQ, iteration, lgPrintDynamics, t_dynamics::lgRecom, t_dynamics::lgStatic_completed, t_rfield::lgTimeVary, t_hextra::lgTurbHeatVaryTime, linint(), t_dynamics::n_initial_relax, t_rfield::nShape, nTime_flux, rfield, t_rfield::time_continuum_scale, t_dynamics::time_elapsed, time_elapsed_time, time_flux_ratio, t_dynamics::timestep, t_hextra::TurbHeat, and t_hextra::TurbHeatSave.

Referenced by IterStart().

Here is the call graph for this function:

void DynaPrtZone ( void )

DynaPrtZone - called to print zone results

Definition at line 2027 of file dynamics.cpp.

References AdvecSpecificEnthalpy, ASSERT, t_dynamics::Cool(), t_thermal::ctot, DEBUG_ENTRY, dense, t_phycon::EnergyBinding, t_phycon::EnergyExcitation, t_phycon::EnergyIonization, EnthalpyDensity, fprintf(), t_dynamics::Heat(), ioQQQ, t_struc::nzlim, nzone, phycon, POW2, t_pressure::PresGasCurr, pressure, scalingDensity(), t_timesc::sound_speed_adiabatic, struc, thermal, timesc, wind, Wind::windv, and t_dense::xMassDensity.

Referenced by PrtZone().

Here is the call graph for this function:

void DynaPunchTimeDep	(	FILE *	ipPnunit,
		const char *	chJob
	)

DynaPunchTimeDep - save info about time dependent solution

Parameters

ipPnunit
*chJob

Definition at line 2057 of file dynamics.cpp.

References cdIonFrac(), cdTemp(), colden, t_colden::colden, column(), cosmology, DEBUG_ENTRY, dense, t_dense::eden, findspecieslocal(), fprintf(), t_dense::gas_phase, ipCOL_HTOT, ipHYDROGEN, t_cosmology::redshift_current, rfield, scalingDensity(), SDIV(), t_rfield::time_continuum_scale, t_dynamics::time_elapsed, t_dynamics::timestep, and TotalInsanity().

Referenced by SaveDo().

Here is the call graph for this function:

void DynaSave	(	FILE *	ipPnunit,
		char	chJob
	)

DynaSave save info related to advection

Parameters

ipPnunit
chJob

Definition at line 2166 of file dynamics.cpp.

References AdvecSpecificEnthalpy, t_dynamics::Cool(), t_dynamics::dCooldT(), DEBUG_ENTRY, t_radius::depth_mid_zone, t_phycon::EnthalpyDensity, fprintf(), t_dynamics::Heat(), t_thermal::htot, ipHYDROGEN, phycon, radius, t_dynamics::Rate, scalingDensity(), t_dynamics::Source, thermal, and TotalInsanity().

Referenced by SaveDo().

Here is the call graph for this function:

void DynaStartZone ( void )

DynaStartZone called at start of iteration when advection is turned on

Definition at line 375 of file dynamics.cpp.

References AdvecSpecificEnthalpy, ASSERT, BIGFLOAT, DEBUG_ENTRY, dense, t_radius::depth, t_dynamics::discretization_error, t_dynamics::dRad, Dyn_dr, t_dynamics::error_scale2, fp_bound(), fprintf(), ioQQQ, ipH_LIKE, iphUpstream, ipHYDROGEN, ipUpstream, ipyUpstream, isnan, iso_sp, iteration, t_dense::lgElmtOn, t_dynamics::lgTracePrint, LIMELM, t_mole_global::list, MAX2, MIN2, mole, mole_global, NISO, nOld_zone, t_mole_global::num_calc, t_iso_sp::numLevels_max, nzone, Old_density, Old_depth, Old_EnthalpyDensity, Old_molecules, Old_StatesElem, Old_xIonDense, POW2, radius, t_dynamics::Rate, scalingDensity(), SMALLFLOAT, t_dynamics::Source, t_mole_local::species, t_dynamics::Upstream_density, Upstream_molecules, UpstreamElem, UpstreamIon, and UpstreamStatesElem.

Referenced by ZoneStart().

Here is the call graph for this function:

void ParseDynaTime ( Parser & p )

ParseDynaTime parse the time command, called from ParseCommands, in dynamics.c

Parameters

*chCard

Definition at line 1669 of file dynamics.cpp.

References advection_set_default(), ASSERT, cdEXIT, DEBUG_ENTRY, EXIT_FAILURE, fprintf(), Parser::getline(), Parser::getNumberCheckAlwaysLog(), Parser::getNumberCheckAlwaysLogLim(), Parser::getNumberDefaultAlwaysLog(), Parser::hasCommand(), ioQQQ, lgPrintDynamics, lgtime_dt_specified, lgtime_Recom, t_dynamics::lgTimeDependentStatic, t_dynamics::lgTracePrint, Parser::m_lgEOF, Parser::nMatch(), NTIME, nTime_flux, Wind::setStatic(), time_dt, time_dt_scale_factor, time_elapsed_time, time_flux_ratio, t_dynamics::timestep, t_dynamics::timestep_init, t_dynamics::timestep_stop, wind, Wind::windv, and Wind::windv0.

Referenced by ParseCommands().

Here is the call graph for this function:

void ParseDynaWind ( Parser & p )

Todo:: 2 Need to include divergence terms in pressure balance if flux index is != 0

Todo:: 2 Need to include divergence terms in pressure balance if flux index is != 0

Definition at line 1816 of file dynamics.cpp.

References advection_set_default(), BIGDOUBLE, cdEXIT, t_dense::chDenseLaw, Wind::comass, DEBUG_ENTRY, dense, EXIT_FAILURE, t_dynamics::FluxCenter, t_dynamics::FluxIndex, t_dynamics::FluxScale, fprintf(), Parser::getNumberCheck(), Parser::getNumberDefault(), Parser::getNumberPlain(), ioQQQ, t_pressure::lgContRadPresOn, Wind::lgDisk, t_dynamics::lgFluxDScale, t_dynamics::lgTracePrint, Wind::lgWindOK, Parser::nMatch(), pow(), pressure, Wind::setBallistic(), Wind::setDefault(), Wind::setStatic(), wind, Wind::windv, and Wind::windv0.

Referenced by ParseCommands().

Here is the call graph for this function:

Data Structures

Functions

Variables

Function Documentation

Variable Documentation