opacity.h

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2022 by Gary J. Ferland and
 * others.  For conditions of distribution and use see copyright notice in license.txt */

#ifndef OPACITY_H_
#define OPACITY_H_

#include "module.h"
#include "vectorize.h"

extern bool lgOpacMalloced;

void OpacityCreateAll(void);

void OpacityAdd1Subshell(
        long int ipOpac, 
        long int ipLowLim, 
        long int ipUpLim, 
        realnum abundance,
        char chStat );

void OpacityAddTotal(void);

void OpacityAdd1Element(
                long int ipZ);

void OpacityZero(void);

void OpacityZeroOld(void);

void OpacityAdd1SubshellInduc(
  long int ipOpac, 
  long int low, 
  long int ihi, 
  double a, 
  double b,
  char chStat );

/* limit to number of opacity cells available in the opacity stack
#define NOPSV   1300000L*/

/*lowest level where hydrogenic nu^-3 photo cross sect used*/
/*#define       NHPLPHOT        400*/

/* limit to number of cells in frequency array */
/* >>chng 00 oct 13, increase from 2000 to 2200 due to increase in FIR resolution */
/* >>chng 01 aug 19, increase from 2000 to 4000 due to increase in mid-range resolution */
/* >>chng 01 sep 24, get rid of NC ELL from code, totally MALLOC */
/*#define       NC ELL  4000L*/

struct t_opac : public module {

        const char *chName() const
        {
                return "opac";
        }
        void zero();
        void comment(t_warnings&) {}
        realnum **TauAbsGeo;
        realnum **TauScatGeo;

        realnum **TauTotalGeo;

        realnum *TauAbsFace, 
          *TauScatFace;

        double *opacity_abs;
         
        double *opacity_sct;

        double *OldOpacSave;

        double *albedo;

        double *opacity_abs_savzon1;
        double *opacity_sct_savzon1;

        double *OpacStatic;

        double *FreeFreeOpacity;
        double *eeFreeFreeOpacity;

        double eeFreeFreeTemp;

        vector< double, allocator_avx<double> > ExpZone;

        realnum *E2TauAbsFace;
        realnum *E2TauAbsTotal;
        realnum *E2TauAbsOut;
        realnum *TauAbsTotal;

        realnum *ExpmTau;

        realnum *tmn;

        bool lgOpacStatic;

        bool lgRedoStatic;

        double *OpacStack;

        realnum taumin;

        realnum tlamin;

        bool lgCaseB;

        bool lgCaseB_HummerStorey;

        bool lgCaseB_no_photo;

        bool lgCaseB_no_pdest;

        realnum telec, 
          thmin;

        bool lgOpacNeg;

        bool lgScatON;

        bool lgNegOpacIO;

        bool lgCompileOpac;

        bool lgUseFileOpac;

        long int nOpacTot;

        long int ipRayScat, 

        iopcom, 

        ippr, 
        ioppr, 

        ipBrems, 

        iphmra, 

        iphmop, 

        ih2pnt[2], 
        ih2pof, 
        /* excited state versions of above H2+ variables */
        ih2pnt_ex[2], 
        ih2pof_ex, 

        iophe1[9], 

        ioptri, 

        ipElement[LIMELM][LIMELM][7][3], 

        in1[3], 

        ipo3exc[3], 
        ipo3exc3[3], 
        ipo1exc[3];

        double o3exc;
        double o3exc3;

        long iopo2d, 

        ipmgex, 
        ipOpMgEx, 

        ica2ex[2], 
        ica2op;

        long int ipCKshell;

        realnum otsmin;

        realnum stimax[2];

        };
extern t_opac opac;



#endif /* OPACITY_H_ */