SatMethods Class Reference

Methods for estimating the saturation vapor pressure of water. More...

#include <ERF_SatMethods.H>

Static Public Member Functions
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real	wv_sat_svp_to_qsat (const amrex::Real &es, const amrex::Real &p)
	Get saturation specific humidity given pressure and SVP. More...
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE void	wv_sat_qsat_water (const amrex::Real &t, const amrex::Real &p, amrex::Real &es, amrex::Real &qs, const int idx=1)
	Calculate SVP over water and saturation specific humidity. More...
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE void	wv_sat_qsat_ice (const amrex::Real &t, const amrex::Real &p, amrex::Real &es, amrex::Real &qs, const int idx=1)
	Calculate SVP over ice and saturation specific humidity. More...
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE void	wv_sat_qsat_trans (const amrex::Real &t, const amrex::Real &p, amrex::Real &es, amrex::Real &qs, const int idx=1)
	Calculate transition SVP and saturation specific humidity. More...
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real	wv_sat_svp_water (const amrex::Real &t, const int idx=1)
	Dispatcher for saturation vapor pressure over water. More...
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real	wv_sat_svp_ice (const amrex::Real &t, const int idx=1)
	Dispatcher for saturation vapor pressure over ice. More...
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real	wv_sat_svp_trans (const amrex::Real &t, const int idx=1)
	Calculate transition saturation vapor pressure. More...
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real	GoffGratch_svp_water (const amrex::Real &t)
	Goff & Gratch (1946) saturation vapor pressure over water. More...
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real	GoffGratch_svp_ice (const amrex::Real &t)
	Goff & Gratch (1946) saturation vapor pressure over ice. More...
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real	MurphyKoop_svp_water (const amrex::Real &t)
	Murphy & Koop (2005) saturation vapor pressure over water. More...
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real	MurphyKoop_svp_ice (const amrex::Real &t)
	Murphy & Koop (2005) saturation vapor pressure over ice. More...
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real	OldGoffGratch_svp_water (const amrex::Real &t)
	Old CAM implementation of Goff & Gratch (1946) SVP over water. More...
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real	OldGoffGratch_svp_ice (const amrex::Real &t)
	Old CAM implementation of Goff & Gratch (1946) SVP over ice. More...
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real	Bolton_svp_water (const amrex::Real &t)
	Bolton (1980) saturation vapor pressure over water. More...

Private Types
enum	Type {   Invalid = -1 , OldGoffGratch = 0 , GoffGratch = 1 , MurphyKoop = 2 ,   Bolton = 3 }

Detailed Description

Methods for estimating the saturation vapor pressure of water.

Member Enumeration Documentation

Type

enum SatMethods::Type

private

Enumerator
Invalid
OldGoffGratch
GoffGratch
MurphyKoop
Bolton

              {
               Invalid = -1,
               OldGoffGratch = 0,
               GoffGratch = 1,
               MurphyKoop = 2,
               Bolton = 3
    };

Member Function Documentation

Bolton_svp_water()

AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real SatMethods::Bolton_svp_water ( const amrex::Real &  t )

inlinestatic

Bolton (1980) saturation vapor pressure over water.

Parameters

[in]

t

Temperature.

Returns

Saturation vapor pressure.

                                                           {
        constexpr auto c1 = amrex::Real(611.2);
        constexpr auto c2 = amrex::Real(17.67);
        constexpr auto c3 = amrex::Real(243.5);
        return c1*exp( (c2*(t - tmelt))/((t - tmelt)+c3) );
    }

Referenced by wv_sat_svp_ice(), and wv_sat_svp_water().

Here is the caller graph for this function:

GoffGratch_svp_ice()

AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real SatMethods::GoffGratch_svp_ice ( const amrex::Real &  t )

inlinestatic

Goff & Gratch (1946) saturation vapor pressure over ice.

Parameters

[in]

t

Temperature.

Returns

Saturation vapor pressure.

                                                             {
        // good down to -100 C
        return pow(amrex::Real(10.), (-amrex::Real(9.09718)*(h2otrip/t-amrex::Real(1))-amrex::Real(3.56654)*
                         log10(h2otrip/t)+amrex::Real(0.876793)*(amrex::Real(1)-t/h2otrip)+
                         log10(amrex::Real(6.1071))))*amrex::Real(100.);
    }

Referenced by wv_sat_svp_ice().

Here is the caller graph for this function:

GoffGratch_svp_water()

AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real SatMethods::GoffGratch_svp_water ( const amrex::Real &  t )

inlinestatic

Goff & Gratch (1946) saturation vapor pressure over water.

Parameters

[in]

t

Temperature.

Returns

Saturation vapor pressure.

                                                               {
        // uncertain below -70 C
        return pow(amrex::Real(10.), (-amrex::Real(7.90298)*(tboil/t-amrex::Real(1))+
                         amrex::Real(5.02808)*std::log10(tboil/t)-
                         amrex::Real(1.3816e-7)*(pow(amrex::Real(10.), (amrex::Real(11.344)*(amrex::Real(1)-t/tboil)))-amrex::Real(1))+
                         amrex::Real(8.1328e-3)*(pow(amrex::Real(10.), (-amrex::Real(3.49149)*(tboil/t-amrex::Real(1))))-amrex::Real(1))+
                         std::log10(amrex::Real(1013.246))))*amrex::Real(100.);
    }

