<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Stuart,<br>
<br>
This one does not seem correct to me:<br>
<br>
'define PRESmsl = PRESmsl*0.01'<br>
<br>
replace with e.g.<br>
<br>
'define PRES = PRESmsl*0.01'<br>
<br>
I am not sure if that helps you to get all time steps, but for sure
on t=2 you would not get correct pressure contours.<br>
<br>
Let me know how it goes, and please show me how you get the ecmwf
grib.<br>
<br>
M.<br>
<br>
<div class="moz-cite-prefix">Dňa 25. 3. 2016 o 13:33 stuart markham
napísal(a):<br>
</div>
<blockquote
cite="mid:AM2PR09MB0339700BDCDD434793D24938F5830@AM2PR09MB0339.eurprd09.prod.outlook.com"
type="cite">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
<div id="divtagdefaultwrapper"
style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Hi Marian, </p>
<p><br>
</p>
<p>I've attached the full script for you, many thanks for your
help. </p>
<p><br>
</p>
<p>regards, </p>
<p><br>
</p>
<p>Stuart.</p>
<br>
<br>
<div style="color: rgb(0, 0, 0);">
<hr style="display:inline-block; width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt"
color="#000000" face="Calibri, sans-serif"><b>From:</b>
<a class="moz-txt-link-abbreviated" href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>
<a class="moz-txt-link-rfc2396E" href="mailto:gradsusr-bounces@gradsusr.org"><gradsusr-bounces@gradsusr.org></a> on behalf of AMS
EmsiWx Bytča <a class="moz-txt-link-rfc2396E" href="mailto:emsiwx@pocasie-bytca.sk"><emsiwx@pocasie-bytca.sk></a><br>
<b>Sent:</b> 25 March 2016 11:56<br>
<b>To:</b> GrADS Users Forum<br>
<b>Subject:</b> Re: [gradsusr] ECMWF Data Issues</font>
<div> </div>
</div>
<div>Hi Stuart,<br>
<br>
in your script you are opening a GFS control file... Replace
with (I suppose) ecmwf_grib.ctl . And I can't see where you
change the t (time), Can you please show me the rest of your
script?<br>
<br>
Regards,<br>
<br>
Marian<br>
<br>
<br>
<div class="moz-cite-prefix">Dňa 25. 3. 2016 o 12:11 stuart
markham napísal(a):<br>
</div>
<blockquote type="cite">
<div id="divtagdefaultwrapper" style="font-size:12pt;
color:#000000; background-color:#FFFFFF;
font-family:Calibri,Arial,Helvetica,sans-serif">
<p style="">Hi all,</p>
<p style=""><br>
</p>
<p style="">I'm currently developing charts using ecmwf
free data but I'm having an issue with the time
steps. The first image is fine but the rest of the
time steps use the same data as the first one. The
charts are not running in sequence. Can anybody advise
where my issues are. I've attached parts of my control
file and script to assist. </p>
<p style=""><br>
</p>
<p style=""><br>
</p>
<div>
<div>dset ^ecmwf_grib.grib2</div>
<div>index ^ecmwf_grib.grib2.idx</div>
<div>undef 9.999E+20</div>
<div>title ecmwf_grib.grib2</div>
<div>* produced by g2ctl v0.0.8.8</div>
<div>* command line options: -verf ecmwf_grib.grib2</div>
<div>* griddef=1:0:(720 x
361):grid_template=0:winds(N/S): lat-lon grid:(720 x
361) units 1e-06 input WE:NS output WE:SN res 48 lat
90.000000 to -90.000000 by 0.500000 lon 0.000000 to
359.500000 by 0.500000 #points=259920:winds(N/S)</div>
<div><br>
</div>
<div>dtype grib2</div>
<div>ydef 361 linear -90.000000 0.5</div>
<div>xdef 720 linear 0.000000 0.500000</div>
<div>tdef 11 linear 00Z25mar2016 1dy</div>
<div>zdef 1 linear 1 1</div>
<div>vars 2</div>
<div>HGT500mb 0,100,50000 0,3,5 ** 500 mb
Geopotential Height [gpm]</div>
<div>PRESmsl 0,101 0,3,0 ** mean sea level Pressure
[Pa]</div>
<div>ENDVARS</div>
<div><br>
</div>
</div>
<div><br>
</div>
----------------------------------------------------------------------------------------------------------------------------------------------------------------
<p style=""><br>
</p>
<p style=""><br>
</p>
<p style=""><br>
</p>
<div>'reinit'</div>
<div>'open gfs_grib.ctl'</div>
<div><br>
</div>
<div>* Get info from the descriptor file</div>
<div>'q ctlinfo'</div>
<div>_ctl = result</div>
<div>_undef = getctl(undef)</div>
<div>_tdef = getctl(tdef)</div>
<div>_zdef = getctl(zdef)</div>
<div><br>
</div>
<div>maps = 25</div>
<div><br>
</div>
<div>* Get the Time axis info</div>
<div>tsize = subwrd(_tdef,2) </div>
<div>_t1 = 1 ;* 2nd half of timeseries</div>
<div>_t2 = 45</div>
<div>tsize = _t2 - _t1 + 1</div>
<div>'set t '_t1' '_t2</div>
<div>'q dims'</div>
<div>times = sublin(result,5)</div>
<div>_time1 = subwrd(times,6) </div>
<div>_time2 = subwrd(times,8)</div>
<div>_tdim = _time1' '_time2</div>
<div><br>
</div>
<div>tincr = subwrd(_tdef,5)</div>
<div>_tdef = 'tdef 'tsize' linear '_time1' 'tincr</div>
<div><br>
</div>
<div>huh = subwrd(_tdef,4)</div>
<div><br>
</div>
<div>Best regards,</div>
<div><br>
</div>
<div>Stuart.</div>
<br>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre>_______________________________________________
gradsusr mailing list
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>
<a moz-do-not-send="true" class="moz-txt-link-freetext" href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a>
</pre>
</blockquote>
<br>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
gradsusr mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>
<a class="moz-txt-link-freetext" href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a>
</pre>
</blockquote>
<br>
</body>
</html>