rt.h

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2013 by Gary J. Ferland and
 * others.  For conditions of distribution and use see copyright notice in license.txt */
 
#ifndef RT_H_
#define RT_H_
 
void RT_line_all( void );
 
double RT_line_driving(void);
 
/* in following two the logical variable says whether to do the
 * escape probabilities too (true) or just the destruction probabilities (false) */
 
void RT_line_one(const TransitionProxy& t, 
                 bool lgShield_this_zone,
                 realnum pestrk,
                 realnum DopplerWidth);
 
double RT_continuum_shield_fcn( const TransitionProxy &t );
 
void RT_diffuse(void);
 
void RT_continuum(void);
 
void RT_OTS(void);
 
void RT_OTS_AddLine(double ots, 
  long int ip );
 
void RT_OTS_Update(double* SumOTS);
 
void RT_OTS_Zero( void );
 
void RT_OTS_ChkSum(
        long int ipPnt);
 
void RT_line_one_tauinc(
        const TransitionProxy & t ,
        long int mas_species,
        long int mas_ion,
        long int mas_hi,
        long int mas_lo,
        realnum DopplerWidth);
 
void RT_tau_init(void);
 
class TransitionProxy;
void RT_line_one_tau_reset(
        const TransitionProxy & t );
 
void RT_tau_reset(void);
 
void RT_tau_inc(void);
 
void RT_OTS_PrtRate(
          double weak ,
          int chFlag );
 
double esc_PRD_1side(double tau, 
  double a);
 
double esc_CRDwing_1side(double tau, 
  double a );
 
double RTesc_lya(
        double *esin, 
        double *dest, 
        double abund, 
        const TransitionProxy& t, 
        realnum DopplerWidth);
 
double esc_CRDwing(double tau, 
  double tout, 
  double damp);
 
double esc_CRDcore(double tau, 
  double tout);
 
double esc_PRD(double tau, 
  double tout, 
  double damp );
 
double esca0k2(double taume);
 
double esccon(double tau, 
  double hnukt);
 
double RT_DestProb(
        double abund, 
          double crsec, 
          long int ipanu, 
          double widl, 
          double escp, 
          int nCore);
 
#if 0
 
double RT_LyaWidth(
                         double tauin, 
                         double tauout, 
                         double a, 
                         double vth);
#endif
 
double RT_LineWidth(
        const TransitionProxy & t,
        realnum DopplerWidth);
 
double RT_DestHummer(
        double beta); 
 
double RT_recom_effic(
        long int ip);
 
void RT_stark(void);
 
/* #define DEST0 1e-8 */
#define DEST0 SMALLFLOAT
 
struct t_rt {
 
        realnum wayin, 
 
          wayout; 
 
        realnum fracin;
 
        realnum TAddHLya, 
          TAddHeI;
 
        realnum DoubleTau;
 
        long int ipxry;
 
        realnum tauxry;
 
        bool lgFstOn;
 
        int nLineContShield;
 
        bool lgElecScatEscape;
 
        realnum dTauMase;
 
        bool lgMaserSetDR;
 
        bool lgMaserCapHit;
 
        long int mas_species , mas_ion , mas_hi , mas_lo;
 
        bool lgStarkON;
 
};
 
extern t_rt rt;
 
void RT_iso_integrate_RRC( const long ipISO, const long nelem, const bool lgUpdateContinuum );
 
#define LINE_CONT_SHIELD_PESC   1
#define LINE_CONT_SHIELD_FEDERMAN       2
#define LINE_CONT_SHIELD_FERLAND        3
 
#endif /* RT_H_ */