clear
mesh=load('mesh.dat');
segment_mesh=load('boundary_segment_mesh.dat');
u=load('u.dat');
vertices=load('vertices.dat');

s=size(vertices);
num_verts=s(1);

for i=1:num_verts
    x(i)=vertices(i,1);
    y(i)=vertices(i,2);
    z(i)=u(i);
end

figure
s=size(segment_mesh);
num_segments=s(1);
for s=1:num_segments
    xs(1)=x(segment_mesh(s,1));
    xs(2)=x(segment_mesh(s,2));
    
    ys(1)=y(segment_mesh(s,1));
    ys(2)=y(segment_mesh(s,2));
    
    zs(1)=z(segment_mesh(s,1));
    zs(2)=z(segment_mesh(s,2));
    
    plot3(xs,ys,zs,'k--','LineWidth',8);
    hold on
end

trisurf(mesh,x,y,z,'EdgeColor','none');

%plot3(xs,ys,zs);