equivalentPotentialTemperature {aiRthermo} | R Documentation |
This function calculates the equivalent potential temperature (K), following the techniques used in Davies-Jones (2009).
equivalentPotentialTemperature(P, Temp, w, TLCL, consts = export_constants())
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. |
This function returns the value of the equivalent potential temp (K).
Davies-Jones, R. (2009). On formulas for equivalent potential temperature. Monthly Weather Review, 137(9), 3137-3148.
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)