diff --git a/DiffractionGratings/grating3D.pro b/DiffractionGratings/grating3D.pro
index dbf6bcc6b387832e17936a0048eeb10ebe7bbea0..0bd27d548318b1ad4027acc153964a151d9e2119 100644
--- a/DiffractionGratings/grating3D.pro
+++ b/DiffractionGratings/grating3D.pro
@@ -268,9 +268,9 @@ FunctionSpace {
       { Name sn; NameOfCoef un; Function BF_Edge; Support Region[{Omega}]; Entity EdgesOf[All]; }
       { Name sn2; NameOfCoef un2; Function BF_Edge_2E;Support Region[{Omega}]; Entity EdgesOf[All]; }
       If(oi==2)
-        { Name sn3; NameOfCoef un3; Function BF_Edge_3F_b; Support Region[Omega]; Entity FacetsOf[Omega, Not SurfBloch]; }
-        { Name sn4; NameOfCoef un4; Function BF_Edge_3F_c; Support Region[Omega]; Entity FacetsOf[Omega, Not SurfBloch]; }
-        { Name sn5; NameOfCoef un5; Function BF_Edge_4E  ; Support Region[Omega]; Entity  EdgesOf[Omega, Not SurfBloch]; }
+        { Name sn3; NameOfCoef un3; Function BF_Edge_3F_b; Support Region[Omega]; Entity FacetsOf[Omega]; }
+        { Name sn4; NameOfCoef un4; Function BF_Edge_3F_c; Support Region[Omega]; Entity FacetsOf[Omega]; }
+        { Name sn5; NameOfCoef un5; Function BF_Edge_4E  ; Support Region[Omega]; Entity  EdgesOf[Omega]; }
       EndIf
     }
     Constraint {
diff --git a/DiffractionGratings/grating3D_data_2Dlamellar.geo b/DiffractionGratings/grating3D_data_2Dlamellar.geo
index 3ba4dce16ab5c30c278b3afd93b388cc317c65bb..7924d9ddf6b9184d46ec2d56246e1afdd3ec2753 100644
--- a/DiffractionGratings/grating3D_data_2Dlamellar.geo
+++ b/DiffractionGratings/grating3D_data_2Dlamellar.geo
@@ -10,8 +10,8 @@ DefineConstant[
     thetadeg      = {30    , Name StrCat[pp1,"/2theta0 [deg]"]},
     phideg        = {90     , Name StrCat[pp1,"/3phi0 [deg]"]},
     psideg        = {0     , Name StrCat[pp1,"/4psi0 [deg]"]},
-    period_x      = {100    , Name StrCat[pp2,"/1X period [nm]"]},
-    period_y      = {100  , Name StrCat[pp2,"/2Y period [nm]"]},
+    period_x      = {10    , Name StrCat[pp2,"/1X period [nm]"]},
+    period_y      = {1000  , Name StrCat[pp2,"/2Y period [nm]"]},
     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]"]},
@@ -49,7 +49,7 @@ DefineConstant[
     og            = {0          , Name StrCat[pp5,"/0geometrical order [-]"]  , Choices {0="1",1="2"} },
     oi            = {1          , Name StrCat[pp5,"/0interpolation order [-]"], Choices {0="1",1="2"} },
     paramaille    = {0.1         , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
-    lc_scat       = {500         , Name StrCat[pp5,"/2Scatterer absolute mesh size [nm]"]},
+    lc_scat       = {30         , Name StrCat[pp5,"/2Scatterer absolute mesh size [nm]"]},
     PML_top       = {100  , Name StrCat[pp5,"/4PML top thickness [nm]"]},
     PML_bot       = {100  , Name StrCat[pp5,"/5PML bot thickness [nm]"]},
     Nmax          = {1          , Name StrCat[pp5,"/6Number of non specular order to output [-]"]},