[gradsusr] ECMWF Data Issues

[AMS EmsiWx Bytča]

Stuart,

This one does not seem correct to me:

'define PRESmsl = PRESmsl*0.01'

replace with e.g.

'define PRES = PRESmsl*0.01'

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.

Let me know how it goes, and please show me how you get the ecmwf grib.

M.

Dňa 25. 3. 2016 o 13:33 stuart markham napísal(a):
>
> Hi Marian,
>
>
> I've attached the full script for you, many thanks for your help.
>
>
> regards,
>
>
> Stuart.
>
>
>
> ------------------------------------------------------------------------
> *From:* gradsusr-bounces at gradsusr.org <gradsusr-bounces at gradsusr.org> 
> on behalf of AMS EmsiWx Bytča <emsiwx at pocasie-bytca.sk>
> *Sent:* 25 March 2016 11:56
> *To:* GrADS Users Forum
> *Subject:* Re: [gradsusr] ECMWF Data Issues
> 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
>
>
>
> _______________________________________________
> 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/70b396ef/attachment.html