diff --git a/DiffractionGratings/grating3D.geo b/DiffractionGratings/grating3D.geo
index a7b2a98e3e04de273b21eba0d6543b4653d18952..84c6cdacedde9c1d933582da9883af005ae86851 100644
--- a/DiffractionGratings/grating3D.geo
+++ b/DiffractionGratings/grating3D.geo
@@ -201,7 +201,11 @@ list_lc(8) = lc_scat;
 
 // This helps meshing: The default behavior of the PointsOf techinque
 // overides points belonging to several domains
-meshing_sequence() = {1,8,2,3,5,6,7,4,9};
+If (tag_geom==7)
+  meshing_sequence() = {1,8,2,3,5,6,7,4,9};
+Else
+    meshing_sequence() = {1,8,2,3,4,6,7,5,9};
+EndIf
 
 // Start with coarsest
 Characteristic Length{:} = lc_PML;
diff --git a/DiffractionGratings/grating3D_runall.sh b/DiffractionGratings/grating3D_runall.sh
index 823b749e86a01a79b9f55ce4847d43ccfc002870..6968b90f6fb65dbd2e8c4083fcaf87c69ea75ae5 100644
--- a/DiffractionGratings/grating3D_runall.sh
+++ b/DiffractionGratings/grating3D_runall.sh
@@ -1,5 +1,5 @@
-for t in bisin checker halfellipsoid hole pyramid torus
+for t in bisin checker halfellipsoid hole pyramid torus 2Dlamellar
 do
-    gmsh grating3D.pro -setstring test_case $t
+    gmsh grating3D.pro -setstring test_case $t -
     mv res3D res3D_$t
 done
\ No newline at end of file