;*************************************************
; NCL Graphics: ice_3.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
;************************************************
  f = addfile("iceh_mavg.0014-02.nc","r")
;************************************************
; read in ice coverage
;************************************************
  hi         = f->hi

  lat2d      = f->TLAT                          ; read in 2D lat/lon
  lon2d      = f->TLON
  hi@lon2d   = lon2d                            ; assign required attributes
  hi@lat2d   = lat2d
;************************************************
; create plot
;************************************************
  wks = gsn_open_wks("ps","ice")                ; open a ncgm file
  gsn_define_colormap(wks,"BlAqGrYeOrReVi200")  ; choose colormap

  res                      = True               ; plot mods desired
  res@cnFillOn             = True               ; turn on color
  res@gsnPolarNH           = True               ; specify the hemisphere
  res@mpMinLatF            = 65                 ; specify min lat

  res@gsnSpreadColors     = True                ; use full colormap
  res@gsnSpreadColorStart = 20                  ; color to start           
  res@gsnSpreadColorEnd   = 193                 ; color to end

  plot = gsn_csm_contour_map_polar(wks,hi(0,:,:),res)

  end