diff --git a/geiger/grids/mortar/mortar_geiger.geo b/geiger/grids/mortar/mortar_geiger.geo
new file mode 100644
index 0000000000000000000000000000000000000000..b4ebdaeff6656841688530049bfb298557814be8
--- /dev/null
+++ b/geiger/grids/mortar/mortar_geiger.geo
@@ -0,0 +1,166 @@
+// Gmsh project created on Mon Feb 13 16:48:53 2017
+
+
+Mesh.Algorithm = 6;
+
+T = 1.0;
+
+L = 0.07*T;
+XL = L*1.3;
+
+
+S = L*.4;
+M = L*0.6;
+
+
+Point(010) = {0, 0, 0, L};
+
+Point(020) = {0.5*T, 0, 0, L};
+Point(021) = {0.5*T, 0, 0, L};
+
+Point(031) = {T, 0, 0, XL};
+
+Point(040) = {0, 0.5*T, 0, M};
+Point(042) = {0, 0.5*T, 0, M};
+
+Point(050) = {0.5*T, 0.5*T, 0, S};
+Point(051) = {0.5*T, 0.5*T, 0, M};
+Point(052) = {0.5*T, 0.5*T, 0, S};
+Point(053) = {0.5*T, 0.5*T, 0, S};
+
+Point(061) = {0.625*T, 0.5*T, 0, M};
+Point(063) = {0.625*T, 0.5*T, 0, M};
+Point(064) = {0.625*T, 0.5*T, 0, S};
+
+Point(071) = {0.75*T, 0.5*T, 0, M};
+Point(074) = {0.75*T, 0.5*T, 0, S};
+Point(077) = {0.75*T, 0.5*T, 0, M};
+
+Point(081) = {T, 0.5*T, 0, L};
+Point(087) = {T, 0.5*T, 0, L};
+
+Point(092) = {0.5*T, 0.625*T, 0, S};
+Point(093) = {0.5*T, 0.625*T, 0, S};
+Point(095) = {0.5*T, 0.625*T, 0, S};
+
+Point(103) = {0.625*T, 0.625*T, 0, S};
+Point(104) = {0.625*T, 0.625*T, 0, S};
+Point(105) = {0.625*T, 0.625*T, 0, S};
+Point(106) = {0.625*T, 0.625*T, 0, S};
+
+Point(114) = {0.75*T, 0.625*T, 0, S};
+Point(116) = {0.75*T, 0.625*T, 0, S};
+Point(117) = {0.75*T, 0.625*T, 0, M};
+
+Point(122) = {0.5*T, 0.75*T, 0, M};
+Point(125) = {0.5*T, 0.75*T, 0, S};
+Point(128) = {0.5*T, 0.75*T, 0, M};
+
+Point(135) = {0.625*T, 0.75, 0, S};
+Point(136) = {0.625*T, 0.75, 0, S};
+Point(138) = {0.625*T, 0.75, 0, M};
+
+Point(146) = {0.75*T, 0.75*T, 0, S};
+Point(147) = {0.75*T, 0.75*T, 0, M};
+Point(148) = {0.75*T, 0.75*T, 0, M};
+Point(149) = {0.75*T, 0.75*T, 0, M};
+
+Point(157) = {T, 0.75*T, 0, L};
+Point(159) = {T, 0.75*T, 0, L};
+
+Point(162) = {0, T, 0, L};
+
+Point(172) = {0.5*T, T, 0, L};
+Point(178) = {0.5*T, T, 0, M};
+
+Point(188) = {0.75*T, T, 0, L};
+Point(189) = {0.75*T, T, 0, L};
+
+Point(199) = {T, T, 0, L};
+
+// Points added for the virtual fractures. 
+
+
+
+
+Line(010) = {010, 020};
+Line(020) = {020, 050};
+Line(030) = {050, 040};
+Line(040) = {040, 010};
+
+Line(051) = {021, 031};
+Line(061) = {031, 081};
+Line(071) = {081, 071};
+Line(081) = {071, 061};
+Line(091) = {061, 051};
+Line(021) = {051, 021};
+
+Line(032) = {042, 052};
+Line(102) = {052, 092};
+Line(112) = {092, 122};
+Line(122) = {122, 172};
+Line(132) = {172, 162};
+Line(142) = {162, 042};
+
+Line(093) = {053, 063};
+Line(153) = {063, 103};
+Line(163) = {103, 093};
+Line(103) = {093, 053};
+
+Line(084) = {064, 074};
+Line(174) = {074, 114};
+Line(184) = {114, 104};
+Line(154) = {104, 064};
+
+Line(165) = {095, 105};
+Line(195) = {105, 135};
+Line(205) = {135, 125};
+Line(115) = {125, 095};
+
+Line(186) = {106, 116};
+Line(216) = {116, 146};
+Line(226) = {146, 136};
+Line(196) = {136, 106};
+
+Line(077) = {077, 087};
+Line(237) = {087, 157};
+Line(247) = {157, 147};
+Line(217) = {147, 117};
+Line(177) = {117, 077};
+
+Line(208) = {128, 138};
+Line(228) = {138, 148};
+Line(258) = {148, 188};
+Line(268) = {188, 178};
+Line(128) = {178, 128};
+
+Line(249) = {149, 159};
+Line(279) = {159, 199};
+Line(289) = {199, 189};
+Line(259) = {189, 149};
+
+
+
+
+
+Line Loop(10) = {010, 020, 030, 040};
+Line Loop(11) = {051, 061, 071, 081, 091, 021};
+Line Loop(12) = {032, 102, 112, 122, 132, 142};
+Line Loop(13) = {093, 153, 163, 103};
+Line Loop(14) = {084, 174, 184, 154};
+Line Loop(15) = {165, 195, 205, 115};
+Line Loop(16) = {186, 216, 226, 196};
+Line Loop(17) = {077, 237, 247, 217, 177};
+Line Loop(18) = {208, 228, 258, 268, 128};
+Line Loop(19) = {249, 279, 289, 259};
+
+Plane Surface(20) = {10};
+Plane Surface(21) = {11};
+Plane Surface(22) = {12};
+Plane Surface(23) = {13};
+Plane Surface(24) = {14};
+Plane Surface(25) = {15};
+Plane Surface(26) = {16};
+Plane Surface(27) = {17};
+Plane Surface(28) = {18};
+Plane Surface(29) = {19};