opacity.h

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2023 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 lgOpacAllocated;

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 );

struct t_opac : public module {

        const char *chName() const
        {
                return "opac";
        }
        void zero();
        void comment(t_warnings&) {}
        vector<realnum> TauAbsGeo[2];
        vector<realnum> TauScatGeo[2];

        vector<realnum> TauTotalGeo[2];

        vector<realnum> TauAbsFace, TauScatFace;

        vector<double> opacity_abs;
         
        vector<double> opacity_sct;

        vector<double> OldOpacSave;

        vector<double> albedo;

        vector<double> opacity_abs_savzon1;
        vector<double> opacity_sct_savzon1;

        vector<double> OpacStatic;

        vector<double> FreeFreeOpacity;
        vector<double> eeFreeFreeOpacity;

        double eeFreeFreeTemp;

        vector_avx<double> ExpZone;

        vector<realnum> E2TauAbsFace;
        vector<realnum> E2TauAbsTotal;
        vector<realnum> E2TauAbsOut;
        vector<realnum> TauAbsTotal;

        vector<realnum> ExpmTau;

        vector<realnum> tmn;

        bool lgOpacStatic;

        bool lgRedoStatic;

        vector<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, 

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