ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
Functions
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_pressure_ice (MultiFab &a_mf_sat_pressure, const MultiFab &a_mf_temperature)
 
AMREX_GPU_HOST void compute_saturation_vapfrac_null (MultiFab &, const 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_vapfrac_ice (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_pressure_ice (amrex::MultiFab &, const amrex::MultiFab &)
 Compute saturation pressure for icy 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...
 
AMREX_GPU_HOST void compute_saturation_vapfrac_ice (amrex::MultiFab &, const amrex::MultiFab &, const amrex::MultiFab &)
 Compute saturation vapour fraction for icy 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 
)

Compute saturation pressure for H2O.

Parameters
[out]a_mf_sat_pressureMultiFab to store saturation pressure.
[in]a_mf_temperatureMultiFab containing temperature.
22  {
23  const auto& gvec = a_mf_sat_pressure.nGrowVect();
24  for (MFIter mfi(a_mf_sat_pressure, TilingIfNotGPU()); mfi.isValid(); ++mfi) {
25  Box bx = mfi.tilebox();
26  bx.grow(gvec);
27  const Array4<Real>& psat_arr = a_mf_sat_pressure.array(mfi);
28  const Array4<Real const>& temperature_arr = a_mf_temperature.array(mfi);
29 
30  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
31  { psat_arr(i,j,k,0) = erf_esatw(temperature_arr(i,j,k,0))*Real(100); } );
32  // formula gives pressure in hPa; we will save it in Pa.
33  }
34  }
erf_esatw
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real erf_esatw(amrex::Real t, bool use_empirical=false)
Definition: ERF_MicrophysicsUtils.H:159
ParallelFor
ParallelFor(fab_box, [=] AMREX_GPU_DEVICE(int i, int j, int k) { qrcuten_arr(i, j, k)=Real(0);qscuten_arr(i, j, k)=Real(0);qicuten_arr(i, j, k)=Real(0);})
Real
amrex::Real Real
Definition: ERF_ShocInterface.H:19

Referenced by MaterialProperties::setProperties_H2O().

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

◆ compute_saturation_pressure_ice() [1/2]

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

Compute saturation pressure for icy air.

◆ compute_saturation_pressure_ice() [2/2]

AMREX_GPU_HOST void saturation_funcs::compute_saturation_pressure_ice ( MultiFab &  a_mf_sat_pressure,
const MultiFab &  a_mf_temperature 
)
39  {
40  const auto& gvec = a_mf_sat_pressure.nGrowVect();
41  for (MFIter mfi(a_mf_sat_pressure, TilingIfNotGPU()); mfi.isValid(); ++mfi) {
42  Box bx = mfi.tilebox();
43  bx.grow(gvec);
44  const Array4<Real>& psat_arr = a_mf_sat_pressure.array(mfi);
45  const Array4<Real const>& temperature_arr = a_mf_temperature.array(mfi);
46 
47  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
48  { psat_arr(i,j,k,0) = erf_esati(temperature_arr(i,j,k,0))*Real(100); } );
49  // formula gives pressure in hPa; we will save it in Pa.
50  }
51  }
erf_esati
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real erf_esati(amrex::Real t)
Definition: ERF_MicrophysicsUtils.H:67

Referenced by MaterialProperties::setProperties_ice().

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 &   
)
11 { }

◆ 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 
)

Compute saturation vapor fraction for H2O.

Parameters
[out]a_mf_sat_vapfracMultiFab to store saturation vapor fraction.
[in]a_mf_temperatureMultiFab containing temperature.
[in]a_mf_pressureMultiFab containing pressure.
67  {
68  const auto& gvec = a_mf_sat_vapfrac.nGrowVect();
69  for (MFIter mfi(a_mf_sat_vapfrac, TilingIfNotGPU()); mfi.isValid(); ++mfi) {
70  Box bx = mfi.tilebox();
71  bx.grow(gvec);
72  const Array4<Real>& qsat_arr = a_mf_sat_vapfrac.array(mfi);
73  const Array4<Real const>& temperature_arr = a_mf_temperature.array(mfi);
74  const Array4<Real const>& pressure_arr = a_mf_pressure.array(mfi);
75 
76  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
77  {
78  // pressure is in Pa; formula takes pressure in hPa
79  erf_qsatw( temperature_arr(i,j,k,0),
80  pressure_arr(i,j,k,0)/Real(100.0),
81  qsat_arr(i,j,k,0) );
82  } );
83  }
84  }
erf_qsatw
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void erf_qsatw(amrex::Real t, amrex::Real p, amrex::Real &qsatw)
Definition: ERF_MicrophysicsUtils.H:264

Referenced by MaterialProperties::setProperties_H2O().

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

◆ compute_saturation_vapfrac_ice() [1/2]

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

Compute saturation vapour fraction for icy air.

◆ compute_saturation_vapfrac_ice() [2/2]

AMREX_GPU_HOST void saturation_funcs::compute_saturation_vapfrac_ice ( MultiFab &  a_mf_sat_vapfrac,
const MultiFab &  a_mf_temperature,
const MultiFab &  a_mf_pressure 
)
90  {
91  const auto& gvec = a_mf_sat_vapfrac.nGrowVect();
92  for (MFIter mfi(a_mf_sat_vapfrac, TilingIfNotGPU()); mfi.isValid(); ++mfi) {
93  Box bx = mfi.tilebox();
94  bx.grow(gvec);
95  const Array4<Real>& qsat_arr = a_mf_sat_vapfrac.array(mfi);
96  const Array4<Real const>& temperature_arr = a_mf_temperature.array(mfi);
97  const Array4<Real const>& pressure_arr = a_mf_pressure.array(mfi);
98 
99  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
100  {
101  // pressure is in Pa; formula takes pressure in hPa
102  erf_qsati( temperature_arr(i,j,k,0),
103  pressure_arr(i,j,k,0)/amrex::Real(100.0),
104  qsat_arr(i,j,k,0) );
105  } );
106  }
107  }
erf_qsati
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void erf_qsati(amrex::Real t, amrex::Real p, amrex::Real &qsati)
Definition: ERF_MicrophysicsUtils.H:254

Referenced by MaterialProperties::setProperties_ice().

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 &  ,
const MultiFab &   
)
54 { }