addpath([pwd '/core/misc']);
%fs=load('L:\array\tcf\rescaled\z001_t001.mat');
fs=load('/Volumes/home/test/4C_interpolant_spill2/tracks/z001_t001.mat').tracks;

%%
figure(1);
clf;

% hold on;
 
gca;
colormap jet
scale=0.01;
u = [fs.u];
x = [fs.x];
u_norm = sqrt(sum(u.^2,1));
% quiverC3D(x(1,:), x(2,:), x(3,:), u(1,:)*scale, u(2,:), u(3,:), 'scale', 0);
% scatter3(x(1,:),x(2,:),x(3,:),ones(1,length(x)),u_norm);
quiver3(x(1,:), x(2,:), x(3,:), u(1,:), u(2,:), u(3,:), 'k-');
axis equal tight xy