<!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>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 <STRONG>5-11d ave:<BR></STRONG> 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 <STRONG>96-264hr
ave:<BR></STRONG> 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> </DIV>
<DIV><BR></DIV></FONT>
<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></BODY></HTML>