/home66/gary/public_html/cloudy/c08_branch/source/atmdat.h

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2008 by Gary J. Ferland and
 * others.  For conditions of distribution and use see copyright notice in license.txt */

#ifndef _ATMDAT_H_
#define _ATMDAT_H_


double atmdat_2phot_shapefunction( double EbyE2nu, long ipISO, long nelem );

void atmdat_2phot_rate(
        long ipZ , 
        long ipISO );

void atmdat_readin(void);

void atmdat_Chianti_readin( long intNS );

void atmdat_lamda_readin( long intNS );


void atmdat_outer_shell(
  long int iz, 
  long int in,
  long int *imax,
  long int *ig0, 
  long int *ig1);

void ChargTranEval( double *O_HIonRate );

double ChargTranSumHeat(void);

/*ChargTranPun punch charge transfer rate coefficients */
void ChargTranPun( FILE* ipPnunit , char* chPunch );

double atmdat_dielrec_fe(long int ion, double t);

void atmdat_H_phot_cs(void);

void atmdat_3body(void);

double atmdat_HS_caseB( 
        long int iHi, long int iLo, long int iZ, double TempIn,
        double DenIn,
        char chCase 
        );

void atmdat_DielSupres(void);

/* these are the vectors that store the original Hummer and Storey case B
 * line data for H and He - the declaration for the interpolator follows */
#define NHSDIM 15 
#define NLINEHS 300  
#define HS_NZ 8 
#define NHCSTE  8 
EXTERN struct t_atmdat {

        double HCharExcIonOf[LIMELM][LIMELM+1], //(cm3 s-1)
          HCharExcRecTo[LIMELM][LIMELM+1],              //(cm3 s-1)
          HCharHeatMax, 
          HCharCoolMax, 
          HCharHeatOn;

        /* rate coefficient (cm3 s-1) for N+(3P) + H+ -> N(2D) + H+ charge transfer*/
        double HCharExcRecTo_N0_2D;

        double HCharExcIonTotal,
          HCharExcRecTotal;

        double HeCharExcIonTotal,
          HeCharExcRecTotal;

        double HIonFrac;

        double HIonFracMax;

        double HeCharExcIonOf[LIMELM][LIMELM], 
          HeCharExcRecTo[LIMELM][LIMELM];

        double HCTAlex;

        bool lgCTOn;

        double Density[2][HS_NZ][NHSDIM], 
                ElecTemp[2][HS_NZ][NHSDIM],
                Emiss[2][HS_NZ][NHSDIM][NHSDIM][NLINEHS];

        long int nDensity[2][HS_NZ] , ntemp[2][HS_NZ] , ncut[2][HS_NZ];

        bool lgHCaseBOK[2][HS_NZ];

        long int nsbig;

        bool lgCollIonOn;

        realnum WaveLengthCaseB[8][25][24];

        realnum *CaseBWlHeI;
        long int nCaseBHeI;

        }       atmdat;

typedef enum { PHFIT_UNDEF, PHFIT95, PHFIT96 } phfit_version;

class t_ADfA : public Singleton<t_ADfA>
{
        friend class Singleton<t_ADfA>;
protected: 
        t_ADfA();
private:
        phfit_version version;
        /* phfit.dat */
        long int L[7];
        long int NINN[30];
        long int NTOT[30];
        realnum PH1[7][30][30][6];
        realnum PH2[30][30][7];
        /* hpfit.dat */
        realnum PHH[NHYDRO_MAX_LEVEL][5];
        /* rec_lines.dat */
        realnum P[8][110];
        realnum ST[9][405];
        /* rad_rec.dat */
        realnum rrec[30][30][2];
        realnum rnew[30][30][4];
        realnum fe[13][3];
        /* h_rad_rec */
        realnum HRF[NHYDRO_MAX_LEVEL][9];
        /* h_phot_cs.dat */
        realnum STH[NHYDRO_MAX_LEVEL];
        /* coll_ion.dat */
        double CF[30][30][5];
        /* h_coll_str.dat */
        /* >>refer      H1      cs      Anderson, H., Ballance, C.P., Badnell, N.R., & Summers, H.P., 
         * >>refercon   2000, J Phys B, 33, 1255; erratum, 2002 */
        double HCS[14][10][8];
public:
        void set_version(phfit_version val) { version = val; }

        phfit_version get_version() const { return version; }

        realnum ph1(int i, int j, int k, int l) const { return PH1[i][j][k][l]; }

        realnum sth(int i) const { return STH[i]; }

        double phfit(long int nz, long int ne, long int is, double e);

        double hpfit(long int iz, long int n, double e);

        double hpfit_rel(long int iz, long int n, double e);

        void rec_lines(double t, realnum r[][471]);

        double rad_rec(long int iz, long int in, double t);

        double H_rad_rec(long int iz, long int n, double t);

        double coll_ion(long int iz, long int in, double t);

        realnum h_coll_str( long ipLo, long ipHi, long ipTe );
};

#endif /* _ATMDAT_H_ */

---