/home66/gary/public_html/cloudy/c13_branch/source/cont_createmesh.cpp File Reference

#include "cddefines.h"
#include "rfield.h"
#include "iterations.h"
#include "physconst.h"
#include "dense.h"
#include "trace.h"
#include "opacity.h"
#include "ipoint.h"
#include "geometry.h"
#include "continuum.h"

Include dependency graph for cont_createmesh.cpp:

Go to the source code of this file.

Functions
STATIC void	read_continuum_mesh (void)
STATIC void	fill (double fenlo, double fenhi, double resolv, long int *n0, long int *ipnt, bool lgCount)
STATIC void	rfield_opac_malloc (void)
STATIC void	ChckFill (void)
void	ContCreateMesh (void)
void	rfield_opac_zero (long lo, long ihi)

---

Function Documentation

STATIC void ChckFill

(

void

)

Definition at line 361 of file cont_createmesh.cpp.

References t_rfield::anu, ASSERT, cdEXIT, continuum, DEBUG_ENTRY, t_rfield::egamry, t_rfield::emm, EXIT_FAILURE, t_continuum::filbnd, ioQQQ, ipoint(), t_continuum::nrange, t_rfield::nupper, rfield, and t_rfield::widflx.

Referenced by ContCreateMesh().

Here is the call graph for this function:

void ContCreateMesh

(

)

set up continuum energy mesh if first call, otherwise reset to original mesh

Definition at line 38 of file cont_createmesh.cpp.

References t_rfield::anu, t_rfield::anu2, t_rfield::AnuOrg, t_rfield::anusqr, cdEXIT, ChckFill(), t_rfield::ConEmitLocal, t_rfield::ConSourceFcnLocal, continuum, t_rfield::csigc, t_rfield::csigh, DEBUG_ENTRY, t_rfield::egamry, t_rfield::emm, EXIT_FAILURE, fill(), t_rfield::fine_anu, t_rfield::fine_ener_hi, t_rfield::fine_ener_lo, geometry, ioQQQ, t_rfield::ipnt_coarse_2_fine, t_iterations::iter_malloc, iterations, t_trace::lgTrace, MALLOC, MAX2, NCELL, t_geometry::nend, t_geometry::nend_max, t_rfield::nfine_malloc, t_rfield::nflux, t_continuum::nrange, t_continuum::nStoredBands, t_rfield::nupper, read_continuum_mesh(), t_rfield::resetCoarseTransCoef(), rfield, rfield_opac_malloc(), rfield_opac_zero(), t_continuum::StoredEnergy, t_continuum::StoredResolution, trace, and t_rfield::widflx.

Referenced by cloudy(), GridGatherInCloudy(), and ParseCompile().

Here is the call graph for this function:

STATIC void fill	(	double	fenlo,
		double	fenhi,
		double	resolv,
		long int *	n0,
		long int *	ipnt,
		bool	lgCount
	)

Definition at line 254 of file cont_createmesh.cpp.

References t_rfield::anu, ASSERT, cdEXIT, continuum, DEBUG_ENTRY, EXIT_FAILURE, t_continuum::filbnd, t_continuum::fildel, t_continuum::filres, t_continuum::ifill0, ioQQQ, t_trace::lgConBug, t_trace::lgPtrace, t_trace::lgTrace, MAX2, t_continuum::nrange, t_continuum::nStoredBands, t_rfield::nupper, pow(), rfield, ShowMe(), trace, and t_rfield::widflx.

Referenced by ContCreateMesh(), and zero().

Here is the call graph for this function:

STATIC void read_continuum_mesh

(

void

)

Definition at line 657 of file cont_createmesh.cpp.

References cdEXIT, continuum, DEBUG_ENTRY, EXIT_FAILURE, FFmtRead(), t_continuum::filbnd, t_continuum::fildel, t_continuum::filres, t_continuum::ifill0, INPUT_LINE_LENGTH, ioQQQ, t_trace::lgTrace, MALLOC, t_continuum::nStoredBands, open_data(), read_whole_line(), t_continuum::ResolutionScaleFactor, t_continuum::StoredEnergy, t_continuum::StoredResolution, and trace.

Referenced by ContCreateMesh().

Here is the call graph for this function:

STATIC void rfield_opac_malloc

(

void

)

Todo:

3 consider making the fine opacity array a double. with a float, the opacity itself often becomes a denormalized number, it then becomes significant when multiplied by dr - can cause numerical noise. this is why the coarse opacity array is a double

Todo:

1 set temp of 1e4K will be too coarse a line for PDRs where H2 line overlap is very important

Definition at line 438 of file cont_createmesh.cpp.

