[RegCNET] Urgent help need on calculte the precipitation value

[Satyaban Bishoyi Ratna]

Hi Rajesh,
You are doing right thing but the easy way is to make the time average at first and than area average.
Ex:
rain = ave(tpr, t=1, t=240)
rain1=  aave(rain, lon=90, lon=93, lat=20, lat=23)

This is your desired average rainfall. You can multiply with 240 in the first line to get the total rainfall.

Thanks.
Satyaban

On Thu, 26 Aug 2010 15:55:06 +0530  wrote
>Hello All,
 
I have divided India into 31 parts like 20-23N & 90-93E so on... I need pricipitation (tpr) over all 31 parts.
Below are the log that I am using on grads to find out precipitation for all 31 parts.
After wrinting this command we got 240 value and then I sum the all 240 value to get monthly precipitation value.
 
ga-> open SRF.2009090100.ctl 
Scanning description file:  SRF.2009090100.ctl
Data file SRF.2009090100 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: 2009:9:1:3 2009:9:1:3 
E set to 1 1 
ga-> set t 1 240
Time values set: 2009:9:1:3 2009:10:1:0 
ga-> d tloop(aave(tpr,lon=90,lon=93,lat=20,lat=23))
 
Here my question is , whether I am using right path to calculate the precipitation value.
 
Please reply as soon as possible
  
 
-- 
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


-------------------------------------------------------------------------------
Dr. Satyaban Bishoyi Ratna
Computational Atmospheric Sciences Team
Scientific & Engineering Computing Group
Centre for Development of Advanced Computing (C-DAC)
Pune University Campus, Ganeshkhind
Pune - 411 007, Maharashtra, India
E-mail: satyabanb at cdac.in
Tel (O): +91-020-2570 4226,  (Cell): +91-094223 56326
--------------------------------------------------------------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20100826/0bc2351d/attachment.html>