0001 function h=pslice2(x,y,u,cax)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017 ind=find(~isnan(u));
0018 if(nargin<4),
0019 cax(1)=min(u(ind));
0020 cax(2)=max(u(ind));
0021
0022 if(cax(2)==cax(1))
0023 cax(1) = cax(1)-0.01;
0024 cax(2) = cax(2)+0.01;
0025 end
0026
0027
0028 else
0029 if(isnan(cax(1))),cax(1)=min(u(ind));end;
0030 if(isnan(cax(2))),cax(2)=max(u(ind));end;
0031 end
0032 u(ind)=max(u(ind),cax(1));
0033 u(ind)=min(u(ind),cax(2));
0034 [m,n]=size(x);
0035 if(min(m,n)>1),
0036 x=.25*(x(1:m-1,1:n-1)+x(2:m,1:n-1)+x(1:m-1,2:n)+x(2:m,2:n));
0037 y=.25*(y(1:m-1,1:n-1)+y(2:m,1:n-1)+y(1:m-1,2:n)+y(2:m,2:n));
0038 u(1,:)=[];
0039 u(:,1)=[];
0040 end
0041 h=pcolor(x,y,u);
0042 colormap('jet');...
0043 shading('flat');
0044 caxis(cax);
0045
0046 matlab_version = version;
0047 if ( matlab_version(1) == '5' )
0048 set ( gca, 'DataAspectRatio', [1 1 1] );
0049 else
0050 set(gca,'aspectratio',[NaN 1])
0051 end