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	edhill	1.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