============================================================

routines:

RTDiffuse in rtdiffuse.c
determines local thermal emission
adds energy to outward beams
called by cloudy

RT_OTS, in rtots.c, continuously updates ots rates, called during solution by 
ConvIonizeOpacityDo

ContCreatePointers ets incontinuum resolution and binning by calling fill
the arguemnts to fill can be changed to change the continuum resolution

ContSetIntensity - sets intensity of final continuum

====================================================================

built in continua are stored in predcont.h, set in zerologic.c
#include "predcont.h"
The array name is EnrPredCont and there are NPREDCONT values
lines entered into em line stack in lines_continuum.c

====================================================================

otscon is updated in AddOTSCon,
which is called by HydroOTS, HeMetalOTS, 

final ots continuum is updated, averaged and set in SumContinuum

====================================================================

unit integration
rfield.ConInterOut[rfield.nflux] is the point

rtdiffuse adds 1e-10 to these two
	rfield.ConEmitLocal[rfield.nflux] = 1.e-10f * Dilution /* */ ;
	rfield.ThrowOut[rfield.nflux] = 1.e-10f * Dilution /* */;

check is done in prtcomment

====================================================================

ConLocInter = local continuum that does interact, thrown into outward beam

====================================================================

rfield.nflux - number of continuum points defined in current continuum
originally set to rfield.nupper in zero and contcreatemesh
set to realistic value in consetintensity
it is RESIZED in convinitemp 

====================================================================

incident continua
rfield.FluxSave[j] - the original incident continuum
rfield.flux - outward attenuated incident continuum
	incident continuum is initally defined in ContSetIntensity

====================================================================

reflected continua
inward or reflected part of continuum, only computed in open geometry 
rfield.ConRefIncid[j] reflected incident continuum
rfield.ConRefDiff[j]) reflected diffuse continua
rfield.reflin[j] reflected lines 

====================================================================

diffuse local thermal continuum 
rfield.ConEmitLocal

====================================================================

outward continua
rfield.ThrowOut - outward beam of interactive continua,
	emitted in this zone, set in RTDiffuse, contains almost all of ConEmitLocal

rfield.ConInterOut - integrated outward beam, interacts

parts of brems and h lyman continua added to ConEmitLocal but not ConInterOut

====================================================================

rfield.SummedCon[i] is sum of interactive fields

contains 
rfield.flux[i] 
rfield.otslin[i] 
rfield.otscon[i] 

====================================================================

rfield.outlin[i] 
outward lines, included in ionizing radiation

rfield.outcon[i] 
outward portion of continuum, incremented in metdif

====================================================================

each cell has associate with it character descriptions
of lines of continua there:  these are:
rfield.chLineLabel[j] ,
rfield.chContLabel[j] );

=================================================================

print diffuse continuum sets flag PrnLine.lgPrnDiff
prints the continuum at the points entered in the vector 
EnrPredCont[NPREDCONT], as set in zerologic.c
added to stack in lineset1

within the continuum array the labels for lines and continua are
LineLabl.chLineLabel[j] ,
/* cont label*/
LineLabl.chContLabel[j] );

====================================================================

continuous optical depths
/* following only optical depth to illuminated face */
opac.depabs[i] += (float)(opac.opac[i]*radius.dReff);
opac.depsct[i] += (float)(opac.scatop[i]*radius.dReff);

/* these are total in inward direction, large if sherical */
opac.tausct[0][i] += (float)(opac.scatop[i]*radius.dReff);
opac.TauTotal[0][i] = opac.tauabs[0][i] + opac.tausct[0][i];

these are incremeneted in radinc.c

dtrans.chDffTrns =  "OTS"
or
dtrans.chDffTrns, "OUn"
default is "OU2" set in zero.c

============================================================

hydrogenic species are discussed in hydrogen.txt

====================================================================

bounds are stored in bounds.h
bounds.emm is low energy limit

============================================================
