[RegCNET] Re : Fwd: Re : Fwd: Suggestion
Bamba Sylla
syllabamba at yahoo.fr
Mon Oct 25 12:37:47 CEST 2010
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/5f4fee63/attachment.html>
More information about the RegCNET
mailing list