
REstar=0;  Bo=1000;

figure(1),clf
figure(2),clf
for t=1003.5:0.5:1021.0,
    %for t=1011.0:1011.0, %0.5:1020.5,
  index=return_basilisk_index(t);

[Xs,Ys,Xraw,Yraw] = read_processed_surface_extract(9,REstar,Bo,t);

ns = length(Xs);
C = [1:ns]'/ns;
%nn = length(Xn);
%Cn = [1:nn]'/nn;
    
figure(1)
hold on
     scatter(Xs,Ys+(t-1002)/5,4,C,'d')
    drawnow
    pause(0.02)
    %    plot(Spc.Location(:,1),Spc.Location(:,2),'kd')

    figure(2)
    hold on
    %for t=1011.0:1011.0, %0.5:1020.5,

nraw = length(Xraw);
Cr = [1:nraw]'/nraw;
%nn = length(Xn);
%Cn = [1:nn]'/nn;
    
figure(2)
hold on
     scatter(Xraw,Yraw+(t-1002)/5,4,Cr,'d')
    drawnow
    pause(0.02)
     
end;
    
    
  %[F,x,y,t]  = read_binary_file(9,REstar,Bo,index,'f');
%[V,x,y,t]  = read_binary_file(9,REstar,Bo,index,'z');
%[P,x,y,t]  = read_binary_file(9,REstar,Bo,index,'p');
