adiabatic_mod_charnock Struct Reference

#include <ERF_MOSTStress.H>

Collaboration diagram for adiabatic_mod_charnock:

Public Member Functions
	adiabatic_mod_charnock (amrex::Real zref, amrex::Real Tflux, amrex::Real Qvflux, amrex::Real depth)
AMREX_GPU_DEVICE AMREX_FORCE_INLINE void	iterate_flux (const int &i, const int &j, const int &k, const int &max_iters, const amrex::Array4< amrex::Real > &z0_arr, const amrex::Array4< const amrex::Real > &umm_arr, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< const amrex::Real > &, const amrex::Array4< amrex::Real > &u_star_arr, const amrex::Array4< amrex::Real > &, const amrex::Array4< amrex::Real > &t_star_arr, const amrex::Array4< amrex::Real > &q_star_arr, const amrex::Array4< amrex::Real > &, const amrex::Array4< amrex::Real > &, const amrex::Array4< amrex::Real > &olen_arr, const amrex::Array4< amrex::Real > &, const amrex::Array4< amrex::Real > &, const amrex::Array4< amrex::Real > &, const amrex::Array4< amrex::Real > &) const

Private Attributes
most_data	mdata
similarity_funs	sfuns
const amrex::Real	tol_z = 1.0e-10
const amrex::Real	WSMIN = 0.1

Detailed Description

Adiabatic with modified charnock roughness

Constructor & Destructor Documentation

adiabatic_mod_charnock()

adiabatic_mod_charnock::adiabatic_mod_charnock ( amrex::Real  zref, amrex::Real  Tflux, amrex::Real  Qvflux, amrex::Real  depth  )

inline

    {
        mdata.zref = zref;
        mdata.surf_temp_flux  = Tflux;
        mdata.surf_moist_flux = Qvflux;
        mdata.Cnk_d = depth;
        mdata.Cnk_b = mdata.Cnk_b1 * std::log(mdata.Cnk_b2 / mdata.Cnk_d);
    }

Member Function Documentation

iterate_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void adiabatic_mod_charnock::iterate_flux ( const int &  i, const int &  j, const int &  k, const int &  max_iters, const amrex::Array4< amrex::Real > &  z0_arr, const amrex::Array4< const amrex::Real > &  umm_arr, const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &  , const amrex::Array4< const amrex::Real > &  , const amrex::Array4< amrex::Real > &  u_star_arr, const amrex::Array4< amrex::Real > &  , const amrex::Array4< amrex::Real > &  t_star_arr, const amrex::Array4< amrex::Real > &  q_star_arr, const amrex::Array4< amrex::Real > &  , const amrex::Array4< amrex::Real > &  , const amrex::Array4< amrex::Real > &  olen_arr, const amrex::Array4< amrex::Real > &  , const amrex::Array4< amrex::Real > &  , const amrex::Array4< amrex::Real > &  , const amrex::Array4< amrex::Real > &    ) const

inline

    {
        amrex::Real ustar  = 0.0;
        amrex::Real z0     = z0_arr(i,j,k);
        amrex::Real z0_old = z0;
        amrex::Real psi_m  = 0.0;
        amrex::Real umm    = std::max(umm_arr(i,j,k), WSMIN);
        amrex::Real C      = std::log(mdata.zref / z0);
 
        // Fixed point iteration on roughness
        int iter_z = 0;
        do {
            // Transfer curr to old
            z0_old = z0;
 
            // Update
            C = std::log(mdata.zref / z0_old);
            ustar = mdata.kappa * umm / (C - psi_m);
            z0    = std::exp( (2.7*ustar - 1.8/mdata.Cnk_b) / (ustar + 0.17/mdata.Cnk_b) );
 
            ++iter_z;
        } while ( (std::abs(z0 - z0_old) > tol_z) && (iter_z <= max_iters) );
        AMREX_ALWAYS_ASSERT_WITH_MESSAGE(iter_z < max_iters,
                                         "Maximum number of MOST roughness iterations reached.");
        C = std::log(mdata.zref / z0);
 
        // Populate the stored MOST arrays
        z0_arr(i,j,k)     = z0;
        olen_arr(i,j,k)   = 1.0e16;
        u_star_arr(i,j,k) = mdata.kappa * umm / (C - psi_m);
        t_star_arr(i,j,k) = 0.0;
        q_star_arr(i,j,k) = 0.0;
    }

Member Data Documentation

mdata

most_data adiabatic_mod_charnock::mdata

private

Referenced by adiabatic_mod_charnock(), and iterate_flux().

sfuns

similarity_funs adiabatic_mod_charnock::sfuns

private

tol_z

const amrex::Real adiabatic_mod_charnock::tol_z = 1.0e-10

private

Referenced by iterate_flux().

WSMIN

const amrex::Real adiabatic_mod_charnock::WSMIN = 0.1

private

Referenced by iterate_flux().

---

The documentation for this struct was generated from the following file:

• Source/BoundaryConditions/ERF_MOSTStress.H