hi Waheed,<br><br>1. aphrodite daily precipitation data can be downloaded in a netcdf file format via: <a href="http://www.chikyu.ac.jp/precip/cgi-bin/aphrodite/script/aphrodite_cgi.cgi/register">http://www.chikyu.ac.jp/precip/cgi-bin/aphrodite/script/aphrodite_cgi.cgi/register</a><br>
2. I have attached a script which I used to produce a monthly average precipitation. you can modify it depending on your need. for example: you have to change "ave" with "sum" which in this case, can produce a monthly total rainfall. Note: I did not use looping to accommodate leap years. it works in GrADS 2.0 or later version because of sdfwrite command which can able to write your output to a netcdf file format.<br>
<br>hope it helps.<br><br>Mar<br><br> <br><br><div class="gmail_quote">On Tue, May 25, 2010 at 11:17 AM, Waheed <span dir="ltr"><<a href="mailto:waheed.met@gmail.com">waheed.met@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div dir="ltr"><br><br>
<div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Waheed</b> <span dir="ltr"><<a href="mailto:waheed.met@gmail.com" target="_blank">waheed.met@gmail.com</a>></span><br>
Date: Mon, May 24, 2010 at 5:18 PM<br>
Subject: how to find monthly sum: when using lats4d<br>To: <a href="mailto:dasilva@alum.mit.edu" target="_blank">dasilva@alum.mit.edu</a><br><br><br>
<div dir="ltr">
<div class="gmail_quote"><br>
<div dir="ltr">
<p>hello<br> <br>dear i am converting daily precipitaion data (aphrodite) form 1961 to 2007 from ctl type to nc(netcdf) type , i have done it but now i want to find the monthly total from 1961-2007 using lats4d<br> but i am not understading how to do?<br>
the script that i have wrtten is as :</p>
<p>'reinit'<br>'open data-aphd\APHRO_PAK_005deg_V0902.1961-2007/APHRO_PAK_005deg_V0902.ctl'<br>ti=1<br>tf=ti+30<br>y=1961<br>while (y<=2007)<br>while (ti<=16832)<br>'set t 'ti<br>'lats4d -v -i <br>
data-aphd\APHRO_PAK_005deg_V0902.1961-2007/APHRO_PAK_005deg_V0902.ctl -ftype ctl -o waheed/JAN'y' -func sum(precip,t='ti',t='tf')'<br>ti=ti+365<br>'set t 'ti<br>'q time'<br>say result<br>
res=subwrd(result,3) <br>yf=substr(res,9,4)<br>if (y=yf)<br>ti=ti+1<br>'set t 'ti<br>endif<br>y=y+1<br>endwhile<br>endwhile</p>
<p> <br> <br> <br>any help <br> <br> <br>thanx in advance <br> <br>Regards<br></p></div><br></div><br><br clear="all"><br>-- <br>waheed iqbal<br>(Meteorologist)<br>Pakistan Meteorological Department<br>R & D islamabad<br>
+923334547814<br></div></div><br><br clear="all"><br>-- <br>waheed iqbal<br>(Meteorologist)<br>Pakistan Meteorological Department<br>R & D islamabad<br>+923334547814<br></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br>