2D lamellar much better with linkfacets enable

DiffractionGratings/grating3D_data_2Dlamellar.geo

@@ -15,14 +15,14 @@ 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     = {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)"]},
+    thick_L_4     = {50     , Name StrCat[pp2,"/6thickness layer 4 [nm]"]},
+    thick_L_5     = {50     , Name StrCat[pp2,"/7thickness layer 5 [nm]"]},
+    thick_L_6     = {50    , 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"}},
+    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"}},
     rx            = { 0     , Name StrCat[pp3,"/1rx"]},
-    ry            = {500        , Name StrCat[pp3,"/2ry"]},
-    rz            = {1000     , Name StrCat[pp3,"/3rz"]},
+    ry            = {500    , Name StrCat[pp3,"/2ry"]},
+    rz            = {1000   , Name StrCat[pp3,"/3rz"]},
     flag_mat_scat = { 0       , Name StrCat[pp3,"/4Scatterer permittivity model"], Choices {0="Custom (Value Below)",1="SiO2",2="Ag (palik)",3="Al (palik)",4="Au (johnson)",5="Nb2O5",6="ZnSe",7="MgF2",8="TiO2",9="PMMA",10="Si",11="ITO",12="Cu (palik)"} },
     eps_re_Scat   = {-44.9757  , Name StrCat[pp3,"/7eps_re_Scat"]},
     eps_im_Scat   = {2.9524  , Name StrCat[pp3,"/8eps_im_Scat"]},
@@ -48,8 +48,8 @@ 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    = {15         , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
-    lc_scat       = {15         , Name StrCat[pp5,"/2Scatterer absolute mesh size [nm]"]},
+    paramaille    = {10         , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
+    lc_scat       = {20         , Name StrCat[pp5,"/2Scatterer absolute mesh size [nm]"]},
     PML_top       = {lambda0*1  , Name StrCat[pp5,"/4PML top thickness [nm]"]},
     PML_bot       = {lambda0/5  , Name StrCat[pp5,"/5PML bot thickness [nm]"]},
     Nmax          = {1          , Name StrCat[pp5,"/6Number of non specular order to output [-]"]},