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.
20  {
21  const auto& gvec = a_mf_sat_pressure.nGrowVect();
22  for (MFIter mfi(a_mf_sat_pressure, TilingIfNotGPU()); mfi.isValid(); ++mfi) {
23  Box bx = mfi.tilebox();
24  bx.grow(gvec);
25  const Array4<Real>& psat_arr = a_mf_sat_pressure.array(mfi);
26  const Array4<Real const>& temperature_arr = a_mf_temperature.array(mfi);
27 
28  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
29  { psat_arr(i,j,k,0) = erf_esatw(temperature_arr(i,j,k,0))*Real(100); } );
30  // formula gives pressure in hPa; we will save it in Pa.
31  }
32  }
erf_esatw
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real erf_esatw(amrex::Real t, bool use_empirical=false)
Definition: ERF_MicrophysicsUtils.H:123
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 
)
37  {
38  const auto& gvec = a_mf_sat_pressure.nGrowVect();
39  for (MFIter mfi(a_mf_sat_pressure, TilingIfNotGPU()); mfi.isValid(); ++mfi) {
40  Box bx = mfi.tilebox();
41  bx.grow(gvec);
42  const Array4<Real>& psat_arr = a_mf_sat_pressure.array(mfi);
43  const Array4<Real const>& temperature_arr = a_mf_temperature.array(mfi);
44 
45  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
46  { psat_arr(i,j,k,0) = erf_esati(temperature_arr(i,j,k,0))*Real(100); } );
47  // formula gives pressure in hPa; we will save it in Pa.
48  }
49  }
erf_esati
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real erf_esati(amrex::Real t)
Definition: ERF_MicrophysicsUtils.H:31

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

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.
65  {
66  const auto& gvec = a_mf_sat_vapfrac.nGrowVect();
67  for (MFIter mfi(a_mf_sat_vapfrac, TilingIfNotGPU()); mfi.isValid(); ++mfi) {
68  Box bx = mfi.tilebox();
69  bx.grow(gvec);
70  const Array4<Real>& qsat_arr = a_mf_sat_vapfrac.array(mfi);
71  const Array4<Real const>& temperature_arr = a_mf_temperature.array(mfi);
72  const Array4<Real const>& pressure_arr = a_mf_pressure.array(mfi);
73 
74  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
75  {
76  // pressure is in Pa; formula takes pressure in hPa
77  erf_qsatw( temperature_arr(i,j,k,0),
78  pressure_arr(i,j,k,0)/Real(100.0),
79  qsat_arr(i,j,k,0) );
80  } );
81  }
82  }
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:228

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 
)
88  {
89  const auto& gvec = a_mf_sat_vapfrac.nGrowVect();
90  for (MFIter mfi(a_mf_sat_vapfrac, TilingIfNotGPU()); mfi.isValid(); ++mfi) {
91  Box bx = mfi.tilebox();
92  bx.grow(gvec);
93  const Array4<Real>& qsat_arr = a_mf_sat_vapfrac.array(mfi);
94  const Array4<Real const>& temperature_arr = a_mf_temperature.array(mfi);
95  const Array4<Real const>& pressure_arr = a_mf_pressure.array(mfi);
96 
97  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
98  {
99  // pressure is in Pa; formula takes pressure in hPa
100  erf_qsati( temperature_arr(i,j,k,0),
101  pressure_arr(i,j,k,0)/amrex::Real(100.0),
102  qsat_arr(i,j,k,0) );
103  } );
104  }
105  }
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:218

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