;*****************************************************
;
; These files are loaded by default in NCL V6.2.0 and newer
; load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
; load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;*****************************************************
begin
  filename = "archv.0021_016_00_3zu.nc"
  in       = addfile(filename,"r")
  u        = in->u
  lat      = in->Latitude
;*****************************************************
;   Create plot
;*****************************************************
 wks  = gsn_open_wks("png","hycom")                 ; send graphics to PNG file

 res                             = True             ; plot mods desired

 res@cnFillOn                    = True             ; turn on color 
 res@cnFillPalette               = "BlWhRe"         ; set color map
 res@cnLinesOn                   = False

 res@cnLevelSelectionMode        = "ManualLevels"   ; manual contour levels
 res@cnMinLevelValF              = -0.2             ; set min contour level
 res@cnMaxLevelValF              =  0.2             ; set max contour level
 res@cnLevelSpacingF             =  0.05            ; contour interval

 res@trYReverse                  = True             ; reverse y axis
 res@tiMainString                = "Latitude: "+lat({30}) 

 plot  = gsn_csm_contour(wks,u(0,:,{30},:),res)

end