adiabatic_wave_coupled Struct Reference

#include <ERF_MOSTStress.H>

Collaboration diagram for adiabatic_wave_coupled:

Public Member Functions
	adiabatic_wave_coupled (amrex::Real zref, amrex::Real Tflux, amrex::Real Qvflux)
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 > &Hwave_arr, const amrex::Array4< amrex::Real > &Lwave_arr, const amrex::Array4< amrex::Real > &eta_arr) const

Private Attributes
most_data	mdata
similarity_funs	sfuns
const amrex::Real	tol = 1.0e-5
const amrex::Real	eps = 1e-15
const amrex::Real	z0_eps = 1.0e-6
const amrex::Real	z0_max = 0.1
const amrex::Real	WSMIN = 0.1

Detailed Description

Adiabatic with wave-coupled roughness

Constructor & Destructor Documentation

adiabatic_wave_coupled()

adiabatic_wave_coupled::adiabatic_wave_coupled ( amrex::Real  zref, amrex::Real  Tflux, amrex::Real  Qvflux  )

inline

    {
        mdata.zref = zref;
        mdata.surf_temp_flux  = Tflux;
        mdata.surf_moist_flux = Qvflux;
    }

Member Function Documentation

iterate_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void adiabatic_wave_coupled::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 > &  Hwave_arr, const amrex::Array4< amrex::Real > &  Lwave_arr, const amrex::Array4< amrex::Real > &  eta_arr  ) const

inline

    {
        int iter = 0;
        amrex::Real umm   = std::max(umm_arr(i,j,k), WSMIN);
        amrex::Real ustar = 0.0;
        amrex::Real z0    = 0.0;
        int ie, je;
        ie = i  < lbound(eta_arr).x ? lbound(eta_arr).x : i;
        je = j  < lbound(eta_arr).y ? lbound(eta_arr).y : j;
        ie = ie > ubound(eta_arr).x ? ubound(eta_arr).x : ie;
        je = je > ubound(eta_arr).y ? ubound(eta_arr).y : je;
        u_star_arr(i,j,k) = mdata.kappa * umm / std::log(mdata.zref / z0_arr(i,j,k));
        do {
            ustar = u_star_arr(i,j,k);
            z0    = std::min( std::max(1200.0 * Hwave_arr(i,j,k) * std::pow( Hwave_arr(i,j,k)/(Lwave_arr(i,j,k)+eps), 4.5 )
                                      + 0.11 * eta_arr(ie,je,k,EddyDiff::Mom_v) / ustar, z0_eps), z0_max );
            u_star_arr(i,j,k) = mdata.kappa * umm / std::log(mdata.zref / z0);
            ++iter;
        } while ((std::abs(u_star_arr(i,j,k) - ustar) > tol) && iter <= max_iters);
        AMREX_ALWAYS_ASSERT_WITH_MESSAGE(iter < max_iters,
                                         "Maximum number of MOST iterations reached.");
 
        // Populate the stored MOST arrays
        z0_arr(i,j,k)     = z0;
        olen_arr(i,j,k)   = 1.0e16;
        t_star_arr(i,j,k) = 0.0;
        q_star_arr(i,j,k) = 0.0;
    }

Member Data Documentation

eps

const amrex::Real adiabatic_wave_coupled::eps = 1e-15

private

Referenced by iterate_flux().

mdata

most_data adiabatic_wave_coupled::mdata

private

Referenced by adiabatic_wave_coupled(), and iterate_flux().

sfuns

similarity_funs adiabatic_wave_coupled::sfuns

private

tol

const amrex::Real adiabatic_wave_coupled::tol = 1.0e-5

private

Referenced by iterate_flux().

WSMIN

const amrex::Real adiabatic_wave_coupled::WSMIN = 0.1

private

Referenced by iterate_flux().

z0_eps

const amrex::Real adiabatic_wave_coupled::z0_eps = 1.0e-6

private

Referenced by iterate_flux().

z0_max

const amrex::Real adiabatic_wave_coupled::z0_max = 0.1

private

Referenced by iterate_flux().

---

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

• Source/BoundaryConditions/ERF_MOSTStress.H