Plot multiple files with one control file

[Wang Lizhi]

Dear Users,

I download the NCEP reanalysis-2 dataset from the site http://nomad3.ncep.noaa.gov/pub/raid1/reanalysis-2/month/  and use the following step to convert the GRIB data
grib2ctl.pl datfile > datfile.ctl
gribmap -i datfile.ctl -0

then every month the file as following: datfile, datfile.ctl and datfile.idx

I want to display multiple file with one control file, so I have the whole control file named ncep.ctl as :

dset ^flx.ft06.%y4%m2.avrg.grib
index ^flx.ft06.%ym%m2.avrg.grib.idx
options template
.....
*  dataset from 197901 to 200709
tdef 345 linear 00Z01jan1979 1mo
vars 54
.....
endvars

when I display the file, the following error is displayed:
Open Error:  Can't open Station/Index map file flx.ft06.%y4%m2.avrg.grib.idx
  --> The invalid description file record is:
  --> index ^flx.ft06.%y4%m2.avrg.grib.idx

How can I read the multiple index files with grads?

Grads:  v1.8SL11 32-bit   OS: SuSE 10 with kernel 2.6.16.27-0.9-smp x86_64

thanks

lizhi


----- Original Message -----
From: "Arlindo da Silva" <dasilva at ALUM.MIT.EDU>
To: <GRADSUSR at LIST.CINECA.IT>
Sent: Wednesday, October 10, 2007 5:12 AM
Subject: Re: Plot multiple files with one control file


> On 10/8/07, ashish no routray <ashishroutray at rediffmail.com> wrote:
>>
>> Dear Users
>> I need your kind help to plot the multile binary files in using one .ctl
>> file.
>> I downloaded TRMM 03 hourly data for a period and also .ctl file from TRMM
>> site
>>
>
> It is always a good idea to give the URL so that someone could try to
> download the same files and reproduce the problem... Also, provide the
> output of "q config" and "!uname -a" so that we can know which version you
> are using...
>
> . we are trying to plot 24hrs accumulated rainfall from all 9
> files.Filesare look like:
>>
>> 3B42.020601.00z.6.precipitation.bin
>> 3B42.020601.03z.6.precipitation.bin
>> ......
>> ......
>> 3B42.020602.00z.6.precipitation.bin
>>
>> Here I am giving the .ctl file as follow:
>>
>
> I am assuming that the individual ctl provided with data works fine.
>
> ***************************************************
>> dset ./3B42.%y2%m2%d2.%h2z.6.precipitation.bin
>> options template byteswapped
>>
> The option "byteswapped" is deprecated. Use either "big_endian" or
> "little_endian", depending on the endianess of the machine TRMM used to
> write the data.
>
> title TRMM 3B42 V6 three hourly TRMM rainfall
>> undef -9999.0
>> xdef 1440 linear -179.875 0.25000000
>> ydef 400  linear -49.8750000 0.25000000
>> zdef 1 levels 1000
>>
> So far, so good...
>
> tdef 9 linear 00:00Z01jan1998 3hr
>>
>
> Why "9"? You need the enter the total number of time steps from 1998 to
> 2007.
>
>
>> * end_time 21:00Z31aug2007 (this_is_comment_line)
>> vars 1
>> r        0  99 Hourly Rain Rate (mm/hr)
>> endvars
>>
>> **************************************************
>>
>> We are running the .ctl in IBM(P5) machine.
>> It showing Gridded are undefined. I hope it could not open the files. when
>> I plot individual files alone it is ploting properly.
>>
>
> Can you plot t=1 through t=9 with the template version of the  ctl?
>
> --
> Arlindo da Silva
> dasilva at alum.mit.edu
>