m=mphload('amb_new_L1')
% c=mphgetcoords(m,'geom1','boundary',1)
% c=mphgetcoords(m,'geom1','boundary',2)
% c=mphgetcoords(m,'geom1','boundary',3)
% c=mphgetcoords(m,'geom1','boundary',4)
[meshstats,meshdata]=mphmeshstats(m)
vtx=meshdata.vertex
tri=meshdata.elem{1}
tri=meshdata.elem{2}
tri=meshdata.elem{3}


% save getdata.mat