tstart=0;
tint=60;
tend=480;
nt=size([tstart:tint:tend]);
 
figure

for i=tstart:tint:tend

subplot(4,3,(i+120)/60);
fnam=sprintf('T.%10.10d',i);

cd /DATA/d31/helen/PRM/MITgcm-r1p8/verification/chimney_coarse/run
tc=rdmds(fnam);
cd /DATA/d31/helen/PRM/MITgcm-r1p8/verification/nplumesIII/1010
tp=rdmds(fnam)

imagesc(squeeze(tc(10,:,:))');caxis([min(min(tp)) 20]);
axis square;colorbar('horiz');
set(gca,'XTickLabels',[]);set(gca,'YTickLabels',[])
hours=i/60;
title([num2str(hours) ' hours']);
end