[RegCNET] Re : Fwd: Re : Fwd: Suggestion
bixq
bixq at ictp.it
Mon Oct 25 12:44:38 CEST 2010
On Mon, 25 Oct 2010, Bamba Sylla wrote:
> 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?
Yes, I think.
My suggestion is that we should also do the comparison of precipitation
with observed with the same unit mm/day. You know, most of the observed
data are actually using this unit (like CRU, CMAP).
>
> 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
>>
>
>
>
>
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
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
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
More information about the RegCNET
mailing list