Hi Bamba Sylla,
I am working for your suggestion d sum(tpr,t=1,t=240) and found for model and observed data are 2983 mm and 897 mm for the same Grid and same time..
Observed Works....
ga-> open rf_2007.ctl
Scanning description file: rf_2007.ctl
Data file rf_2007.grd is open as file 1
LON set to 66.5 100.5
LAT set to 6.5 38.5
LEV set to 1 1
Time values set: 2007:1:1:0 2007:1:1:0
E set to 1 1
ga-> deffine r=sum(rf,t=152,t=273)
Unknown command: deffine
ga-> define r=sum(rf,t=152,t=273)
SUMing. dim = 3, start = 152, end = 273
Define memory allocation size = 9240 bytes
ga-> d aave(r,lon=73,lon=81,lat=15,lat=31)
Result value = 896.992
Model works....
ga-> open /home/praveen/singh/RegCM/2007/output/SRF.2007060100.ctl
Scanning description file: /home/praveen/singh/RegCM/2007/output/SRF.2007060100
.ctl
Data file /home/praveen/singh/RegCM/2007/output/SRF.2007060100 is open as file 1
LON set to 49.4242 111.019
LAT set to -10.48 39.77
LEV set to 1008.43 1008.43
Time values set: 2007:6:1:3 2007:6:1:3
E set to 1 1
ga-> define r1=sum(tpr,t=1,t=240)
SUMing. dim = 3, start = 1, end = 240
Define memory allocation size = 133584 bytes
Name already DEFINEd: r1. Will be deleted and replaced.
ga-> close 1
File 1 has been closed
ga-> open /home/praveen/singh/RegCM/2007/output/SRF.2007070100.ctl
Scanning description file: /home/praveen/singh/RegCM/2007/output/SRF.2007070100
.ctl
Data file /home/praveen/singh/RegCM/2007/output/SRF.2007070100 is open as file 1
LON set to 49.4242 111.019
LAT set to -10.48 39.77
LEV set to 1008.43 1008.43
Time values set: 2007:7:1:3 2007:7:1:3
E set to 1 1
ga-> define r2=sum(tpr,t=1,t=248)
SUMing. dim = 3, start = 1, end = 248
Define memory allocation size = 133584 bytes
Name already DEFINEd: r2. Will be deleted and replaced.
ga-> close 1
File 1 has been closed
ga-> open /home/praveen/singh/RegCM/2007/output/SRF.2007080100.ctl
Scanning description file: /home/praveen/singh/RegCM/2007/output/SRF.2007080100
.ctl
Data file /home/praveen/singh/RegCM/2007/output/SRF.2007080100 is open as file 1
LON set to 49.4242 111.019
LAT set to -10.48 39.77
LEV set to 1008.43 1008.43
Time values set: 2007:8:1:3 2007:8:1:3
E set to 1 1
ga-> define r3=sum(tpr,t=1,t=248)
SUMing. dim = 3, start = 1, end = 248
Define memory allocation size = 133584 bytes
Name already DEFINEd: r3. Will be deleted and replaced.
ga-> close 1
File 1 has been closed
ga-> open /home/praveen/singh/RegCM/2007/output/SRF.2007090100.ctl
Scanning description file: /home/praveen/singh/RegCM/2007/output/SRF.2007090100
.ctl
Data file /home/praveen/singh/RegCM/2007/output/SRF.2007090100 is open as file 1
LON set to 49.4242 111.019
LAT set to -10.48 39.77
LEV set to 1008.43 1008.43
Time values set: 2007:9:1:3 2007:9:1:3
E set to 1 1
ga-> define r4=sum(tpr,t=1,t=240)
SUMing. dim = 3, start = 1, end = 240
Define memory allocation size = 133584 bytes
Name already DEFINEd: r4. Will be deleted and replaced.
ga-> define r=r1+r2+r3+r4
Define memory allocation size = 133584 bytes
Name already DEFINEd: r. Will be deleted and replaced.
ga-> d aave(r,lon=73,lon=81,lat=15,lat=31)
Result value = 2983.43
Please help me.....
Thanks
Rajesh Kumar Maurya
M.Sc(Tech)- Geophysics(Meteorology); Junior Research Fellow & Research Scholar
Department of Geophysics
Banaras Hindu University
Varanasi 221005
---------- Forwarded message ----------
From:
Bamba Sylla <syllabamba@yahoo.fr>Date: Mon, Oct 25, 2010 at 1:31 PM
Subject: Re : [RegCNET] Fwd: Suggestion
To: rajesh maurya <
rksmaurya@gmail.com>, regcnet <
regcnet@lists.ictp.it>
Hi Rajesh,
If you have some doubts just sum up values for all timesteps and you will have the cumulative rainfall for the month.
d sum(tpr,t=1,t=240)
cheers !
***************************************************
Mouhamadou Bamba SYLLA, PhD
Post-Doc Fellow, Earth System Physics
International Centre for Theoretical Physics
Strada Costiera 11, 34100 Trieste, Italy
Office tel: 0039 040 2240 385
Cell: 0039 34 71 94 06 55
emails:
msylla@ictp.it,
syllabamba@yahoo.fr ***************************************************
De : rajesh maurya <rksmaurya@gmail.com>
À : regcnet <regcnet@lists.ictp.it>
Envoyé le : Lun 25 octobre 2010, 9h 54min 28s
Objet : [RegCNET] Fwd: Suggestion
Hello All,
I am plotting the contour for precipitation , in that it is showing 3 hourly data( In SRF file ). Here I am plotting monthly precipitation the command i am using is
ga-> d ave(tpr,t=1,t=240)*30 (mm) and found the value 542, 479,455,512,513 mm ( year- 2000, 2001,2002, 2003, 2004) but observed data (IMD Delhi, India) same time and same grid are found 799, 797, 658, 854, 815 mm respectively.
again using Mr. Md. Mizanur Rahman suggestion using command ga-> d ave(tpr,t=1,t=240)*8*30 (mm) and found the values 4340, 3832, 3643, 4099, 4099, 4103 mm respectively.
Please let me know If i am missing something.
Thanks
Rajesh Kumar Maurya
M.Sc(Tech)- Geophysics(Meteorology); Junior Research Fellow & Research Scholar
Department of Geophysics
Banaras Hindu University
Varanasi 221005
---------- Forwarded message ----------
From:
mizanur rahman <mrahman426@yahoo.com>
Date: Mon, Oct 25, 2010 at 11:11 AM
Subject: Re: [RegCNET] Suggestion
To: rajesh maurya <
rksmaurya@gmail.com>
hi rajesh
I think your formula is wrong. Because you are multiplying only by 30, it is not correct. It should be multiplied by 24*30, because you are doing average (t=1, t=240), after averging you will get one hour rainfall after that you will multiply 24 to average rainfall then you will get 24 hours rainfall. i, e complete formula for one month will be ave(t=1, t=240)*24*30 then you will get one month rainfall (mm).
regards
Md. Mizanur Rahman
Research Officer
SAARC Meteorological Research Centre (SMRC)
Plot No# E-4/C, Agargaon
Sher-e-Bangla Nagar
Dhaka-1207, BANGLADESH
Cell Phone: +88-01726652541
Phone: +88-02-9144374 (office)
Fax: +88-02-9144373
From: rajesh maurya <rksmaurya@gmail.com>
To: regcnet <regcnet@lists.ictp.it>
Sent: Fri, October 22, 2010 4:51:04 PM
Subject: [RegCNET] Suggestion
Hello All,
I am plotting the contour for precipitation , in that it is showing 3 hourly data( In SRF file ). Here I am plotting monthly precipitation the command i am using is ga-> d ave(tpr,t=1,t=240)*30 (mm) My question is about the command whether it is right one to get correct value or else .
Please let me know If i am missing something.
--
Rajesh Kumar Maurya
M.Sc(Tech)- Geophysics(Meteorology); Junior Research Fellow & Research Scholar
Department of Geophysics
Banaras Hindu University
Varanasi 221005