<div dir="ltr"><div>Daniel,<br>I used arbitrary names for the height field. Use 'q file' to check to see what the actual name is for the height field in your data (if it even exists). Also, check for syntax errors. It looks like I missed a ')' at the end of the command.<br><br></div>Jeff<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Nov 14, 2014 at 2:05 AM, Daniel Nairn <span dir="ltr"><<a href="mailto:dannairn@gmail.com" target="_blank">dannairn@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Thanks Jeff, I wasn't convinced that my original formula was going to be that accurate. I've updated my script, as you suggested, for the 700-500 lapse rates however the images generated are now blank? The script completes without any errors. If I revert back to my original script then the images are produced as expected. <br><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On 13 November 2014 22:59, Jeff Duda <span dir="ltr"><<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>It would be more accurate if you used the actual height fields at the pressure levels instead of assuming a fixed height difference. In other words, try<br><br></div>'define LR75 = (tmpprs(lev=500)-tmpprs(lev=700)) / (hgtprs(lev=500)-hgtprs(lev=700)'<br><br></div>for 700-500 mb lapse rates. Note how I flipped the order of the subtraction from your formula. You'll want to use delta_T = T_upper - T_lower when computing lapse rates.<br><br></div>0-3 km lapse rates can be difficult if you only have isobaric fields. There are scripts in the script library to interpolate isobaric fields to constant height fields. If you don't need high accuracy, subbing in a pressure level that is close to 3 km works. Depending on the time of year, usually 700 mb is pretty close, but keep in mind that people usually use 0-3 km ABOVE GROUND LEVEL lapse rates, so in areas of high elevation, 700 mb will be far less than 3 km above ground.<br><br></div>Jeff Duda<br></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Thu, Nov 13, 2014 at 3:53 PM, Daniel Nairn <span dir="ltr"><<a href="mailto:dannairn@gmail.com" target="_blank">dannairn@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Hi <br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Could somebody please help verify my two GrAds scripts for displaying Lapse Rates? I have created two individual scripts, one for displaying 700-500mb lapse rates and the other for displaying 0-3km lapse rates, they are as follows:-<br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">0-3 Script<br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">'define temp = ((tmpprs(lev=1000)-tmpprs(lev=700))/2.5908)'<br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">700-500 Script<br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">'define temp = ((tmpprs(lev=700)-tmpprs(lev=500))/2.5908)'<br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">They work and seem to verify against other sources but they have been pieced together form many Google searches and I'm a little unsure if they are accurate/correct :-)<br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Thank you<br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div></div>
<br></div></div>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><span><font color="#888888"><br><br clear="all"><span class="HOEnZb"><font color="#888888"><br>-- <br><div>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br></div>
</font></span></font></span></div><span class="HOEnZb"><font color="#888888">
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></font></span></blockquote></div><br></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature">Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br></div>
</div>