[gradsusr] ECMWF Data Issues

[stuart markham]

Hi all,


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.



dset ^ecmwf_grib.grib2
index ^ecmwf_grib.grib2.idx
undef 9.999E+20
title ecmwf_grib.grib2
* produced by g2ctl v0.0.8.8
* command line options: -verf ecmwf_grib.grib2
* 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)

dtype grib2
ydef 361 linear -90.000000 0.5
xdef 720 linear 0.000000 0.500000
tdef 11 linear 00Z25mar2016 1dy
zdef 1 linear 1 1
vars 2
HGT500mb   0,100,50000   0,3,5 ** 500 mb Geopotential Height [gpm]
PRESmsl  0,101   0,3,0 ** mean sea level Pressure [Pa]
ENDVARS


----------------------------------------------------------------------------------------------------------------------------------------------------------------




'reinit'
'open gfs_grib.ctl'

* Get info from the descriptor file
'q ctlinfo'
_ctl = result
_undef = getctl(undef)
_tdef = getctl(tdef)
_zdef = getctl(zdef)

maps = 25

* Get the Time axis info
tsize = subwrd(_tdef,2)
_t1 = 1       ;* 2nd half of timeseries
_t2 = 45
tsize = _t2 - _t1 + 1
'set t '_t1' '_t2
'q dims'
times  = sublin(result,5)
_time1 = subwrd(times,6)
_time2 = subwrd(times,8)
_tdim = _time1' '_time2

tincr = subwrd(_tdef,5)
_tdef = 'tdef 'tsize' linear '_time1' 'tincr

huh = subwrd(_tdef,4)

Best regards,

Stuart.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160325/922eaac2/attachment.html