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	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')