thermal.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 THERMAL_H_
#define THERMAL_H_

/* thermal.h */

#include "module.h"
#include "depth_table.h"

#define NCOLNT  10000L

void TempChange(
                         double TempNew ,
                         /* option to force update of all variables */
                         bool lgForceUpdate);

void TempChange(
                         double TempNew );

class t_thermal : public module
{
public:
        const char *chName() const
        {
                return "thermal";
        }

        void zero();
        void comment(t_warnings&) {}

        bool lgPredNextTe;

        bool lgTemperatureConstant;

        bool lgTemperatureConstantCommandParsed;

        realnum ConstTemp;

        realnum ConstGrainTemp;

        long int nUnstable;
        bool lgUnstable;

        realnum thist, 
          tlowst;

        bool lgTeHigh;

        bool lgEdnGTcm;

        bool lgTLaw;

        bool lgTeBD96;
        realnum T0BD96,
                SigmaBD96;

        realnum T0SN99;
        bool lgTeSN99;

        DepthTable tlaw;
        bool lgTeTLaw;

        realnum collam[NCOLNT];
        double cooling[NCOLNT], 
          heatnt[NCOLNT];
        long int ncltot;
#       define  NCOLNT_LAB_LEN  15
        char chClntLab[NCOLNT][NCOLNT_LAB_LEN+1];

        /* element coolants, the last one is molecular */
        double elementcool[LIMELM + 1];
        /* heavy element collisional ionization cooling
         * zero with CoolHeavy.colmet, add to element cooling
         * at cool_save.cpp */ 
        double heavycollcool[LIMELM];

        /* cooling due to level 2 lines */
        double dima;

        bool lgCNegChk;

        realnum CoolHeatMax;
        realnum wlCoolHeatMax;
        char chCoolHeatMax[NCOLNT_LAB_LEN+1];

        double totcol, 
          ctot, 
          heatl,
          coolheat;

        double dCooldT;

        double tsq1, 
          halfte;

        double te_update;

        bool lgCExtraOn;
        realnum CoolExtra, 
          cextpw;

        bool lgCoolEvalOK;

        realnum GBarMax;
        long int ipMaxExtra;

        double char_tran_heat , char_tran_cool;
        
        double  htot,

                power, 

                dHeatdT;

        double FreeFreeTotHeat;

        realnum HeatLineMax;
private:
        double m_heating[LIMELM][LIMELM];
public:
        double heating( long nelem, long ion )
        {
                return m_heating[nelem][ion];
        }
        void setHeating( long nelem, long ion, double heating )
        {
                m_heating[nelem][ion] = heating;
        }
        void AddHeating( long nelem, long ion, double heating )
        {
                m_heating[nelem][ion] += heating;
        }
};
extern t_thermal thermal;

struct t_phoHeat
{
        double HeatNet;
        double HeatLowEnr;
        double HeatHiEnr;
};

#endif /* THERMAL_H_ */