timour_matlab/mscripts/plottt.m

clear 
clear path 
whos
global Nx Ny Nz
global lat long dz dm mdep
global delt_su su_its t_su delt
global descriptor this_path
global f deltaf Q beta r_expt r_heat H
global time rots it
global g Cp rho_bar alpha
global u v t w 
global iterations

for iii=1:4,
getuvwt;

%i=find(t<20.1);
%t(i)=NaN;
t(1,:,:)=NaN;
t(Nx,:,:)=NaN;
t(:,:,3)=30;
t(:,:,Nz)=NaN;

        colormap gray
        cmap=flipud(colormap);
        colormap(cmap);


cmin=20;
cmax=30;
V=[cmin cmax];
% figure
caxis('manual')

subplot(4,2,2*iii-1)
%        title='in-situ temperature at z=5';
        imagesc(lat,long,squeeze(t(:,:,5))');shading flat;axis image;caxis(V);colorbar
        set(gca,'ydir','norm','XTickLabel','|||','YTickLabel',' |||')
        xlabel('x','FontSize',15)
        ylabel('y','FontSize',15,'Rotation',[0])

        if iii==1
        text(-30,110,'a)','FontSize',15)
        end

        if iii==2
        text(-30,110,'b)','FontSize',15)
        end

        if iii==3
        text(-30,110,'c)','FontSize',15)
        end

        if iii==4
        text(-30,110,'d)','FontSize',15)
        end

%        text(0,-50,descriptor);
%        text(0,140,title);
%        text(0,-30,'timestep');text(80,-30,num2str(it))
subplot(4,2,2*iii)
        [c,h]=contourf(flipud(squeeze(t(:,Ny/2,3:30))'),7);shading flat;caxis(V);box on;
%       set(h(:),'LineWidth',[0.001])
%        imagesc(lat,long,squeeze(t(:,Ny/2,3:30))');shading flat;axis image;caxis(V);
        set(gca,'DataAspectRatio',[4,1,1],'XTickLabel','|||','YTickLabel',' |||')
        xlabel('x','FontSize',15)
        ylabel('z','FontSize',15,'Rotation',[0])
%        text(0,-2,'in-situ east-west temperature at mid-tank')

end
1	clear
2	clear path
3	whos
4	global Nx Ny Nz
5	global lat long dz dm mdep
6	global delt_su su_its t_su delt
7	global descriptor this_path
8	global f deltaf Q beta r_expt r_heat H
9	global time rots it
10	global g Cp rho_bar alpha
11	global u v t w
12	global iterations
13
14	for iii=1:4,
15	getuvwt;
16
17	%i=find(t<20.1);
18	%t(i)=NaN;
19	t(1,:,:)=NaN;
20	t(Nx,:,:)=NaN;
21	t(:,:,3)=30;
22	t(:,:,Nz)=NaN;
23
24	colormap gray
25	cmap=flipud(colormap);
26	colormap(cmap);
27
28
29	cmin=20;
30	cmax=30;
31	V=[cmin cmax];
32	% figure
33	caxis('manual')
34
35	subplot(4,2,2*iii-1)
36	% title='in-situ temperature at z=5';
37	imagesc(lat,long,squeeze(t(:,:,5))');shading flat;axis image;caxis(V);colorbar
38	set(gca,'ydir','norm','XTickLabel','\|\|\|','YTickLabel',' \|\|\|')
39	xlabel('x','FontSize',15)
40	ylabel('y','FontSize',15,'Rotation',[0])
41
42	if iii==1
43	text(-30,110,'a)','FontSize',15)
44	end
45
46	if iii==2
47	text(-30,110,'b)','FontSize',15)
48	end
49
50	if iii==3
51	text(-30,110,'c)','FontSize',15)
52	end
53
54	if iii==4
55	text(-30,110,'d)','FontSize',15)
56	end
57
58	% text(0,-50,descriptor);
59	% text(0,140,title);
60	% text(0,-30,'timestep');text(80,-30,num2str(it))
61	subplot(4,2,2*iii)
62	[c,h]=contourf(flipud(squeeze(t(:,Ny/2,3:30))'),7);shading flat;caxis(V);box on;
63	% set(h(:),'LineWidth',[0.001])
64	% imagesc(lat,long,squeeze(t(:,Ny/2,3:30))');shading flat;axis image;caxis(V);
65	set(gca,'DataAspectRatio',[4,1,1],'XTickLabel','\|\|\|','YTickLabel',' \|\|\|')
66	xlabel('x','FontSize',15)
67	ylabel('z','FontSize',15,'Rotation',[0])
68	% text(0,-2,'in-situ east-west temperature at mid-tank')
69
70	end