//  Discussion:
//
//    This example reads a mesh created by mesh_write.edp.
//
//  Location:
//
//    http://people.sc.fsu.edu/~jburkardt/freefem_src/mesh_read/mesh_read.edp
//
//  Licensing:
//
//    This code is distributed under the MIT license.
//
//  Modified:
//
//    21 June 2015
//
//  Author:
//
//    John Burkardt
//
//  Reference:
//
//    Frederic Hecht,
//    Freefem++,
//    Third Edition, version 3.22
//
cout << "\n";
cout << "mesh_read\n";
cout << "  FreeFem++ version\n";
cout << "  Read a file defining the mesh to be used.\n";
//
//  Define Th, the triangulation of the "left" side of the boundaries,
//  by reading in the mesh created by the "mesh_write" program.
//
mesh Th = readmesh ( "mesh_write_mesh.msh" );
//
//  Plot the mesh.
//
plot ( Th, ps = "mesh_read.ps" );
//
//  When reading a mesh file back into FreeFem++, the "names" of the
//  borders are lost, but the label values are available.
//
//  Terminate.
//
cout << "\n";
cout << "mesh_read\n";
cout << "  Normal end of execution.\n";

exit ( 0 );