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_ */