% create ncgeodataset 
rtv=ncgeodataset('http://hfrnet.ucsd.edu/thredds/dodsC/HFRNet/USWC/1km/hourly/RTV');

% use geovariable, this grabs all available data
U=rtv.geovariable('u');
V=rtv.geovariable('v');

% subset lat and lon
%s.lat=[33.47 33.50];  %Orig Location
%s.lon=[-118.56 -118.52];  %Orig Location

s.lat=[33.46 33.49];  %Alt Location
s.lon=[-118.53 -118.49];  %Alt Location


% subset time by time
tstart=datenum(2011,1,1,0,0,0);
tstart = datenum(2011,1,1,0,0,0);
tend=now;
s.time=[tstart tend];

% use geosubset
Us=U.geosubset(s);
Vs=V.geosubset(s);

% clear full records
clear U V;

% assign variables
times=Us.grid.time;
lat=Us.grid.lat;
lon=Us.grid.lon;
Udata=Us.data;
Vdata=Vs.data;