atomfeii.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 ATOMFEII_H_
#define ATOMFEII_H_
 
void AssertFeIIDep( double *pred , double *BigError , double *StdDev );
 
void FeIICreate(void);
void FeIIPrint(void);
 
void FeIILevelPops(void);
 
double FeIISumBand(realnum wl1, 
                                   realnum wl2,
                                   double *SumBandInward );
 
void FeII_RT_TauInc(void);
 
void FeII_RT_tau_reset(void);
 
void FeIIPoint(void);
 
 
void FeIIAccel(double *fe2drive);
 
void FeII_RT_Make( void );
 
void FeIIAddLines( void );
 
void FeIIPunchLevels( FILE * ioPUN );
 
void FeIIPunchColden(
  FILE *ioPUN );
 
void FeII_Colden( const char *chLabel );
 
 
void FeIISaveLines( FILE * ioPUN );
 
void FeIIPunchOpticalDepth( FILE * ioPUN );
 
void FeII_LineZero(void);
 
void FeIIIntenZero(void);
 
double FeIIRadPress(void);
 
double FeII_InterEnergy(void);
 
void FeIIPunDepart(FILE* ioPUN , 
        bool lgDoAll );
 
void PunFeII( FILE * io );
 
void FeIIPun1Depart(FILE * ioPUN , long int nPUN );
void FeIIPunData(
        FILE* ioPUN ,
        bool lgDoAll );
 
void FeIIPunPop(FILE* ioPUN , 
        bool lgPunchRange ,
        long int ipRangeLo ,
        long int ipRangeHi ,
        bool lgPunchDensity );
 
 
void FeIIPunchLineStuff( FILE * io , realnum xLimit  , long index);
 
#if 0
 
void FeIIPun1Pop( 
        FILE * ioPUN , 
        long int nPUN );
#endif
 
void FeIIZero(void);
 
void FeIIReset(void);
 
void FeII_OTS(void);
 
void FeII_RT_Out(void);
 
class Parser;
 
void ParseAtomFeII( Parser &p );
 
#define NFE2LEVN        371
 
extern bool lgFeIIMalloc;
 
struct t_FeII {
 
        long int nFeIILevel_local;
        long int nFeIILevel_malloc;
 
        bool lgFeIILargeOn;
 
        bool lgSlow;
 
        bool lgPrint;
 
        double FeIINRGs[NFE2LEVN];
 
        double FeIISTWT[NFE2LEVN];
 
        double FeIIAul[NFE2LEVN][NFE2LEVN];
 
        double FeIIColl[NFE2LEVN][NFE2LEVN][8];
 
        bool lgSimulate;
 
        char chFeIIAtom[7];
 
        bool lgShortFe2;
 
        bool lgLyaPumpOn;
 
        realnum fe2ener[2];
 
        realnum fe2thresh;
 
        realnum feconwlLo , feconwlHi;
        long int nfe2con;
 
        int ipRedisFcnResonance;
        int ipRedisFcnSubordinate;
 
        double Fe2_large_cool, 
          ddT_Fe2_large_cool, 
          Fe2_large_heat;
 
        double 
                ddT_Fe2_UVsimp_cool ,
                Fe2_UVsimp_cool,
                for7;
 
        };
extern t_FeII FeII;
 
        /* this info used to estimate destruction of Lya by overlap with FeII
 * in case where large atom is not turned on */
 
const int NFEII = 373;
const int NFE2PR = 61;
 
class t_fe2ovr_la : public Singleton<t_fe2ovr_la>
{
        friend class Singleton<t_fe2ovr_la>;
protected:
        t_fe2ovr_la();
private:
        realnum fe2lam[NFEII];
        realnum fe2osc[NFEII];
        realnum fe2enr[NFEII];
        realnum fe2gs[NFEII];
 
        long int ipfe2[NFEII];
 
        realnum feopc[NFEII];
        realnum Fe2TauLte[NFEII];
        realnum Fe2PopLte[NFEII];
 
        double fe2pt[NFE2PR];
        double fe2pf[NFE2PR];
 
        double fe2par(double te);
public:
        void zero_opacity();
 
        void init_pointers();
 
        void tau_inc();
 
        void atoms_fe2ovr(void);
};
 
#endif /* ATOMFEII_H_ */