cd ../uncompressed
fn='calus200_flt.nc'

%open float file, get coordinates of trajectories
% nc=netcdf(fn);
% x=nc{'x'}(2:end,:);  %last index for for the number of floats released
% y=nc{'y'}(2:end,:);
% close(nc);

x_=nc_varget(fn,'x');
x=x_(2:end,:);

y_=nc_varget(fn,'y');
y=y_(2:end,:);

%open his.nc file to get coordinates of computational grid

% nc=netcdf('tmpnc34.nc');
% Xgrd=nc{x_rho}(:,:)
% Ygrd=nc{y_rho}(:,:)
% 
% close(nc)

Xgrd=nc_varget(nc,'x_rho');
Ygrd=nc(y_rho);




cd ../analout

figure
%blank grid
pcolor(Xgrd,Ygrd,Xgr*NaN);
hold on
%axis equal
for i=1:length(x(:,1)),
plot(x(i,:),y(i,:),'.b')
pause(0.2)
end


