load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "./WRFUserARW.ncl"

begin

; Open up the workstation so we can make a plot
wks  = gsn_open_wks ("X11","t2_plot")

; Open up the T2 max/min file that we just made
a = addfile("t2_maxmin_daily_jja_1990_rtty_1990.nc","r")
; Read in t2max, lat and lon
t2max = a->T2MAX
printVarSummary(t2max)

; Options for our data contour
opts = True
opts@cnFillOn = True  ; Setting the contour color fill
pltres = True
mpres = True
mpres@mpGeophysicalLineColor = "black"
mpres@mpGridSpacingF = 15

; wrf_contour draw contours for the t2max zoomed data
contour = wrf_contour(a,wks,t2,opts)

; Here we are plotting the t2max contoured data over a map using the function wrf_map_overlays
plot = wrf_map_overlays(a,wks,contour,pltres,mpres)

; End of program
end