David Small wrote:
> Hi,
>
> Has anyone developed a function that will calculate equivalent potential
> temperature that doesn't assume saturation? Has anyone converted the
> wrf function wrf_eth.ncl to work for input that was not produced by the
> wrf model?
>
> Thanks,
>
> Dave Small
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
The following is fortran code provided by PROFS.
I am sure it could be implemented in NCL.
real function ept(t,td,p)
c include 'lib_dev:[gudoc]edfvaxbox.for/list'
c baker, schlatter 17-may-1982 original version.
c this function returns the equivalent potential temperature ept
c (celsius) for a parcel of air initially at temperature t (celsius),
c dew point td (celsius) and pressure p (millibars). the formula used
c is eq.(43) in bolton, david, 1980: "the computation of equivalent
c potential temperature," monthly weather review, vol. 108, no. 7
c (july), pp. 1046-1053. the maximum error in ept in 0.3c. in most
c cases the error is less than 0.1c.
c
c compute the mixing ratio (grams of water vapor per kilogram of
c dry air).
w = wmr(p,td)
c compute the temperature (celsius) at the lifting condensation level.
tlcl = tcon(t,td)
tk = t+273.15
tl = tlcl+273.15
pt = tk*(1000./p)**(0.2854*(1.-0.00028*w))
eptk = pt*exp((3.376/tl-0.00254)*w*(1.+0.00081*w))
ept= eptk-273.15
return
end
-- ====================================================== Dennis J. Shea tel: 303-497-1361 | P.O. Box 3000 fax: 303-497-1333 | Climate Analysis Section | Climate & Global Dynamics Div. | National Center for Atmospheric Research | Boulder, CO 80307 | USA email: shea 'at' ucar.edu | ====================================================== _______________________________________________ ncl-talk mailing list List instructions, subscriber options, unsubscribe: http://mailman.ucar.edu/mailman/listinfo/ncl-talkReceived on Mon May 18 2009 - 12:55:27 MDT
This archive was generated by hypermail 2.2.0 : Tue May 19 2009 - 09:17:29 MDT