2D lamellar much better with linkfacets enable

DiffractionGratings/grating3D.pro

@@ -31,7 +31,7 @@ Group {
 
   
 	SurfDirichlet  = Region[{401,402}];
-  SurfBloch      = Region[{SurfBlochXm,SurfBlochXp,SurfBlochYm,SurfBlochYp}];
+  SurfBloch      = Region[{SurfBlochYm,SurfBlochYp}];
 
   L_1 = Region[{L_1_temp,SurfIntTop}];
   L_6 = Region[{L_6_temp,SurfIntBot}];
@@ -214,8 +214,8 @@ Constraint {
   }
   { Name BlochX;
     Case {
-    { Region SurfBlochXp; Type LinkCplx ; RegionRef SurfBlochXm;
-      Coefficient dephX[]; Function Vector[$X-period_x,$Y,$Z] ; }
+    // { Region SurfBlochXp; Type LinkCplx ; RegionRef SurfBlochXm;
+    //   Coefficient dephX[]; Function Vector[$X-period_x,$Y,$Z] ; }
     }
   }
   { Name BlochY;

DiffractionGratings/grating3D_data_2Dlamellar.geo

@@ -15,8 +15,8 @@ DefineConstant[
     thick_L_1     = {200    , Name StrCat[pp2,"/3thickness layer 1 [nm] (superstrate)"]},
     thick_L_2     = {200    , Name StrCat[pp2,"/4thickness layer 2 [nm]"]},
     thick_L_3     = {1050   , Name StrCat[pp2,"/5thickness layer 3 [nm]"]},
-    thick_L_4     = {150     , Name StrCat[pp2,"/6thickness layer 4 [nm]"]},
-    thick_L_5     = {150     , Name StrCat[pp2,"/7thickness layer 5 [nm]"]},
+    thick_L_4     = {100     , Name StrCat[pp2,"/6thickness layer 4 [nm]"]},
+    thick_L_5     = {100     , Name StrCat[pp2,"/7thickness layer 5 [nm]"]},
     thick_L_6     = {100    , Name StrCat[pp2,"/8thickness layer 6 [nm] (substrate)"]},
 
     tag_geom      = {  7      , Name StrCat[pp3,"/0Shape"], Choices {1="Pyramid",2="Cylindrical Hole",3="Torus",4="HalfEllipspoid",5="Checkerboard",6="bi-sinusoidal",7="2D lamellar"}},