--- MITgcm/verification/lab_sea/matlab/mypcolor.m	2002/12/05 08:06:18	1.1
+++ MITgcm/verification/lab_sea/matlab/mypcolor.m	2002/12/05 08:43:03	1.2
@@ -0,0 +1,20 @@
+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')