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/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"

begin

fils = systemfunc("ls *.nc")

  f  = addfiles(fils, "r")                               
printVarSummary(f)

;ListSetType(f,"join")                            

;---Read data
  sit      = f[:]->SST
printVarSummary(sit)

sit!0="time"
sit!1="lat"
sit!2="lon"
TIME=sit!0
j=sit!1
i=sit!2   

sit@lat2d = f[:]->J
sit@lon2d = f[:]->I

printVarSummary(sit)
  src_lat  = f[:]->J
  src_lon  = f[:]->I
  src_mask = where(.not.ismissing(sit(0,:,:)),1,0)

;---Remove these because they are not the coordinates we want to use
 ;delete(sit&i)
 ; delete(sit&j)

;---Set up regridding options
  Opt                   = True

  Opt@InterpMethod      = "bilinear"        ; "bilinear","patch","conserve","neareststod"
  Opt@SrcGridLat        = src_lat           ; source grid
  Opt@SrcGridLon        = src_lon
  Opt@SrcMask2D         = src_mask

  Opt@DstGridType       = "2deg"            ; destination grid

  Opt@ForceOverwrite    = True
  Opt@PrintTimings      = True
  Opt@Debug             = True

  if(Opt@InterpMethod.eq."conserve") then
    Opt@SrcGridCornerLat = f[:]->lat_vertices    ; corners are needed for"conserve" 
    Opt@SrcGridCornerLon = f[:]->lon_vertices    ; method
  end if

  sit_regrid = ESMF_regrid(sit,Opt)  ; Do the regridding
  printVarSummary(sit_regrid)

;---Write to file
  newfile = "sit1.nc"
  system("rm -f "+newfile)
  file_create = addfile(newfile, "c")    ;criando arquivo.
  file_create->tos= sit_regrid

end