up pyramid

DiffractionGratings/grating3D.geo

@@ -116,16 +116,18 @@ If (tag_geom==1)
   Line(97) = {25, 65};
   Line(98) = {65, 27};
   Line(99) = {65, 29};
-  Line(100) = {31, 65};2
-  Curve Loop(49) = {38, -982, -97};
-  Plane Surface(49) = {49};2
-  Curve Loop(50) = {98, 48,2 100};
-  Plane Surface(50) = {50};2
-  Curve Loop(51) = {99, 42,2 100};
+  Line(100) = {31, 65};
+  Curve Loop(49) = {97, 99, -46};
+  Plane Surface(49) = {49};
+  Curve Loop(50) = {99, 42, 100};
+  Plane Surface(50) = {50};
+  Curve Loop(51) = {98, -38, 97};
   Plane Surface(51) = {51};
-  Curve Loop(52) = {99, -46, 97};
+  Curve Loop(52) = {98, 48, 100};
   Plane Surface(52) = {52};
-  Surface Loop(9) = {52, 51, 50, 49, 24};
+  Curve Loop(53) = {38, 48, -42, -46};
+  Plane Surface(53) = {53};
+  Surface Loop(9) = {49, 51, 52, 50, 53};
   Volume(9) = {9};
 EndIf
 

DiffractionGratings/grating3D_data_pyramid.geo

@@ -48,7 +48,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    = {8     , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
+    paramaille    = {10     , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
     lc_scat       = {lambda0/(1.5*paramaille)    , Name StrCat[pp5,"/2Scatterer absolute mesh size [nm]"]},
     PML_top       = {lambda0  , Name StrCat[pp5,"/4PML top thickness [nm]"]},
     PML_bot       = {lambda0  , Name StrCat[pp5,"/5PML bot thickness [nm]"]},