39-the-python…

acurnier_2D_2.geo
View Options

			cl1 = 2;
			cl2 = -0.15;
			cl3 = -0.5;
			Point(1) = {-1, cl3, 0, cl1};
			Point(2) = {-1, 0.5, 0, cl1};
			Point(3) = {1, 0.5, 0, cl1};
			Point(4) = {1, cl3, 0, cl1};
			Point(5) = {0.35, 0.7+cl2, 0, cl1};
			Point(6) = {0.35, 1.4+cl2, 0, cl1};
			Point(7) = {-0.35, 1.4+cl2, 0, cl1};
			Point(8) = {-0.35, 0.7+cl2, 0, cl1};
			Line(1) = {1, 4};
			Line(2) = {4, 3};
			Line(3) = {3, 2};
			Line(4) = {2, 1};
			Line(5) = {8, 5};
			Line(6) = {5, 6};
			Line(7) = {6, 8};
			Line(8) = {7, 8};
			Line(9) = {6, 7};
			Line Loop(10) = {5, 6, 7};
			Plane Surface(11) = {10};
			Line Loop(12) = {-7, 9, 8};
			Plane Surface(13) = {12};
			Line Loop(14) = {1, 2, 3, 4};
			Plane Surface(15) = {14};
			Physical Line("Contact") = {3};
			Physical Line("Top") = {9};
			Physical Line("Bottom") = {1};
			Physical Surface("Top_body") = {11, 13};
			Physical Surface("Bottom_body") = {15};