equivalentPotentialTemperature: Equivalent Potential Temperature

In aiRthermo: Atmospheric Thermodynamics and Visualization

Description

This function calculates the equivalent potential temperature (K), following the techniques used in Davies-Jones (2009).

Usage

equivalentPotentialTemperature(P, Temp, w, TLCL, consts = export_constants())

Arguments

P	The pressure (Pa) of the air parcel.
Temp	The temperature (K) of the parcel.
w	The mixing ratio (kg/kg) of the parcel.
TLCL	The temperature (K) at the Lifting Condensation Level (LCL).
consts	The constants defined in aiRthermoConstants data are necessary.

Value

This function returns the value of the equivalent potential temp (K).

References

Davies-Jones, R. (2009). On formulas for equivalent potential temperature. Monthly Weather Review, 137(9), 3137-3148.

See Also

PT2Theta and moistCp

Examples

data(RadiosondeA)
aPs<-RadiosondeA[,1]*100
aP0<-aPs[1]
aT0<-C2K(RadiosondeA[1,3])
aw0<-RadiosondeA[1,6]/1000
deltaP=1
Na=length(aPs)
Ptop=aPs[Na]
fndlcl=find_lcl(Ptop,aP0,aT0,aw0,deltaP)
TLCL=fndlcl$Tlcl
equivalentPotentialTemperature(aP0,aT0,aw0,TLCL)

Example output

[1] 312.2493

aiRthermo documentation built on May 1, 2019, 9:24 p.m.