adiabatic_mod_charnock Struct Reference

#include <ERF_MOSTStress.H>

Collaboration diagram for adiabatic_mod_charnock:

Public Member Functions
	adiabatic_mod_charnock (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< const amrex::Real > &zref_arr, 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 = amrex::Real(1.0e-10)
const amrex::Real	WSMIN = amrex::Real(0.1)

Detailed Description

Adiabatic with modified charnock roughness

Constructor & Destructor Documentation

adiabatic_mod_charnock()

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

inline

Construct adiabatic modified-Charnock MOST data.

Parameters

[in]	Tflux	prescribed heat flux
[in]	Qvflux	prescribed moisture flux
[in]	depth	water depth for modified Charnock roughness

    {
        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< const amrex::Real > &  zref_arr, 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

Iterate MOST state for one surface point.

Calling sequence

i, j, k, max_iters, zref_arr, z0_arr, umm_arr, tm_arr, tvm_arr, qvm_arr, u_star_arr, w_star_arr, t_star_arr, q_star_arr, t_surf_arr, q_surf_arr, olen_arr, pblh_arr, Hwave_arr, Lwave_arr, eta_arr.

    {
        amrex::Real ustar  = zero;
        amrex::Real zref   = zref_arr(i,j,k);
        amrex::Real z0     = z0_arr(i,j,k);
        amrex::Real z0_old = z0;
        amrex::Real psi_m  = zero;
        amrex::Real umm    = std::max(umm_arr(i,j,k), WSMIN);
        amrex::Real C      = std::log(zref / z0);
 
        // Fixed point iteration on roughness
        int iter_z = 0;
        do {
            // Transfer curr to old
            z0_old = z0;
 
            // Update
            C = std::log(zref / z0_old);
            ustar = mdata.kappa * umm / (C - psi_m);
            z0    = std::exp( (amrex::Real(2.7)*ustar - amrex::Real(1.8)/mdata.Cnk_b) / (ustar + amrex::Real(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(zref / z0);
 
        // Populate the stored MOST arrays
        z0_arr(i,j,k)     = z0;
        olen_arr(i,j,k)   = bogus_large_value;
        u_star_arr(i,j,k) = mdata.kappa * umm / (C - psi_m);
        t_star_arr(i,j,k) = zero;
        q_star_arr(i,j,k) = zero;
    }

Here is the call graph for this function:

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 = amrex::Real(1.0e-10)

private

Referenced by iterate_flux().

WSMIN

const amrex::Real adiabatic_mod_charnock::WSMIN = amrex::Real(0.1)

private

Referenced by iterate_flux().

---

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

• Source/BoundaryConditions/ERF_MOSTStress.H