ERF_Kessler.H File Reference

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

Include dependency graph for ERF_Kessler.H:

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

Go to the source code of this file.

Classes
class	Kessler

Namespaces
	MicVar_Kess

Enumerations
enum	{   MicVar_Kess::rho =0 , MicVar_Kess::theta , MicVar_Kess::tabs , MicVar_Kess::pres ,   MicVar_Kess::qt , MicVar_Kess::qv , MicVar_Kess::qcl , MicVar_Kess::qp ,   MicVar_Kess::rain_accum , MicVar_Kess::NumVars }

Functions
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void	kessler_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 > &qv, const amrex::Array4< amrex::Real > &qc, const amrex::Array4< amrex::Real > &qp, const amrex::Array4< amrex::Real > &qt, const amrex::Array4< amrex::Real > &tabs, const amrex::Array4< amrex::Real > &pres, const amrex::Real rdOcp, const bool use_anelastic_reference_pressure, const int i, const int j, const int k) noexcept

Function Documentation

kessler_copy_state_to_micro_cell()

AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void kessler_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 > &  qv, const amrex::Array4< amrex::Real > &  qc, const amrex::Array4< amrex::Real > &  qp, const amrex::Array4< amrex::Real > &  qt, const amrex::Array4< amrex::Real > &  tabs, const amrex::Array4< amrex::Real > &  pres, 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);
    const amrex::Real qv_value = states(i,j,k,RhoQ1_comp) / rho_value;
 
    rho(i,j,k) = rho_value;
    theta(i,j,k) = rho_theta / rho_value;
    qv(i,j,k) = qv_value;
    qc(i,j,k) = states(i,j,k,RhoQ2_comp) / rho_value;
    qp(i,j,k) = states(i,j,k,RhoQ3_comp) / rho_value;
    qt(i,j,k) = qv(i,j,k) + qc(i,j,k);
 
    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_value, rdOcp,
        use_anelastic_reference_pressure, p0);
    tabs(i,j,k) = thermo.temperature;
    // Kessler's saturation helpers use mbar / hPa.
    pres(i,j,k) = amrex::Real(0.01) * thermo.pressure_pa;
}

Here is the call graph for this function: