timour_matlab/mscripts/DIAG1a.m

x=1:Nx-2;
y=1:Ny-2;

subplot(2,1,1)
imagesc(abs(squeeze(sumD1(4:Nx-2,3:Ny-1,1)))')
%contourf(x,y,squeeze(sumD1(1:Nx-2,1:Ny-2,1))',8)
 set(gca,'DataAspectRatio',[2,2,2]) 
        set(gca,'ydir','norm')
xlabel('X ','FontSize',17)
  ylabel('Y ','FontSize',17)  
 set(gca,'XtickLabel','||')
 set(gca,'YtickLabel','||')
text(Nx2/2-25,Ny2+5,'\int Q D_p dz','Fontsize',17);
text(-25,Ny2-5,'(a)','Fontsize',17);

subplot(2,1,2)
imagesc(abs(squeeze(sumM1(4:Nx-2,3:Ny-1,1)))')
%contourf(x,y,-squeeze(sumM1(1:Nx-2,1:Ny-2,1))',8)
 set(gca,'DataAspectRatio',[2,2,2]) 
        set(gca,'ydir','norm')
xlabel('X ','FontSize',17)
ylabel('Y ','FontSize',17)  
set(gca,'XtickLabel','|||')
set(gca,'YtickLabel','|||')
text(Nx2/2-25,Ny2+5,'\int Q M_p dz','Fontsize',17);
text(-25,Ny2-5,'(b)','Fontsize',17);

colormap gray   
cmap=colormap;  
cmap=flipud(cmap);
%cmap=sqrt(cmap);
cmap=0.5*(1-cos(3.1415*cmap));
colormap(cmap) 

%figure

%subplot(2,1,1)
%contourf(x,y,squeeze(sumD1(1:Nx-2,1:Ny-2,1))',8)
% set(gca,'DataAspectRatio',[2,2,2])
%        set(gca,'ydir','norm')
%xlabel('X ','FontSize',17)
%  ylabel('Y ','FontSize',17)
% set(gca,'XtickLabel','||')
% set(gca,'YtickLabel','||')
%text(Nx2/2-25,Ny2+5,'\int Q D_p dz','Fontsize',17);
%text(-25,Ny2-5,'(a)','Fontsize',17);

%subplot(2,1,2)
%contourf(x,y,-squeeze(sumM1(1:Nx-2,1:Ny-2,1))',8)
% set(gca,'DataAspectRatio',[2,2,2])
%        set(gca,'ydir','norm')
%xlabel('X ','FontSize',17)
%ylabel('Y ','FontSize',17)
%set(gca,'XtickLabel','|||')
%set(gca,'YtickLabel','|||')
%text(Nx2/2-25,Ny2+5,'\int Q M_p dz','Fontsize',17);
%text(-25,Ny2-5,'(b)','Fontsize',17);

% cmap=flipud(cmap);
%colormap(cmap) 
1	edhill	1.1	x=1:Nx-2;
2			y=1:Ny-2;
3
4			subplot(2,1,1)
5			imagesc(abs(squeeze(sumD1(4:Nx-2,3:Ny-1,1)))')
6			%contourf(x,y,squeeze(sumD1(1:Nx-2,1:Ny-2,1))',8)
7			set(gca,'DataAspectRatio',[2,2,2])
8			set(gca,'ydir','norm')
9			xlabel('X ','FontSize',17)
10			ylabel('Y ','FontSize',17)
11			set(gca,'XtickLabel','\|\|')
12			set(gca,'YtickLabel','\|\|')
13			text(Nx2/2-25,Ny2+5,'\int Q D_p dz','Fontsize',17);
14			text(-25,Ny2-5,'(a)','Fontsize',17);
15
16			subplot(2,1,2)
17			imagesc(abs(squeeze(sumM1(4:Nx-2,3:Ny-1,1)))')
18			%contourf(x,y,-squeeze(sumM1(1:Nx-2,1:Ny-2,1))',8)
19			set(gca,'DataAspectRatio',[2,2,2])
20			set(gca,'ydir','norm')
21			xlabel('X ','FontSize',17)
22			ylabel('Y ','FontSize',17)
23			set(gca,'XtickLabel','\|\|\|')
24			set(gca,'YtickLabel','\|\|\|')
25			text(Nx2/2-25,Ny2+5,'\int Q M_p dz','Fontsize',17);
26			text(-25,Ny2-5,'(b)','Fontsize',17);
27
28			colormap gray
29			cmap=colormap;
30			cmap=flipud(cmap);
31			%cmap=sqrt(cmap);
32			cmap=0.5(1-cos(3.1415cmap));
33			colormap(cmap)
34
35			%figure
36
37			%subplot(2,1,1)
38			%contourf(x,y,squeeze(sumD1(1:Nx-2,1:Ny-2,1))',8)
39			% set(gca,'DataAspectRatio',[2,2,2])
40			% set(gca,'ydir','norm')
41			%xlabel('X ','FontSize',17)
42			% ylabel('Y ','FontSize',17)
43			% set(gca,'XtickLabel','\|\|')
44			% set(gca,'YtickLabel','\|\|')
45			%text(Nx2/2-25,Ny2+5,'\int Q D_p dz','Fontsize',17);
46			%text(-25,Ny2-5,'(a)','Fontsize',17);
47
48			%subplot(2,1,2)
49			%contourf(x,y,-squeeze(sumM1(1:Nx-2,1:Ny-2,1))',8)
50			% set(gca,'DataAspectRatio',[2,2,2])
51			% set(gca,'ydir','norm')
52			%xlabel('X ','FontSize',17)
53			%ylabel('Y ','FontSize',17)
54			%set(gca,'XtickLabel','\|\|\|')
55			%set(gca,'YtickLabel','\|\|\|')
56			%text(Nx2/2-25,Ny2+5,'\int Q M_p dz','Fontsize',17);
57			%text(-25,Ny2-5,'(b)','Fontsize',17);
58
59			% cmap=flipud(cmap);
60			%colormap(cmap)