close all
clear all

[x, y] = meshgrid(1:7, 1:10);
z = peaks(10);
z = z(:, 1:7);
z(3, 1) = NaN;
z(5, 1) = NaN;

figure
pcolor(x, y, z)
colorbar
title('pcolor(x, y, z)')

figure
pcolor(x, y, z)
shading interp
colorbar
grid on
set(gca, 'layer', 'top')
title('pcolor(x, y, z), shading interp')

figure
contourf(x, y, z)
colorbar
grid on
set(gca, 'layer', 'top')
title('contourf(x, y, z)')

figure
flat_pcolor(x, y, z)
shading faceted
colorbar
title('flat_pcolor(x, y, z)')

x = [112.5000  112.6000  112.7000  112.8000  112.9000
  112.5000  112.6000  112.7000  112.8000  112.9000
  112.5000  112.6000  112.7000  112.8000  112.9000
  112.5000  112.6000  112.7000  112.8000  112.9000
  112.5000  112.6000  112.7000  112.8000  112.9000
  112.5000  112.6000  112.7000  112.8000  112.9000
  112.5000  112.6000  112.7000  112.8000  112.9000
  112.5000  112.6000  112.7000  112.8000  112.9000];
y = [-21.0000  -21.0000  -21.0000  -21.0000  -21.0000
  -20.9000  -20.9000  -20.9000  -20.9000  -20.9000
  -20.8000  -20.8000  -20.8000  -20.8000  -20.8000
  -20.7000  -20.7000  -20.7000  -20.7000  -20.7000
  -20.6000  -20.6000  -20.6000  -20.6000  -20.6000
  -20.5000  -20.5000  -20.5000  -20.5000  -20.5000
  -20.4000  -20.4000  -20.4000  -20.4000  -20.4000
  -20.3000  -20.3000  -20.3000  -20.3000  -20.3000];
z = [NaN   NaN     0     0     0
   NaN   NaN     0     0     0
   NaN     0     0     0     0
   NaN     0     0     0     0
     0     0     0     0     0
   NaN   NaN     0     0     0
   NaN   NaN   NaN   NaN     0
   NaN   NaN   NaN   NaN   NaN];
figure
flat_pcolor(x, y, z)
shading faceted
colorbar
title('flat_pcolor(x, y, z)')
