Plot multiple files with one control file

[Sergey Varlamov]

Dear Wang Lizhi,

Grads do not support templates in index file names, but:
1) please change "index" line in your ctl file to some fixed name, like

index ^flx.ft06.avrg.grib.idx

then create this 'single' index file:

gribmap -v -i ncep.ctl -0

It must work, -v is for the verbose output:
if work goes, you will see a lot of MATCH strings,
else - look on your ctl file again.

Best regards,

Sergey Varlamov
vsm at riam.kyushu-u.ac.jp


Wang Lizhi wrote:

>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
>>
>>
>>
>
>
>
>