ECMWF ensemble forecasts problems: perhaps tdef?
José Paulo Bonatti
bonatti at CPTEC.INPE.BR
Thu Nov 24 11:42:30 EST 2005
Hi,
May be this work
gribmap -0 -e -v -i myfile.ctl
-0 : ignore the forecast time int the mapping; only use the base time
-e : ignore junk bytes at end of ECMWF files because of blocking
Bonatti
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
>
>***************************************************************************
>
>
>
--
*****************************************************************
Jose Paulo Bonatti
Dr. em Meteorologia - Pesquisador Titular
Instituto Nacional de Pesquisas Espaciais - INPE
Centro de Previsao de Tempo e Estudos Climaticos - CPTEC
Divisao de Modelagem e Desenvolvimento - DMD
(http://www.cptec.inpe.br/)
Rodovia Presidente Dutra, km 40 SPRJ
12630-000 Cachoeira Paulista - SP - Brasil
Tel: (+55) (12) 3186.8421 Fax: (+55) (12) 3101.2835
*****************************************************************
More information about the gradsusr
mailing list