Referenced by wv_sat_svp_water().

Here is the caller graph for this function:

MurphyKoop_svp_ice()

AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real SatMethods::MurphyKoop_svp_ice ( const amrex::Real &  t )

inlinestatic

Murphy & Koop (2005) saturation vapor pressure over ice.

Parameters

[in]

t

Temperature.

Returns

Saturation vapor pressure.

                                                             {
        // (good down to 110 K)
        return exp(amrex::Real(9.550426) - (amrex::Real(5723.265) / t) + (amrex::Real(3.53068) * log(t))
                   - (amrex::Real(0.00728332) * t));
    }

Referenced by wv_sat_svp_ice().

Here is the caller graph for this function:

MurphyKoop_svp_water()

AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real SatMethods::MurphyKoop_svp_water ( const amrex::Real &  t )

inlinestatic

Murphy & Koop (2005) saturation vapor pressure over water.

Parameters

[in]

t

Temperature.

Returns

Saturation vapor pressure.

                                                               {
        // (good for 123 < T < 332 K)
        return exp(amrex::Real(54.842763) - (amrex::Real(6763.22) / t) - (amrex::Real(4.210) * log(t)) +
                   (amrex::Real(0.000367) * t) + (tanh(amrex::Real(0.0415) * (t - amrex::Real(218.8))) *
                                     (amrex::Real(53.878) - (amrex::Real(1331.22) / t) - (amrex::Real(9.44523) * log(t)) +
                                      amrex::Real(0.014025) * t)));
    }

Referenced by wv_sat_svp_water().

Here is the caller graph for this function:

OldGoffGratch_svp_ice()

AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real SatMethods::OldGoffGratch_svp_ice ( const amrex::Real &  t )

inlinestatic

Old CAM implementation of Goff & Gratch (1946) SVP over ice.

Parameters

[in]

t

Temperature.

Returns

Saturation vapor pressure.

                                                                {
        auto term1 = amrex::Real(2.01889049)/(tmelt/t);
        auto term2 = amrex::Real(3.56654)*log(tmelt/t);
        auto term3 = amrex::Real(20.947031)*(tmelt/t);
        return amrex::Real(575.185606e10)*exp(-(term1 + term2 + term3));
    }

Referenced by wv_sat_svp_ice().

Here is the caller graph for this function:

OldGoffGratch_svp_water()

AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real SatMethods::OldGoffGratch_svp_water ( const amrex::Real &  t )

inlinestatic

Old CAM implementation of Goff & Gratch (1946) SVP over water.

Parameters

[in]

t

Temperature.

Returns

Saturation vapor pressure.

                                                                  {
        auto ps = amrex::Real(1013.246);
        auto e1 = amrex::Real(11.344)*(amrex::Real(1) - t/tboil);
        auto e2 = -amrex::Real(3.49149)*(tboil/t - amrex::Real(1));
        auto f1 = -amrex::Real(7.90298)*(tboil/t - amrex::Real(1));
        auto f2 = amrex::Real(5.02808)*log10(tboil/t);
        auto f3 = -amrex::Real(1.3816)*(pow(amrex::Real(10.0), e1) - amrex::Real(1))/amrex::Real(10000000.0);
        auto f4 = amrex::Real(8.1328)*(pow(amrex::Real(10.0), e2) - amrex::Real(1))/amrex::Real(1000.0);
        auto f5 = log10(ps);
        auto f  = f1 + f2 + f3 + f4 + f5;
        return (pow(amrex::Real(10.0), f))*amrex::Real(100.0);
    }

Referenced by wv_sat_svp_water().

Here is the caller graph for this function:

wv_sat_qsat_ice()

AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE void SatMethods::wv_sat_qsat_ice ( const amrex::Real &  t, const amrex::Real &  p, amrex::Real &  es, amrex::Real &  qs, const int  idx = 1  )

inlinestatic

Calculate SVP over ice and saturation specific humidity.

Parameters

[in]	t	Temperature.
[in]	p	Total pressure.
[out]	es	Saturation vapor pressure.
[out]	qs	Saturation specific humidity.
[in]	idx	Scheme index.

                                                                                  {
        // Purpose:
        //   Calculate SVP over ice at a given temperature, and then
        //   calculate and return saturation specific humidity.
 
        es = wv_sat_svp_ice(t, idx);
 
        qs = wv_sat_svp_to_qsat(es, p);
 
        // Ensures returned es is consistent with limiters on qs.
        es = std::min(es, p);
    }

Here is the call graph for this function:

wv_sat_qsat_trans()

AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE void SatMethods::wv_sat_qsat_trans ( const amrex::Real &  t, const amrex::Real &  p, amrex::Real &  es, amrex::Real &  qs, const int  idx = 1  )

inlinestatic

Calculate transition SVP and saturation specific humidity.

Parameters

