radius_next.cpp File Reference

#include "cddefines.h"
#include "iso.h"
#include "geometry.h"
#include "h2.h"
#include "hyperfine.h"
#include "opacity.h"
#include "dense.h"
#include "heavy.h"
#include "grainvar.h"
#include "elementnames.h"
#include "dynamics.h"
#include "thermal.h"
#include "hmi.h"
#include "coolheavy.h"
#include "timesc.h"
#include "stopcalc.h"
#include "colden.h"
#include "phycon.h"
#include "rt.h"
#include "trace.h"
#include "wind.h"
#include "save.h"
#include "pressure.h"
#include "iterations.h"
#include "struc.h"
#include "radius.h"
#include "dark_matter.h"
#include "mole.h"
#include "rfield.h"
#include "freebound.h"
#include "lines_service.h"
#include "cosmology.h"

Include dependency graph for radius_next.cpp:

Classes
class	drChoiceItem
class	drList

Functions
STATIC void	ContRate (double *freqm, double *opacm)
STATIC void	GrainRateDr (double *grfreqm, double *gropacm)
void	radius_next ()

Function Documentation

ContRate()

STATIC void ContRate	(	double *	freqm,
		double *	opacm )

Todo

2 - this routine is very important since it sets the pace for the calculation, and directly affects the convergence of the code. Most of the logic is very old and messy.
1) make sure all test cases have save dr 2) cat all these reasons together into one file and sort on the reason 3) discover what logic is the main pacesetter for the code 4) which are never triggered and so can be removed

References ASSERT, CoolHeavy, DEBUG_ENTRY, dense, fprintf(), gv, Heavy, ioQQQ, ipCARBON, ipH1s, ipH2s, ipH_LIKE, ipHYDROGEN, iso_sp, MAX2, MIN2, opac, rfield, safe_div(), SDIV(), SMALLFLOAT, STATIC, and thermal.

Referenced by radius_next().

Here is the call graph for this function:

GrainRateDr()

STATIC void GrainRateDr	(	double *	grfreqm,
		double *	gropacm )

References CoolHeavy, DEBUG_ENTRY, dense, Heavy, ipCARBON, ipH2s, ipH_LIKE, ipHELIUM, ipHYDROGEN, iso_sp, MAX2, opac, rfield, safe_div(), STATIC, and thermal.

Referenced by radius_next().

Here is the call graph for this function:

radius_next()

void radius_next

(

)

radius_next use adaptive logic to find next zone thickness return 0 if ok, 1 for abort

References abund, ASSERT, TransitionProxy::associated(), drList::begin(), cdEXIT, chLineLbl(), drList::clear(), colden, TransitionProxy::Coll(), column(), molezone::column, ContRate(), cosmology, EmissionProxy::damp(), EmissionProxy::dampXvel(), dark, DEBUG_ENTRY, molezone::den, dense, dense_fabden(), dense_parametric_wind(), diatoms, dynamics, elementnames, TransitionProxy::Emis(), drList::end(), EXIT_FAILURE, exp10(), findspecieslocal(), findspecieslocal_validate(), fixit, FndLineHt(), fprintf(), geometry, GetHubbleFactor(), GrainRateDr(), CollisionProxy::heat(), hmi, hyperfine, drList::insert(), ioQQQ, ipCARBON, ipCOL_HTOT, ipH_LIKE, ipHELIUM, ipHYDROGEN, iso_sp, iteration, iterations, LIMELM, MAX2, MIN2, mole, mole_global, nzone, opac, EmissionProxy::opacity(), EmissionProxy::Pesc(), phycon, EmissionProxy::PopOpc(), pressure, EmissionProxy::pump(), radius, rfield, rt, safe_div(), SAFETY, save, SDIV(), drList::size(), SMALLFLOAT, StopCalc, struc, EmissionProxy::TauIn(), thermal, timesc, TotalInsanity(), trace, EmissionProxy::VoigtLineCen(), and wind.

Referenced by cloudy().

Here is the call graph for this function: