rad_rcslw.h

Go to the documentation of this file.

 
 
#pragma once
 
#include <string>
#include <vector> 
#include "rad.h"
 
 
class rad_rcslw : public rad {
 
 
 
    protected:
 
        double       P;                   
        double       Tref;                
        double       Tb;                  
 
        std::vector<double> P_table;      
        std::vector<double> C_table;      
        std::vector<double> Tg_table;     
        std::vector<double> Tb_table;     
        std::vector<double> xH2O_table;   
 
        std::vector<double> F_pts;        
        std::vector<double> Ft_pts;       
 
        std::vector<double> Falbdf_CO2;   
        std::vector<double> Falbdf_CO;    
        std::vector<double> Falbdf_H2O;   
 
        int nP;                           
        int nC;                           
        int nTg;                          
        int nTb;                          
        int ny_H2O;                       
 
        double Cmin;                      
        double Cmax;                      
        double Fmin;                      
        double Fmax;                      
 
 
    public:
 
        void get_k_a(std::vector<double> &kabs,     // absorption coefficient (1/m)
                     std::vector<double> &awts,     // gas weight (sum to one)
                     const double   T,              // gas temperature (K)
                     const double   P_not_used,     // gas pressure (Pa); (interface only, set by constructor)
                     const double   fvsoot,         // volume fraction soot (not ppmv)
                     const double   xH2O,           // mole fraction H2O
                     const double   xCO2,           // mole fraction CO2
                     const double   xCO,            // mole fraction CO
                     const double   xCH4_not_used); // mole fraction CH4 (interface only, not used)
 
        void get_k_a_oneband(double         &kabs,            // absorption coefficient (1/m)
                             double         &awts,            // gas weight (sum to one)
                             const int      iband,            // band/gas to get k_a for
                             const double   T,                // gas temperature (K)
                             const double   P_not_used,       // gas pressure (Pa)
                             const double   fvsoot,           // volume fraction soot (not ppmv)
                             const double   xH2O,             // mole fraction h2o
                             const double   xCO2,             // mole fraction co2
                             const double   xCO,              // mole fraction co
                             const double   xCH4_not_used);   // mole fraction ch4
 
    private:
 
        double get_F_albdf(const double C, double Tg, double Tb, double xCO2, double xCO, double xH2O,
                           double fvsoot);
        double get_FI_albdf(const double F, const double Tg, const double Tb, 
                            const double xCO2, const double xCO, const double xH2O,
                            const double fvsoot);
        void get_FI_albdf_tables(const std::string Ptable_file_name, 
                                 const int nx, const int ny, const int nz,
                                 std::vector<double> &myarray);
        void get_FI_albdf_tables(const std::string Ptable_file_name, 
                                 const int nx, const int ny, const int nz, const int nw,
                                 std::vector<double> &myarray);
        void set_Fpts(void);
        void set_Falbdf_CO2_CO_H2O_at_P();
        double F_albdf_soot(const double C, const double Tg, const double Tb, const double fvsoot);
        
    
    public: 
    
        rad_rcslw(const int    p_nGG,
                  const double TbTref,
                  const double p_P,
                  const double fvsoot,
                  const double xH2O,
                  const double xCO2,
                  const double xCO);
 
        virtual ~rad_rcslw(){}
 
};