Dear all users<br>this is my gs file; the  ctl file is given below. i want to file u-wind anomaly fro the grid 0-90N and 0-357.5<br>at 1000 t0 10 hPa level for a single day(time= 1feb2008)<br>The grd file has shows the output at a single level<br>
The problem is output result didnot vary at different levels <br><br>&#39;reinit&#39;<br>&#39;sdfopen /media/resmi-1/windows/ncep/uwnd/200h_uwnd/<a href="http://uwnd.2000_10.nc">uwnd.2000_10.nc</a>&#39;<br>&#39;set lat 0 90&#39;<br>
&#39;set lon 0 357.5&#39;<br>*&#39;set lev 850&#39;<br>&#39;set time 1feb2008&#39;<br><br>&#39;set gxout fwrite&#39;<br>&#39;set fwrite stj850_26j.grd&#39;<br>n=1<br>while(n&lt;=17)<br><br>&#39;set z &#39;n <br>&#39;define tt=uwnd(time=1feb2000)&#39;<br>
&#39;define tt1=uwnd(time=1feb2001)&#39;<br>&#39;define tt2=uwnd(time=1feb2002)&#39;<br>&#39;define tt3=uwnd(time=1feb2003)&#39;<br>&#39;define tt4=uwnd(time=1feb2004)&#39;<br>&#39;define tt5=uwnd(time=1feb2005)&#39;<br>&#39;define tt6=uwnd(time=1feb2006)&#39;<br>
&#39;define tt7=uwnd(time=1feb2007)&#39;<br>&#39;define tt8=uwnd(time=1feb2008)&#39;<br>&#39;define tt9=uwnd(time=1feb2009)&#39;<br>&#39;define tt10=uwnd(time=1feb2010)&#39;<br><br>&#39;define daymean=(tt+tt1+tt2+tt3+tt4+tt5+tt6+tt7+tt8+tt9+tt10)/11&#39;<br>
*&#39;d daymean&#39;<br>&#39;define dayanom=tt9-daymean&#39;<br>*&#39;d dayanom&#39;<br><br>n=n+1<br>endwhile<br>&#39;disable fwrite&#39;<br><br>Control file<br><br> DSET  /media/resmi-1/windows/ncep/uwnd/200h_uwnd/stj850_26j.grd<br>
 undef -999.99<br> xdef 144 linear 0 2.5<br> ydef 37 linear 0 2.5<br> zdef 17 levels 1000 925 850 700 600 500 400 300 200 250 150 100 70 50 30 20 10<br> tdef 1 linear 01feb2008 1dy<br> TITLE  stj850_26j.grd  144x37 grd<br>
 vars 1<br> uwnd  0 -999 uwnd temperature<br> endvars<br><br>