[gradsusr] Trouble Calculating ThetaE

[Corby Dickerson - NOAA Federal]

Andrew,

Thanks for your reply! I knew it had to be forgetting something simple.
I'll give it a try and see if that solves my problem. I'll letcha know!

Again, thank you for the help!

corBy

*C. Corby Dickerson IV*
-General Forecaster-
Weather Story Program Manager
*406-329-4840*
[image: Inline image 1] <https://twitter.com/NWSMissoula>   [image: Inline
image 2]<https://www.facebook.com/#!/US.NationalWeatherService.Missoula.gov>
 [image: Inline image 3] <https://www.youtube.com/user/NWSMissoula>   [image:
Inline image 4] <http://www.wrh.noaa.gov/mso/>
<http://www.nssl.noaa.gov/projects/ping/>


On Wed, May 7, 2014 at 6:28 PM, Andrew Friedman <andfried at berkeley.edu>wrote:

> Hi Corby,
> Where in your script do you set the vertical levels? It should be done
> prior to the calculation of level-varying variables.
> Andrew
>
> On May 7, 2014, at 3:57 PM, Corby Dickerson - NOAA Federal <
> corby.dickerson at noaa.gov> wrote:
>
> > Greetings!
> >
> > I have written the following script to calculate ThetaE from NARR data
> using a new formula referenced in Volume 137 of AMS's Monthly Weather
> Review. The script runs successfully and displays data which makes sense.
> However it plots the same values of ThetaE no matter what level I choose. I
> am hoping someone out there will be able to provide me with some direction
> as to what I might be doing wrong.
> >
> > Thanks in advance for the help!
> >
> >
> ******************************************************************************************************************************************
> > *Script to calculate Equivalent Potential Temperature referencing
> equation 6.5 from
> > *Davies-Jones' September 2009 Monthly Weather Review findings On
> Formulas for Equivalent
> > 'open narr.ctl'
> > 'set lat 30 55'
> > 'set lon -140 -100'
> > *First calculate 2m Dew Point Temperature
> > *
> > 'define tc2m = TMP2m - 273.15'
> > 'define B = (log(RH2m/100)+((17.27*tc2m)/(237.3+tc2m)))/17.27'
> > 'define TdC2m = (237.3*B)/(1-B)'
> > 'define dewp2m = TdC2m + 273.15'
> > *
> > *Then we can calculate Temperature at LCL [K]
> > 'define templcl=(((1/(1/(dewp2m-56)+log((tmp2m/dewp2m))/800))+56))'
> > *
> > *And calculate the Pressure at LCL
> > 'define preslcl=(pressfc*pow(((templcl+273.16)/tmp2m),(7/2)))/1000'
> > *
> > *With this data in hand we can calculate Potential Temperature at LCL
> *(K)
> > *Saturation Vapor Pressure (mb)=6.112 * exp[17.67(Tk-C)/(Tk-C+243.5)]
> > *
> > 'define eslcl =
> 6.112*exp((17.67*(templcl-273.15))/(templcl-273.15+243.5))'
> > *
> > 'define thetaLCL = templcl*pow(1000/(lev-eslcl),0.2854)'
> > *
> > *Calculate Equivalent Potential Temperature (ThetaE)
> > *Define constants
> > *
> > 'define L0 = 2501000'
> > 'define L1 = 1754'
> > 'define L3 = L0 - L1'
> > 'define K2 = 1137000'
> > 'define cpd = 1105.7'
> > 'define C = 273.15'
> > 'define mythetae =
> thetaLCL*exp((L3*(tmpprs-C)+(K2*ICMRprs))*ICMRprs/(cpd*templcl))'
> >
> > C. Corby Dickerson IV
> > -General Forecaster-
> > Weather Story Program Manager
> > 406-329-4840
> > <image.png>   <image.png>   <image.png>   <image.png>
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140511/d88844e4/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 3091 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20140511/d88844e4/attachment.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: image/png
Size: 947 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20140511/d88844e4/attachment-0001.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 2165 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20140511/d88844e4/attachment-0002.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 1147 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20140511/d88844e4/attachment-0003.png