[RegCNET] regarding 6 hourly data files
Dr. Md. Nazrul Islam
mnislam at phy.buet.ac.bd
Fri Jul 13 11:57:12 CEST 2007
Dear Dr. Bi,
I am not using this script. I extract rainfall for each record using the
script below. Then convert data using Fortran code to get text values for
each record. After that I use Excel for daily or monthly average.
'reinit'
'open gridded.ctl'
'set fwrite out.dat'
'set gxout fwrite'
'set t 1 120'
'd aave(tpr, lon=88.8, lon=92, lat=21.8, lat=25.2)'
'disable fwrite'
'close 1'
return
The suggested commands in my last mail are in GrADS manual for gathering
data in a certain interval. That is why I asked Mishra to check with his
outputs.
Dr. Md. Nazrul Islam
Associate Professor
Department of Physics
Bangladesh University of Engineering and Technology
Dhaka-1000, Bangladesh
http://www.buet.ac.bd/phy/nazrul.htm
Phone: 880 2 966 5650 (PABX)
FAX: 880 2 861 3026, 880 2 861 3046
---------- Original Message -----------
From: XUNQIANG BI <bixq at ictp.it>
To: "Dr. Md. Nazrul Islam" <mnislam at phy.buet.ac.bd>
Cc: regcnet at lists.ictp.it
Sent: Fri, 13 Jul 2007 10:30:03 +0200 (CEST)
Subject: Re: [RegCNET] regarding 6 hourly data files
> Hi, Dr. Md. Nazrul Islam:
>
> On Fri, 13 Jul 2007, Dr. Md. Nazrul Islam wrote:
>
> > Hi Mishra,
> >
> > I think you can compare your daily average RF with the same extracted
using
> > the following in GrADS script:
> >
> > 'define 00RF=ave(trp,t=1,t=120,4)'
> > 'define 06RF=ave(trp,t=1,t=120,5)'
> > 'define 12RF=ave(trp,t=1,t=120,6)'
> > 'define 18RF=ave(trp,t=1,t=120,7)'
>
> Did you use the above script ? Based on my knowledge, it's totally
> wrong.
>
> The right one should be
>
> 'define 00RF=ave(tpr,t=1,t=120,4)'
> 'define 06RF=ave(tpr,t=2,t=120,4)'
> 'define 12RF=ave(tpr,t=3,t=120,4)'
> 'define 18RF=ave(tpr,t=4,t=120,4)'
> >
> > 'd (00RF+06RF+12RF+18RF)/4'
> >
> > This will display daily rainfall in a monthly file (6 hourly ATM file)
> > To add for other months you may use set dfile2 etc.
> >
>
> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
> Dr. Xunqiang Bi email:bixq at ictp.it
> Earth System Physics Group
> The Abdus Salam ICTP
> Strada Costiera, 11
> P.O. BOX 586, 34100 Trieste, ITALY
> Tel: +39-040-2240302 Fax: +39-040-2240449
> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
------- End of Original Message -------
More information about the RegCNET
mailing list