<div dir="ltr">Stuart,<div>Use whatever timeseries (or timestep range) fits your needs. What I sent you was just an example with six timesteps. The point was that you need to change your timestep through a time loop, either by specifying the timestep (or time/date) as an argument for the variable, or by setting the timestep (or time/date) for each increment of time.</div><div><br></div><div>Stephen Mc</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Apr 27, 2015 at 9:06 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:27-04-2015%2013" value="+12704201513">27-04-2015 13</a>:33, Stephen McMillan wrote:<br>
> For the heights to change, you need to change the timestep. The<br>
> example you sent has t=5. If, for example, you want it to vary from<br>
> timestep 5 to 10, you need t=5 on the first, t=6 on the second, and so<br>
> forth. Here's a simple loop, using your sample script:<br>
><br>
> 'set lev 500'<br>
><br>
> t1=5;t2=10<br>
> ts=t1<br>
> while(ts<=t2);#start time loop'set gxout shaded'<br>
> 'set csmooth on'<br>
> 'set clevs 4920 4960 5000 5040 5080 5120 5160 5200 5240 5280 5320 5360<br>
> 5400 5440 5480 5520 5560 5600 5640 5680 5720 5760 5800 5840 5880 5920<br>
> 5960 6000 6040 6080 6120'<br>
> 'set ccols 61 63 65 67 69 71 73 75 77 79 81 83 <br>
> 85 87 89 18 20 22 24 26 28 30 32 34 36 38 <br>
> 41 43 45 47 49 51'<br>
><br>
> 'd hgtprs(t='ts')'<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>
> ts=ts+1<br>
> endwhile;#end time loop<br>
> 'cbarn'<br>
><br>
> Stephen Mc<br>
><br>
> On Mon, Apr 27, 2015 at 8:08 AM, Chorley Weather<br>
> <<a href="mailto:weatherstu@chorleyweather.com">weatherstu@chorleyweather.com</a>> wrote:<br>
><br>
</div></div><div><div class="h5">>> On 26-04-2015 20 [1]:28, Stephen McMillan wrote:<br>
>><br>
>> Stuart,<br>
>> It would help if you shared the portion of your script that you<br>
>> 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<br>
>> to<br>
>> do?<br>
>><br>
>> Whether you need to define them depends on how you're using<br>
>> 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">weatherstu@chorleyweather.com</a>> wrote:<br>
>><br>
>> 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<br>
>> 500mb<br>
>> height chart. Ive scripted all the colour and shading, including<br>
>> 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<br>
>> 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">gradsusr@gradsusr.org</a><br>
</div></div>>> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a> [2] [1]<br>
>><br>
>> Links:<br>
>> ------<br>
>> [1] <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a> [2]<br>
<span class="">>><br>
>> _______________________________________________<br>
>> gradsusr mailing list<br>
>> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
</span>>> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a> [2]<br>
<div><div class="h5">> Hi Stephen,<br>
><br>
> Thank you very much for your help. The sample worked fine and I have<br>
> manged to create a chart including a pressure overlay. However, the<br>
> issue I now have is the heights don't change when I run the gfs data.<br>
> I still get the same heights but the pressure works fine. I've<br>
> attached some of my script for you to view. Not sure what I'm missing<br>
> to make the heights run with the gfs data. All my other charts work<br>
> 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<br>
> 5360 5400 5440 5480 5520 5560 5600 5640 5680 5720 5760 5800 5840 5880<br>
> 5920 5960 6000 6040 6080 6120'<br>
> 'set ccols 61 63 65 67 69 71 73 75 77 79 81 <br>
> 83 85 87 89 18 20 22 24 26 28 30 32 34 36 <br>
> 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<br>
><br>
> --<br>
> Chief Forecaster & Owner<br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
</div></div>> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a> [2]<br>
><br>
><br>
><br>
> Links:<br>
> ------<br>
> [1] tel:26-04-2015%2020<br>
> [2] <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<span class="">><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>
Hi Stephen,<br>
<br>
</span>Again many thanks for your assistance. The time series I have been using<br>
for all my charts is below. Would I need to change this?<br>
<br>
function main(args)<br>
<br>
* Parse the arguments: date, hour<br>
if (args = '')<br>
prompt 'Enter forecast date (example, 20110717) --> '<br>
pull date<br>
prompt 'Enter forecast hour (example, 00 or 06 or 12 or 18) --> '<br>
pull hour<br>
else<br>
date = subwrd(args,1)<br>
hour = subwrd(args,2)<br>
endif<br>
<br>
'reinit'<br>
'sdfopen<br>
<a href="http://nomads.ncep.noaa.gov:9090/dods/gfs_0p25/gfs'date'/gfs_0p25_'hour'z" target="_blank">http://nomads.ncep.noaa.gov:9090/dods/gfs_0p25/gfs'date'/gfs_0p25_'hour'z</a>'<br>
<br>
* Get info from the descriptor file<br>
'q ctlinfo'<br>
_ctl = result<br>
_undef = getctl(undef)<br>
_tdef = getctl(tdef)<br>
_zdef = getctl(zdef)<br>
<br>
maps = 81<br>
<br>
* Get the Time axis info<br>
tsize = subwrd(_tdef,2)<br>
_t1 = 1 ;* 2nd half of timeseries<br>
_t2 = 45<br>
tsize = _t2 - _t1 + 1<br>
'set t '_t1' '_t2<br>
'q dims'<br>
times = sublin(result,5)<br>
_time1 = subwrd(times,6)<br>
_time2 = subwrd(times,8)<br>
_tdim = _time1' '_time2<br>
<br>
tincr = subwrd(_tdef,5)<br>
_tdef = 'tdef 'tsize' linear '_time1' 'tincr<br>
<br>
huh = subwrd(_tdef,4)<br>
<br>
Regards,<br>
<div class="HOEnZb"><div class="h5">Stuart.<br>
<br>
<br>
--<br>
Chief Forecaster & Owner<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>
</div></div></blockquote></div><br></div>