<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=iso-8859-1">
<META content="MSHTML 6.00.2800.1607" name=GENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT face=Arial size=2>Yes Christian, I used gribmap.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>After your suggestion I tried:</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>gribmap -ie
EMEMO1500_20080313_011_032.ctl</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>gribmap: Scanning binary GRIB
file(s):<BR>gribmap: Opening GRIB file
EMEMO1500_20080313_011_032.grb<BR>gribmap: Reached EOF<BR>gribmap:
Writing the map...</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>But I had no success:</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>grads -cl "open
EMEMO1500_20080313_011_032.ctl"</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>ga-> d zprs<BR>Cannot contour grid - all
undefined values<BR></FONT></DIV>
<DIV><FONT face=Arial size=2>Thank you anyway.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Daniele</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<BLOCKQUOTE
style="PADDING-RIGHT: 0px; PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #000000 2px solid; MARGIN-RIGHT: 0px">
<DIV style="FONT: 10pt arial">----- Original Message ----- </DIV>
<DIV
style="BACKGROUND: #e4e4e4; FONT: 10pt arial; font-color: black"><B>From:</B>
<A title=marquardt.christian@GMAIL.COM
href="mailto:marquardt.christian@GMAIL.COM">Christian Marquardt</A> </DIV>
<DIV style="FONT: 10pt arial"><B>To:</B> <A title=GRADSUSR@LIST.CINECA.IT
href="mailto:GRADSUSR@LIST.CINECA.IT">GRADSUSR@LIST.CINECA.IT</A> </DIV>
<DIV style="FONT: 10pt arial"><B>Sent:</B> Wednesday, March 19, 2008 11:23
AM</DIV>
<DIV style="FONT: 10pt arial"><B>Subject:</B> Re: Problems with new ECMWF
monthly forecast files</DIV>
<DIV><BR></DIV>If I remember correctly, there used to be a -e option to
gribmap which takes care of some filler bits ECMWF tends to add to their grib
records. It won't work without that option.<BR><BR>BTW, you did use gribmap,
did you?<BR><BR>Hope this helps,<BR><BR> Christian.<BR><BR><BR><BR>
<DIV class=gmail_quote>On Wed, Mar 19, 2008 at 10:39 AM, Daniele Gandini
<<A
href="mailto:Daniele.Gandini@arpa.piemonte.it">Daniele.Gandini@arpa.piemonte.it</A>>
wrote:<BR>
<BLOCKQUOTE class=gmail_quote
style="PADDING-LEFT: 1ex; MARGIN: 0pt 0pt 0pt 0.8ex; BORDER-LEFT: rgb(204,204,204) 1px solid">
<DIV bgcolor="#ffffff">
<DIV><FONT face=Arial size=2>Dear GrADS users,</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>I run GrADS 1.8sl11 on Sun Solaris Blade-1000
5.9</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Last week there was a change in ECMWF monthly
forecast files</FONT></DIV>
<DIV><FONT face=Arial size=2>and now I get the message "Entire grid
undefined" when I try</FONT></DIV>
<DIV><FONT face=Arial size=2>to display a variable.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Here is an example of wgrib output for an old
file which successfully works:</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>wgrib -V EMEMO1500_20080306_011_032.grb |
more<BR>rec 1:0:date 2008030600 Z kpds5=129 kpds6=100 kpds7=500
levels=(1,244) grid=255<BR>500 mb ECext=23 class=operations type=Ensemble
means stream=1206 Version=0001 <B>5-11d ave:<BR></B> Z=Geopotential
[m**2 s**-2]<BR> timerange 3 P1 5 P2 11 TimeU 2 nx 62 ny 28 GDS
grid 0 num_in_ave 14 missing 0<BR> center 98 subcenter 0 process 130
Table 128 scan: WE:NS winds(N/S)<BR> latlon: lat 70.500000 to
30.000000 by 1.500000 nxny
1736<BR> long
-51.000000 to 40.500000 by 1.500000, (62 x 28) scan 0 mode 128 bdsgrid
1<BR> min/max data 49641.9 56994.9 num bits 16 BDS_Ref
49641.9 DecScale 0 BinScale -3</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>And now wgrib output for new troublesome
file:</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>wgrib -V EMEMO1500_20080313_011_032.grb |
more</FONT></DIV>
<DIV><FONT face=Arial size=2>rec 1:0:date 2008031300 Z kpds5=129 kpds6=100
kpds7=500 levels=(1,244) grid=255<BR>500 mb ECext=30 class=operations
type=Ensemble means stream=ensemble forecasts Version=0001 <B>96-264hr
ave:<BR></B> Z=Geopotential [m**2 s**-2]<BR> timerange 3 P1 16
P2 44 TimeU 11 nx 62 ny 28 GDS grid 0 num_in_ave 14 missing
0<BR> center 98 subcenter 0 process 131 Table 128 scan: WE:NS
winds(N/S)<BR> latlon: lat 70.500000 to 30.000000 by
1.500000 nxny
1736<BR> long
-51.000000 to 40.500000 by 1.500000, (62 x 28) scan 0 mode 128 bdsgrid
1<BR> min/max data 49871.8 56654.3 num bits 16 BDS_Ref
49871.8 DecScale 0 BinScale -3</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>As you can see, time definition is different,
in hours instead of months.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>I create ctl file with grb2ctl.pl:</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>perl $HOME/bin/grib2ctl.pl -verf
EMEMO1500_20080313_011_032.grb EMEMO1500_20080313_011_032.gmp >
EMEMO1500_20080313_011_032.ctl</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>more
EMEMO1500_20080313_011_032.ctl</FONT></DIV>
<DIV><FONT face=Arial size=2><BR>dset
^EMEMO1500_20080313_011_032.grb<BR>index
^EMEMO1500_20080313_011_032.gmp<BR>undef 9.999E+20<BR>title
EMEMO1500_20080313_011_032.grb<BR>* analysis 2008031300 1stforecast 96-264
hr<BR>* produced by grib2ctl v0.9.12.5p10<BR>dtype grib 255<BR>options
yrev<BR>ydef 28 linear 30.000000 1.5<BR>xdef 62 linear -51.000000
1.500000<BR>tdef 4 linear 00Z24mar2008 7dy<BR>zdef 2 levels 850 500<BR>vars
6<BR>t2m 0 167,1,0 ** 2 metre temperature K<BR>mslmsl 0
151,1,0 ** Mean sea level pressure Pa<BR>probprec 0
228,1,0 ** Pot. evaporation [kg/m^2]<BR>probsnow 0 144,1,0
** Volumetric soil moisture [fraction]<BR>tprs 2 130,100,0 ** Temperature
K<BR>zprs 2 129,100,0 ** Geopotential (at the surface = orography) m**2
s**-2<BR>ENDVARS</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>I run GrADS:</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>grads -cl "open
EMEMO1500_20080313_011_032.ctl"</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>ga-> d zprs<BR>Cannot contour grid - all
undefined values</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>In attachment you have both files.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Thank you in advance for every
suggestions.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Daniele Gandini</FONT></DIV>
<DIV><FONT face=Arial size=2>ARPA Piemonte</FONT></DIV>
<DIV><FONT face=Arial size=2>Servizio Meteorologico</FONT></DIV>
<DIV><FONT face=Arial size=2>Turin - Italy</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><BR></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial
size=2></FONT> </DIV></DIV></BLOCKQUOTE></DIV><BR></BLOCKQUOTE></BODY></HTML>