/home66/gary/public_html/cloudy/c10_branch/source/zero.cpp File Reference

#include "cddefines.h"
#include "physconst.h"
#include "iterations.h"
#include "hydrogenic.h"
#include "oxy.h"
#include "doppvel.h"
#include "dense.h"
#include "hextra.h"
#include "grains.h"
#include "magnetic.h"
#include "state.h"
#include "rt.h"
#include "he.h"
#include "struc.h"
#include "h2.h"
#include "coolheavy.h"
#include "lines.h"
#include "dynamics.h"
#include "carb.h"
#include "mean.h"
#include "atomfeii.h"
#include "iso.h"
#include "conv.h"
#include "geometry.h"
#include "timesc.h"
#include "peimbt.h"
#include "ionbal.h"
#include "continuum.h"
#include "atmdat.h"
#include "mole.h"
#include "ca.h"
#include "input.h"
#include "atoms.h"
#include "pressure.h"
#include "numderiv.h"
#include "colden.h"
#include "yield.h"
#include "hmi.h"
#include "rfield.h"
#include "abund.h"
#include "radius.h"
#include "opacity.h"
#include "secondaries.h"
#include "called.h"
#include "phycon.h"
#include "warnings.h"
#include "thermal.h"
#include "cooling.h"
#include "fe.h"
#include "hyperfine.h"
#include "init.h"
#include "dark_matter.h"

Go to the source code of this file.

Functions

void zero (void)

Function Documentation

void zero ( void   ) 

zero actively zero out or initialize variables needed for model calculation this is the old one and should be removed - its vars moved into those above

Todo:
2 change this to CS_new
Todo:
2 should iso.ipLyaRedist[ipHE_LIKE] use ipLY_A as does H-like?
Todo:
2 file opacities are disabled for now - reinstate this when arrays settle down

option to set proton elimination rates to zero >>refer CO chemistry Huntress, W. T., 1977, ApJS, 33, 495 By default, this is false - changed with set chemistry command

option to not include neutrals in the non-equilibrium scheme >> refer Federman, S. R. & Zsargo, J. 2003, ApJ, 589, 319 By default, this is false - changed with set chemistry command

HeatH2Dexc_used is heating due to collisional deexcitation of vib-excited H2 actually used

these are derivative wrt temp for collisional processes within X

the Solomon process rate H2 dissociates into X continuum - actually used

set_NaN( hmi.H2_Solomon_dissoc_rate_used );

H2 + hnu => 2H from TH85

H2 + hnu => 2H actually used

Definition at line 72 of file zero.cpp.

