ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
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-3
 
const amrex::Real tol_z = 1.0e-10
 
const amrex::Real alpha = 0.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
1204  {
1205  mdata.zref = zref;
1206  mdata.surf_temp_flux = Tflux;
1207  mdata.surf_moist_flux = Qvflux;
1208  mdata.Cnk_a = cnk_a;
1209  mdata.visc = cnk_visc;
1210  spec_qflux = cons_qflux;
1211  }
amrex::Real surf_moist_flux
Moisture flux.
Definition: ERF_MOSTStress.H:20
amrex::Real zref
Reference height (m)
Definition: ERF_MOSTStress.H:15
amrex::Real Cnk_a
Standard Charnock constant https://doi.org/10.1175/JAMC-D-17-0137.1.
Definition: ERF_MOSTStress.H:22
bool visc
Use viscous Charnock formulation.
Definition: ERF_MOSTStress.H:27
amrex::Real surf_temp_flux
Heat flux.
Definition: ERF_MOSTStress.H:19
most_data mdata
Definition: ERF_MOSTStress.H:1349
bool spec_qflux
Definition: ERF_MOSTStress.H:1350

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
1236  {
1237  amrex::Real Rib = 0.0;
1238  amrex::Real zeta = 0.0;
1239  amrex::Real zeta_old = 0.0;
1240  amrex::Real psi_m = 0.0;
1241  amrex::Real psi_h = 0.0;
1242  amrex::Real num = 0.0;
1243  amrex::Real den = 0.0;
1244  amrex::Real z0 = z0_arr(i,j,k);
1245  amrex::Real z0_old = z0;
1246  amrex::Real ustar = 0.0;
1247  amrex::Real umm = std::max(umm_arr(i,j,k), WSMIN);
1248  amrex::Real C = std::log(mdata.zref / z0);
1249 
1250  // First iteration we assume neutral (L -> inf)
1251  if (u_star_arr(i,j,k) == 1.E34) { olen_arr(i,j,k) = 1.0e3; }
1252  zeta = mdata.zref / olen_arr(i,j,k);
1253 
1254  // Water vapor in atmos and surface
1255  amrex::Real qv_s, qv_a;
1256  if (q_surf_arr(i,j,k) > 0.) {
1257  qv_s = q_surf_arr(i,j,k);
1258  } else {
1259  // First iteration and no qv_surf was specified
1260  // Use mean since there will be no flux
1261  qv_s = qvm_arr(i,j,k);
1262  }
1263  qv_a = qvm_arr(i,j,k);
1264 
1265  // update w* and Umagmean from Beljaars (1995)
1266  if (w_star_arr) {
1267  // NOTE: Thv flux is lagged, similar to WRF
1268  psi_m = sfuns.calc_psi_m2(zeta);
1269  psi_h = sfuns.calc_psi_h2(zeta);
1270  ustar = mdata.kappa * umm / (C - psi_m);
1271  amrex::Real tstar = mdata.kappa * (tm_arr(i,j,k) - t_surf_arr(i,j,k)) / (C - psi_h);
1273  -ustar * mdata.kappa * (qv_a - qv_s) / (C - psi_h);
1274  amrex::Real tflux = -ustar*tstar*(1 + 0.61*qvm_arr(i,j,k)) + 0.61*tm_arr(i,j,k)*qflux;
1275  w_star_arr(i,j,k) = calc_wstar(tflux, pblh_arr(i,j,k), tvm_arr(i,j,k));
1276  amrex::Real wstar = mdata.Bjr_beta * w_star_arr(i,j,k);
1277  umm = std::sqrt(umm_arr(i,j,k)*umm_arr(i,j,k) + wstar*wstar);
1278  umm = std::max(umm, WSMIN);
1279  }
1280 
1281  // Bulk Richardson number w/ moisture
1282  amrex::Real thv_s = t_surf_arr(i,j,k) * (1.0 + 0.61*qv_s);
1283  amrex::Real thv_a = tm_arr(i,j,k) * (1.0 + 0.61*qv_a);
1284  Rib = ( (mdata.gravity * mdata.zref) / tm_arr(i,j,k) ) *
1285  ( (thv_a - thv_s) / (umm * umm) );
1286  Rib = std::min(std::max(Rib,-4.0),4.0);
1287 
1288  // Fixed point iteration on zeta
1289  int iter = 0;
1290  do {
1291  // Transfer curr to old
1292  zeta_old = zeta;
1293 
1294  // Stability functions
1295  psi_m = sfuns.calc_psi_m2(zeta_old);
1296  psi_h = sfuns.calc_psi_h2(zeta_old);
1297 
1298  // Fixed point iteration on roughness
1299  int iter_z = 0;
1300  do {
1301  // Transfer curr to old
1302  z0_old = z0;
1303 
1304  // Update
1305  C = std::log(mdata.zref / z0_old);
1306  ustar = mdata.kappa * umm / (C - psi_m);
1307  if (mdata.Cnk_a > 0) {
1308  z0 = (mdata.Cnk_a / mdata.gravity) * ustar * ustar;
1309  if (mdata.visc) {
1310  z0 += air_viscosity(tm_arr(i,j,k)) / std::max(ustar, 0.05);
1311  }
1312  } else {
1313  z0 = COARE3_roughness(mdata.zref, umm, ustar);
1314  }
1315 
1316  ++iter_z;
1317  } while ( (std::abs(z0 - z0_old) > tol_z) && (iter_z <= max_iters) );
1318  AMREX_ALWAYS_ASSERT_WITH_MESSAGE(iter_z < max_iters,
1319  "Maximum number of MOST roughness iterations reached.");
1320  C = std::log(mdata.zref / z0);
1321 
1322  // Limiting
1323  num = std::max(C - psi_m, 1.0);
1324  den = std::max(C - psi_h, 1.0);
1325 
1326  // Update with under relaxation
1327  zeta = (1.0 - alpha) * zeta_old + alpha * Rib * num * num / den;
1328 
1329  ++iter;
1330  } while ( (std::abs(zeta - zeta_old) > tol) && (iter <= max_iters) );
1331  AMREX_ALWAYS_ASSERT_WITH_MESSAGE(iter < max_iters,
1332  "Maximum number of MOST iterations reached.");
1333 
1334  // Populate the stored MOST arrays
1335  z0_arr(i,j,k) = z0;
1336  olen_arr(i,j,k) = mdata.zref / zeta;
1337  u_star_arr(i,j,k) = mdata.kappa * umm / (C - psi_m);
1338  t_star_arr(i,j,k) = mdata.kappa * (tm_arr(i,j,k) - t_surf_arr(i,j,k)) / (C - psi_h);
1339  if (spec_qflux) {
1340  q_surf_arr(i,j,k) = mdata.surf_moist_flux * (C - psi_h) /
1341  (u_star_arr(i,j,k) * mdata.kappa) + qvm_arr(i,j,k);
1342  q_star_arr(i,j,k) = -mdata.surf_moist_flux / u_star_arr(i,j,k);
1343  } else {
1344  q_star_arr(i,j,k) = mdata.kappa * (qvm_arr(i,j,k) - q_surf_arr(i,j,k)) / (C - psi_h);
1345  }
1346  }
@ num
Definition: ERF_DataStruct.H:22
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real COARE3_roughness(amrex::Real zref, amrex::Real umm, amrex::Real ustar)
Definition: ERF_MOSTRoughness.H:9
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real air_viscosity(amrex::Real T_degK)
Definition: ERF_MOSTStress.H:126
amrex::Real Real
Definition: ERF_ShocInterface.H:16
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
amrex::Real kappa
von Karman constant
Definition: ERF_MOSTStress.H:17
amrex::Real gravity
Acceleration due to gravity (m/s^2)
Definition: ERF_MOSTStress.H:18
const amrex::Real Bjr_beta
Definition: ERF_MOSTStress.H:29
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real calc_psi_h2(amrex::Real zeta) const
Definition: ERF_MOSTStress.H:64
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real calc_psi_m2(amrex::Real zeta) const
Definition: ERF_MOSTStress.H:44
const amrex::Real tol_z
Definition: ERF_MOSTStress.H:1353
const amrex::Real alpha
Definition: ERF_MOSTStress.H:1354
const amrex::Real WSMIN
Definition: ERF_MOSTStress.H:1355
const amrex::Real tol
Definition: ERF_MOSTStress.H:1352
similarity_funs sfuns
Definition: ERF_MOSTStress.H:1351
Here is the call graph for this function:

Member Data Documentation

◆ alpha

const amrex::Real surface_temp_charnock::alpha = 0.5
private

Referenced by iterate_flux().

◆ mdata

most_data surface_temp_charnock::mdata
private

◆ sfuns

similarity_funs surface_temp_charnock::sfuns
private

Referenced by iterate_flux().

◆ spec_qflux

bool surface_temp_charnock::spec_qflux
private

◆ tol

const amrex::Real surface_temp_charnock::tol = 1.0e-3
private

Referenced by iterate_flux().

◆ tol_z

const amrex::Real surface_temp_charnock::tol_z = 1.0e-10
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: