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

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

begin
    dfn1 = "norcpm1_tropac_analysis_ens_TREFHT_trend.nc"
    dfn2 = "norcpm-cmip6_pacemaker_ens_TREFHT_trend.nc"
    vn  = "TREFHT"
    ofig= "dif_pac-atl_TREFHT_ens_trend_ANN.ps"

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

    df2 = addfile(dfn2,"r")
    var2 = df2->$vn$(0,:,:)
    var2 = var2*120 ;; 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, 9)
        ;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 TroPac - Atl DA"
        res@gsnRightString = "ANN"
        
    wks = gsn_open_wks("ps",ofig)
    plot = gsn_csm_contour_map(wks,dvar,res)
end