a   = addfile("test.siconc_SIday_BCC-CSM2-MR.nc","r")
  startDay= 0
  endDay  = 364
  numDays = endDay-startDay+1

  latS    = 50.0

  sic = a->siconc(startDay:endDay,:,:) 
  printVarSummary(sic)
  printMinMax(sic,0)
  print("---")

;;ymdh = cd_calendar(sic&time, -3)
;;print(ymdh)
  ymd  = cd_calendar(sic&time, -2)
  print(ymd)
  
  
  wks = gsn_open_wks("png","siconc_polar")  

  res            = True                          ; plot mods desired
  res@gsnPolar   = "NH"                          ; specify the hemisphere
  res@cnFillOn   = True
  res@cnFillPalette = "WhiteBlueGreenYellowRed"

  res@cnLevelSelectionMode = "ManualLevels"     ; set manual contour levels
  res@cnMinLevelValF       =  5.0               ; set min contour level
  res@cnMaxLevelValF       = 95.0               ; set max contour level
  res@cnLevelSpacingF      =  5.0               ; set contour spacing

  res@cnLevelSelectionMode = "ExplicitLevels"              
  res@cnLevels             = (/5,10,20,30,40,50,60,70,80,90,92,94,96,98,99/) ; %

  res@mpMinLatF  = latS                         ; minimum lat to plot
  res@mpFillOn   = True    ; True->Light Gray;  False-> Background color

  nt   = 364
  res@tiMainString = ymd(nt) 
  plot = gsn_csm_contour_map_polar(wks,sic(nt,:,:),res)    ; create the plot