;*************************************************
; NCL Graphics: polar_1.ncl
;************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" 
;************************************************
begin
;************************************************
; read in netCDF file
;************************************************
  a = addfile("uv300.nc","r")
;************************************************
; read in zonal winds
;************************************************
  u = a->U(1,:,:)                                ; read July zonal winds
;************************************************
; create plot
;************************************************
  wks = gsn_open_wks("ps" ,"polar")               ; open a ps file
  res                      = True                 ; plot mods desired
  res@gsnPolar             = "NH"                 ; specify the hemisphere
  res@gsnPolarLabelSpacing = 90                   ; how frequently to label

  plot = gsn_csm_contour_map_polar(wks,u,res)     ; create the plot

end