All,<br><br>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 <a href="http://plotskew.gs">plotskew.gs</a> program, but would like to show the heights on a plan-view plot. Has anyone accomplished this? I dug pretty deep into the <a href="http://plotskew.gs">plotskew.gs</a> 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...<br>
<br>'define tempk=tmpsfc+273.15'<br>'define dewpk=dptprs(lev=1013)+273.15'<br>'define Parta=1/(dewpk-56)'<br>'define Partb=log(tempk/dewpk)/800'<br>'define Tlcl=1/(Parta+Partb)+56'<br>
<br><br>'define Tlclk=Tlcl'<br>'define tempk=tmpprs+273.15'<br>'define theta=tempk*pow(1000/990,0.286)'<br>'define plcl=1000*pow(Tlclk/theta,3.48)'<br>'d plcl'<br><br>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.<br>
<br>Any thoughts on how I might be able to accomplish this?<br><br>Thanks,<br>Dan<br>