;;begin      ; not needed for MAIN script

  varname="pr_2a23_hfreez2A23"                                
 ;f   = addfile("/media/stefane/Seagate Expansion Drive/TRMM-PF/199801/1Z09.19980102.00560.7.HDF", "r")
  f   = addfile("./1Z09.19980102.00560.7.HDF", "r")

      ; type short but no unpacking parameter [eg: scale, offset]. Really, no need to unpack
  var = short2flt_hdf(f->$varname$)     ; (:,:) produto 2A23 altura do nível de congelamento em metros
  var@_FillValue = -9999.0
  var@units = "???"
  var!0     = "north_south"
  var!1     = "east_weat"
  printVarSummary(var)
  printMinMax(var,0)
  print("-----")

  lat_dbz   = (f->pr_lat)*0.01            ; (:,:)
  lat_dbz@long_name = "Latitude"
  lat_dbz@units     = "degrees_north"
  lat_dbz!0 = var!0        
  lat_dbz!1 = var!1        
  printVarSummary(lat_dbz)
  printMinMax(lat_dbz,0)
  print("-----")

  lon_dbz   = (f->pr_lon)*0.01            ; (:,:)
  lon_dbz@long_name = "Longitude"
  lon_dbz@units     = "degrees_east"
  lon_dbz!0 = var!0        
  lon_dbz!1 = var!1        
  printVarSummary(lon_dbz)
  printMinMax(lon_dbz,0)
  print("-----")

;---------------------------------------------------------------------------------------------------------
;Subscript para selecionar dados apenas dentro de área de interesse
;
;*****Caixa na Amazõnia coordenadas*************

  latS  = -11.                  
  latN  =  2.
  lonW  = -72.
  lonE  = -45.
  ji    = region_ind(lat_dbz, lon_dbz, latS, latN, lonW, lonE)

  jStrt = ji(0)      ; lat start
  jLast = ji(1)      ; lat last
  iStrt = ji(2)      ; lon start
  iLast = ji(3)      ; lon last

  lat2d = lat_dbz(jStrt:jLast,iStrt:iLast)
  lat2d@long_name = "Latitude Subset"
  lat2d@units="degrees_north"
  printVarSummary(lat2d)
  printMinMax(lat2d, True)
  print("-----")

  lon2d = lon_dbz(jStrt:jLast,iStrt:iLast)
  lon2d@long_name = "Longitude Subset"
  lon2d@units="degrees_east"
  printVarSummary(lon2d)
  printMinMax(lon2d, True)
  print("-----")

  sub = var(jStrt:jLast,iStrt:iLast)
  sub@long_name = "subset: "+sub@long_name 
  printVarSummary(sub)
  printMinMax(sub,0)
  print("-----")

;;sub@lat2d = lat2d                     ; @lat2d only needed for plotting
;;sub@lon2d = lon2d                     ; @lon2d  "     "     "     "

;;end   ; not needed