[Z,R] = arcgridread('montcanhead.grd');
mapshow(Z,R,'DisplayType','surface');
xlabel('x (easting in meters)'); ylabel('y (northing in meters)')
colormap(demcmap(Z))

% View the terrain in 3D
axis normal; view(3); axis equal; grid on
zlabel('elevation in feet')