load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin
  fname     = "19980601orcau.nc"
  f         = addfile(fname,"r")          ; Open file
  ssu       = f->sozocrtx(0,:,:)          ; Read first timestep
  ssu@lon2d = f->lon
  ssu@lat2d = f->lat

  wks = gsn_open_wks("x11","orca")

  res                     = True
  res@gsnMaximize         = True               ; Maximize plot in frame

  res@trGridType          = "TriangularMesh"   ; Define the mesh type. This may not be needed.

  res@cnFillOn            = True               ; Turn on color fill
  res@cnLinesOn           = False              ; Turn lines off
  res@cnLineLabelsOn      = False              ; Turn labels off

  res@tiMainString       =  "ORCA Grid - Smooth fill"

  map = gsn_csm_contour_map(wks,ssu,res)

  res@tiMainString       = "ORCA Grid - Raster fill"
  res@cnFillMode         = "RasterFill"

  map = gsn_csm_contour_map(wks,ssu,res)

end