begin
a = addfile("/scale_wlg_persistent/filesets/project/uoa02450/2017_Test/wrfout_d03_2017-01-01_00:00:00_Time.nc","r")
type = "pdf"
wks = gsn_open_wks(type,"Average 55 m Wind Speed - Jan_TEST 2017")
res = True
res@mpProjection  = "Mercator"    ; The default
res@mpDataBaseVersion = "HighRes"

times = wrf_user_getvar(a,"times",-1) 
ntimes = dimsizes(times)       

u = wrf_user_getvar(a,"ua",-1) 
v = wrf_user_getvar(a,"va",-1)   
z = wrf_user_getvar(a,"z",-1)
ter = wrf_user_getvar(a,"ter",-1) 

nheight = conform(z,ter,(/0,2,3/))
height = z - nheight

u_plane = wrf_user_intrp3d(u,height,"h",55.,0.,False)
v_plane = wrf_user_intrp3d(v,height,"h",55.,0.,False)

spd = (u_plane*u_plane + v_plane*v_plane)^(0.5)
spd@description = "Wind Speed"
spd@units = "m s-1" 
     
spd_avg = dim_avg_n(spd(0:31,:,:,:),0)
spd_avg@description = "Average Wind Speed"
spd_avg@units = "m s-1" 
              
res = True
opts =  res
opts@cnFillOn = True  
opts@cnLinesOn = True
opts@ContourParameters = (/ 0., 12., 0.5/)
opts@gsnSpreadColorEnd = -3  
contour_spd_avg = wrf_contour(a,wks,spd_avg(0,:,:),opts)
delete(opts)

plot = wrf_map_overlays(a,wks,(/contour_spd_avg/),res,res)
      
end