visibility bottom surface

DiffractionGratings/grating3D.geo

@@ -238,6 +238,7 @@ Physical Surface("STOP",301 ) = Surface In BoundingBox{-period_x/2-E,-period_y/2
 Physical Surface("SBOT",302 ) = Surface In BoundingBox{-period_x/2-E,-period_y/2-E, hh_L_6-e,period_x/2+E    , period_y/2+E, hh_L_6+e};
 Physical Surface("SPMLTOP",401 ) = Surface In BoundingBox{-period_x/2-E,-period_y/2-E, PML_top_hh+PML_top-e,period_x/2+E, period_y/2+E, PML_top_hh+PML_top+e};
 Physical Surface("SPMLBOT",402 ) = Surface In BoundingBox{-period_x/2-E,-period_y/2-E, PML_bot_hh-e,period_x/2+E, period_y/2+E, PML_bot_hh+e};
+Physical Surface("SVIS",500 )    = Surface In BoundingBox{-period_x/2-E,-period_y/2-E, hh_L_3-e,period_x/2+E    , period_y/2+E, hh_L_3+e};
 
 pts_PMLBOT()   = PointsOf{ Physical Volume{1}; };
 pts_LAYER_L6() = PointsOf{ Physical Volume{2}; };

DiffractionGratings/grating3D_data_hole.geo

@@ -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    = {4.         , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
+    paramaille    = {5.         , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
     lc_scat       = {lambda0/(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]"]},

DiffractionGratings/gratings_clean_all.sh

@@ -5,5 +5,3 @@ rm *.pre
 rm *.db
 rm *.pdf
 rm *.npz
-
-