[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