[gradsusr] ECMWF Data Issues

[AMS EmsiWx Bytča]

Hi Stuart,

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?

Regards,

Marian


Dňa 25. 3. 2016 o 12:11 stuart markham napísal(a):
>
> 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.
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

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