adiabatic_wave_coupled Struct Reference

#include <ERF_MOSTStress.H>

Collaboration diagram for adiabatic_wave_coupled:

Public Member Functions
	adiabatic_wave_coupled (amrex::Real Tflux, amrex::Real Qvflux, 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 > &Hwave_arr, const amrex::Array4< amrex::Real > &Lwave_arr, 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	eps = amrex::Real(1e-15)
const amrex::Real	WSMIN = most_min_wind_speed

Detailed Description

Adiabatic with wave-coupled roughness

Constructor & Destructor Documentation

adiabatic_wave_coupled()

adiabatic_wave_coupled::adiabatic_wave_coupled ( amrex::Real  Tflux, amrex::Real  Qvflux, bool  use_visc  )

inline

Construct adiabatic wave-coupled roughness MOST data.

Parameters

[in]	Tflux	prescribed heat flux
[in]	Qvflux	prescribed moisture flux

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

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< 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 > &  Hwave_arr, const amrex::Array4< amrex::Real > &  Lwave_arr, 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);
            // NOTE: air_viscosity should take T and not Th!
            z0    = WaveCoupled_roughness(Hwave_arr(i,j,k), Lwave_arr(i,j,k), eps, 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

eps

const amrex::Real adiabatic_wave_coupled::eps = amrex::Real(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_z

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

private

Referenced by iterate_flux().

WSMIN

const amrex::Real adiabatic_wave_coupled::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