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

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

begin
    season = getenv("SEASON")
    case1 = getenv("CASE1")
    case2 = getenv("CASE2")
    difstr = getenv("DIFSTR")
    difstrfn = getenv("DIFSTRFN")
    dfn1 = case1+"_ens_TREFHT_"+season+"_trend.nc"
    dfn2 = case2+"_ens_TREFHT_"+season+"_trend.nc"
    vn  = "TREFHT"
    ofig= "dif_"+difstrfn+"_TREFHT_ens_trend_"+season+".ps"

    df1 = addfile(dfn1,"r")
    var1 = df1->$vn$(0,:,:)
    var1 = var1*10 ;; monthly trend to annual trend

    df2 = addfile(dfn2,"r")
    var2 = df2->$vn$(0,:,:)
    var2 = var2*10 ;; monthly trend to annual trend

    dvar = var1
    dvar = var1-var2
    res = True
        res@cnFillOn = True
        res@cnLinesOn = False
        res@cnFillMode = "CellFill"
        res@mpCenterLonF = 210.
        res@cnFillPalette = "BlueDarkRed18"
        res@cnLevelSelectionMode = "ExplicitLevels"
        res@cnLevels = fspan(-.8,.8, 17)
        ;res@cnFillColors = (/(/225,239,223,255/), (/255,255,255,255/), (/248,246,226,255/), (/244,218,140,255/), (/241,192,45,255/), (/168,129,32,255/)/)
        ;res@cnFillColors := res@cnFillColors / 255.
        
        res@gsnLeftString = "Decadal trend t2m "+difstr
        res@gsnRightString = season
        
    wks = gsn_open_wks("ps",ofig)
    plot = gsn_csm_contour_map(wks,dvar,res)
end