[gradsusr] area average of meridional over 40o-55oE at the equator from surface to 700hpa

[syairah atiqah atiq]

Hi, I really need help with my script. I want to find thearea average of meridional wind over 40o-55oE at the equator from surface to700hpa, and I write a scripts as below. But the result was error. Can I knowwhich line was error.


 
Thank you


 
Scripts:


 
 'reinit'

*

 'sdfopen D:\atiqah\nlev.uv.nc'

*

******************************************************

*

   'cJa=ave(u,time=jun1979,time=dec2013,12)'

   'cJla=ave(u,time=jul1979,time=dec2013,12)'

   'cAa=ave(u,time=aug1979,time=dec2013,12)'

*

   'au=(cJa+cJla+cAa)/3'

*

   'cJb=ave(v,time=jun1979,time=dec2013,12)'

   'cJlb=ave(v,time=jul1979,time=dec2013,12)'

   'cAb=ave(v,time=aug1979,time=dec2013,12)'

*

   'av=(cJb+cJlb+cAb)/3'

*

  y.1 =1979

  y.2 =1980

  y.3 =1981

  y.4 =1982

  y.5 =1983

  y.6 =1984

  y.7 =1985

  y.8 =1986

  y.9 =1987

  y.10= 1988

  y.11= 1989

  y.12= 1990

  y.13= 1991

  y.14= 1992

  y.15= 1993

  y.16= 1994

  y.17= 1995

  y.18= 1996

  y.19= 1997

  y.20= 1998

  y.21= 1999

  y.22= 2000

  y.23= 2001

  y.24= 2002

  y.25= 2003

  y.26= 2004

  y.27= 2005

  y.28= 2006

  y.29= 2007

  y.30= 2008

  y.31= 2009

  y.32= 2010

  y.33= 2011

  y.34= 2012

  y.35= 2013

*

   n=1

   while (n <= 35)

*

   'set lat 0' 

   'set y 1'

   'set z 1'

*

   'mi=('y.n'- 1979)*12'

    'dmi'

   tid=subwrd(result,4)

   ti=6+tid

   tf=8+tid 

*

    'auw'y.n'=aave(ave(u,t='ti',t='tf'),lon=40,lon=55,lev=850,lev=1000)'

    'avw'y.n'=aave(ave(v,t='ti',t='tf'),lon=40,lon=55,lev=850,lev=1000)'

*

     'anomuw'y.n'=auw'y.n'- au'

     'anomvw'y.n'=avw'y.n'- av'

      n = n + 1

*

  endwhile

*

********************************************************************************************************************************

* 

 'reset'

*

  'setgxout fwrite'

  'setfwrite D:/atiqah/desktop/Latest/wind/35.8sj.dat'

*

  n=1

  while (n <= 35)

     'c'

     'set lat 0'

     'set y 1'

     'set z 1'

     'd anomuw'y.n

     'd anomvw'y.n

     n=n+1

*

  endwhile

*

 'disable fwrite'

*


 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20161006/528f3ba1/attachment-0001.html