function [f, s, x, z] = plotMe(fileName)
% fileName = '/Volumes/publicSonarData/SRTM15/com/lame.xyz.surface.grdsample.grd'
D = grdinfo(fileName);
grd = grdread(fileName);

x0 = D(1);
x1 = D(2);
y0 = D(3);
y1 = D(4);
zMin = D(5);
zMax = D(6);
pixelRegistered = D(7)
dx = D(8);
dy = D(9);

if(pixelRegistered)
    x = linspace(x0+dx/2, x1-dx/2, int32((x1-x0)/dx));
    y = linspace(y0+dy/2, y1-dy/2, int32((y1-y0)/dy));
else
    x = linspace(x0, x1, 1+int32((x1-x0)/dx));
    y = linspace(y0, y1, 1+int32((y1-y0)/dy));
end


figure
img = imagesc(x, y, (grd));
% fix the freaking lat axis...
set(imgca,'YDir','normal');

colorbar;
caxis([zMin zMax])
% caxis([-11 13])
title('Topography (m)');
xlabel('Longitude (deg)');
ylabel('Latitude (deg)');


% figure
% img = surf(x, y, grd)
% % set(imgca,'YDir','normal');
% colorbar;
% % colormap(haxby(32)); % n <=38 for some reason?
% caxis([zMin zMax])
% caxis([-11 13])
