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@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
> _______________________________________________
> RegCNET mailing list
>
RegCNET@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@ictp.it*************************************************
_______________________________________________
RegCNET mailing list
RegCNET@lists.ictp.ithttps://lists.ictp.it/mailman/listinfo/regcnet