[RegCNET] Re : Re : Fwd: Re : Fwd: Suggestion
Bamba Sylla
syllabamba at yahoo.fr
Mon Oct 25 12:52:12 CEST 2010
Yes Dogar,
I understand that he wants cumulative monthly rainfall but from raw RegCM output
that are actually 3-hourly not daily. So summing will not work, you will
overestimate 8 times.
It only works for the observations he has if they are daily of course.
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 : mubashar dogar <mubashardogar at yahoo.com>
À : rajesh maurya <rksmaurya at gmail.com>; Coppola Erika <coppolae at ictp.it>; Bi X
Q <bixq at ictp.it>; Bamba Sylla <syllabamba at yahoo.fr>
Cc : regcnet <regcnet at lists.ictp.it>
Envoyé le : Lun 25 octobre 2010, 12h 45min 15s
Objet : Re: Re : [RegCNET] Fwd: Re : Fwd: Suggestion
Dear Bamba,
Do you understand that he wants monthly accumulated precipitation or the rate of
precipitation?
If he wants accumulated precipitation then making average and multiply by 30,
31, 28 or 29 is same if he just adds all the daily precipitation. Just do it
manually and you will get the same answer.
Cheers :)
--- On Mon, 25/10/10, Bamba Sylla <syllabamba at yahoo.fr> wrote:
>From: Bamba Sylla <syllabamba at yahoo.fr>
>Subject: Re : [RegCNET] Fwd: Re : Fwd: Suggestion
>To: "mubashar dogar" <mubashardogar at yahoo.com>, "rajesh maurya"
><rksmaurya at gmail.com>, "Coppola Erika" <coppolae at ictp.it>, "Bi X Q"
><bixq at ictp.it>
>Cc: "regcnet" <regcnet at lists.ictp.it>
>Date: Monday, 25 October, 2010, 3:37 PM
>
>
>Hi again
>
>
>It seems that it each of us has its own perception but just to clarify:
>The rainfall unit in RegCM is mm/day (even for the 3-hourly). Then, just take
>few seconds and think about: to get the monthly cumulative rainfall (in mm or
>mm/month) you need to average (forget about summing) for all timesteps and
>multiply by 30, 31, 28 or 29 depending on the month your consider. Bi, am I
>right?
>
>
>Again, Rajesh if your observed data is daily (I guess it is), just sum up all
>timesteps and you will have your cumulative rainfall in mm.
>
>
>Hope it clear now.
>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 : mubashar dogar <mubashardogar at yahoo.com>
>À : rajesh maurya <rksmaurya at gmail.com>; Coppola Erika <coppolae at ictp.it>
>Cc : regcnet <regcnet at lists.ictp.it>
>Envoyé le : Lun 25 octobre 2010, 12h 20min 47s
>Objet : Re: [RegCNET] Fwd: Re : Fwd: Suggestion
>
>
>Dear all,
>
>To convert model into mm/month, you need only to sum up the daily precipitation
>or hourly precipitation in each month and this will automatically give mm/month.
>
>
>Regards,
>
>Mubashar Dogar
>
>
>
>--- On Mon, 25/10/10, Coppola Erika <coppolae at ictp.it> wrote:
>
>
>>From: Coppola Erika <coppolae at ictp.it>
>>Subject: Re: [RegCNET] Fwd: Re : Fwd: Suggestion
>>To: "rajesh maurya" <rksmaurya at gmail.com>
>>Cc: "regcnet" <regcnet at lists.ictp.it>
>>Date: Monday, 25 October, 2010, 2:39 PM
>>
>>
>>Dear Rajesh,
>>I think the best way to sort your problem out is to work in term of unit
>>you have in the model and in the observation.
>>The model SRF output is in mm/day, therefore if you want to have the mm/
>>month you need to average the precipitation in that month and multiply by
>>the number of day in the month.
>>
>>This assume you have your observation in mm/month.
>>Best wishes
>>Erika
>>
>>> 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 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
>>> _______________________________________________
>>> RegCNET mailing list
>>> RegCNET at lists.ictp.it
>>> https://lists.ictp.it/mailman/listinfo/regcnet
>>
>>
>>--
>>*************************************************
>>Coppola Erika , Ph.D
>>Earth System Physics Section
>>The Abdus Salam International Centre for Theoretical Physics
>>Strada Costiera 11
>>34100 Trieste, ITALY
>>Phone: + 39 040 2240 407
>>Fax: + 39 040 2240 449
>>
>>email: coppolae at ictp.it
>>*************************************************
>>
>>
>>_______________________________________________
>>RegCNET mailing list
>>RegCNET at lists.ictp.it
>>https://lists.ictp.it/mailman/listinfo/regcnet
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20101025/ec1bcf89/attachment.html>
More information about the RegCNET
mailing list