timour_matlab/mscripts/DIAG5a.m



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)
contourf(x,y,abs(squeeze(sumD1(4:Nx-2,1:Ny,1))'),-0.00:0.0007:0.01)
 set(gca,'DataAspectRatio',[1,2,2]) 
        set(gca,'ydir','norm')

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

xlabel('X ','FontSize',17)
  ylabel('Y ','FontSize',17,'Rotation',[0])  
 set(gca,'XtickLabel','||')
 set(gca,'YtickLabel','||')
text(Nx2/2-15,Ny+5,'\int D_p dz','Fontsize',17);
text(-25,Ny-5,'(a)','Fontsize',17);

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