adiabatic_charnock Struct Reference

#include <ERF_MOSTStress.H>

Collaboration diagram for adiabatic_charnock:

Public Member Functions
	adiabatic_charnock (amrex::Real Tflux, amrex::Real Qvflux, amrex::Real cnk_a, bool use_visc)
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 > &tm_arr, 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 = most_min_wind_speed

Detailed Description

Adiabatic with charnock roughness

Constructor & Destructor Documentation

adiabatic_charnock()

adiabatic_charnock::adiabatic_charnock ( amrex::Real  Tflux, amrex::Real  Qvflux, amrex::Real  cnk_a, bool  use_visc  )

inline

Construct adiabatic Charnock-roughness MOST data.

Parameters

[in]	Tflux	prescribed heat flux
[in]	Qvflux	prescribed moisture flux
[in]	cnk_a	Charnock parameter
[in]	cnk_visc	whether to include viscous roughness

    {
        mdata.surf_temp_flux  = Tflux;
        mdata.surf_moist_flux = Qvflux;
        mdata.Cnk_a = cnk_a;
        mdata.visc  = use_visc;
    }

Member Function Documentation

iterate_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void adiabatic_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 > &  tm_arr, 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);
            if (mdata.Cnk_a > 0) {
                // NOTE: air_viscosity should take T and not Th!
                z0 = Charnock_roughness(ustar, mdata.Cnk_a,
                                        air_viscosity(tm_arr(i,j,k)), mdata.visc);
            } else {
                // NOTE: air_viscosity should take T and not Th!
                z0 = COARE3_roughness(zref, umm, ustar,
                                      air_viscosity(tm_arr(i,j,k)), mdata.visc);
            }
 
            ++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_charnock::mdata

private

Referenced by adiabatic_charnock(), and iterate_flux().

sfuns

similarity_funs adiabatic_charnock::sfuns

private

tol_z

const amrex::Real adiabatic_charnock::tol_z = amrex::Real(1.0e-10)

private

Referenced by iterate_flux().

WSMIN

const amrex::Real adiabatic_charnock::WSMIN = most_min_wind_speed

private

Referenced by iterate_flux().

---

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

• Source/BoundaryConditions/ERF_MOSTStress.H