References abund, AbundancesZero(), t_colden::ajmmin, atmdat, ATOMIC_MASS_UNIT, atoms, t_conv::autocv, t_conv::AverEdenError, t_conv::AverHeatCoolError, t_conv::AverPressError, t_timesc::BigCOMoleForm, t_conv::BigEdenError, t_hmi::BiggestH2, t_conv::BigHeatCoolError, t_phycon::BigJumpCO, t_phycon::BigJumpH2, t_phycon::BigJumpne, t_phycon::BigJumpTe, t_conv::BigPressError, t_CoolHeavy::brems_cool_h, t_CoolHeavy::brems_cool_net, t_co::C12_C13_isotope_ratio, t_colden::C1Colden, t_colden::C1Pops, t_colden::C2Colden, t_colden::C2Pops, t_colden::C3Colden, t_colden::C3Pops, ca, t_ca::Ca2RmLya, t_ca::Ca3d, t_ca::Ca4p, called, carb, t_iso::CaseBCheck, cdEXIT, t_thermal::char_tran_cool, t_thermal::char_tran_heat, t_conv::chConvEden, t_conv::chConvIoniz, t_hmi::chGrainFormPump, t_hmi::chH2_small_model_type, t_hydro::chHTopType, t_hmi::chJura, t_LineSave::chNormLab, t_conv::chNotConverged, t_conv::chSolverEden, t_conv::chSolverTemp, t_input::chTitle, t_iso::chTypeAtomUsed, t_hydro::cintot, t_timesc::CloudAgeSet, co, t_co::codfrc, t_co::CODissHeat, t_co::codtot, t_colden::colden, colden, t_colden::coldenH2_ov_vel, t_ionbal::CollIonRate_Ground, t_CoolHeavy::colmet, t_ionbal::CompRecoilHeatRate, t_ionbal::CompRecoilHeatRateSave, t_ionbal::CompRecoilIonRate, t_ionbal::CompRecoilIonRateSave, continuum, conv, t_thermal::CoolExtra, t_hmi::CoolH2DexcMax, t_thermal::coolheat, t_thermal::CoolHeatMax, CoolHeavy, CoolZero(), t_ionbal::CotaRate, t_geometry::covgeo, t_geometry::covrt, t_hextra::CrsSecNeutron, t_secondaries::csupra, t_secondaries::csupra_effic, t_thermal::ctot, t_radius::CylindHigh, t_hydro::D2H_ratio, t_oxy::d4363, t_oxy::d5007r, t_oxy::d5007t, t_oxy::d6300, t_hydro::DampOnFac, t_radius::darea_x_fillfac, dark, t_conv::dCmHdT, DEBUG_ENTRY, dense, t_abund::Depletion, t_abund::depset, t_radius::depth, t_radius::depth_mid_zone, DEPTH_OFFSET, t_radius::depth_x_fillfac, t_hmi::deriv_HeatH2Dexc_BD96, t_hmi::deriv_HeatH2Dexc_BHT90, t_hmi::deriv_HeatH2Dexc_BigH2, t_hmi::deriv_HeatH2Dexc_ELWERT, t_hmi::deriv_HeatH2Dexc_TH85, t_hmi::deriv_HeatH2Dexc_used, t_rfield::DiffPumpOn, t_DoppVel::DispScale, t_radius::distance, DoppVel, t_rt::DoubleTau, t_ionbal::DR_Badnell_rate_coef, t_ionbal::DR_Badnell_rate_coef_mean_ion, t_ionbal::DR_old_rate_coef, t_ionbal::DR_rate_coef_used, t_radius::drad, t_radius::drad_mid_zone, t_radius::drad_x_fillfac, t_radius::drad_x_fillfac_mean, t_radius::dRadSign, t_radius::dRNeff, t_radius::drNext, t_ca::dstCala, t_rt::dTauMase, t_radius::dVeffAper, t_radius::dVeffVol, dynamics, DynaZero(), t_dense::eden, t_conv::EdenErrorAllowed, t_dense::EdenExtra, t_dense::EdenSet, t_hextra::effneu, t_continuum::EnergyKshell, t_conv::failmx, fe, t_fe::Fe7CoolTot, FeIIZero(), t_geometry::FillFac, t_geometry::filpow, t_rfield::fine_opac_nelem, t_rfield::fine_opac_nresolv, t_geometry::fiscal, fixit(), fnzone, t_he::frac_he0dest_23S, t_he::frac_he0dest_23S_photo, t_hextra::frcneu, t_thermal::FreeFreeTotHeat, t_dense::gas_phase, geometry, GrainZero(), t_pressure::gravity_symmetry, t_ionbal::guess_noise, t_colden::H0_21cm_lower, t_colden::H0_21cm_upper, h2, t_hmi::H2_formation_scale, t_hmi::H2_frac_abund_set, H2_Init(), t_hmi::H2_photodissoc_BHT90, t_hmi::H2_photodissoc_BigH2_H2g, t_hmi::H2_photodissoc_BigH2_H2s, t_hmi::H2_photodissoc_ELWERT_H2g, t_hmi::H2_photodissoc_ELWERT_H2s, t_hmi::H2_photodissoc_TH85, t_hmi::H2_photodissoc_used_H2g, t_hmi::H2_photodissoc_used_H2s, t_hmi::H2_Solomon_dissoc_rate_BD96_H2g, t_hmi::H2_Solomon_dissoc_rate_BD96_H2s, t_hmi::H2_Solomon_dissoc_rate_BHT90_H2g, t_hmi::H2_Solomon_dissoc_rate_BHT90_H2s, t_hmi::H2_Solomon_dissoc_rate_BigH2_H2g, t_hmi::H2_Solomon_dissoc_rate_BigH2_H2s, t_hmi::H2_Solomon_dissoc_rate_ELWERT_H2g, t_hmi::H2_Solomon_dissoc_rate_ELWERT_H2s, t_hmi::H2_Solomon_dissoc_rate_TH85_H2g, t_hmi::H2_Solomon_dissoc_rate_TH85_H2s, t_hmi::H2_Solomon_dissoc_rate_used_H2g, t_hmi::H2_Solomon_dissoc_rate_used_H2s, t_hmi::H2_total, H2_Zero(), t_hmi::h2dep, t_hmi::H2Opacity, t_hmi::h2pdep, t_hmi::h2plus_heat, t_hmi::h3pdep, t_dense::H_sum_in_CO, t_atmdat::HCharCoolMax, t_atmdat::HCharExcIonOf, t_atmdat::HCharExcIonTotal, t_atmdat::HCharExcRecTo, t_atmdat::HCharExcRecTotal, t_atmdat::HCharHeatMax, t_atmdat::HCharHeatOn, t_hydro::HCollIonMax, t_dense::HCorrFac, t_atmdat::HCTAlex, he, t_colden::He123S, t_conv::HeatCoolRelErrorAllowed, t_secondaries::HeatEfficPrimary, t_hmi::HeatH2Dexc_BD96, t_hmi::HeatH2Dexc_BHT90, t_hmi::HeatH2Dexc_BigH2, t_hmi::HeatH2Dexc_ELWERT, t_hmi::HeatH2Dexc_TH85, t_hmi::HeatH2Dexc_used, t_hmi::HeatH2DexcMax, t_hmi::HeatH2Dish_BD96, t_hmi::HeatH2Dish_BHT90, t_hmi::HeatH2Dish_BigH2, t_hmi::HeatH2Dish_ELWERT, t_hmi::HeatH2Dish_TH85, t_hmi::HeatH2Dish_used, t_thermal::heatl, HeatZero(), t_DoppVel::Heiles_Troland_F, hextra, t_atmdat::HIonFrac, t_atmdat::HIonFracMax, hmi, t_hmi::hmidep, t_hmi::hmihet, t_hmi::hmitot, hmole_init(), t_co::hmole_mass, t_hmi::Hmolec, t_thermal::htot, hydro, hyperfine, t_geometry::iEmissPower, t_ionbal::ifail, t_ionbal::ihthn, t_ionbal::ilt, t_ionbal::ilthn, t_ionbal::iltln, input, Singleton< t_yield >::Inst(), Singleton< t_fe2ovr_la >::Inst(), Singleton< T >::Inst(), t_pressure::IntegRhoGravity, ionbal, ioPrnErr, ioQQQ, t_ionbal::ipCompRecoil, ipCRD, ipCRDW, ipH2p, ipH_LIKE, ipHe2p1P, ipHE_LIKE, ipHELIUM, ipHYDROGEN, ipIRON, ipLY_A, t_iso::ipLyaRedist, t_LineSave::ipNormWavL, t_pressure::ipPradMax_nzone, ipPRD, t_iso::ipResoRedist, t_iso::ipSubRedist, iso, t_iterations::iter_malloc, iteration, iterations, t_iterations::itermx, t_iterations::IterPrnt, t_ionbal::lg_guess_coef, lgAbort, t_conv::lgAutoIt, t_rfield::lgBlockHIon, t_thermal::lgCExtraOn, t_atmdat::lgChiantiOn, t_thermal::lgCNegChk, t_iso::lgColl_excite, t_iso::lgColl_ionize, t_iso::lgColl_l_mixing, t_atmdat::lgCollIonOn, t_iso::lgCollStrenThermAver, t_rfield::lgCompileGauntFF, t_opac::lgCompileOpac, t_iso::lgCompileRecomb, t_ionbal::lgCompRecoil, t_continuum::lgCon0, t_iso::lgContinuumLoweringEnabled, t_conv::lgConvEden, t_iterations::lgConverge_set, t_conv::lgConvPres, t_conv::lgConvTemp, t_continuum::lgCoStarInterpolationCaution, t_iso::lgCritDensLMix, t_iso::lgCS_None, t_iso::lgCS_therm_ave, t_iso::lgCS_Vriens, t_iso::lgCS_Vrinceanu, t_secondaries::lgCSetOn, t_atmdat::lgCTOn, t_radius::lgCylnOn, t_abund::lgDepln, t_iso::lgDielRecom, t_rfield::lgDoLineTrans, t_radius::lgdR2Small, t_ionbal::lgDR_Badnell_rate_coef_exist, t_ionbal::lgDR_recom_Badnell_use, t_radius::lgDrMinUsed, t_radius::lgDrMnOn, t_rt::lgElecScatEscape, t_geometry::lgEndDflt, t_iso::lgErrGenDone, t_co::lgFederman, lgFirstCall, t_CoolHeavy::lgFreeOn, t_iso::lgFSM, t_rt::lgFstOn, t_state::lgGet_state, t_mole::lgGrain_mole_deplete, t_ionbal::lgGrainIonRecom, t_co::lgH2Ozer, t_rfield::lgHabing, t_rfield::lgHeIIOTS, t_hydro::lgHiPop2, t_ionbal::lgHO_ct_chem, t_iso::lgInd2nu_On, t_ionbal::lgInnerShell_Gu06, t_ionbal::lgInnerShell_Kisielius, t_ionbal::lgInnerShellLine_on, t_rfield::lgIonizReevaluate, t_rfield::lgKillOTSLine, t_rfield::lgKillOutCont, t_rfield::lgKillOutLine, t_atmdat::lgLamdaOn, t_hmi::lgLeiden_Keep_ipMH2s, t_hmi::lgLeidenCRHack, t_iso::lgLevelsEverLowered, t_iso::lgLevelsLowered, t_pressure::lgLineRadPresOn, t_hyperfine::lgLya_pump_21cm, t_rfield::lgLyaOTS, t_hydro::lgLymanPumping, t_conv::lgMap, t_rt::lgMaserCapHit, t_rt::lgMaserSetDR, t_rfield::lgMustBlockHIon, t_iso::lgMustReeval, t_opac::lgNegOpacIO, t_co::lgNeutrals, t_hextra::lgNeutrnHeatOn, t_dark_matter::lgNFW_Set, t_co::lgNoCOMole, t_ionbal::lgNoCota, t_hmi::lgNoH2Mole, t_co::lgNonEquilChem, t_iso::lgNoRecombInterp, t_LineSave::lgNormSet, t_NumDeriv::lgNumDeriv, t_rfield::lgOpacityReevaluate, t_opac::lgOpacNeg, t_opac::lgOpacStatic, t_iso::lgPessimisticErrors, t_ionbal::lgPhotoIoniz_On, t_phycon::lgPhysOK, t_pressure::lgPradCap, t_pressure::lgPradDen, t_iso::lgPrintNumberOfLevels, lgPrnErr, t_co::lgProtElim, t_iso::lgPrtDepartCoef, t_iso::lgPrtLevelPops, t_state::lgPut_state, t_radius::lgRadiusKnown, t_pressure::lgRadPresAbortOK, t_iso::lgRandErrGen, t_ionbal::lgRecom_Badnell_print, t_ionbal::lgRR_Badnell_rate_coef_exist, t_ionbal::lgRR_recom_Badnell_use, t_rfield::lgSaveOpacityFine, t_opac::lgScatON, t_radius::lgSdrmaxRel, t_radius::lgSdrminRel, t_secondaries::lgSecOFF, t_radius::lgSMinON, t_pressure::lgSonicPoint, t_pressure::lgSonicPointAbortOK, t_geometry::lgSphere, t_state::lgState_print, t_geometry::lgStatic, t_geometry::lgStaticNoIt, t_pressure::lgStrongDLimbo, t_ionbal::lgSupDie, t_called::lgTalk, t_called::lgTalkSave, lgTestCodeCalled, lgTestCodeEnabled, t_iso::lgTopoff, t_ionbal::lgTrimhiOn, t_DoppVel::lgTurb_pressure, t_DoppVel::lgTurbEquiMag, t_DoppVel::lgTurbLawOn, t_co::lgUMISTrates, t_opac::lgUseFileOpac, t_geometry::lgZoneSet, t_geometry::lgZoneTrp, LIMELM, t_conv::LimFail, t_conv::limPres2Ioniz, LineSave, Magnetic_init(), MALLOC, mean, t_mean::MeanZero(), mole, N_H_MOLEC, N_OI_LEVELS, NCOLD, t_iso::nCS_new, t_ionbal::nDR_S_guess, t_hmi::nElectron, t_geometry::nend, t_geometry::nEndDflt, t_conv::nGrainFail, t_conv::nIonFail, NISO, t_hydro::nLyaHot, t_iso::nLyaLevel, t_iso::nLyman, t_iso::nLyman_malloc, t_conv::nNeFail, t_atoms::nNegOI, t_conv::nPopFail, t_conv::nPreFail, t_conv::nPres2Ioniz, t_geometry::nprint, t_hmi::nProton, t_atmdat::nsbig, NSHELLS, t_conv::nTeFail, t_conv::nTotalFailures, t_conv::nTotalIoniz, NumDeriv, nzone, t_he::nzone, t_colden::O1Colden, t_colden::O1Pops, opac, t_h2::ortho_colden, t_opac::otsmin, oxy, t_oxy::p1401, t_oxy::p1666, t_carb::p1909, t_carb::p2326, t_h2::para_colden, peimbt, PHFIT96, t_ionbal::PhotoRate_Shell, phycon, t_oxy::poiexc, t_oxy::poiii2, t_oxy::poiii3, t_hydro::pop2mx, t_ca::popca2ex, t_atoms::popmg2, t_atoms::popoi, t_thermal::power, t_pressure::pres_radiation_lines_curr, t_pressure::PresMax, t_pressure::PresRamCurr, pressure, t_conv::PressureErrorAllowed, t_pressure::PresTotlCurr, t_pressure::PresTotlInit, t_radius::r1r0sq, t_dark_matter::r_200, t_dark_matter::r_s, t_pressure::RadBetaMax, t_radius::Radius, radius, t_radius::Radius_mid_zone, t_iso::RadRec_caseB, t_ionbal::RateIoniz, t_ionbal::RateRecomTot, t_radius::rdfalt, t_yield::reset_yield(), t_continuum::ResolutionScaleFactor, rfield, t_pressure::RhoGravity, t_pressure::RhoGravity_dark, t_pressure::RhoGravity_external, t_pressure::RhoGravity_self, t_radius::rinner, t_colden::rjnmin, t_ionbal::RR_Badnell_rate_coef, t_ionbal::RR_rate_coef_used, t_ionbal::RR_Verner_rate_coef, rt, t_hmi::ScaleJura, t_abund::ScaleMetals, t_LineSave::ScaleNormLine, t_radius::sdrmax, t_radius::sdrmin, t_secondaries::sec2total, t_secondaries::SecExcitLya2PrimaryErg, t_secondaries::SecHIonMax, t_secondaries::SecIon2PrimaryErg, secondaries, t_pressure::self_mass_factor, set_NaN(), t_secondaries::SetCsupra, ShowMe(), t_colden::Si2Colden, t_colden::Si2Pops, t_LineSave::sig_figs, t_mole::sink, t_iso::SmallA, SMALLFLOAT, t_abund::solar, t_timesc::sound_speed_isothermal, t_mole::source, state, t_opac::stimax, t_radius::StopRadius, t_radius::StopThickness, t_hmi::Tad, TempChange(), t_hydro::TexcLya, thermal, t_rfield::time_continuum_scale, t_timesc::time_H2_Dest_here, t_timesc::time_H2_Dest_longest, t_timesc::time_H2_Form_here, t_timesc::time_H2_Form_longest, t_timesc::TimeH21cm, timesc, t_hydro::TLyaMax, t_colden::tmas, t_thermal::totcol, t_colden::TotMassColl, t_hextra::totneu, t_ionbal::trimhi, t_ionbal::trimlo, t_peimbt::tsqden, t_DoppVel::TurbVel, t_DoppVel::TurbVelLaw, t_iso::TwoNu_induc_dn_max, t_dynamics::Upstream_hden, t_ionbal::UTA_heat_rate, t_ionbal::UTA_ionize_rate, t_hmi::UV_Cont_rel2_Draine_DB96_depth, t_hmi::UV_Cont_rel2_Draine_DB96_face, t_hmi::UV_Cont_rel2_Habing_TH85_depth, t_hmi::UV_Cont_rel2_Habing_TH85_face, t_LineSave::WavLNorm, wcnint(), t_thermal::wlCoolHeatMax, t_colden::wmas, t_secondaries::x12tot, t_dense::xIonDense, t_hydro::xLymanPumpingScaleFactor, t_dense::xMassDensity0, t_dense::xMassTotal, t_mole::xMoleChTrRate, t_dense::xMolecules, t_dense::xNucleiTotal, and t_fe2ovr_la::zero_opacity().

Referenced by cdInit(), hypho(), and optimize_func().

Here is the call graph for this function:

 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines

Generated on 15 Nov 2012 for cloudy by  doxygen 1.6.1