ERF_WDM6.H File Reference

#include <memory>
#include <string>
#include <vector>
#include <AMReX_FArrayBox.H>
#include <AMReX_Geometry.H>
#include <AMReX_MultiFabUtil.H>
#include "ERF_Constants.H"
#include "ERF_MicrophysicsConstants.H"
#include "ERF_DataStruct.H"
#include "ERF_IndexDefines.H"
#include "ERF_MicrophysicsColumns.H"
#include "ERF_MicrophysicsUtils.H"
#include "ERF_NullMoist.H"
#include "ERF_SurfacePrecipitation.H"

Include dependency graph for ERF_WDM6.H:

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

Go to the source code of this file.

Classes
class	WDM6

Namespaces
	MicVar_WDM6
	WDM6Sed

Enumerations
enum	{   MicVar_WDM6::rho = 0 , MicVar_WDM6::theta , MicVar_WDM6::tabs , MicVar_WDM6::pres ,   MicVar_WDM6::qv , MicVar_WDM6::qc , MicVar_WDM6::qi , MicVar_WDM6::qr ,   MicVar_WDM6::qs , MicVar_WDM6::qg , MicVar_WDM6::nn , MicVar_WDM6::nc ,   MicVar_WDM6::nr , MicVar_WDM6::rain_accum , MicVar_WDM6::snow_accum , MicVar_WDM6::graupel_accum ,   MicVar_WDM6::precip_accum , MicVar_WDM6::NumVars }
enum	{   WDM6Sed::qi = 0 , WDM6Sed::qr , WDM6Sed::qs , WDM6Sed::qg ,   WDM6Sed::nr , WDM6Sed::NumComps }

Functions
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void	wdm6_copy_state_to_micro_cell (const amrex::Array4< const amrex::Real > &states, const amrex::Array4< const amrex::Real > &base, const amrex::Array4< amrex::Real > &rho, const amrex::Array4< amrex::Real > &theta, const amrex::Array4< amrex::Real > &tabs, const amrex::Array4< amrex::Real > &pres, const amrex::Array4< amrex::Real > &qv, const amrex::Array4< amrex::Real > &qc, const amrex::Array4< amrex::Real > &qi, const amrex::Array4< amrex::Real > &qr, const amrex::Array4< amrex::Real > &qs, const amrex::Array4< amrex::Real > &qg, const amrex::Real rdOcp, const bool use_anelastic_reference_pressure, const int i, const int j, const int k) noexcept
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real	wdm6_theta_from_temperature_and_pressure (const amrex::Real temperature, const amrex::Real pressure, const amrex::Real configured_rdOcp, const bool use_anelastic_reference_pressure) noexcept
constexpr amrex::Real	wdm6_literal (double d)

Function Documentation

wdm6_copy_state_to_micro_cell()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void wdm6_copy_state_to_micro_cell ( const amrex::Array4< const amrex::Real > &  states, const amrex::Array4< const amrex::Real > &  base, const amrex::Array4< amrex::Real > &  rho, const amrex::Array4< amrex::Real > &  theta, const amrex::Array4< amrex::Real > &  tabs, const amrex::Array4< amrex::Real > &  pres, const amrex::Array4< amrex::Real > &  qv, const amrex::Array4< amrex::Real > &  qc, const amrex::Array4< amrex::Real > &  qi, const amrex::Array4< amrex::Real > &  qr, const amrex::Array4< amrex::Real > &  qs, const amrex::Array4< amrex::Real > &  qg, const amrex::Real  rdOcp, const bool  use_anelastic_reference_pressure, const int  i, const int  j, const int  k  )

noexcept

{
    const amrex::Real rho_value = states(i,j,k,Rho_comp);
    const amrex::Real rho_theta = states(i,j,k,RhoTheta_comp);
    rho(i,j,k) = rho_value;
    theta(i,j,k) = rho_theta / rho_value;
    qv(i,j,k) = amrex::max(amrex::Real(0.0), states(i,j,k,RhoQ1_comp) / rho_value);
    qc(i,j,k) = amrex::max(amrex::Real(0.0), states(i,j,k,RhoQ2_comp) / rho_value);
    qi(i,j,k) = amrex::max(amrex::Real(0.0), states(i,j,k,RhoQ3_comp) / rho_value);
    qr(i,j,k) = amrex::max(amrex::Real(0.0), states(i,j,k,RhoQ4_comp) / rho_value);
    qs(i,j,k) = amrex::max(amrex::Real(0.0), states(i,j,k,RhoQ5_comp) / rho_value);
    qg(i,j,k) = amrex::max(amrex::Real(0.0), states(i,j,k,RhoQ6_comp) / rho_value);
 
    const amrex::Real p0 = use_anelastic_reference_pressure
        ? base(i,j,k,BaseState::p0_comp) : amrex::Real(0.0);
    const MicrophysicsThermoState thermo = diagnose_microphysics_thermo_state(
        rho_value, rho_theta, qv(i,j,k), rdOcp,
        use_anelastic_reference_pressure, p0);
    tabs(i,j,k) = thermo.temperature;
    pres(i,j,k) = thermo.pressure_pa;
}

Here is the call graph for this function:

wdm6_literal()

constexpr amrex::Real wdm6_literal ( double  d )

constexpr

{
#ifdef ERF_WDM6_F32_LITERALS
    // Match the Fortran as compiled today: its literals carry no kind suffix,
    // so they are evaluated in single precision and then widened.
    return amrex::Real(static_cast<float>(d));
#else
    // True doubles. Correct when the Fortran is built -fdefault-real-8.
    return amrex::Real(d);
#endif
}

Referenced by WDM6::AdvanceColumns(), WDM6::initialize_coeffs(), wdm6_diffus(), wdm6_lamdac_exact(), wdm6_lamdar(), wdm6_nislfv_rain_plm6_column(), wdm6_slope_rain_cell(), wdm6_venfac(), and wdm6_viscos().

Here is the caller graph for this function:

wdm6_theta_from_temperature_and_pressure()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real wdm6_theta_from_temperature_and_pressure ( const amrex::Real  temperature, const amrex::Real  pressure, const amrex::Real  configured_rdOcp, const bool  use_anelastic_reference_pressure  )

noexcept

{
    const amrex::Real theta_rdOcp = use_anelastic_reference_pressure
        ? configured_rdOcp : RdoCp;
    return getThgivenTandP(temperature, pressure, theta_rdOcp);
}

Referenced by WDM6::AdvanceColumns().

Here is the call graph for this function:

Here is the caller graph for this function: