function surf_movie
g=loadgrid('MESH/stw2');
x=dir('CRUISE_Q5.3/STF*icq5')
n=length(x);
e=g.e';
load survey.1.loc.dat;
xx=survey(:,1);
yy=survey(:,2);
yd=survey(:,3);
t1=1:12;
t2=13:26;
t3=27:40;
t4=41:54;
t5=55:69;
t6=70:82;
for k=1:n
    s=read_icq4_nomex(['CRUISE_Q5.3/',x(k).name]);
    yd_now=datenum(s.year,s.month,s.day)+s.curr_seconds/(24*2600)-12.4206/24/2-datenum(2004,1,1);
    
    j=find(yd<yd_now);

    if ~isempty(j)
       u=hslice(s.ZMID,s.UZMID,-1);
       v=hslice(s.ZMID,s.VZMID,-1);
       X=scattered_interpm(g.x,g.y,[u,v],g.e,xx(j),yy(j),0);
       xx(j)=xx(j)+12.4206*3600*X(:,1);
       yy(j)=yy(j)+12.4206*3600*X(:,2);
    end
     
    clf;
    t=s.SALMID(:,21);
    patch(g.x(e),g.y(e),t(e));
    colormap('jet(20)');
    axis equal;shading interp;caxis([29,33]);
    hold on;
%    plot(xx,yy,'k.');
    plot(...
        xx(t1),yy(t1),'k.-',...
        xx(t2),yy(t2),'kx-',...
        xx(t3),yy(t3),'ko-',...
        xx(t4),yy(t4),'kd-',...
        xx(t5),yy(t5),'k^-',...
        xx(t6),yy(t6),'k+-'...
    );
    colorbar 'h';
    str=sprintf('%3i',k);str(find(str==' '))='0';
    str
    set(gca,'xtick',[],'ytick',[]);
    set(gca,'xticklabel',[],'yticklabel',[]);
    title([int2str(s.month),'/',int2str(s.day),'/',int2str(s.year),': ',int2str(s.curr_seconds)]);
    eval(['print -djpeg movie/frame_',str,'.jpg']);
end
