Unexpected…

square.geo
View Options

			// CUBE ON CUBE (2 BODIES) 3D:

			// REMARKS:
			// Physical Surfaces are defined so that load_mesh_msh can load
			// directly the surface element structure. When creating Plane Surface,
			// the surface normal has to point to the inside of the body. Otherwise the
			// solvecontact3d algorithm won't work. The Physical Surface number
			// corresponds to the face number. The face numbering has to start with 1!
			// Physical Volume defines the material of each body. This can be read
			// by the load_mesh_msh function in adlib.

			// Mesh size
			h = 0.02; // Top cube
			hy = 0.5; // Top cube

			// Dimensions of top cube
			Lx = 1;
			Ly = 1;


			// ------------------------------------------
			// Geometry
			// ------------------------------------------

			// Base Cube
			Point(101) = { 0.0, 0.0, 0.0, h}; // Bottom Face
			Point(102) = { Lx, 0.0, 0.0, h}; // Bottom Face
			Point(103) = { Lx, Ly, 0.0, h}; // Bottom Face
			Point(104) = { 0.0, Ly, 0.0, h}; // Bottom Face

			// Base Cube
			Line(101) = {101,102}; // Bottom Face
			Line(102) = {102,103}; // Bottom Face
			Line(103) = {103,104}; // Bottom Face
			Line(104) = {104,101}; // Bottom Face

			// Base Cube
			Line Loop(101) = {101:104};

			// Base Cube
			Plane Surface(101) = {101};
			Physical Line("Edges") = {101, 102, 103, 104};
			Physical Surface(6) = {101};