mesh.entity =load('mesh.entity.dat');
mesh.elem2node =load('mesh.elem2node.dat');
mesh.node2coord =load('mesh.node2coord.dat');