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

 begin

; Make a list of all files we are interested in
  DATADir = "./"
  FILES = systemfunc (" ls -1 " + DATADir + "wrfout_d04* ")

  f     = addfiles(FILES+".nc","r")
  times = wrf_user_list_times(f)
  ntimes = dimsizes(times)
lat = f[0]->XLAT(0,:,:)
  lon = f[0]->XLONG(0,:,:)

  do it = 0,ntimes-1,1          ; TIME LOOP
;    print("Working on time: " + it + ") = " + times(it))
    u = wrf_user_getvar(f,"ua",it)      ; u averaged to mass points
    v = wrf_user_getvar(f,"va",it)      ; v averaged to mass points
;    printVarSummary(u)
;   printMinMax(u,0)
;   printMinMax(v,0)

    loc = wrf_user_ll_to_ij(f,36.85693,2.53115,True)
    mx = loc(0) - 1
    ny = loc(1) - 1
    lat1  = lat(ny,mx)    ; the is the nearest grid point
    lon1  = lon(ny,mx)
  ;  print("lat1/lon1 = " + lat1 + "/" + lon1)
     res               = True
    res@tiMainString  = "WRF: All Times: grid point ["+sprintf("%5.2f",lat1)+" , "+sprintf("%5.2f",lon1) +"]"
    wspd1= ndtooned( sqrt(u(2,ny,mx)^2 + v(2,ny,mx)^2) )
    print("Working on time: " + times(it) + " wind_speed="+wspd1)
   end do
end