ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Pages
Public Member Functions | Private Attributes | List of all members
surface_temp_donelan Struct Reference

#include <ERF_MOSTStress.H>

Collaboration diagram for surface_temp_donelan:
Collaboration graph
[legend]

Public Member Functions

 surface_temp_donelan (amrex::Real zref, amrex::Real Tflux, amrex::Real Qvflux, 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 donelan roughness

Constructor & Destructor Documentation

◆ surface_temp_donelan()

surface_temp_donelan::surface_temp_donelan ( amrex::Real  zref,
amrex::Real  Tflux,
amrex::Real  Qvflux,
bool  cons_qflux 
)
inline
1338  {
1339  mdata.zref = zref;
1340  mdata.surf_temp_flux = Tflux;
1341  mdata.surf_moist_flux = Qvflux;
1342  spec_qflux = cons_qflux;
1343  }
most_data::surf_moist_flux
amrex::Real surf_moist_flux
Moisture flux.
Definition: ERF_MOSTStress.H:20
most_data::zref
amrex::Real zref
Reference height (m)
Definition: ERF_MOSTStress.H:15
most_data::surf_temp_flux
amrex::Real surf_temp_flux
Heat flux.
Definition: ERF_MOSTStress.H:19
surface_temp_donelan::mdata
most_data mdata
Definition: ERF_MOSTStress.H:1437
surface_temp_donelan::spec_qflux
bool spec_qflux
Definition: ERF_MOSTStress.H:1438

Member Function Documentation

◆ iterate_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void surface_temp_donelan::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
1368  {
1369  int iter = 0;
1370  amrex::Real ustar = 0.0;
1371  amrex::Real wstar = 0.0;
1372  amrex::Real z0 = 0.0;
1373  amrex::Real tflux = 0.0;
1374  amrex::Real qflux = 0.0;
1375  amrex::Real zeta = 0.0;
1376  amrex::Real psi_m = 0.0;
1377  amrex::Real psi_h = 0.0;
1378  amrex::Real Olen = 0.0;
1379  amrex::Real Oleno = 0.0;
1380  amrex::Real umm = std::max(umm_arr(i,j,k), WSMIN);
1381  if (u_star_arr(i,j,k) == 1.E34) {
1382  u_star_arr(i,j,k) = mdata.kappa * umm / std::log(mdata.zref / z0_arr(i,j,k));
1383  } else {
1384  Olen = olen_arr(i,j,k);
1385  Oleno = Olen;
1386  zeta = mdata.zref / Olen;
1387  psi_m = sfuns.calc_psi_m(zeta);
1388  psi_h = sfuns.calc_psi_h(zeta);
1389  }
1390  do {
1391  ustar = u_star_arr(i,j,k);
1392  z0 = Donelan_roughness(ustar);
1393  tflux = -(tm_arr(i,j,k) - t_surf_arr(i,j,k)) * ustar * mdata.kappa /
1394  (std::log(mdata.zref / z0) - psi_h); // <w'T'>
1395  tflux *= (1 + 0.61*qvm_arr(i,j,k));
1396  qflux = (spec_qflux) ? mdata.surf_moist_flux :
1397  -(qvm_arr(i,j,k) - q_surf_arr(i,j,k)) * ustar * mdata.kappa /
1398  (std::log(mdata.zref / z0_arr(i,j,k)) - psi_h); // <w'Qv'>
1399  tflux += 0.61*tm_arr(i,j,k) * qflux; // ~= <w'Tv'>
1400  if (w_star_arr) {
1401  // update w* and Umagmean
1402  w_star_arr(i,j,k) = calc_wstar(tflux, pblh_arr(i,j,k), tvm_arr(i,j,k));
1403  wstar = mdata.Bjr_beta * w_star_arr(i,j,k);
1404  umm = std::sqrt(umm_arr(i,j,k)*umm_arr(i,j,k) + wstar*wstar);
1405  umm = std::max(umm, WSMIN);
1406  }
1407  Olen = -ustar * ustar * ustar * tvm_arr(i,j,k) / (mdata.kappa * mdata.gravity * tflux);
1408  if ( (((Olen >= 0.0) && (Oleno <= 0.0)) ||
1409  ((Olen <= 0.0) && (Oleno >= 0.0))) &&
1410  std::fabs(Olen) + std::fabs(Oleno) < 1.) {
1411  Olen = 0.5 * (Olen + Oleno);
1412  }
1413  Oleno = Olen;
1414  zeta = mdata.zref / Olen;
1415  psi_m = sfuns.calc_psi_m(zeta);
1416  psi_h = sfuns.calc_psi_h(zeta);
1417  u_star_arr(i,j,k) = mdata.kappa * umm / (std::log(mdata.zref / z0) - psi_m);
1418  ++iter;
1419  } while ((std::abs(u_star_arr(i,j,k) - ustar) > tol) && iter <= max_iters);
1420  AMREX_ASSERT_WITH_MESSAGE(iter < max_iters, "Maximum number of MOST iterations reached.");
1421 
1422  t_star_arr(i,j,k) = mdata.kappa * (tm_arr(i,j,k) - t_surf_arr(i,j,k)) /
1423  (std::log(mdata.zref / z0) - psi_h);
1424  olen_arr(i,j,k) = Olen;
1425  z0_arr(i,j,k) = z0;
1426  if (spec_qflux) {
1427  q_surf_arr(i,j,k) = mdata.surf_moist_flux * (std::log(mdata.zref / z0_arr(i,j,k)) - psi_h) /
1428  (u_star_arr(i,j,k) * mdata.kappa) + qvm_arr(i,j,k);
1429  q_star_arr(i,j,k) = -mdata.surf_moist_flux / u_star_arr(i,j,k);
1430  } else {
1431  q_star_arr(i,j,k) = mdata.kappa * (qvm_arr(i,j,k) - q_surf_arr(i,j,k)) /
1432  (std::log(mdata.zref / z0_arr(i,j,k)) - psi_h);
1433  }
1434  }
Donelan_roughness
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real Donelan_roughness(amrex::Real ustar)
Definition: ERF_MOSTRoughness.H:29
calc_wstar
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real calc_wstar(const amrex::Real &ust, const amrex::Real &tst, const amrex::Real &qst, const amrex::Real &pblh, const amrex::Real &th, const amrex::Real &thv, const amrex::Real &qv=0.0)
Definition: ERF_Wstar.H:13
most_data::kappa
amrex::Real kappa
von Karman constant
Definition: ERF_MOSTStress.H:17
most_data::gravity
amrex::Real gravity
Acceleration due to gravity (m/s^2)
Definition: ERF_MOSTStress.H:18
most_data::Bjr_beta
const amrex::Real Bjr_beta
Definition: ERF_MOSTStress.H:29
similarity_funs::calc_psi_m
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real calc_psi_m(amrex::Real zeta) const
Definition: ERF_MOSTStress.H:41
similarity_funs::calc_psi_h
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real calc_psi_h(amrex::Real zeta) const
Definition: ERF_MOSTStress.H:55
surface_temp_donelan::sfuns
similarity_funs sfuns
Definition: ERF_MOSTStress.H:1439
surface_temp_donelan::tol
const amrex::Real tol
Definition: ERF_MOSTStress.H:1440
surface_temp_donelan::WSMIN
const amrex::Real WSMIN
Definition: ERF_MOSTStress.H:1441
Here is the call graph for this function:

Member Data Documentation

◆ mdata

most_data surface_temp_donelan::mdata
private

Referenced by iterate_flux(), and surface_temp_donelan().

◆ sfuns

similarity_funs surface_temp_donelan::sfuns
private

Referenced by iterate_flux().

◆ spec_qflux

bool surface_temp_donelan::spec_qflux
private

Referenced by iterate_flux(), and surface_temp_donelan().

◆ tol

const amrex::Real surface_temp_donelan::tol = 1.0e-5
private

Referenced by iterate_flux().

◆ WSMIN

const amrex::Real surface_temp_donelan::WSMIN = 0.1
private

Referenced by iterate_flux().

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