Calculating LCL on plan-view plot

Shawn Delaney sd19surf at GMAIL.COM
Thu May 1 08:06:27 EDT 2008


       Do you run your own WRF model or are you obtaining the GRIB files
from else where.   The reason that I ask is there is a variable called
"hgtadcl" that can be enabled and it would solve your LCL caluculation
problems quickly.  I am currently working on something similar and here's
what I am using.

'set gxout stat'
'set stat on'
'd tmp2m-273.16'
'd dpt2m-273.16'
'd pressfc/100'
lclp = Preslcl(temp,dewp,pres)
lclp = math_format(fmt,lclp)

Where Preslcl is the function from Bob Hart's

function Preslcl(temp,dewp,pres)

* Calculate press of LCL given temp & dewp in C and pressure



I am using the pressfc/100 and this gets me within ~10mb of the BUFKIT LCL

This is rough code and hasn't been rehacked. I apologize for the slop in


On Wed, Apr 30, 2008 at 10:33 PM, Dan Leins <theedge981 at> wrote:

> All,
> I am trying to calculate and plot LCL heights (in millibars) across an
> entire WRF model domain. I am already accomplishing this on a point-by-point
> basis using Bob Hart's program, but would like to show the
> heights on a plan-view plot.  Has anyone accomplished this?  I dug pretty
> deep into the code to determine how Bob calculated LCL heights
> and I was able to adapt a portion of his script for my needs. Here is what
> I've got so far...
> 'define tempk=tmpsfc+273.15'
> 'define dewpk=dptprs(lev=1013)+273.15'
> 'define Parta=1/(dewpk-56)'
> 'define Partb=log(tempk/dewpk)/800'
> 'define Tlcl=1/(Parta+Partb)+56'
> 'define Tlclk=Tlcl'
> 'define tempk=tmpprs+273.15'
> 'define theta=tempk*pow(1000/990,0.286)'
> 'define plcl=1000*pow(Tlclk/theta,3.48)'
> 'd plcl'
> This seems to match up "poorly" with the LCL heights indicated on my
> skew-T images, and is just plain horrendous when compared to the LCL heights
> shown in Bufkit. Not that Bufkit is considered the correct solution, but I
> expect to be within the ballpark with my calculations. In some spots I was
> off by over 200mb.
> Any thoughts on how I might be able to accomplish this?
> Thanks,
> Dan
-------------- next part --------------
An HTML attachment was scrubbed...

More information about the gradsusr mailing list