ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
saturation_funcs Namespace Reference

Functions

AMREX_GPU_HOST void compute_saturation_pressure_null (MultiFab &, const MultiFab &)
 
AMREX_GPU_HOST void compute_saturation_pressure_H2O (MultiFab &a_mf_sat_pressure, const MultiFab &a_mf_temperature)
 
AMREX_GPU_HOST void compute_saturation_vapfrac_null (MultiFab &, const MultiFab &)
 
AMREX_GPU_HOST void compute_saturation_vapfrac_H2O (MultiFab &a_mf_sat_vapfrac, const MultiFab &a_mf_temperature, const MultiFab &a_mf_pressure)
 
AMREX_GPU_HOST void compute_saturation_pressure_null (amrex::MultiFab &, const amrex::MultiFab &)
 Null function for saturation pressure. More...
 
AMREX_GPU_HOST void compute_saturation_pressure_H2O (amrex::MultiFab &, const amrex::MultiFab &)
 Compute saturation pressure for moist air. More...
 
AMREX_GPU_HOST void compute_saturation_vapfrac_null (amrex::MultiFab &, const amrex::MultiFab &, const amrex::MultiFab &)
 Null function for saturation vapour fraction. More...
 
AMREX_GPU_HOST void compute_saturation_vapfrac_H2O (amrex::MultiFab &, const amrex::MultiFab &, const amrex::MultiFab &)
 Compute saturation vapour fraction for moist air. More...
 

Detailed Description

Functions for computing saturation pressure and vapour fraction

Function Documentation

◆ compute_saturation_pressure_H2O() [1/2]

AMREX_GPU_HOST void saturation_funcs::compute_saturation_pressure_H2O ( amrex::MultiFab &  ,
const amrex::MultiFab &   
)

Compute saturation pressure for moist air.

◆ compute_saturation_pressure_H2O() [2/2]

AMREX_GPU_HOST void saturation_funcs::compute_saturation_pressure_H2O ( MultiFab &  a_mf_sat_pressure,
const MultiFab &  a_mf_temperature 
)
14  {
15  const auto& gvec = a_mf_sat_pressure.nGrowVect();
16  for (MFIter mfi(a_mf_sat_pressure, TilingIfNotGPU()); mfi.isValid(); ++mfi) {
17  Box bx = mfi.tilebox();
18  bx.grow(gvec);
19  const Array4<Real>& psat_arr = a_mf_sat_pressure.array(mfi);
20  const Array4<Real const>& temperature_arr = a_mf_temperature.array(mfi);
21 
22  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
23  { psat_arr(i,j,k,0) = erf_esatw(temperature_arr(i,j,k,0))*100; } );
24  // formula gives pressure in hPa; we will save it in Pa.
25  }
26  }
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real erf_esatw(amrex::Real t)
Definition: ERF_MicrophysicsUtils.H:68

Referenced by MaterialProperties::setProperties_H2O().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ compute_saturation_pressure_null() [1/2]

AMREX_GPU_HOST void saturation_funcs::compute_saturation_pressure_null ( amrex::MultiFab &  ,
const amrex::MultiFab &   
)

Null function for saturation pressure.

◆ compute_saturation_pressure_null() [2/2]

AMREX_GPU_HOST void saturation_funcs::compute_saturation_pressure_null ( MultiFab &  ,
const MultiFab &   
)
9 { }

◆ compute_saturation_vapfrac_H2O() [1/2]

AMREX_GPU_HOST void saturation_funcs::compute_saturation_vapfrac_H2O ( amrex::MultiFab &  ,
const amrex::MultiFab &  ,
const amrex::MultiFab &   
)

Compute saturation vapour fraction for moist air.

◆ compute_saturation_vapfrac_H2O() [2/2]

AMREX_GPU_HOST void saturation_funcs::compute_saturation_vapfrac_H2O ( MultiFab &  a_mf_sat_vapfrac,
const MultiFab &  a_mf_temperature,
const MultiFab &  a_mf_pressure 
)
35  {
36  const auto& gvec = a_mf_sat_vapfrac.nGrowVect();
37  for (MFIter mfi(a_mf_sat_vapfrac, TilingIfNotGPU()); mfi.isValid(); ++mfi) {
38  Box bx = mfi.tilebox();
39  bx.grow(gvec);
40  const Array4<Real>& qsat_arr = a_mf_sat_vapfrac.array(mfi);
41  const Array4<Real const>& temperature_arr = a_mf_temperature.array(mfi);
42  const Array4<Real const>& pressure_arr = a_mf_pressure.array(mfi);
43 
44  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
45  {
46  // pressure is in Pa; formula takes pressure in hPa
47  erf_qsatw( temperature_arr(i,j,k,0),
48  pressure_arr(i,j,k,0)/100.0,
49  qsat_arr(i,j,k,0) );
50  } );
51  }
52  }
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void erf_qsatw(amrex::Real t, amrex::Real p, amrex::Real &qsatw)
Definition: ERF_MicrophysicsUtils.H:166

Referenced by MaterialProperties::setProperties_H2O().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ compute_saturation_vapfrac_null() [1/2]

AMREX_GPU_HOST void saturation_funcs::compute_saturation_vapfrac_null ( amrex::MultiFab &  ,
const amrex::MultiFab &  ,
const amrex::MultiFab &   
)

Null function for saturation vapour fraction.

◆ compute_saturation_vapfrac_null() [2/2]

AMREX_GPU_HOST void saturation_funcs::compute_saturation_vapfrac_null ( MultiFab &  ,
const MultiFab &   
)
29 { }