load /nas/edl/OFES_UV_KOE/matlib/KOE_GRD.mat

%% Load velocity data
disp('Loading velocity data ...');

tstart=datenum(YEAR1, MONTH1, DAY1);
tend=datenum(YEAR2, MONTH2, DAY2);

istr=find(ctl.time < tstart);
iend=find(ctl.time > tend);
in=istr:iend;

U=rnt_loadvar(ctl,in, 'u');
V=rnt_loadvar(ctl,in, 'v');
UV_TIME=ctl.time(in);


U=U/100; % convert to m/s from cm/s
V=V/100;
U(isnan(U))=0;
V(isnan(V))=0;