;---------------------------------------------------------------------------------------------------
begin
  a = addfile("/home/ew/NCL/bin/indices/level850/sst/test.nc","r")

  s  = a->sst

 ;printMinMax(s, True)
 ;print(s)
 ;printVarSummary(s)


;---Graphics
  wks = gsn_open_wks("png","local")    ; Open a workstation


  ;cmap = read_colormap_file("MPL_RdBu")   ; read the color map

;--------------
  res              = True

  res@gsnDraw      = True                   ; draw the plot
  res@gsnFrame     = True                   ; advance the frame 
  res@gsnAddCyclic = False
  res@gsnMaximize  = True

  ;res@cnFillPalette = cmap(::-1,:)          ; draw the reversed color map
  ;res@cnFillPalette = cmap(10:100,:)        ; subset color map

  ;res@cnFillPalette        = "SVG_es_landscape_79"
  ;res@cnFillPalette         = "MPL_afmhot"
  ;res@cnFillPalette        = "MPL_YlOrRd"
  ;res@cnFillPalette        = "BlAqGrYeOrRe"
  res@cnFillPalette        = "MPL_Reds"
  ;res@cnFillPalette        = "WhiteBlueGreenYellowRed"
  ;res@cnFillOpacityF       = 0.8

  res@mpLimitMode         = "LatLon"
  res@mpMinLatF           = -5
  res@mpMinLonF           = -135
  res@mpMaxLatF           = 5
  res@mpMaxLonF           = -120  

  res@cnLineLabelsOn       = False             ; do use line labels
  res@cnFillOn             = True              ; color fill
  res@cnLinesOn            = False             ; do draw contour lines                

  ;res@cnInfoLabelOn        = True

  res@lbOrientation         =   "vertical"          ; vertical label bar

  ;res@cnLevelSelectionMode = "ExplicitLevels"
  ;res@cnLevels = ispan(75800,78800,200)

  ;res@cnLevels = (/0.44,0.48,0.52,0.56,0.6,0.64,0.68,0.72,0.76,0.8,0.84,0.88/)


  res@gsnCenterString     = " "
  res@gsnLeftString       = "sst "
  res@gsnRightString      = " deg.K"   

  ;res@txFontHeightF   = 0.025

  ;res@gsnRightStringFontHeightF   = 0.025

  res@gsnCenterStringOrthogonalPosF = 0.018
  res@gsnRightStringOrthogonalPosF  = 0.015
  res@gsnLeftStringOrthogonalPosF   = 0.020

  ;res@gsnRightStringParallelPosF   = 0.020

  res@tmXTOn               = False       ; bottom off
  res@tmYROn               = False       ; right off

  res@tiMainString      = "The East Pacific Trade Wind Index  ~C~       ~Z90~ (ERA5 Timeseries of 2008)"
  res@tiMainFont        = "times-roman"
  res@tiMainFontHeightF =  0.018        ;-- decrease the size of the title
  res@tiMainOffsetYF    =  0.05

  ;res@lbTitlePosition          = "Bottom"     ; put it below the plot
  ;res@lbTitleString            = "Potential Energy (m^2/s^2)" 
  ;res@lbTitleFontHeightF       =  0.012 
  ;res@lbTitleFont              = "times-roman"
  ;res@lbLabelFontHeightF       = 0.013
  ;res@lbLabelAngleF            =    65                 ; angle labels
  ;res@lbLabelAlignment = "InteriorEdges"
  ;res@pmLabelBarOrthogonalPosF = 0.10           ; move whole thing down

;---create the plot
  plot = gsn_csm_contour_map(wks, s(0,:,:), res)

end