MOData.H

Go to the documentation of this file.

#ifndef MODATA_H
#define MODATA_H
 
#include <cstdint>
#include "src/utilities/trig_ops.H"
#include "AMReX_REAL.H"
 
using namespace amrex::literals;
 
namespace kynema_sgf {

struct MOData
{
    enum class ThetaCalcType : std::uint8_t {
        HEAT_FLUX = 0,      
        SURFACE_TEMPERATURE 
    };
 
    amrex::Real zref{0.0_rt};           
    amrex::Real z0{0.1_rt};             
    amrex::Real z0t{z0};                
    amrex::Real utau;                   
    amrex::Real kappa{0.41_rt};         
    amrex::Real gravity{9.81_rt};       
    amrex::Real obukhov_len{1.0e16_rt}; 
 
    amrex::RealArray vel_mean; 
    amrex::Real vmag_mean;     
    amrex::Real Su_mean;       
    amrex::Real Sv_mean;       
    amrex::Real theta_mean;    
 
    amrex::Real surf_temp_flux{0.0_rt}; 
    amrex::Real surf_temp;              
 
    amrex::Real gamma_m{5.0_rt};
    amrex::Real gamma_h{5.0_rt};
    amrex::Real beta_m{16.0_rt};
    amrex::Real beta_h{16.0_rt};
 
    ThetaCalcType alg_type{ThetaCalcType::HEAT_FLUX};
 
    [[nodiscard]] amrex::Real phi_m() const
    {
        return std::log(zref / z0) - calc_psi_m(zref / obukhov_len);
    }
 
    [[nodiscard]] amrex::Real phi_m(amrex::Real z) const
    {
        return std::log(z / z0) - calc_psi_m(z / obukhov_len);
    }
 
    [[nodiscard]] amrex::Real phi_h() const
    {
        return std::log(zref / z0t) - calc_psi_h(zref / obukhov_len);
    }
 
    [[nodiscard]] amrex::Real phi_h(amrex::Real z) const
    {
        return std::log(z / z0t) - calc_psi_h(z / obukhov_len);
    }
 
    [[nodiscard]] amrex::Real calc_psi_m(amrex::Real zeta) const;
    static amrex::Real
    calc_psi_m(amrex::Real zeta, amrex::Real beta_m, amrex::Real gamma_m);
    [[nodiscard]] amrex::Real calc_psi_h(amrex::Real zeta) const;
    static amrex::Real
    calc_psi_h(amrex::Real zeta, amrex::Real beta_h, amrex::Real gamma_h);
    void update_fluxes(int max_iters = 25);
};
 
} // namespace kynema_sgf
 
#endif /* MODATA_H */