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	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