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

begin
  orca = addfile("ctorca.nc","r")
  sst  = orca->sosstsst(0,:147,1:)      ; Read data
  lon  = orca->nav_lon
  lat  = orca->nav_lat

;
; Open PS workstation and change color map.
;
  wks = gsn_open_wks("ps","orca")
  gsn_define_colormap(wks,"gui_default")

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

  res@trGridType          = "TriangularMesh"    ; Define the mesh type.

  res@sfXArray            = lon
  res@sfYArray            = lat
  res@sfMissingValueV     = sst@missing_value

; Define colors for filled contours

  res@cnFillOn            = True
  res@cnFillColors        = (/ 2, 3, 4, 5, 7, 8, 9, 11, 12, 13, 15, 16, 17, \
                               19, 20, 21, 23/)

  res@cnLinesOn           = False         ; Turn lines off
  res@cnLineLabelsOn      = False         ; Turn labels off

  res@mpProjection         = "Orthographic"
  res@mpCenterLatF         =  50
  res@mpGridLineColor      =  -1       ; Turn off grid lat/lon lines
  res@mpGridAndLimbOn      =  True

  res@cnFillMode           = "AreaFill"       ; The default.
  res@tiMainString         = "Orca Grid - Default fill mode"

  map = gsn_csm_contour_map(wks,sst,res)

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

  map = gsn_csm_contour_map(wks,sst,res)

end