surface_temp_charnock Struct Reference

#include <ERF_MOSTStress.H>

Collaboration diagram for surface_temp_charnock:

Public Member Functions
	surface_temp_charnock (amrex::Real zref, amrex::Real Tflux, amrex::Real Qvflux, amrex::Real cnk_a, bool cnk_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< 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 = 1.0e-5
const amrex::Real	WSMIN = 0.1

Detailed Description

Surface temperature with charnock roughness

Constructor & Destructor Documentation

surface_temp_charnock()

surface_temp_charnock::surface_temp_charnock ( amrex::Real  zref, amrex::Real  Tflux, amrex::Real  Qvflux, amrex::Real  cnk_a, bool  cnk_visc, bool  cons_qflux  )

inline

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

Member Function Documentation

iterate_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void surface_temp_charnock::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 > &  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

    {
        int iter = 0;
        amrex::Real ustar = 0.0;
        amrex::Real wstar = 0.0;
        amrex::Real z0    = 0.0;
        amrex::Real tflux = 0.0;
        amrex::Real qflux = 0.0;
        amrex::Real zeta  = 0.0;
        amrex::Real psi_m = 0.0;
        amrex::Real psi_h = 0.0;
        amrex::Real Olen  = 0.0;
        amrex::Real Oleno = 0.0;
        amrex::Real umm   = std::max(umm_arr(i,j,k), WSMIN);
        if (u_star_arr(i,j,k) == 1.E34) {
            u_star_arr(i,j,k) = mdata.kappa * umm / std::log(mdata.zref / z0_arr(i,j,k));
        } else {
            Olen  = olen_arr(i,j,k);
            Oleno = Olen;
            zeta  = mdata.zref / Olen;
            psi_m = sfuns.calc_psi_m(zeta);
            psi_h = sfuns.calc_psi_h(zeta);
        }
        do {
            ustar = u_star_arr(i,j,k);
            if (mdata.Cnk_a > 0) {
                z0 = (mdata.Cnk_a / mdata.gravity) * ustar * ustar;
                if (mdata.visc) {
                    z0 += air_viscosity(tm_arr(i,j,k)) / std::max(ustar, 0.05);
                }
            } else {
                z0 = COARE3_roughness(mdata.zref, umm, ustar);
            }
            tflux = -(tm_arr(i,j,k) - t_surf_arr(i,j,k)) * ustar * mdata.kappa /
                     (std::log(mdata.zref / z0) - psi_h); // <w'T'>
            tflux *= (1 + 0.61*qvm_arr(i,j,k));
            qflux = (spec_qflux) ? mdata.surf_moist_flux :
                                 -(qvm_arr(i,j,k) - q_surf_arr(i,j,k)) * ustar * mdata.kappa /
                                  (std::log(mdata.zref / z0_arr(i,j,k)) - psi_h); // <w'Qv'>
            tflux += 0.61*tm_arr(i,j,k) * qflux; // ~= <w'Tv'>
            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);
            if ( (((Olen >= 0.0) && (Oleno <= 0.0)) ||
                  ((Olen <= 0.0) && (Oleno >= 0.0))) &&
                 std::fabs(Olen) + std::fabs(Oleno) < 1.) {
                Olen = 0.5 * (Olen + Oleno);
            }
            Oleno = Olen;
            zeta  = mdata.zref / Olen;
            psi_m = sfuns.calc_psi_m(zeta);
            psi_h = sfuns.calc_psi_h(zeta);
            u_star_arr(i,j,k) = mdata.kappa * umm / (std::log(mdata.zref / z0) - psi_m);
            ++iter;
        } while ((std::abs(u_star_arr(i,j,k) - ustar) > tol) && iter <= max_iters);
        AMREX_ASSERT_WITH_MESSAGE(iter < max_iters, "Maximum number of MOST iterations reached.");
 
        t_star_arr(i,j,k) = mdata.kappa * (tm_arr(i,j,k) - t_surf_arr(i,j,k)) /
                            (std::log(mdata.zref / z0) - psi_h);
          olen_arr(i,j,k) = Olen;
            z0_arr(i,j,k) = z0;
        if (spec_qflux) {
            q_surf_arr(i,j,k) = mdata.surf_moist_flux * (std::log(mdata.zref / z0_arr(i,j,k)) - psi_h) /
                                (u_star_arr(i,j,k) * mdata.kappa) + qvm_arr(i,j,k);
            q_star_arr(i,j,k) = -mdata.surf_moist_flux / u_star_arr(i,j,k);
        } else {
            q_star_arr(i,j,k) = mdata.kappa * (qvm_arr(i,j,k) - q_surf_arr(i,j,k)) /
                                (std::log(mdata.zref / z0_arr(i,j,k)) - psi_h);
        }
    }

Here is the call graph for this function:

Member Data Documentation

mdata

most_data surface_temp_charnock::mdata

private

Referenced by iterate_flux(), and surface_temp_charnock().

sfuns

similarity_funs surface_temp_charnock::sfuns

private

Referenced by iterate_flux().

spec_qflux

bool surface_temp_charnock::spec_qflux

private

Referenced by iterate_flux(), and surface_temp_charnock().

tol

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

private

Referenced by iterate_flux().

WSMIN

const amrex::Real surface_temp_charnock::WSMIN = 0.1

private

Referenced by iterate_flux().

---

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

• Source/BoundaryConditions/ERF_MOSTStress.H