up

DiffractionGratings/grating3D.pro

@@ -59,7 +59,7 @@ Function{
 	om0         = 2*Pi*cel/lambda0;
 	k0          = 2.*Pi/lambda0;
 	Ae          = 1;
-	Pinc        =  0.5*Ae^2*Sqrt[eps_re_L_1*ep0/mu0] * Cos[theta0];
+	Pinc        = 0.5*Ae^2*Sqrt[eps_re_L_1*ep0/mu0] * Cos[theta0];
 
   // permittivities
   For i In {1:6}

DiffractionGratings/grating3D_data_hole.geo

@@ -10,12 +10,12 @@ DefineConstant[
     thetadeg      = {0     , Name StrCat[pp1,"/2theta0 [deg]"]},
     phideg        = {0     , Name StrCat[pp1,"/3phi0 [deg]"]},
     psideg        = {45    , Name StrCat[pp1,"/4psi0 [deg]"]},
-    period_x      = {999.9 , Name StrCat[pp2,"/1X period [nm]"]},
-    period_y      = {999.9 , Name StrCat[pp2,"/2Y period [nm]"]},
+    period_x      = {1000.001 , Name StrCat[pp2,"/1X period [nm]"]},
+    period_y      = {1000.001 , Name StrCat[pp2,"/2Y period [nm]"]},
     thick_L_1     = {100   , Name StrCat[pp2,"/3thickness layer 1 [nm] (superstrate)"]},
-    thick_L_2     = {100   , Name StrCat[pp2,"/4thickness layer 2 [nm]"]},
+    thick_L_2     = {50   , Name StrCat[pp2,"/4thickness layer 2 [nm]"]},
     thick_L_3     = {50    , Name StrCat[pp2,"/5thickness layer 3 [nm]"]},
-    thick_L_4     = {100   , Name StrCat[pp2,"/6thickness layer 4 [nm]"]},
+    thick_L_4     = {50   , 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)"]},
 
@@ -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    = {5     , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
-    lc_scat       = {lambda0/(3*paramaille)    , Name StrCat[pp5,"/2Scatterer absolute mesh size [nm]"]},
+    paramaille    = {4.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]"]},
     Nmax          = {2      , Name StrCat[pp5,"/6Number of non specular order to output [-]"]},