References t_opac::albedo, t_rfield::anu, t_rfield::anu2, t_rfield::anu3, t_rfield::anulog, t_rfield::AnuOrg, t_rfield::anusqr, ASSERT, ATOMIC_MASS_UNIT, t_dense::AtomicWeight, BOLTZMANN, t_rfield::chContLabel, t_rfield::chLineLabel, t_rfield::comdn, t_rfield::comup, t_rfield::ConEmitOut, t_rfield::ConEmitReflec, t_rfield::ConInterOut, t_rfield::ConOTS_local_OTS_rate, t_rfield::ConOTS_local_photons, t_rfield::ConRefIncid, t_rfield::ContBoltz, t_rfield::convoc, t_rfield::csigc, t_rfield::csigh, DEBUG_ENTRY, dense, t_rfield::DiffuseEscape, t_rfield::DiffuseLineEmission, t_opac::E2TauAbsFace, t_opac::E2TauAbsOut, t_opac::E2TauAbsTotal, t_rfield::emm, t_opac::ExpmTau, t_opac::ExpZone, t_rfield::ExtinguishFactor, t_rfield::fine_anu, t_rfield::fine_ener_hi, t_rfield::fine_ener_lo, t_rfield::fine_opac_nelem, t_rfield::fine_opac_nresolv, t_rfield::fine_opac_velocity_width, t_rfield::fine_opac_zone, t_rfield::fine_opt_depth, t_rfield::fine_resol, t_rfield::flux, t_rfield::flux_accum, t_rfield::flux_beam_const, t_rfield::flux_beam_const_save, t_rfield::flux_beam_time, t_rfield::flux_isotropic, t_rfield::flux_isotropic_save, t_rfield::flux_time_beam_save, t_rfield::flux_total_incident, t_opac::FreeFreeOpacity, t_rfield::gff, t_rfield::ipFineConVelShift, t_rfield::ipnt_coarse_2_fine, lgRfieldMalloced, LIMELM, t_rfield::line_count, MALLOC, NCELL, t_rfield::nfine, t_rfield::nfine_malloc, t_rfield::nupper, t_rfield::OccNumbBremsCont, t_rfield::OccNumbContEmitOut, t_rfield::OccNumbDiffCont, t_rfield::OccNumbIncidCont, t_opac::OldOpacSave, opac, t_opac::opacity_abs, t_opac::opacity_abs_savzon1, t_opac::opacity_sct, t_opac::opacity_sct_savzon1, t_opac::OpacStatic, t_rfield::otscon, t_rfield::otslin, t_rfield::otssav, t_rfield::outlin, t_rfield::outlin_noplot, pow(), t_rfield::reflin, rfield, t_rfield::setCoarseTransCoefPtr(), SPEEDLIGHT, t_rfield::SummedCon, t_rfield::SummedDif, t_rfield::SummedDifSave, t_rfield::SummedOcc, t_opac::TauAbsFace, t_opac::TauAbsGeo, t_opac::TauAbsTotal, t_opac::TauScatFace, t_opac::TauScatGeo, t_opac::TauTotalGeo, t_opac::tmn, TotalInsanity(), t_rfield::TotDiff2Pht, and t_rfield::widflx.

Referenced by ContCreateMesh().

Here is the call graph for this function:

void rfield_opac_zero	(	long	lo,
		long	ihi
	)

zero out rfield arrays between certain limits, code in zero.c

Definition at line 801 of file cont_createmesh.cpp.

References t_opac::albedo, t_rfield::ConEmitOut, t_rfield::ConEmitReflec, t_rfield::ConInterOut, t_rfield::ConOTS_local_OTS_rate, t_rfield::ConOTS_local_photons, t_rfield::ConRefIncid, t_rfield::ContBoltz, t_rfield::convoc, t_opac::E2TauAbsFace, t_opac::E2TauAbsOut, t_opac::E2TauAbsTotal, t_opac::ExpmTau, t_opac::ExpZone, t_rfield::fine_opac_zone, t_rfield::fine_opt_depth, t_rfield::flux, t_rfield::flux_accum, t_rfield::flux_beam_const_save, t_rfield::flux_isotropic_save, t_rfield::flux_time_beam_save, t_rfield::flux_total_incident, t_opac::FreeFreeOpacity, lgRfieldMalloced, t_rfield::nfine_malloc, t_rfield::OccNumbBremsCont, t_rfield::OccNumbContEmitOut, t_rfield::OccNumbDiffCont, t_opac::OldOpacSave, opac, t_opac::opacity_abs, t_opac::opacity_sct, t_opac::OpacStatic, t_rfield::otscon, t_rfield::otslin, t_rfield::outlin, t_rfield::outlin_noplot, t_rfield::reflin, rfield, t_rfield::SummedCon, t_rfield::SummedDif, t_rfield::SummedOcc, t_opac::TauAbsGeo, t_opac::TauAbsTotal, t_opac::taumin, t_opac::TauScatGeo, t_opac::TauTotalGeo, and t_opac::tmn.

Referenced by ContCreateMesh(), ConvInitSolution(), and InitSimPostparse().