% Plot the interpolated Merimbula bathymetry data

%To read saved files use open command

load building_boundary_Hobart_58890_triangles.pts
pts = building_boundary_Hobart_58890_triangles;
x = pts(:,1);
y = pts(:,2);

load building_boundary_Hobart_58890_triangles.tri
tri = building_boundary_Hobart_58890_triangles + 1;

load building_boundary_Hobart_58890_triangles_pt01.bed
z = building_boundary_Hobart_58890_triangles_pt01;

x_max = ceil(max(x)/1000)*1000;
x_min = floor(min(x)/1000)*1000;

y_max = ceil(max(y)/1000)*1000;
y_min = floor(min(y)/1000)*1000;

z_max = ceil(max(z)/100)*100;
z_min = floor(min(z)/100)*100;
cmax = z_max
cmin = z_min

figure(1)

axis([x_min x_max y_min y_max z_min z_max])
trisurf(tri,x,y,z)
shading interp
colormap(jet)
caxis([cmin cmax])
view(0,90)
colorbar
title('\alpha = 0.01')
axis off
save building_boundary_Hobart_58890_triangles_pt01.fig
