diri="./"
  sfile = "chathu.nc" 
  f = addfile(diri+sfile,"r")

  tauvo_ssp           = f->tauvo
  printVarSummary(tauvo_ssp)
  printMinMax(tauvo_ssp,0)
  print("---------------")

  nZero = num(tauvo_ssp.eq.0.0)  ; count the number of exact 0.0 values
  print("nZero="+nZero)
  print("---------------")
  
  LAT2D = f->lat
  LON2D = f->lon
                               ; explore coordinates
;;printVarSummary(LAT2D)
  printMinMax(LAT2D,0)
  print("---------------")
;;printVarSummary(LON2D)
  printMinMax(LON2D,0)
  print("---------------")
  
; print all longitudes at arbitrarily specified latitude indices
  
  NT = 120      ; arbitrary time
  JJ =  75  
  print("JJ="+JJ+"  "+LAT2D(JJ,:)+"   "+LON2D(JJ,:)+"   "+tauvo_ssp(NT,JJ,:))
  print("-----")

  JJ = 125
  print("JJ="+JJ+"  "+LAT2D(JJ,:)+"   "+LON2D(JJ,:)+"   "+tauvo_ssp(NT,JJ,:))
  print("-----")