1 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
  2 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
  3 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
  4 
  5 begin
  6 
  7 do i = 0,0  ; 19 
  8 
  9 filename   =  "/data1/share/CSM/mom4/"+ (21100101+50000*i)  + ".ice_day.nc"
 10 varname    =  "HI"
 11 
 12 fi         =  addfile(filename,"r")
 13 data       =  fi->$varname$
 14 lon        =  fi->GEOLON
 15 lat        =  fi->GEOLAT
 16 print(min(lat) + " -- " + max(lat))
 17 print(min(lon) + " -- " + max(lat))
 18 printVarSummary(data)
 19 printVarSummary(lon)
 20 printVarSummary(lat)
 21 
 22 x          =  dim_avg_Wrap(data(yt|:,xt|:,time|:))
 23 
 24 end do
 25 
 26 
 27 ; plot 
 28 
 29   wks = gsn_open_wks("eps","test")
 30   ;gsn_define_colormap(wks,"BlWhRe") 
 31   gsn_define_colormap(wks,"rainbow") 
 32 
 33   res  =  True
 34  ;res@gsnDraw  =  False
 35  ;res@gsnFrame =  False
 36 
 37   
 38   res@gsnMaximize           = True
 39   
 40   res@trGridType          = "TriangularMesh"    ; Define the mesh type.
 41   res@sfXArray              = fi->GEOLON
 42   res@sfYArray              = fi->GEOLAT
 43 
 44   res@cnFillOn              = True
 45   res@cnFillMode            = "RasterFill"
 46   res@gsnSpreadColors       = True
 47   res@gsnSpreadColorStart   = 88
 48   res@gsnSpreadColorEnd     = 8
 49 
 50   res@cnLinesOn             = False
 51   res@cnLineLabelsOn        = False
 52 
 53   res@lbBoxLinesOn          = False
 54   res@lbLabelAutoStride     = True
 55 
 56   res@tiMainString          = "MOM output test " 
 57 
 58 ;symMinMaxPlt(x, 10, False, res)   ; set symmetric plot min/max
 59   res@gsnRightString = ""
 60   res@gsnCenterString = ""
 61   res@cnInfoLabelOn   = False 
 62 ; res@gsnPolar    =  "NH"
 63 ; res@mpMinLatF   =  50
 64 
 65   plot  = gsn_csm_contour_map(wks, x, res)
 66 ; plot  = gsn_csm_contour_map_polar(wks, x ,res)
 67 
 68 
 69 
 70 
 71 end
 72 
 73