/home66/gary/public_html/cloudy/c08_branch/source/ionbal.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 _IONBAL_H_
#define _IONBAL_H_


void ion_recom_calculate( void );

void ion_trim(
        long int nelem );

void ion_zero(long int nelem);

void ion_collis(
        long nelem);

void ion_solver(long int nelem, 
  bool lgPrintIt);

void ion_photo(
        long int nelem , 
        bool lgPrintIt );

void ion_recomb(bool,const double*,const double*,const double[],const double[],
                const double[],const double[],const double[],const double[],long);

void ion_recombAGN( FILE * io );

void Badnell_rec_init( void );

/* routines to do heavy element ionization balance */
void IonAlumi(void);
void IonArgon(void);
void IonBeryl(void);
void IonBoron(void);
void IonCalci(void);
void IonCarbo(void);
void IonChlor(void);
void IonChrom(void);
void IonCobal(void);
void IonCoppe(void);
void IonFluor(void);
void IonHelium(void);
void IonIron(void);
void IonLithi(void);
void IonMagne(void);
void IonManga(void);
void IonNeon(void);
void IonNicke(void);
void IonNitro(void);
void IonOxyge(void);
void IonPhosi(void);
void IonPotas(void);
void IonScand(void);
void IonSilic(void);
void IonSodiu(void);
void IonSulph(void);
void IonTitan(void);
void IonVanad(void);
void IonZinc(void);

#define NSHELLS 7

EXTERN struct t_ionbal {

        double trimhi, 

          trimlo;

        bool lgTrimhiOn;

        /* ==============================================================
         * all following deals with ionization processes */

        double ****PhotoRate_Shell;

        bool lgPhotoIoniz_On;

        bool lgHO_ct_chem;

        double ***CollIonRate_Ground;

        double CosRayIonRate;

        double CosRayHeatNeutralParticles;

        double CosRayHeatThermalElectrons;

        double ExtraHeatRate;

        double xNeutronHeatRate;

        double PairProducPhotoRate[3];

        /* ==============================================================
         * following deal with Compton recoil ionization of bound electrons */
        
        bool lgCompRecoil;

        double CompRecoilHeatLocal;

        long int **ipCompRecoil;

        double **CompRecoilIonRate;

        double **CompRecoilIonRateSave;

        double **CompRecoilHeatRate;

        double **CompRecoilHeatRateSave;

        double **UTA_ionize_rate;
        double **UTA_heat_rate;

        bool lgInnerShellLine_on;
        bool lgInnerShell_Kisielius;
        bool lgInnerShell_Gu06;

        double **RateIonizTot;

        long int  nCompRecoilElec[LIMELM];

        double CompHeating_Max;
        /* ==============================================================
         * end Compton recoil ionization of bound electrons */

        /* ==============================================================
         * all following deals with recombination */

        double **RateRecomTot;

        double **RR_Badnell_rate_coef ,
                   **DR_Badnell_rate_coef,
                   *DR_Badnell_rate_coef_mean_ion;

        int **lgDR_Badnell_rate_coef_exist ,
                **lgRR_Badnell_rate_coef_exist;

        bool lg_use_DR_Badnell_rate_coef_mean_ion;

        bool lgDR_recom_Badnell_use,
                lgRR_recom_Badnell_use,
                lgRecom_Badnell_print;

        double **DR_old_rate_coef;

        double **RR_rate_coef_used,
                   **DR_rate_coef_used;

        double **RR_Verner_rate_coef;

        int nDR_S_guess;
        realnum DR_S_scale[5];

        realnum **GrainCreat,
                **GrainDestr;

        int 
          lgGrainIonRecom;

        realnum DielSupprs[2][LIMELM];

        realnum lg_guess_coef;
        realnum guess_noise;

        bool lgSupDie[2];

        realnum GuessDiel[4];

        bool lgNoCota;

        realnum CotaRate[LIMELM];

        long int ilt, 
          iltln, 
          ilthn, 
          ihthn, 
          ifail;

        }       ionbal;


#endif /* _IONBAL_H_ */

---