if ~exist('uvb')
uvb=load(['../exp_forw_dan/uv_bck_dan.dat']);
end
filelonlat='../../data/data/MODEL_CYCO/uveta_cyco1_2013090112.nc';
lon=ncread(filelonlat,'lon');
lat=ncread(filelonlat,'lat');
[Lon,Lat]=meshgrid(lon,lat);
it=1:2:25
nlon=length(unique(uvb(:,2)));
nlat=length(unique(uvb(:,3)));
ntime=length(unique(uvb(:,1)));
Ub=reshape(uvb(:,end-1),nlat,nlon,ntime);
Vb=reshape(uvb(:,end),nlat,nlon,ntime);

for i=it
Hb=quiver(Lon,Lat,squeeze(Ub(:,:,i)),squeeze(Vb(:,:,i)),'b');
axis([35.5 36 34.5 35.5])
title(['it=' int2str(i)]);
pause(0.2)
end