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