ECMWF ensemble forecasts problems: perhaps tdef?

[Becker, Bernd]

Martynas,

use wgrib or gribex to slice  the input file up into one file per
ensemble member, write one .ctl file per member and take things from
there.
Do you have the Provost set of Cdroms? Good. There are examples on it
how to deal with ensemble data before Grib2 is fully up and compatible
with grads.

Cheers,
Bernd.


On Thu, 2005-11-24 at 14:07, Martynas wrote:
> Hi,
>
> I've got some grib files containing ECMWF monthly forecasts. Each file
> contains 9 members and 6 forecast means (e.g. for forecast started on 1 Feb,
> the forecasts means are for Feb, Mar, Apr, May, Jun, Jul).
> Grib2ctl script seems to generate a .ctl file which doesn't work properly.
>
> My original ctl file looks like:
>
> dtype grib 255
> options yrev
> ydef 8 linear 49.500000 1.5
> xdef 12 linear 16.500000 1.500000
> eas_mm_GEO_500hPa_19900201_cnrm
> index ^seas_mm_GEO_500hPa_19900201_cnrm.idx
> undef 9.999E+20
> title seas_mm_GEO_500hPa_19900201_cnrm
> *  produced by grib2ctl v0.9.12.5p33k
> dtype grib 255
> options yrev
> ydef 8 linear 49.500000 1.5
> xdef 12 linear 16.500000 1.500000
> tdef 1 linear 00Z01jan1990 1mo
> *********(i got some results with tdef 6 linear 00Z01jan1990 720hr and etc.)
> zdef 1 linear 1 1
> vars 1
> MSLMAprs  0 129,100,500  ** Mean sea level pressure (MAPS) [Pa]
> ENDVARS
> def 1 linear 00Z01jan1990 1mo
> zdef 1 linear 1 1
> vars 1
> MSLMAprs  0 129,100,500  ** Mean sea level pressure (MAPS) [Pa]
> ENDVARS
>
>
> To my mind, the problem is the tdef line: There are 9 ensemble variables in
> the file, containing 6 forecast means;
>
> I'm not able to plot them unless i change 1mo to 720hr-but gribmap doesn't
> make all matches for the .idx file (several months are of 744hr and i'm not
> able to plot them)
>
> If i make tdef 6 linear 00Z01jan1990 1mo, i get "time out of bounds" and no
> solution;
>
> I've tried to slice the original grib file "to pieces containing" only
> several records from the original file and the result was the same problem -
> even if 2 forecast means are there in the grib file, "1mo" gives "time out
> of bounds" and the time steb "720hr" works and i was able to plot smth.
>
> But not with the original file. Perhaps the timestep is inncorect? Number of
> timesteps?
>
> I would be glad to get any sugestions on that (hope that the problem is
> described clear enough).
>
> Martynas
>
>
>
> Here's my wgrib output:
>
> *****************Wgrib output***********************
>
> Undefined parameter table (center 98-231 table 190), using NCEP-opn
> 1:0:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=0:P2=24:TimeU=1:500
> mb:24hr fcst:ECMWF stream?:NAve=0
> 2:360:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=0:P2=24:TimeU=1:500
> mb:24hr fcst:ECMWF stream?:NAve=0
> 3:720:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=0:P2=24:TimeU=1:500
> mb:24hr fcst:ECMWF stream?:NAve=0
> 4:1080:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=0:P2=24:TimeU=1:500
> mb:24hr fcst:ECMWF stream?:NAve=0
> 5:1440:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=0:P2=24:TimeU=1:500
> mb:24hr fcst:ECMWF stream?:NAve=0
> 6:1800:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=0:P2=24:TimeU=1:500
> mb:24hr fcst:ECMWF stream?:NAve=0
> 7:2160:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=0:P2=24:TimeU=1:500
> mb:24hr fcst:ECMWF stream?:NAve=0
> 8:2520:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=0:P2=24:TimeU=1:500
> mb:24hr fcst:ECMWF stream?:NAve=0
> 9:2880:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=0:P2=24:TimeU=1:500
> mb:24hr fcst:ECMWF stream?:NAve=0
> 10:3240:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=2:P2=184:TimeU=1:500
> mb:696hr fcst:ECMWF stream?:NAve=0
> 11:3600:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=2:P2=184:TimeU=1:500
> mb:696hr fcst:ECMWF stream?:NAve=0
> 12:3960:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=2:P2=184:TimeU=1:500
> mb:696hr fcst:ECMWF stream?:NAve=0
> 13:4320:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=2:P2=184:TimeU=1:500
> mb:696hr fcst:ECMWF stream?:NAve=0
> 14:4680:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=2:P2=184:TimeU=1:500
> mb:696hr fcst:ECMWF stream?:NAve=0
> 15:5040:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=2:P2=184:TimeU=1:500
> mb:696hr fcst:ECMWF stream?:NAve=0
> 16:5400:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=2:P2=184:TimeU=1:500
> mb:696hr fcst:ECMWF stream?:NAve=0
> 17:5760:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=2:P2=184:TimeU=1:500
> mb:696hr fcst:ECMWF stream?:NAve=0
> 18:6120:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=2:P2=184:TimeU=1:500
> mb:696hr fcst:ECMWF stream?:NAve=0
> 19:6480:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=5:P2=160:TimeU=1:500
> mb:1440hr fcst:ECMWF stream?:NAve=0
> 20:6840:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=5:P2=160:TimeU=1:500
> mb:1440hr fcst:ECMWF stream?:NAve=0
> 21:7200:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=5:P2=160:TimeU=1:500
> mb:1440hr fcst:ECMWF stream?:NAve=0
> 22:7560:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=5:P2=160:TimeU=1:500
> mb:1440hr fcst:ECMWF stream?:NAve=0
> 23:7920:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=5:P2=160:TimeU=1:500
> mb:1440hr fcst:ECMWF stream?:NAve=0
> 24:8280:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=5:P2=160:TimeU=1:500
> mb:1440hr fcst:ECMWF stream?:NAve=0
> 25:8640:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=5:P2=160:TimeU=1:500
> mb:1440hr fcst:ECMWF stream?:NAve=0
> 26:9000:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=5:P2=160:TimeU=1:500
> mb:1440hr fcst:ECMWF stream?:NAve=0
> 27:9360:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=5:P2=160:TimeU=1:500
> mb:1440hr fcst:ECMWF stream?:NAve=0
> 28:9720:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=8:P2=112:TimeU=1:500
> mb:2160hr fcst:ECMWF stream?:NAve=0
> 29:10080:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=8:P2=112:TimeU=1:500
> mb:2160hr fcst:ECMWF stream?:NAve=0
> 30:10440:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=8:P2=112:TimeU=1:500
> mb:2160hr fcst:ECMWF stream?:NAve=0
> 31:10800:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=8:P2=112:TimeU=1:500
> mb:2160hr fcst:ECMWF stream?:NAve=0
> 32:11160:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=8:P2=112:TimeU=1:500
> mb:2160hr fcst:ECMWF stream?:NAve=0
> 33:11520:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=8:P2=112:TimeU=1:500
> mb:2160hr fcst:ECMWF stream?:NAve=0
> 34:11880:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=8:P2=112:TimeU=1:500
> mb:2160hr fcst:ECMWF stream?:NAve=0
> 35:12240:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=8:P2=112:TimeU=1:500
> mb:2160hr fcst:ECMWF stream?:NAve=0
> 36:12600:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=8:P2=112:TimeU=1:500
> mb:2160hr fcst:ECMWF stream?:NAve=0
> 37:12960:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=11:P2=88:TimeU=1:500
> mb:2904hr fcst:ECMWF stream?:NAve=0
> 38:13320:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=11:P2=88:TimeU=1:500
> mb:2904hr fcst:ECMWF stream?:NAve=0
> 39:13680:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=11:P2=88:TimeU=1:500
> mb:2904hr fcst:ECMWF stream?:NAve=0
> 40:14040:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=11:P2=88:TimeU=1:500
> mb:2904hr fcst:ECMWF stream?:NAve=0
> 41:14400:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=11:P2=88:TimeU=1:500
> mb:2904hr fcst:ECMWF stream?:NAve=0
> 42:14760:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=11:P2=88:TimeU=1:500
> mb:2904hr fcst:ECMWF stream?:NAve=0
> 43:15120:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=11:P2=88:TimeU=1:500
> mb:2904hr fcst:ECMWF stream?:NAve=0
> 44:15480:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=11:P2=88:TimeU=1:500
> mb:2904hr fcst:ECMWF stream?:NAve=0
> 45:15840:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=11:P2=88:TimeU=1:500
> mb:2904hr fcst:ECMWF stream?:NAve=0
> 46:16200:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=14:P2=40:TimeU=1:500
> mb:3624hr fcst:ECMWF stream?:NAve=0
> 47:16560:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=14:P2=40:TimeU=1:500
> mb:3624hr fcst:ECMWF stream?:NAve=0
> 48:16920:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=14:P2=40:TimeU=1:500
> mb:3624hr fcst:ECMWF stream?:NAve=0
> 49:17280:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=14:P2=40:TimeU=1:500
> mb:3624hr fcst:ECMWF stream?:NAve=0
> 50:17640:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=14:P2=40:TimeU=1:500
> mb:3624hr fcst:ECMWF stream?:NAve=0
> 51:18000:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=14:P2=40:TimeU=1:500
> mb:3624hr fcst:ECMWF stream?:NAve=0
> 52:18360:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=14:P2=40:TimeU=1:500
> mb:3624hr fcst:ECMWF stream?:NAve=0
> 53:18720:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=14:P2=40:TimeU=1:500
> mb:3624hr fcst:ECMWF stream?:NAve=0
> 54:19080:d=90020100:MSLMA:kpds5=129:kpds6=100:kpds7=500:TR=10:P1=14:P2=40:TimeU=1:500
> mb:3624hr fcst:ECMWF stream?:NAve=0
>
> ***************************************************************************
--
Bernd Becker   The Monthly Outlook
Met Office  FitzRoy Road  Exeter   Devon EX1 3PB  United Kingdom
Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050
E-mail:bernd.becker at metoffice.gov.uk -  http://www.metoffice.gov.uk