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

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

begin
    dfn = "ens_Z200_trend.nc"
    vn  = "Z200"
    ofig= "Z200_ens_trend_ANN.ps"

    df = addfile(dfn,"r")
    var = df->$vn$(0,:,:)
    var = var*120 ;; monthly trend to annual trend

    res = True
        res@cnFillOn = True
        res@cnLinesOn = False
        res@cnFillMode = "CellFill"
        res@mpCenterLonF = 210.
        res@cnLevelSelectionMode = "ExplicitLevels"
        res@cnLevels = (/-5.,5.,10.,15.,20./)
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 Z200"
        res@gsnRightString = "ANN"
        
    wks = gsn_open_wks("ps",ofig)
    plot = gsn_csm_contour_map(wks,var,res)
end