* To obtain climatological pentad mean 850 hPa winds
* by Ooi See Hai on 12 September 2006
*
 'reinit'
*
  'set gxout fwrite'
  'set fwrite c:\onset\data\o_pentad.dat'
*
*
    'sdfopen g:\cdc\uvw\uwnd.1986.nc'
    'sdfopen g:\cdc\uvw\vwnd.1986.nc'
*
* Note: Initialization can only be done after opening files
*
   i=1
   while (i <= 73)
     'up'i'=0'
     'vp'i'=0'
     i=i+1
   endwhile
*
   ny=1986
   m=1
*
   while (m <= 40)
*
    'set lev 850'
*
    i=1
    ti=1
    tf=17
*
* The extra day in leap year is not taken into account
*
    while (i <= 73)
*
      'up'i'=up'i'+ave(uwnd.'m',t='ti',t='tf',4)'
      'vp'i'=vp'i'+ave(vwnd.'m+1',t='ti',t='tf',4)'
      i=i+1
      ti=ti+20
      tf=tf+20
    endwhile
*
*
    ny=ny+1
    m=m+2
*
    if (m <= 40)
     'sdfopen g:\cdc\uvw\uwnd.'ny'.nc'
     'sdfopen g:\cdc\uvw\vwnd.'ny'.nc'
     'set dfile 'm
    endif
*
   endwhile
*
   i=1
   while (i <= 73)
     'set t 'i
     'upm'i'=up'i'/'m-2
     'vpm'i'=vp'i'/'m-2
     'd ave(upm'i',lon=97,lon=108)'
     'd ave(vpm'i',lon=97,lon=108)'
     i=i+1
   endwhile
*
  'disable fwrite'
*
*