Plot multiple files with one control file

[Wang Lizhi]

Dear Sergey Varlamov,

It works.

Thank you very much.

Best regards,

lizhi

----- Original Message -----
From: "Sergey Varlamov" <vsm at RIAM.KYUSHU-U.AC.JP>
To: <GRADSUSR at LIST.CINECA.IT>
Sent: Friday, October 12, 2007 5:14 PM
Subject: Re: Plot multiple files with one control file


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