surface_flux_charnock Struct Reference

#include <ERF_MOSTStress.H>

Collaboration diagram for surface_flux_charnock:

Public Member Functions
	surface_flux_charnock (amrex::Real Tflux, amrex::Real Qvflux, amrex::Real cnk_a, bool use_visc, bool cons_qflux)
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 > &tvm_arr, const amrex::Array4< const amrex::Real > &qvm_arr, const amrex::Array4< amrex::Real > &u_star_arr, const amrex::Array4< amrex::Real > &w_star_arr, const amrex::Array4< amrex::Real > &t_star_arr, const amrex::Array4< amrex::Real > &q_star_arr, const amrex::Array4< amrex::Real > &t_surf_arr, const amrex::Array4< amrex::Real > &q_surf_arr, const amrex::Array4< amrex::Real > &olen_arr, const amrex::Array4< amrex::Real > &pblh_arr, const amrex::Array4< amrex::Real > &, const amrex::Array4< amrex::Real > &, const amrex::Array4< amrex::Real > &) const

Private Attributes
most_data	mdata
bool	spec_qflux
similarity_funs	sfuns
const amrex::Real	tol = amrex::Real(1.0e-5)
const amrex::Real	alpha = amrex::Real(0.5)
const amrex::Real	WSMIN = most_min_wind_speed

Detailed Description

Surface flux with charnock roughness

Constructor & Destructor Documentation

surface_flux_charnock()

surface_flux_charnock::surface_flux_charnock ( amrex::Real  Tflux, amrex::Real  Qvflux, amrex::Real  cnk_a, bool  use_visc, bool  cons_qflux  )

inline

Construct specified-flux, 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
[in]	cons_qflux	whether the moisture flux is specified

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

Member Function Documentation

iterate_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void surface_flux_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 > &  tvm_arr, const amrex::Array4< const amrex::Real > &  qvm_arr, const amrex::Array4< amrex::Real > &  u_star_arr, const amrex::Array4< amrex::Real > &  w_star_arr, const amrex::Array4< amrex::Real > &  t_star_arr, const amrex::Array4< amrex::Real > &  q_star_arr, const amrex::Array4< amrex::Real > &  t_surf_arr, const amrex::Array4< amrex::Real > &  q_surf_arr, const amrex::Array4< amrex::Real > &  olen_arr, const amrex::Array4< amrex::Real > &  pblh_arr, 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.

    {
        int iter = 0;
        amrex::Real ustar = zero;
        amrex::Real wstar = zero;
        amrex::Real tflux = zero;
        amrex::Real qflux = zero;
        amrex::Real z0    = zero;
        amrex::Real zeta  = zero;
        amrex::Real psi_m = zero;
        amrex::Real psi_h = zero;
        amrex::Real Olen  = zero;
        amrex::Real zref  = zref_arr(i,j,k);
        amrex::Real C     = std::log(zref / z0_arr(i,j,k));
        amrex::Real zeta_fold = sfuns.calc_zeta_fold(C);
        bool at_fold = false;
        amrex::Real umm   = std::max(umm_arr(i,j,k), WSMIN);
        if (u_star_arr(i,j,k) == bogus_large_value) {
            u_star_arr(i,j,k) = mdata.kappa * umm / most_bracket(C, zero);
        } else {
            Olen  = olen_arr(i,j,k);
            zeta  = most_clamp_zeta(zref / Olen, zeta_fold, at_fold);
            psi_m = sfuns.calc_psi_m(zeta);
            psi_h = sfuns.calc_psi_h(zeta);
        }
        most_aitken accel;
        do {
            ustar = u_star_arr(i,j,k);
            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);
            }
            C         = std::log(zref / z0);
            zeta_fold = sfuns.calc_zeta_fold(C);
            qflux = (spec_qflux) ? mdata.surf_moist_flux :
                                 -(qvm_arr(i,j,k) - q_surf_arr(i,j,k)) * ustar * mdata.kappa /
                                  most_bracket(C, psi_h); // <w'Qv'>
            tflux = mdata.surf_temp_flux*(one + epsv*qvm_arr(i,j,k)) + qflux*epsv*tm_arr(i,j,k);
            if (w_star_arr) {
                // update w* and Umagmean
                w_star_arr(i,j,k) = calc_wstar(tflux, pblh_arr(i,j,k), tvm_arr(i,j,k));
                wstar = mdata.Bjr_beta * w_star_arr(i,j,k);
                umm = std::sqrt(umm_arr(i,j,k)*umm_arr(i,j,k) + wstar*wstar);
                umm = std::max(umm, WSMIN);
            }
            Olen  = -ustar * ustar * ustar * tvm_arr(i,j,k) / (mdata.kappa * mdata.gravity * tflux);
            // Stable side bounded at the fold, where the prescribed flux stops having a solution
            zeta  = most_clamp_zeta(zref / Olen, zeta_fold, at_fold);
            psi_m = sfuns.calc_psi_m(zeta);
            psi_h = sfuns.calc_psi_h(zeta);
            {
                amrex::Real ustar_new = mdata.kappa * umm / most_bracket(C, psi_m);
                u_star_arr(i,j,k) = (one - alpha) * ustar + alpha * ustar_new;
            }
            ++iter;
            u_star_arr(i,j,k) = accel.update(u_star_arr(i,j,k), iter);
        } while ((accel.extrapolated || 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; below the fold, apply the flux the wind can carry
        z0_arr(i,j,k)     = z0;
        most_flux_outputs(mdata, spec_qflux, at_fold, zeta_fold, C, psi_h, u_star_arr(i,j,k), zref,
                          tm_arr(i,j,k), tvm_arr(i,j,k), qvm_arr(i,j,k), qflux, Olen,
                          t_surf_arr(i,j,k), t_star_arr(i,j,k), q_surf_arr(i,j,k), q_star_arr(i,j,k));
        olen_arr(i,j,k)   = Olen;
    }

Here is the call graph for this function:

Member Data Documentation

alpha

const amrex::Real surface_flux_charnock::alpha = amrex::Real(0.5)

private

Referenced by iterate_flux().

mdata

most_data surface_flux_charnock::mdata

private

Referenced by iterate_flux(), and surface_flux_charnock().

sfuns

similarity_funs surface_flux_charnock::sfuns

private

Referenced by iterate_flux().

spec_qflux

bool surface_flux_charnock::spec_qflux

private

Referenced by iterate_flux(), and surface_flux_charnock().

tol

const amrex::Real surface_flux_charnock::tol = amrex::Real(1.0e-5)

private

Referenced by iterate_flux().

WSMIN

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