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

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

title = "Sea level pressure trend" 
figfn = "PSL_trend_ensmean"
f = addfile("PSL_trend_ensmean.nc","r")
var = f->PSL
var = var*120 ;; monthly trend to decadal
;var@lat2d = f->TLAT ;; for cice output, T grid
;var@lon2d = f->TLON

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

res@cnLevelSelectionMode = "ExplicitLevels"
res@cnLevels = fspan(- 50., 50.,11)


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