function h=gethT(TT)

tris=getTriNodeIndices1(TT);
b1=sqrt(sum((TT.Nodes(tris(:,2),:)-TT.Nodes(tris(:,1),:)).^2,2));
b2=sqrt(sum((TT.Nodes(tris(:,3),:)-TT.Nodes(tris(:,2),:)).^2,2));
b3=sqrt(sum((TT.Nodes(tris(:,1),:)-TT.Nodes(tris(:,3),:)).^2,2));
h=max(max([b1,b2,b3],[],2));

return
end

