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

/* dense.h density related variables */

double dense_fabden(double radius, 
  double depth);

double dense_tabden(double r0, 
  double depth);

const int LIMTABDLAW = 500;

EXTERN struct t_dense
{
        t_dense() {
                /* list of atomic weights, mass in AMU, used for thermal line widths */
                /* >>refer      all     atomic weight   Coplen, T.B. 2001, J. Phys. Chem REf Data, 30, 701 */
                AtomicWeight[0] = 1.00794f;
                AtomicWeight[1] = 4.0026022f;
                AtomicWeight[2] = 6.9412f;
                AtomicWeight[3] = 9.0121823f;
                AtomicWeight[4] = 10.8117f;
                AtomicWeight[5] = 12.01078f;
                AtomicWeight[6] = 14.00672f;
                AtomicWeight[7] = 15.99943f;
                AtomicWeight[8] = 18.9984032f;
                AtomicWeight[9] = 20.17976f;
                AtomicWeight[10] = 22.989770f;
                AtomicWeight[11] = 24.30506f;
                AtomicWeight[12] = 26.9815382f;
                AtomicWeight[13] = 28.08553f;
                AtomicWeight[14] = 30.9737612f;
                AtomicWeight[15] = 32.0655f;
                AtomicWeight[16] = 35.4532f;
                AtomicWeight[17] = 39.9481f;
                AtomicWeight[18] = 39.09831f;
                AtomicWeight[19] = 40.0784f;
                AtomicWeight[20] = 44.9559108f;
                AtomicWeight[21] = 47.8671f;
                AtomicWeight[22] = 50.94151f;
                AtomicWeight[23] = 51.99616f;
                AtomicWeight[24] = 54.9380499f;
                AtomicWeight[25] = 55.8472f;
                AtomicWeight[26] = 58.9332009f;
                AtomicWeight[27] = 58.69342f;
                AtomicWeight[28] = 63.5463f;
                AtomicWeight[29] = 65.392f;
        }

        realnum gas_phase[LIMELM];

        realnum AtomicWeight[LIMELM];

        realnum xMolecules[LIMELM];

        realnum H_sum_in_CO;

        realnum xMassDensity;
                        
        realnum xMassDensity0;
                
        realnum pden;
                
        realnum wmole;
                
        realnum xNucleiTotal;
                
        realnum xMassTotal;

        realnum HCorrFac;

        long int IonLow[LIMELM+1];
        long int IonHigh[LIMELM+1];

        realnum xIonDense[LIMELM+3][LIMELM+1];

        realnum AbundanceLimit;

        bool lgElmtOn[LIMELM];

        bool lgSetIoniz[LIMELM];

        realnum SetIoniz[LIMELM][LIMELM+1];

        char chDenseLaw[5];

        /* this says keep initial density constant, 
         * so pressure from iter to iter not really const */
        bool lgDenseInitConstant;

        double DensityLaw[10];

        bool lgAsChoose[LIMELM][LIMELM];
        bool lgCSChoose[LIMELM][LIMELM];

        realnum frad[LIMTABDLAW];
        realnum fhden[LIMTABDLAW];

        long int nvals;

        bool lgDLWDepth;

        double eden;

        realnum eden_f;

        double density_low_limit;

        long int nzEdenBad;

        realnum EdenSet;

        realnum EdenExtra;

        double SqrtEden;

        double EdenHCorr;

        double EdenHontoHCorr;

        double EdenTrue;

        double eden_from_metals;

        bool lgEdenBad;

        double edensqte; 

        double cdsqte;

        realnum DensityPower;
        realnum rscale;
        realnum den0;

        bool lgDenFlucOn;

        bool lgDenFlucRadius;

        realnum flong;
        realnum cfirst; 
        realnum csecnd;
        realnum flcPhase;

} dense;


#endif /* _DENSE_H_ */

---