cont_createmesh.cpp File Reference

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

Include dependency graph for cont_createmesh.cpp:

Functions
STATIC void	rfield_opac_alloc ()
void	ContCreateMesh ()
void	rfield_opac_zero (long lo, long ihi)

Function Documentation

void ContCreateMesh

(

)

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

References multi_arr< T, d, ALLOC, lgBC >::alloc(), t_mesh::anu(), ASSERT, t_rfield::ConEmitLocal, t_rfield::ConSourceFcnLocal, t_rfield::csigc, t_rfield::csigh, DEBUG_ENTRY, t_rfield::fine_anu, t_rfield::fine_ener_hi, t_rfield::fine_ener_lo, fprintf(), geometry, ioQQQ, t_rfield::ipnt_coarse_2_fine, t_iterations::iter_alloc, iterations, t_mesh::lgMeshSetUp(), t_trace::lgTrace, MAX2, t_mesh::ncells(), t_iterations::nend, t_geometry::nend_max, t_rfield::nfine, t_rfield::nflux, t_rfield::nflux_with_check, t_rfield::nPositive, opac, t_rfield::resetCoarseTransCoef(), rfield, rfield_opac_alloc(), rfield_opac_zero(), t_opac::TauAbsFace, trace, and vzero().

Referenced by cloudy(), and ParseCompile().

Here is the call graph for this function:

STATIC void rfield_opac_alloc

(

)

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

References t_opac::albedo, multi_arr< T, d, ALLOC, lgBC >::alloc(), ASSERT, t_dense::AtomicWeight, 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::ContBoltzAvg, t_rfield::ContBoltzHelp1, t_rfield::ContBoltzHelp2, 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::eeBremsDif, t_opac::eeFreeFreeOpacity, t_mesh::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::ipFineConVelShift, t_rfield::ipnt_coarse_2_fine, lgRfieldAllocated, t_rfield::line_count, t_rfield::nfine, t_rfield::nflux_with_check, 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, t_rfield::reflin, Spectrum::resize(), rfield, t_rfield::setCoarseTransCoefPtr(), 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, t_rfield::vexp_arg, and vzero().

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

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::ContBoltzAvg, t_rfield::ContBoltzHelp1, t_rfield::ContBoltzHelp2, 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, lgRfieldAllocated, 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, t_opac::tmn, and vzero().

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

Here is the call graph for this function: