;; Empty NCL template
;; run code in recipes, use it carefully

;;DIAG_NORCPM; RUNTHESECODES: print("    no ncl code run here")

title = "Sea ice fraction trend" 
figfn = "aice_trend_ensmean"
f = addfile("aice_trend_ensmean.nc","r")
var = f->aice
var = var*120 ;; monthly trend to decadal

if (True)then
    gridf = addfile("/cluster/shared/noresm/inputdata/ocn/micom/gx1v6/20101119/grid.nc","r")
    plat = gridf->plat
    plon = gridf->plon
    var@lat2d = plat
    var@lon2d = plon
end if

res = True
    res@gsnPolar = "NH"
    res@cnFillOn = True
    res@mpMinLatF = 60.
    res@gsnLeftString = title
    res@gsnRightString = "% / 10yr"
    cmap = read_colormap_file("temp_19lev")
    cmap = cmap(::-1,:)
    res@cnFillPalette = cmap

res@cnLevelSelectionMode = "ExplicitLevels"
res@cnLevels = fspan(-8.,8.,17)


wks = gsn_open_wks("ps",figfn)
plot = gsn_csm_contour_map_polar(wks,var(0,:,:),res)