function mir3d_movie
g=loadgrid('mesh/ga_la');
e=g.e';
xb=g.x(g.bnd');
yb=g.y(g.bnd');
load obs/sst.mir3d
t=unique(sst(:,1));
clf;
ax=[-66207      357918           0      321343];
t=t([1,2,4,5,12,13,15,16,17]);
for k=1:length(t)
    subplot(3,3,k)
    j=find(sst(:,1)==t(k));
    plot(xb,yb,'k');hold on;colorscatter(sst(j,2),sst(j,3),sst(j,5),'.',[6,26]);axis(ax)
    set(gca,'xtick',[],'ytick',[]);
    title([num2str(t(k))]);
end
colorbar 'v';
