dense.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 DENSE_H_
#define DENSE_H_
 
/* dense.h density related variables */
 
double dense_fabden(double radius, 
  double depth);
 
double dense_tabden(double r0, 
  double depth);
 
/*dense_parametric_wind called by dlaw command, returns density for any density law */
double dense_parametric_wind(double rad);
 
const int LIMTABDLAW = 500;
 
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;
                
                zero();
        }
 
        realnum gas_phase[LIMELM];
        void SetGasPhaseDensity( const long nelem, const realnum density );
 
        realnum AtomicWeight[LIMELM];
 
private:
        realnum m_xMolecules[LIMELM];
 
public:
        realnum xMolecules(long nelem)
        {
                return m_xMolecules[nelem];
        }
        void updateXMolecules();
        void zero();
 
        realnum xMassDensity;
                        
        realnum xMassDensity0;
                
        realnum pden;
                
        realnum wmole;
                
        realnum xNucleiTotal;
                
        realnum xMassTotal;
 
        realnum HCorrFac;
 
        long int IonLow[LIMELM+1];
        long int IonHigh[LIMELM+1];
 
        double xIonDense[LIMELM][LIMELM+1];
 
        bool lgIonChiantiOn[LIMELM][LIMELM+1];
 
        bool lgIonStoutOn[LIMELM][LIMELM+1];
 
        double maxWN[LIMELM][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;
 
        // let initial pressure vary with time
        bool lgPressureVaryTime;
 
        //  required number is timescale for time variation
        double PressureVaryTimeTimescale;
        //  optional number is index for time variation
        double PressureVaryTimeIndex;
 
        double DensityLaw[10];
 
        bool lgAsChoose[LIMELM][LIMELM];
        bool lgCSChoose[LIMELM][LIMELM];
 
        realnum frad[LIMTABDLAW];
        realnum fhden[LIMTABDLAW];
 
        long int nvals;
 
        bool lgDLWDepth;
 
        double eden;
 
        double EdenMax , EdenMin;
 
        double density_low_limit;
 
        long int nzEdenBad;
 
        realnum EdenSet;
 
        realnum EdenExtra;
 
        realnum EdenFraction;
 
        double SqrtEden;
 
        double EdenHCorr;
 
        realnum EdenHCorr_f;
 
        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;
 
};
extern t_dense dense;
 
// test whether nuclei conservation holds
bool lgElemsConserved(void);
 
#include "iso.h"
 
// test whether the sum of a set of states agrees with the ion population to the requested tolerance.
void lgStatesConserved ( long nelem, long ionStage, qList states,  long numStates, realnum err_tol, long loop_ion );
 
void SumDensities( void );
 
void ScaleAllDensities( const realnum factor );
void ScaleIonDensities( const long nelem, const realnum factor );
bool AbundChange( );
 
realnum scalingDensity(void);
realnum scalingZoneDensity(long i);
 
#endif /* DENSE_H_ */