[RegCNET] Re : Fwd: Re : Fwd: Suggestion
Bamba Sylla
syllabamba at yahoo.fr
Mon Oct 25 11:29:52 CEST 2010
Hi again,
Always remember that the unit for rainfall (even the 3-hourly) is mm/day.
Now, if you use this method then you should divide by 8 for each month. However,
the easiest way is to always average and multiply by 30 or 31 depending on the
lenght of the month.
Hope this will help.
Cheers!
Ps: you can try the two methods and see if you gonna have the same results.
***************************************************
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 at ictp.it, syllabamba at yahoo.fr
***************************************************
________________________________
De : rajesh maurya <rksmaurya at gmail.com>
À : regcnet <regcnet at lists.ictp.it>
Cc : Bamba Sylla <syllabamba at yahoo.fr>
Envoyé le : Lun 25 octobre 2010, 11h 12min 10s
Objet : Fwd: Re : [RegCNET] Fwd: Suggestion
Hi Bamba Sylla,
I am working for your suggestion d sum(tpr,t=1,t=240) and found for
model and observed data are2983 mm and897 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 at yahoo.fr>
Date: Mon, Oct 25, 2010 at 1:31 PM
Subject: Re : [RegCNET] Fwd: Suggestion
To: rajesh maurya <rksmaurya at gmail.com>, regcnet <regcnet at 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 at ictp.it, syllabamba at yahoo.fr
***************************************************
________________________________
De : rajesh maurya <rksmaurya at gmail.com>
À : regcnet <regcnet at 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 at yahoo.com>
Date: Mon, Oct 25, 2010 at 11:11 AM
Subject: Re: [RegCNET] Suggestion
To: rajesh maurya <rksmaurya at 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 at gmail.com>
To: regcnet <regcnet at 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20101025/f149f2a4/attachment.html>
More information about the RegCNET
mailing list