MITgcm_contrib/grid_bathymetry_matlab/mypcolor.m

function h=mypcolor(x,y,z)
% mypcolor(x,y,z)     simulate pcolor using image (much faster)

if nargin == 1
  h=image(x,'CDataMapping','scaled');
elseif nargin ==3
  dx=diff(x);
  dy=diff(y);
  if min(dx)==max(dx) & min(dy)==max(dy)
    h=image(x,y,z,'CDataMapping','scaled');
  else
    dx=min(x):min(abs(dx)):max(x);
    dy=min(y):min(abs(dy)):max(y);
    dz=interp2(x,y,z,dx,dy');
    h=image(dx,dy,dz,'CDataMapping','scaled');
  end
else
  error('wrong number of arguments')
end 
set(gca,'ydir','normal')
1	heimbach	1.1	function h=mypcolor(x,y,z)
2			% mypcolor(x,y,z) simulate pcolor using image (much faster)
3
4			if nargin == 1
5			h=image(x,'CDataMapping','scaled');
6			elseif nargin ==3
7			dx=diff(x);
8			dy=diff(y);
9			if min(dx)==max(dx) & min(dy)==max(dy)
10			h=image(x,y,z,'CDataMapping','scaled');
11			else
12			dx=min(x):min(abs(dx)):max(x);
13			dy=min(y):min(abs(dy)):max(y);
14			dz=interp2(x,y,z,dx,dy');
15			h=image(dx,dy,dz,'CDataMapping','scaled');
16			end
17			else
18			error('wrong number of arguments')
19			end
20			set(gca,'ydir','normal')