load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl"
begin
  case     = "DISPERSION_CURVE:"  
  diro     = "/media/ANN/"        ; output dir: location of plots 

  latN     = 20
  latS     = -15       

  nDayWin  = 96      ; Wheeler-Kiladis [WK] temporal window length (days)
  nDaySkip = 30      ; time (days) between temporal windows [segments]
                     ; negative means there will be overlapping 
                     ; temporal segments

  opt      = True
  opt@debug= True
  opt@pltType     = "pdf"
  opt@cnLinesOn   = False   ; turn off contour lines
                     ; specify the plot intervals

opt@FillValue= -999

  spd      = 4       ; samples per day (6 hour sampling rate)
  var      = "OMEGA500"  ; user specified 

  diri     = "/media/ANN/"

  fili     = "may_eta.nc"    
 

  f        = addfile (diri+fili, "r") 
  eta    = f->ETAA(:,0,{latS:latN},:)     ; (time,lat,lon)
nMsg = num(ismissing(eta))
                     ; calculation and create plots
  wkSpaceTime (eta, diro, case, var,latN, spd, nDayWin, nDaySkip, opt  )
end