<div dir="ltr">For the heights to change, you need to change the timestep. The example you sent has t=5. If, for example, you want it to vary from timestep 5 to 10, you need t=5 on the first, t=6 on the second, and so forth. Here's a simple loop, using your sample script:<div><br></div><div><span style="font-size:12.8000001907349px">'set lev 500'</span><br style="font-size:12.8000001907349px"><div>t1=5;t2=10</div><div>ts=t1</div><div>while(ts<=t2);#start time loop</div><span style="font-size:12.8000001907349px">'set gxout shaded'</span><br style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">'set csmooth on'</span><br style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">'set clevs 4920 4960 5000 5040 5080 5120 5160 5200 5240 5280 5320 5360 5400 5440 5480 5520 5560 5600 5640 5680 5720 5760 5800 5840 5880 5920 5960 6000 6040 6080 6120'</span><br style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">'set ccols 61 63 65 67 69 71 73 75 77 79 81 83 85 87 89 18 20 22 24 26 28 30 32 34 36 38 41 43 45 47 49 51'</span><br style="font-size:12.8000001907349px"><br style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">'d hgtprs(t='ts')'</span><br style="font-size:12.8000001907349px"><br style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">'define slp = const((prmslmsl*0.01),0,-u)'</span><br style="font-size:12.8000001907349px"><br style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">'set gxout contour'</span><br style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">'set ccolor 0'</span><br style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">'set cthick 10'</span><br style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">'set cint 5'</span><br style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">'set clopts 0'</span><br style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">'set clab masked'</span><br style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">'d slp'</span><br style="font-size:12.8000001907349px">ts=ts+1</div><div>endwhile;#end time loop<br><span style="font-size:12.8000001907349px">'cbarn'</span><br><div><br></div></div><div>Stephen Mc</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Apr 27, 2015 at 8:08 AM, Chorley Weather <span dir="ltr"><<a href="mailto:weatherstu@chorleyweather.com" target="_blank">weatherstu@chorleyweather.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="HOEnZb"><div class="h5">On <a href="tel:26-04-2015%2020" value="+12604201520" target="_blank">26-04-2015 20</a>:28, Stephen McMillan wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
Stuart,<br>
It would help if you shared the portion of your script that you used<br>
to attempt doing the 500mb height charts. What is the snag? What<br>
GrADS error are you getting? What are you trying to get your code to<br>
do?<br>
<br>
Whether you need to define them depends on how you're using them. <br>
Following is a simple sample, assuming I've opened the gfs 0p25:<br>
<br>
'set lev 500'<br>
'set clevs 5160 5280 5400 5520 5640'<br>
'set ccols 9 4 11 3 7'<br>
'd hgtprs(t=5)'<br>
<br>
You can also define it first...<br>
'define hgt500=hgtprs(t=5)'<br>
'd hgt500'<br>
<br>
Stephen Mc<br>
<br>
On Sun, Apr 26, 2015 at 9:53 AM, Chorley Weather<br>
<<a href="mailto:weatherstu@chorleyweather.com" target="_blank">weatherstu@chorleyweather.com</a>> wrote:<br>
<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
Hi guys,<br>
<br>
Still rather new to grads but have managed to create some good<br>
charts running the GFS 0.25, however, I have hit a snag on the 500mb<br>
height chart. Ive scripted all the colour and shading, including the<br>
contour bar but the hgt500 code will not work. Does this need<br>
defining and if so can someone please share the definition. Also I'd<br>
love a jet stream chart, however, I am also struggling with that.<br>
Any info would be much appreciated. I've attached my precip &<br>
rainfall chart to show I'm not completely without hope :)<br>
<br>
Thanks in advance,<br>
Stuart.<br>
<br>
--<br>
Chief Forecaster & Owner<br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
</div></div><a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a> [1]<br>
</blockquote>
<br>
<br>
<br>
Links:<br>
------<br>
[1] <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><span class=""><br>
<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>
</span></blockquote>
Hi Stephen,<br>
<br>
Thank you very much for your help. The sample worked fine and I have manged to create a chart including a pressure overlay. However, the issue I now have is the heights don't change when I run the gfs data. I still get the same heights but the pressure works fine. I've attached some of my script for you to view. Not sure what I'm missing to make the heights run with the gfs data. All my other charts work fine so I know the noaa/dods link is fine.<br>
<br>
'set lev 500'<br>
'set gxout shaded'<br>
'set csmooth on'<br>
'set clevs 4920 4960 5000 5040 5080 5120 5160 5200 5240 5280 5320 5360 5400 5440 5480 5520 5560 5600 5640 5680 5720 5760 5800 5840 5880 5920 5960 6000 6040 6080 6120'<br>
'set ccols 61 63 65 67 69 71 73 75 77 79 81 83 85 87 89 18 20 22 24 26 28 30 32 34 36 38 41 43 45 47 49 51'<br>
<br>
'd hgtprs(t=5)'<br>
<br>
'define slp = const((prmslmsl*0.01),0,-u)'<br>
<br>
'set gxout contour'<br>
'set ccolor 0'<br>
'set cthick 10'<br>
'set cint 5'<br>
'set clopts 0'<br>
'set clab masked'<br>
'd slp'<br>
<br>
'cbarn'<br>
<br>
Regards,<br>
Stuart<div class="HOEnZb"><div class="h5"><br>
-- <br>
Chief Forecaster & Owner</div></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></div>