<div dir="ltr"><div>Dear all</div>
<div> </div>
<div><font style="BACKGROUND-COLOR: #ffffff" color="#ff0000">Grads version: 2.0.a8</font></div>
<div><font style="BACKGROUND-COLOR: #ffffff" color="#ff0000">OS : RHL 5 (32bit).</font></div>
<div> </div>
<div>i am trying to find the monthly average precipitation from yearly nc files  of daily data. i want to get ouput in nc format but monthly average for each month.</div>
<div>dkiktev had sent me a script and i made modifications in it. now i m running it for 10 years from 1951-1960 for the month of May . But i am having trouble with it.</div>
<div> </div>
<div>the script is as follows :</div>
<div>* monthly average for the month on May from 1951-1960</div>
<div> </div>
<div>&#39;sdfopen <a href="http://1951.nc">1951.nc</a>&#39;</div>
<div>&#39;set lat 23.025 37.475&#39;</div>
<div>&#39;set lon 60.025 79.975&#39;</div>
<div>&#39;define p51=ave(precip,t=121,t=151)&#39;</div>
<div>&#39;close 1&#39;</div>
<div>&#39;sdfopen <a href="http://1952.nc">1952.nc</a>&#39;</div>
<div>&#39;set lat 23.025 37.475&#39;</div>
<div>&#39;set lon 60.025 79.975&#39;</div>
<div>&#39;define p52=ave(precip,t=122,t=152)&#39;</div>
<div>&#39;define sum1=p51+p52&#39;</div>
<div>&#39;close 1&#39;</div>
<div>&#39;undefine p51&#39;</div>
<div>&#39;undefine p52&#39;</div>
<div>&#39;sdfopen <a href="http://1953.nc">1953.nc</a>&#39;</div>
<div>&#39;set lat 23.025 37.475&#39;</div>
<div>&#39;set lon 60.025 79.975&#39;</div>
<div>&#39;define p53=ave(precip,t=121,t=151)&#39;</div>
<div>&#39;define sum2=sum1+p53&#39;</div>
<div>&#39;undefine sum1&#39;</div>
<div>&#39;undefine p53&#39;</div>
<div>&#39;close 1&#39;</div>
<div>&#39;sdfopen <a href="http://1954.nc">1954.nc</a>&#39;</div>
<div>&#39;set lat 23.025 37.475&#39;</div>
<div>&#39;set lon 60.025 79.975&#39;</div>
<div>&#39;define p54=ave(precip,t=121,t=151)&#39;</div>
<div>&#39;define sum3=sum2+p54&#39;</div>
<div>&#39;undefine sum2&#39;</div>
<div>&#39;undefine p54&#39;</div>
<div>&#39;close 1&#39;</div>
<div>&#39;sdfopen <a href="http://1955.nc">1955.nc</a>&#39;</div>
<div>&#39;set lat 23.025 37.475&#39;</div>
<div>&#39;set lon 60.025 79.975&#39;</div>
<div>&#39;define p55=ave(precip,t=121,t=151)&#39;</div>
<div>&#39;define sum4=sum3+p55&#39;</div>
<div>&#39;undefine sum3&#39;</div>
<div>&#39;undefine p55&#39;</div>
<div>&#39;close 1&#39;</div>
<div>&#39;sdfopen <a href="http://1956.nc">1956.nc</a>&#39;</div>
<div>&#39;set lat 23.025 37.475&#39;</div>
<div>&#39;set lon 60.025 79.975&#39;</div>
<div>&#39;define p56=ave(precip,t=122,t=152)&#39;</div>
<div>&#39;define sum5=sum4+p56&#39;</div>
<div>&#39;undefine sum4&#39;</div>
<div>&#39;undefine p56&#39;</div>
<div>&#39;close 1&#39;</div>
<div>&#39;sdfopen <a href="http://1957.nc">1957.nc</a>&#39;</div>
<div>&#39;set lat 23.025 37.475&#39;</div>
<div>&#39;set lon 60.025 79.975&#39;</div>
<div>&#39;define p57=ave(precip,t=121,t=151)&#39;</div>
<div>&#39;define sum6=sum5+p57&#39;</div>
<div>&#39;undefine sum5&#39;</div>
<div>&#39;undefine p57&#39;</div>
<div>&#39;close 1&#39;</div>
<div>&#39;sdfopen <a href="http://1958.nc">1958.nc</a>&#39;</div>
<div>&#39;set lat 23.025 37.475&#39;</div>
<div>&#39;set lon 60.025 79.975&#39;</div>
<div>&#39;define p58=ave(precip,t=121,t=151)&#39;</div>
<div>&#39;define sum7=sum6+p58&#39;</div>
<div>&#39;undefine sum6&#39;</div>
<div>&#39;undefine p58&#39;</div>
<div>&#39;close 1&#39;</div>
<div>&#39;sdfopen <a href="http://1959.nc">1959.nc</a>&#39;</div>
<div>&#39;set lat 23.025 37.475&#39;</div>
<div>&#39;set lon 60.025 79.975&#39;</div>
<div>&#39;define p59=ave(precip,t=121,t=151)&#39;</div>
<div>&#39;define sum8=sum7+p59&#39;</div>
<div>&#39;undefine sum7&#39;</div>
<div>&#39;undefine p59&#39;</div>
<div>&#39;close 1&#39;</div>
<div>&#39;sdfopen <a href="http://1960.nc">1960.nc</a>&#39;</div>
<div>&#39;set lat 23.025 37.475&#39;</div>
<div>&#39;set lon 60.025 79.975&#39;</div>
<div>&#39;define p60=ave(precip,t=122,t=152)&#39;</div>
<div>&#39;define sum9=sum8+p60&#39;</div>
<div>&#39;undefine sum8&#39;</div>
<div>&#39;undefine p60&#39;</div>
<div>&#39;close 1&#39;</div>
<div>&#39;define pave=sum9/10&#39;</div>
<div>&#39;set sdfwrite <a href="http://waheed.nc">waheed.nc</a>&#39;</div>
<div>&#39;sdfwrite pave&#39;</div>
<div>&#39;quit&#39;</div>
<div> </div>
<div><font color="#cc33cc">It gives an error &quot; define error : no files open yet&quot;</font></div>
<div> </div>
<div>any suggestions for this script or if there is something better please tell .</div>
<div> </div>
<div>Thanx in advance</div>
<div> </div>
<div> </div>
<div>Best regards</div>
<div>-- </div>
<div>Waheed Iqbal</div>
<div> (Meteorologist)<br>Pakistan Meteorological Department<br>R &amp; D Islamabad,Pakistan<br>+923334547814</div>
<div> </div>
<div> </div></div>