ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
adiabatic_donelan Struct Reference

#include <ERF_MOSTStress.H>

Collaboration diagram for adiabatic_donelan:

Public Member Functions

 adiabatic_donelan (amrex::Real Tflux, amrex::Real Qvflux)
 
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 > &, 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 = 1.0e-5
 
const amrex::Real WSMIN = 0.1
 

Detailed Description

Adiabatic with Donelan roughness

Constructor & Destructor Documentation

◆ adiabatic_donelan()

adiabatic_donelan::adiabatic_donelan ( amrex::Real  Tflux,
amrex::Real  Qvflux 
)
inline
362  {
363  mdata.surf_temp_flux = Tflux;
364  mdata.surf_moist_flux = Qvflux;
365  }
most_data mdata
Definition: ERF_MOSTStress.H:418
amrex::Real surf_moist_flux
Moisture flux.
Definition: ERF_MOSTStress.H:19
amrex::Real surf_temp_flux
Heat flux TODO: decide whether this is <θ'w'> or <θv'w'> under moist conditions.
Definition: ERF_MOSTStress.H:18

Member Function Documentation

◆ iterate_flux()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void adiabatic_donelan::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 > &  ,
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
391  {
392  int iter = 0;
393  amrex::Real umm = std::max(umm_arr(i,j,k), WSMIN);
394  amrex::Real ustar = 0.0;
395  amrex::Real zref = zref_arr(i,j,k);
396  amrex::Real z0 = 0.0;
397  u_star_arr(i,j,k) = mdata.kappa * umm / std::log(zref / z0_arr(i,j,k));
398  if (u_star_arr(i,j,k) == 1.E34) {
399  u_star_arr(i,j,k) = mdata.kappa * umm / std::log(zref / z0_arr(i,j,k));
400  }
401  do {
402  ustar = u_star_arr(i,j,k);
403  z0 = Donelan_roughness(ustar);
404  u_star_arr(i,j,k) = mdata.kappa * umm / std::log(zref / z0);
405  ++iter;
406  } while ((std::abs(u_star_arr(i,j,k) - ustar) > tol) && iter <= max_iters);
407  AMREX_ALWAYS_ASSERT_WITH_MESSAGE(iter < max_iters,
408  "Maximum number of MOST iterations reached.");
409 
410  // Populate the stored MOST arrays
411  z0_arr(i,j,k) = z0;
412  olen_arr(i,j,k) = 1.0e16;
413  t_star_arr(i,j,k) = 0.0;
414  q_star_arr(i,j,k) = 0.0;
415  }
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real Donelan_roughness(amrex::Real ustar)
Definition: ERF_MOSTRoughness.H:29
amrex::Real Real
Definition: ERF_ShocInterface.H:19
const amrex::Real WSMIN
Definition: ERF_MOSTStress.H:421
const amrex::Real tol
Definition: ERF_MOSTStress.H:420
amrex::Real kappa
von Karman constant
Definition: ERF_MOSTStress.H:16
Here is the call graph for this function:

Member Data Documentation

◆ mdata

most_data adiabatic_donelan::mdata
private

Referenced by adiabatic_donelan(), and iterate_flux().

◆ sfuns

similarity_funs adiabatic_donelan::sfuns
private

◆ tol

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

Referenced by iterate_flux().

◆ WSMIN

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

Referenced by iterate_flux().


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