//  square_split.edp
//
//  Discussion:
//
//    This example reads a mesh created by mesh_write.edp.
//
//  Location:
//
//    http://people.sc.fsu.edu/~jburkardtfreefem++/square_split/square_split.edp
//
//  Licensing:
//
//    This code is distributed under the MIT license.
//
//  Modified:
//
//    26 July 2015
//
//  Author:
//
//    John Burkardt
//
//  Reference:
//
//    Frederic Hecht,
//    Freefem++,
//    Third Edition, version 3.22
//
cout << "\n";
cout << "SQUARE_SPLIT\n";
cout << "  FreeFem++ version.\n";
cout << "  Read a mesh created by mesh_write.\n";
//
//  Retrieve mesh data from the MSH files
//    square_split_o.msh
//  and
//    square_split_g.msh
//
int n = 10;
mesh Th = readmesh ( "square_split_o.msh" );
//
//  Plot the mesh.
//
plot ( Th, wait = true, ps = "square_split.ps" );
//
//  Here are some things you can query about a mesh.
//
cout << "\n";
cout << "  Some mesh data:\n";
cout << "\n";
cout << "  Th.area           = " << Th.area << "\n";;
cout << "  Th.nt             = " << Th.nt << "\n";
cout << "  Th.nv             = " << Th.nv << "\n";
//
//  Terminate.
//
cout << "\n";
cout << "square_split\n";
cout << "  Normal end of execution.\n";

exit ( 0 );