cg/inputs/plotchimney.m

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

1	cnh	1.1	tstart=0;
2			tint=60;
3			tend=480;
4			nt=size([tstart:tint:tend]);
5
6			figure
7
8			for i=tstart:tint:tend
9
10			subplot(4,3,(i+120)/60);
11			fnam=sprintf('T.%10.10d',i);
12
13			cd /DATA/d31/helen/PRM/MITgcm-r1p8/verification/chimney_coarse/run
14			tc=rdmds(fnam);
15			cd /DATA/d31/helen/PRM/MITgcm-r1p8/verification/nplumesIII/1010
16			tp=rdmds(fnam)
17
18			imagesc(squeeze(tc(10,:,:))');caxis([min(min(tp)) 20]);
19			axis square;colorbar('horiz');
20			set(gca,'XTickLabels',[]);set(gca,'YTickLabels',[])
21			hours=i/60;
22			title([num2str(hours) ' hours']);
23			end
24