atmdat_adfa.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 ATMDAT_ADFA_H_
#define ATMDAT_ADFA_H_

#include "iso.h"
#include "atmdat.h"

const int NRECCOEFCNO=471;
typedef enum { PHFIT_UNDEF, PHFIT95, PHFIT96 } phfit_version;

class t_ADfA : public Singleton<t_ADfA>
{
        friend class Singleton<t_ADfA>;
protected: 
        t_ADfA();
private:
        phfit_version version;
        /* phfit.dat */
        static const int NSHELLS = 7;
        long int L[NSHELLS];
        long int NINN[LIMELM];
        long int NTOT[LIMELM];
        static const int NFIT_PH1 = 6,
                        NFIT_PH2 = 7;
        realnum PH1[NSHELLS][LIMELM][LIMELM][NFIT_PH1];
        realnum PH2[LIMELM][LIMELM][NFIT_PH2];

        inline long set_vshell_index( long ne, long Z ) const;

        /* hpfit.dat */
        realnum PHH[NHYDRO_MAX_LEVEL][5];
        /* rec_lines.dat */
        realnum P[8][110];
        realnum ST[9][405];
        /* rad_rec.dat */
        realnum rrec[LIMELM][LIMELM][2];
        realnum rnew[LIMELM][LIMELM][4];
        realnum fe[13][3];
        /* h_rad_rec */
        realnum HRF[NHYDRO_MAX_LEVEL][9];
        /* h_phot_cs.dat */
        realnum STH[NHYDRO_MAX_LEVEL];
        /* coll_ion.dat */
        double CF[LIMELM][LIMELM][5];
        /* h_coll_str.dat */
        /* >>refer      H1      cs      Anderson, H., Ballance, C.P., Badnell, N.R., & Summers, H.P., 
         * >>refercon   2000, J Phys B, 33, 1255; erratum, 2002 */
        map<QNPair, array<realnum,NHCSTE>> HCS;
public:
        void set_version(phfit_version val) { version = val; }

        phfit_version get_version() const { return version; }

        realnum ph1(int i, int j, int k, int l) const { return PH1[i][j][k][l]; }

        double getEthresh( long nshell, long nel, long Z ) const;

        realnum sth(int i) const { return STH[i]; }

        double phfit(long int nz, long int ne, long int is, double e);

        double hpfit(long int iz, long int n, double e);

        void rec_lines(double t, realnum r[][NRECCOEFCNO]);

        double rad_rec(long int iz, long int in, double t);

        double H_rad_rec(long int iz, long int n, double t);

        double coll_ion(long int iz, long int in, double t);

        double coll_ion_wrapper(long int z, long int n, double t);

        /* coll_ion_hybrid computes hybrid collisional ionization rates */
        double coll_ion_hybrid( long int z, long int n, double t);

        const realnum* h_coll_str(long nHi, long lHi, long nLo, long lLo) const;
};

#endif