ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
Public Member Functions | Private Attributes | List of all members
adiabatic_wave_coupled Struct Reference

#include <ERF_MOSTStress.H>

Collaboration diagram for adiabatic_wave_coupled:
Collaboration graph
[legend]

Public Member Functions

 adiabatic_wave_coupled (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 > &Hwave_arr, const amrex::Array4< amrex::Real > &Lwave_arr, const amrex::Array4< amrex::Real > &eta_arr) const
 

Private Attributes

most_data mdata
 
similarity_funs sfuns
 
const amrex::Real tol = 1.0e-5
 
const amrex::Real eps = 1e-15
 
const amrex::Real z0_eps = 1.0e-6
 
const amrex::Real z0_max = 0.1
 
const amrex::Real WSMIN = 0.1
 

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 
)
inline
432  {
433  mdata.surf_temp_flux = Tflux;
434  mdata.surf_moist_flux = Qvflux;
435  }
adiabatic_wave_coupled::mdata
most_data mdata
Definition: ERF_MOSTStress.H:491
most_data::surf_moist_flux
amrex::Real surf_moist_flux
Moisture flux.
Definition: ERF_MOSTStress.H:19
most_data::surf_temp_flux
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_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 > &  ,
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 > &  eta_arr 
) const
inline
461  {
462  int iter = 0;
463  amrex::Real umm = std::max(umm_arr(i,j,k), WSMIN);
464  amrex::Real ustar = 0.0;
465  amrex::Real zref = zref_arr(i,j,k);
466  amrex::Real z0 = 0.0;
467  int ie, je;
468  ie = i < lbound(eta_arr).x ? lbound(eta_arr).x : i;
469  je = j < lbound(eta_arr).y ? lbound(eta_arr).y : j;
470  ie = ie > ubound(eta_arr).x ? ubound(eta_arr).x : ie;
471  je = je > ubound(eta_arr).y ? ubound(eta_arr).y : je;
472  u_star_arr(i,j,k) = mdata.kappa * umm / std::log(zref / z0_arr(i,j,k));
473  do {
474  ustar = u_star_arr(i,j,k);
475  z0 = std::min( std::max(1200.0 * Hwave_arr(i,j,k) * std::pow( Hwave_arr(i,j,k)/(Lwave_arr(i,j,k)+eps), 4.5 )
476  + 0.11 * eta_arr(ie,je,k,EddyDiff::Mom_v) / ustar, z0_eps), z0_max );
477  u_star_arr(i,j,k) = mdata.kappa * umm / std::log(zref / z0);
478  ++iter;
479  } while ((std::abs(u_star_arr(i,j,k) - ustar) > tol) && iter <= max_iters);
480  AMREX_ALWAYS_ASSERT_WITH_MESSAGE(iter < max_iters,
481  "Maximum number of MOST iterations reached.");
482 
483  // Populate the stored MOST arrays
484  z0_arr(i,j,k) = z0;
485  olen_arr(i,j,k) = 1.0e16;
486  t_star_arr(i,j,k) = 0.0;
487  q_star_arr(i,j,k) = 0.0;
488  }
Real
amrex::Real Real
Definition: ERF_ShocInterface.H:19
EddyDiff::Mom_v
@ Mom_v
Definition: ERF_IndexDefines.H:175
adiabatic_wave_coupled::eps
const amrex::Real eps
Definition: ERF_MOSTStress.H:494
adiabatic_wave_coupled::z0_eps
const amrex::Real z0_eps
Definition: ERF_MOSTStress.H:495
adiabatic_wave_coupled::WSMIN
const amrex::Real WSMIN
Definition: ERF_MOSTStress.H:497
adiabatic_wave_coupled::tol
const amrex::Real tol
Definition: ERF_MOSTStress.H:493
adiabatic_wave_coupled::z0_max
const amrex::Real z0_max
Definition: ERF_MOSTStress.H:496
most_data::kappa
amrex::Real kappa
von Karman constant
Definition: ERF_MOSTStress.H:16

Member Data Documentation

◆ eps

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

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

Referenced by iterate_flux().

◆ WSMIN

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

Referenced by iterate_flux().

◆ z0_eps

const amrex::Real adiabatic_wave_coupled::z0_eps = 1.0e-6
private

Referenced by iterate_flux().

◆ z0_max

const amrex::Real adiabatic_wave_coupled::z0_max = 0.1
private

Referenced by iterate_flux().

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