[in]	t	Temperature.
[in]	p	Total pressure.
[out]	es	Saturation vapor pressure.
[out]	qs	Saturation specific humidity.
[in]	idx	Scheme index.

                                                                                    {
        // Purpose:
        //   Calculate SVP over ice at a given temperature, and then
        //   calculate and return saturation specific humidity.
 
        es = wv_sat_svp_trans(t, idx);
 
        qs = wv_sat_svp_to_qsat(es, p);
 
        // Ensures returned es is consistent with limiters on qs.
        es = std::min(es, p);
    }

Here is the call graph for this function:

wv_sat_qsat_water()

AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE void SatMethods::wv_sat_qsat_water ( const amrex::Real &  t, const amrex::Real &  p, amrex::Real &  es, amrex::Real &  qs, const int  idx = 1  )

inlinestatic

Calculate SVP over water and saturation specific humidity.

Parameters

[in]	t	Temperature.
[in]	p	Total pressure.
[out]	es	Saturation vapor pressure.
[out]	qs	Saturation specific humidity.
[in]	idx	Scheme index.

                                                                                    {
        // Purpose:
        //   Calculate SVP over water at a given temperature, and then
        //   calculate and return saturation specific humidity.
 
        es = wv_sat_svp_water(t, idx);
 
        qs = wv_sat_svp_to_qsat(es, p);
 
        // Ensures returned es is consistent with limiters on qs.
        es = std::min(es, p);
    }

Here is the call graph for this function:

wv_sat_svp_ice()

AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real SatMethods::wv_sat_svp_ice ( const amrex::Real &  t, const int  idx = 1  )

inlinestatic

Dispatcher for saturation vapor pressure over ice.

Parameters

[in]	t	Temperature.
[in]	idx	Scheme index.

Returns

Saturation vapor pressure over ice.

                                                                            {
        amrex::Real es;
        switch(idx)
            {
            case GoffGratch:
                es = GoffGratch_svp_ice(t); break;
            case MurphyKoop:
                es = MurphyKoop_svp_ice(t); break;
            case OldGoffGratch:
                es = OldGoffGratch_svp_ice(t); break;
            case Bolton:
                es = Bolton_svp_water(t); break;
            }
        return es;
    }

Referenced by wv_sat_qsat_ice(), and wv_sat_svp_trans().

Here is the call graph for this function:

Here is the caller graph for this function:

wv_sat_svp_to_qsat()

AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real SatMethods::wv_sat_svp_to_qsat ( const amrex::Real &  es, const amrex::Real &  p  )

inlinestatic

Get saturation specific humidity given pressure and SVP.

Parameters

[in]	es	Saturation vapor pressure.
[in]	p	Total pressure.

Returns

Saturation specific humidity.

                                                                                  {
        // If pressure is less than SVP, set qs to maximum of amrex::Real(1)
        if ( (p - es) <= amrex::Real(0) )
            return amrex::Real(1);
        else
            return epsilo*es / (p - omeps*es);
    }

Referenced by wv_sat_qsat_ice(), wv_sat_qsat_trans(), and wv_sat_qsat_water().

Here is the caller graph for this function:

wv_sat_svp_trans()

AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real SatMethods::wv_sat_svp_trans ( const amrex::Real &  t, const int  idx = 1  )

inlinestatic

Calculate transition saturation vapor pressure.

Parameters

[in]	t	Temperature.
[in]	idx	Scheme index.

Returns

Transition saturation vapor pressure.

                                                                              {
        amrex::Real esice;      // Saturation vapor pressure over ice
        amrex::Real weight;     // Intermediate scratch variable for es transition
        amrex::Real es;
 
        // Water
        if (t >= (tmelt - ttrice))
            es = wv_sat_svp_water(t,idx);
        else
            es = amrex::Real(0);
 
        // Ice
        if (t < tmelt) {
            esice = wv_sat_svp_ice(t,idx);
            if ( (tmelt - t) > ttrice )
                weight = amrex::Real(1);
            else
                weight = (tmelt - t)/ttrice;
 
            es = weight*esice + (amrex::Real(1) - weight)*es;
        }
        return es;
    }

Referenced by wv_sat_qsat_trans().

Here is the call graph for this function:

Here is the caller graph for this function:

wv_sat_svp_water()

AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real SatMethods::wv_sat_svp_water ( const amrex::Real &  t, const int  idx = 1  )

inlinestatic

Dispatcher for saturation vapor pressure over water.

Parameters

[in]	t	Temperature.
[in]	idx	Scheme index.

Returns

Saturation vapor pressure over water.

                                                                              {
        amrex::Real es;
        switch(idx)
            {
            case GoffGratch:
                es = GoffGratch_svp_water(t); break;
            case MurphyKoop:
                es = MurphyKoop_svp_water(t); break;
            case OldGoffGratch:
                es = OldGoffGratch_svp_water(t); break;
            case Bolton:
                es = Bolton_svp_water(t); break;
            }
        return es;
    }

Referenced by wv_sat_qsat_water(), and wv_sat_svp_trans().

Here is the call graph for this function:

Here is the caller graph for this function:

---

The documentation for this class was generated from the following file:

• Source/Utils/ERF_SatMethods.H