cd ../uncompressed
fn='calus200_fltW.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)

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