Template files

[Isabella Bovolo]

Dear Oscar,

thank you for your reply.

I have the following files to use as an example (currently 5 days but I
would like to be able to repeat the procedure for several months/years of
data):

gpfs1958010100cp_24h.ctl
gpfs1958010100cp_24h.grb
gpfs1958010100cp_24h.grb.idx
gpfs1958010200cp_24h.ctl
gpfs1958010200cp_24h.grb
gpfs1958010200cp_24h.grb.idx
gpfs1958010300cp_24h.ctl
gpfs1958010300cp_24h.grb
gpfs1958010300cp_24h.grb.idx
gpfs1958010400cp_24h.ctl
gpfs1958010400cp_24h.grb
gpfs1958010400cp_24h.grb.idx
gpfs1958010500cp_24h.ctl
gpfs1958010500cp_24h.grb
gpfs1958010500cp_24h.grb.idx

The control file of gpfs1958010100cp_24h.ctl is:
dset ^gpfs1958101000cp_24h.grb
index ^gpfs1958101000cp_24h.grb.idx
undef 9.999E+20
title ^gpfs1958101000cp_24h.grb
* produced by /usr/local/badcdat/bin/badcGrib.py
dtype grib 255
options yrev
ydef 181 linear -90.000000 1
xdef 360 linear 0.000000 1.000000
tdef 1 linear 00Z11oct1958 1dy
zdef 1 linear 1 1
vars 1
cp  0 143,1,0 Convective precipitation (accum.) (m)
ENDVARS


I can open each file individually and plot the contents of each file:

ga-> query file 1
File 1 : ^gpfs1958010200cp_24h.grb
   Descriptor: gpfs1958010200cp_24h.ctl
   Binary: gpfs1958010200cp_24h.grb
   Type = Gridded
   Xsize = 360  Ysize = 181  Zsize = 1  Tsize = 1  Esize = 1
   Number of Variables = 1
      cp  0  143  Convective precipitation (accum.) (m)
ga-> d cp.1
Contouring: 0 to 0.07 interval 0.01

But when I try to open more than one file at once, I cannot seem to be
able to plot any content except that of the first file:

File 1 : ^gpfs1958010200cp_24h.grb
   Descriptor: gpfs1958010200cp_24h.ctl
   Binary: gpfs1958010200cp_24h.grb
   Type = Gridded
   Xsize = 360  Ysize = 181  Zsize = 1  Tsize = 1  Esize = 1
   Number of Variables = 1
      cp  0  143  Convective precipitation (accum.) (m)
ga-> q file 2
File 2 : ^gpfs1958010300cp_24h.grb
   Descriptor: gpfs1958010300cp_24h.ctl
   Binary: gpfs1958010300cp_24h.grb
   Type = Gridded
   Xsize = 360  Ysize = 181  Zsize = 1  Tsize = 1  Esize = 1
   Number of Variables = 1
      cp  0  143  Convective precipitation (accum.) (m)
ga-> d cp.1
Contouring: 0 to 0.07 interval 0.01
ga-> d cp.2
Data Request Warning:  Request is completely outside file limits
   Entire grid contents are set to missing data
   Grid limits of file:     X = 1 360  Y = 1 181  Z = 1 1  T = 1 1  E = 1 1
   Grid limits of request:  X = 1 361  Y = 1 181  Z = 1 1  T = 0 0  E = 1 1
   Warning issued for variable = cp.2
Cannot contour grid - all undefined values

(Note I can plot the second file if I open just that file:)
ga-> open gpfs1958010300cp_24h
Scanning description file:  gpfs1958010300cp_24h
Data file gpfs1958010300cp_24h.grb is open as file 1
LON set to 0 360
LAT set to -90 90
LEV set to 1 1
Time values set: 1958:1:4:0 1958:1:4:0
E set to 1 1
ga-> d cp
Contouring: 0 to 0.07 interval 0.01


Isabella


On Thu, 24 Sep 2009, Oscar Andres Frumento wrote:

> Isabella,
>
> what is the original ctl file?
> how many files do you have to open?
> can you display cp variable opening only 1 file?
>
> oscar
>