ERF_MOSTRoughness.H File Reference

#include <ERF_Constants.H>

Include dependency graph for ERF_MOSTRoughness.H:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	COARE3_roughness (amrex::Real zref, amrex::Real umm, amrex::Real ustar)
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	Donelan_roughness (amrex::Real ustar)

Function Documentation

COARE3_roughness()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real COARE3_roughness	(	amrex::Real	zref,
		amrex::Real	umm,
		amrex::Real	ustar
	)

COARE3.0, following WRF phys/module_sf_mynn.F

{
    amrex::Real ws10m = umm*std::log(10./1e-4) / std::log(zref/1e-4);
    amrex::Real Czc = 0.011 + 0.007*std::min(std::max((ws10m-10.)/8., 0.), 1.0);
    amrex::Real z0 = (Czc / CONST_GRAV) * ustar * ustar;
    // From Davis et al 2008 MWR, but their lower limit is 0.125e-6 -- WRF typo?
    return std::min(std::max(z0, 1.27e-7), 2.85e-3);
}

Referenced by adiabatic_charnock::iterate_flux(), surface_flux_charnock::iterate_flux(), and surface_temp_charnock::iterate_flux().

Here is the caller graph for this function:

Donelan_roughness()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real Donelan_roughness

(

amrex::Real

ustar

)

Updated Donelan 2004 as found in AHW

Matches implementations in WRF phys/physics_mmm/sf_sfclayrev.F90 and phys/module_sf_mynn.F

{
    constexpr amrex::Real ozo{1.59e-5};
    amrex::Real zw  = std::min(std::pow(ustar/1.06, 0.3),1.0);
    amrex::Real zn1 = 0.011*ustar*ustar/CONST_GRAV + ozo;
    amrex::Real zn2 = 10.*std::exp(-9.5*std::pow(ustar, -.3333))
                    + 0.11*1.5e-5/std::max(ustar,0.01);
    amrex::Real z0 = (1.0-zw) * zn1 + zw * zn2;
    // From Davis et al 2008 MWR, but their lower limit is 0.125e-6 -- WRF typo?
    return std::min(std::max(z0, 1.27e-7), 2.85e-3);
}

Referenced by adiabatic_donelan::iterate_flux(), surface_flux_donelan::iterate_flux(), and surface_temp_donelan::iterate_flux().

Here is the caller graph for this function: