;;fils = systemfunc("ls [1-2][8-9-0-1]*/*nc")
  fils = systemfunc("ls AVH*nc")
  print(fils)

   f   = addfiles (fils, "r")
  ;ListSetType (f, "cat")
   prc   = short2flt(f[:]->NDVI)

   printVarSummary(prc)
   time = f[:]->time
   yyyymm = cd_calendar(time, -1)
   ntim = dimsizes(time)

  months = (/"January", "February", "March", "April" \
            ,"May", "June", "July", "August"         \
            ,"September", "October", "November"      \
            ,"December" /)
;************************************************
; Compute the climatology using a function in contributed.ncl
;************************************************
  opt = True
  opt@nval_crit = 12   ; require at least 12 values for the "avg" is calculated.

 xM = calculate_monthly_values (prc, "avg", 0,opt)
 printVarSummary(xM)