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

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

   printVarSummary(x)
;************************************************
; 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 (x, "avg", 0,opt)
 printVarSummary(xM)