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

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

procedure plotice(n)
begin
    if(n.eq.1)then
        vn = "aice"
        title = "Sea ice fraction trend dif. TroPac-Atl DA" 
        unit = "%/10yr"
        cnlevels = fspan(-8.,8.,9)
    end if
    if(n.eq.2)then
        vn = "hi"
        title = "Sea ice thickness trend dif. TroPac-Atl DA" 
        unit = "m/10yr"
        cnlevels = fspan(-.5,.5,9)
    end if

    figfn = "dif_"+vn+"_trend_ensmean"
    fn1 = "norcpm1_tropac_analysis_"+vn+"_trend_ensmean.nc"
    fn2 = "norcpm-cmip6_pacemaker_"+vn+"_trend_ensmean.nc"

    f = addfile(fn1,"r")
    var1 = f->$vn$
    var1 = var1*120 ;; monthly trend to decadal

    f = addfile(fn2,"r")
    var2 = f->$vn$
    var2 = var2*120 ;; monthly trend to decadal

    dvar = var1
    dvar = var1-var2

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

    res = True
        res@gsnPolar = "NH"
        res@cnFillOn = True
        res@mpMinLatF = 40.
        res@gsnLeftString = title
        res@gsnRightString = unit
        res@cnFillPalette = "BlueDarkRed18"

        res@cnLevelSelectionMode = "ExplicitLevels"
        res@cnLevels = cnlevels


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

begin
    plotice(1)
    plotice(2)
end