#include "cddefines.h"
#include "dynamics.h"
#include "trace.h"
#include "elementnames.h"
#include "save.h"
#include "phycon.h"
#include "secondaries.h"
#include "stopcalc.h"
#include "grainvar.h"
#include "highen.h"
#include "hmi.h"
#include "pressure.h"
#include "taulines.h"
#include "rt.h"
#include "grains.h"
#include "atmdat.h"
#include "ionbal.h"
#include "ion_trim.h"
#include "opacity.h"
#include "cooling.h"
#include "thermal.h"
#include "iso.h"
#include "conv.h"
#include "h2.h"
#include "deuterium.h"
#include "mole.h"
#include "rfield.h"
#include "dense.h"

Include dependency graph for conv_base.cpp:

Functions
STATIC bool	lgNetEdenSrcSmall (void)

void	UpdateUTAs (void)

int	ConvBase (long loopi)

Function Documentation

int ConvBase ( long loopi )

ConvBase main routine to drive ionization solution for all species, find total opacity called by ConvIoniz/ lgConverg check whether ionization of element nelem has converged

Parameters

loopi this tells how many times ConvBase has been called by ConvIoniz while trying to converge electron density == 0 on first call - allows logic in ConvBase to check for ots oscillations

Todo:: 0 PvH prefers test against err/10

Definition at line 188 of file conv_base.cpp.

References ASSERT, atmdat, atmdat_3body(), t_thermal::char_tran_cool, t_thermal::char_tran_heat, t_atmdat::CharExcIonOf, t_atmdat::CharExcRecTo, ChargTranEval(), ChargTranSumHeat(), t_conv::chConvIoniz(), t_elementnames::chElementSym, t_save::chHashString, conv, CoolEvaluate(), CoolSave(), t_secondaries::csupra, t_thermal::ctot, DEBUG_ENTRY, molezone::den, dense, deut, diatoms, DumpLine(), DynaIonize(), dynamics, t_dense::eden, eden_sum(), t_conv::EdenErrorAllowed, t_dense::EdenTrue, t_mole_local::elec, elementnames, findspecies(), findspecieslocal(), fnzone, fp_equal(), fprintf(), t_deuterium::gas_phase, t_dense::gas_phase, GrainDrive(), gv, t_hmi::H2_total, t_conv::HeatCoolRelErrorAllowed, HeatSum(), HeatZero(), highen(), hmi, t_thermal::htot, INPUT_LINE_LENGTH, ion_recom_calculate(), ion_trim(), ion_trim2(), ion_trim_from_set(), ion_trim_validate(), ion_widen(), ion_wrapper(), ionbal, t_dense::IonHigh, t_conv::IonizErrorAllowed, t_dense::IonLow, ioQQQ, ipH_LIKE, ipHELIUM, ipHYDROGEN, ipOXYGEN, t_save::ipTraceConvergeBase, iso_collapsed_update(), iso_multiplet_opacities(), iso_renorm(), iso_solve(), iso_sp, iso_update_rates(), iteration, lgAbort, t_dynamics::lgAdvection, lgConvBaseHeatTest, t_conv::lgConvIoniz(), lgElemsConserved(), t_deuterium::lgElmtOn, t_dense::lgElmtOn, t_conv::lgFirstSweepThisZone, GrainVar::lgGrainElectrons, t_rfield::lgIonizReevaluate, t_conv::lgIonStageTrimed, lgNetEdenSrcSmall(), t_ionbal::lgNewTrim, t_rfield::lgOpacityReevaluate, t_opac::lgOpacStatic, t_conv::lgOscilOTS, t_trace::lgOTSBug, t_opac::lgRedoStatic, t_conv::lgSearch, lgStatesConserved(), t_thermal::lgTemperatureConstant, t_dynamics::lgTimeDependentStatic, t_trace::lgTrace, t_save::lgTraceConvergeBase, t_save::lgTraceConvergeBaseHash, t_conv::lgUpdateCouplings, LIMELM, t_conv::limPres2Ioniz, t_mole_global::list, MAX2, MIN2, mole, mole_dominant_rates(), mole_drive(), mole_global, mole_save(), mole_update_sources(), t_atmdat::NCX, NISO, t_conv::nPres2Ioniz, t_conv::nTotalIoniz, t_StopCalc::nTotalIonizStop, t_trace::nTrConvg, t_mole_global::num_calc, nzone, opac, OpacityAddTotal(), phycon, PresTotCurrent(), t_ionbal::RateIonizTot(), t_conv::register_(), t_conv::resetConvIoniz(), rfield, RT_line_all_escape(), RT_OTS(), RT_OTS_PrtRate(), RT_OTS_Update(), save, SaveHeat(), SDIV(), t_secondaries::sec2total, SecIoniz(), secondaries, t_conv::setConvIonizFail(), SetDeuteriumIonization(), SMALLFLOAT, t_mole_local::species, t_iso_sp::st, StopCalc, t_phycon::te, t_thermal::te_update, thermal, TorF(), GrainVar::TotalEden, trace, UpdateUTAs(), x0, x1, x2, t_deuterium::xIonDense, t_dense::xIonDense, and t_dense::xNucleiTotal.

Referenced by ConvInitSolution(), ConvIoniz(), and ion_recombAGN().

STATIC bool lgNetEdenSrcSmall ( void )

Definition at line 1324 of file conv_base.cpp.

References conv, DEBUG_ENTRY, dense, t_dense::EdenTrue, t_ionbal::elecsnk, t_ionbal::elecsrc, findspecies(), fixit, fprintf(), GrainVar::GrainChTrRate, gv, molecule::index, ionbal, ioQQQ, GrainVar::lgDustOn(), t_dense::lgElmtOn, t_conv::lgSearch, LIMELM, mole, nzone, SDIV(), t_conv::setConvIonizFail(), t_mole_local::species, and t_dense::xIonDense.

Referenced by ConvBase().

Here is the call graph for this function:

void UpdateUTAs ( void )

Definition at line 1280 of file conv_base.cpp.

References conv, DEBUG_ENTRY, dense, TransitionList::Emis(), fprintf(), ionbal, ioQQQ, t_conv::lgFirstSweepThisZone, LIMELM, TransitionList::size(), t_ionbal::UTA_heat_rate, t_ionbal::UTA_ionize_rate, UTALines, and t_dense::xIonDense.

Referenced by ConvBase().

Here is the call graph for this function:

Functions

Function Documentation