* Script 3
* no 'set lon 90' and 'set lat -5' before x'k'
* 'close '1 at the end of the script
* output is perfect

   'reinit'
*
*
     'sdfopen h:/cdc/monthlymean/olr.day.mean.nc'
*
     j=11872
     k=1
     'xsum=0'
     while (j <= 11902)
       'set t 'j
       'x'k'=aave(olr,lon=100,lon=110,lat=0,lat=5)'
       'xsum=xsum+x'k
       j=j+1
       k=k+1
     endwhile
     'xbar=xsum/('k-1')'
*
*
      'y=x15-xbar'
      'y1=pow(y,2)'
*
      i=1
      while (i <= 31)
         'yd'i'=x'i'-xbar'
         i=i+1
      endwhile
*
      j=1
      'yp=0'
      while (j <= 2)
        'yp=yp+(pow(yd'j',2))'
        j=j+1
      endwhile
*
     'close '1
