b = addfile("temp_seltime.nc","r")      ; source curvilinear grid
  t = rm_single_dims(b->TSOI(:,0,:,:))
  printVarSummary(t)                      ; [y | 218] x [x | 305]
  printMinMax(t,0)
  print("=====")

  lat2d = b->lat
  lon2d = b->lon
  printVarSummary(lat2d)
  printMinMax(lat2d,0)
  print("=====")
  printVarSummary(lon2d)
  printMinMax(lon2d,0)
  print("=====")

; Entire grid 

  opt_sd = True
  opt_sd@PrintStat = True
  stat_tsoil = stat_dispersion(t, opt_sd )

; Specified Region

  latS   = 5                ; region
  latN   = 40
  lonL   = 65
  lonR   = 100

; Extract indices that 'surround' the specified region
; http://www.ncl.ucar.edu/Document/Functions/Contributed/region_ind.shtml

  ji     = region_ind (lat2d,lon2d, latS, latN, lonL, lonR)

  jStrt  = ji(0)      ; lat start
  jLast  = ji(1)      ; lat last
  iStrt  = ji(2)      ; lon start
  iLast  = ji(3)      ; lon last
  t_region = t(jStrt:jLast,iStrt:iLast)          ; explicitly extract for convenience and clarity

  printVarSummary(t_region)
  printMinMax(t_region,0)
  print("==============")

  stat_region = stat_dispersion(t_region, opt_sd)
  print("==============")

  lat2d_region = lat2d(jStrt:jLast,iStrt:iLast)  ; explicitly extract for convenience and clarity
  lon2d_region = lon2d(jStrt:jLast,iStrt:iLast)

  printVarSummary(lat2d_region)
  printMinMax(lat2d_region,0)
  print("=====")
  printVarSummary(lon2d_region)
  printMinMax(lon2d_region,0)
  print("=====")