timour_matlab/mscripts/DIAG1b.m

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


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

sumD1(4:Nx-2,Ny,1)=0;
sumD1(4:Nx-2,Ny-1,1)=0;
sumM1(4:Nx-2,Ny-1,1)=0;
sumM1(4,1:Ny-1,1)=0;
subplot(2,1,1)
contour(x,y,abs(squeeze(sumD1(4:Nx-2,1:Ny,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)
contour(x,y,abs(squeeze(sumM1(4:Nx-2,1:Ny,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=sqrt(cmap);
cmap=0.5*(1-cos(3.1415*cmap));
cmap=flipud(cmap)
cmap=[0 0 0];
colormap(cmap) 

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