add order controls

DiffractionGratings/grating3D_data_2Dlamellar.geo

@@ -46,11 +46,13 @@ DefineConstant[
     eps_re_L_6    = {-44.9757 , Name StrCat[pp4,"/layer 6: real part of relative permittivity"]},
     eps_im_L_6    = {2.9524   , Name StrCat[pp4,"/layer 6: imag part of relative permittivity"]},
     
+    og            = {1          , Name StrCat[pp5,"/0geometrical order [-]"]},
+    oi            = {2          , Name StrCat[pp5,"/0interpolation order [-]"]},
     paramaille    = {10         , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
-    lc_scat       = {30        , Name StrCat[pp5,"/2Scatterer absolute mesh size [nm]"]},
-    PML_top       = {lambda0   , 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 [-]"]},
+    lc_scat       = {30         , Name StrCat[pp5,"/2Scatterer absolute mesh size [nm]"]},
+    PML_top       = {lambda0    , 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 [-]"]},
     refine_mesh_L_1= {1         , Name StrCat[pp5,"/7refine layers/1refine mesh layer 1 [-]"]},
     refine_mesh_L_2= {1         , Name StrCat[pp5,"/7refine layers/2refine mesh layer 2 [-]"]},
     refine_mesh_L_3= {1         , Name StrCat[pp5,"/7refine layers/3refine mesh layer 3 [-]"]},

DiffractionGratings/grating3D_data_bisin.geo

@@ -46,6 +46,8 @@ DefineConstant[
     eps_re_L_6    = {4     , Name StrCat[pp4,"/layer 6: real part of relative permittivity"]},
     eps_im_L_6    = {0     , Name StrCat[pp4,"/layer 6: imag part of relative permittivity"]},
     
+    og            = {1          , Name StrCat[pp5,"/0geometrical order [-]"]},
+    oi            = {2          , Name StrCat[pp5,"/0interpolation order [-]"]},
     paramaille    = {7     , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
     lc_scat       = {lambda0/(2*paramaille)    , Name StrCat[pp5,"/2Scatterer absolute mesh size [nm]"]},
     PML_top       = {lambda0, Name StrCat[pp5,"/4PML top thickness [nm]"]},

DiffractionGratings/grating3D_data_checker.geo

@@ -46,6 +46,8 @@ DefineConstant[
     eps_re_L_6    = {1     , Name StrCat[pp4,"/layer 6: real part of relative permittivity"]},
     eps_im_L_6    = {0     , Name StrCat[pp4,"/layer 6: imag part of relative permittivity"]},
     
+    og            = {1          , Name StrCat[pp5,"/0geometrical order [-]"]},
+    oi            = {2          , Name StrCat[pp5,"/0interpolation order [-]"]},
     paramaille    = {6     , 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]"]},

DiffractionGratings/grating3D_data_halfellipsoid.geo

@@ -46,6 +46,8 @@ DefineConstant[
     eps_re_L_6    = {4    , Name StrCat[pp4,"/layer 6: real part of relative permittivity"]},
     eps_im_L_6    = {0    , Name StrCat[pp4,"/layer 6: imag part of relative permittivity"]},
     
+    og            = {1          , Name StrCat[pp5,"/0geometrical order [-]"]},
+    oi            = {2          , Name StrCat[pp5,"/0interpolation order [-]"]},
     paramaille    = {8    , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
     lc_scat       = {10   , Name StrCat[pp5,"/2metal mesh size [nm]"]},
     PML_top       = {lambda0 , Name StrCat[pp5,"/4PML top thickness [nm]"]},

DiffractionGratings/grating3D_data_hole.geo

@@ -46,6 +46,8 @@ DefineConstant[
     eps_re_L_6    = {2.25  , Name StrCat[pp4,"/layer 6: real part of relative permittivity"]},
     eps_im_L_6    = {0     , Name StrCat[pp4,"/layer 6: imag part of relative permittivity"]},
     
+    og            = {1          , Name StrCat[pp5,"/0geometrical order [-]"]},
+    oi            = {2          , Name StrCat[pp5,"/0interpolation order [-]"]},
     paramaille    = {5     , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
     lc_scat       = {lambda0/(3*paramaille)    , Name StrCat[pp5,"/2Scatterer absolute mesh size [nm]"]},
     PML_top       = {lambda0, Name StrCat[pp5,"/4PML top thickness [nm]"]},

DiffractionGratings/grating3D_data_pyramid.geo

@@ -46,6 +46,8 @@ DefineConstant[
     eps_re_L_6    = {2.25  , Name StrCat[pp4,"/layer 6: real part of relative permittivity"]},
     eps_im_L_6    = {0     , Name StrCat[pp4,"/layer 6: imag part of relative permittivity"]},
 
+    og            = {1          , Name StrCat[pp5,"/0geometrical order [-]"]},
+    oi            = {2          , Name StrCat[pp5,"/0interpolation order [-]"]},
     paramaille    = {8     , 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]"]},

DiffractionGratings/grating3D_data_torus.geo

@@ -46,6 +46,8 @@ DefineConstant[
     eps_re_L_6    = {2.25  , Name StrCat[pp4,"/7Custom Values/layer 6: real part of relative permittivity"]},
     eps_im_L_6    = {0     , Name StrCat[pp4,"/7Custom Values/layer 6: imag part of relative permittivity"]},
     
+    og            = {1          , Name StrCat[pp5,"/0geometrical order [-]"]},
+    oi            = {2          , Name StrCat[pp5,"/0interpolation order [-]"]},
     paramaille    = {12      , Name StrCat[pp5,"/1Number of mesh elements per wavelength [-]"]},
     lc_scat       = {10       , Name StrCat[pp5,"/2Scatterer absolute mesh size [nm]"]},
     PML_top       = {lambda0*1.5 , Name StrCat[pp5,"/4PML top thickness [nm]"]},

DiffractionGratings/grating3D_postplot.py

 import numpy as np
-myDir = 'res3D/'
+import sys
+myDir = sys.argv[1]
+# myDir = res3D
 Rnm = np.loadtxt(myDir+'eff_r.txt')[:,1]
 Tnm = np.loadtxt(myDir+'eff_t.txt')[:,1]
 Q = [np.loadtxt(myDir+'temp-Q_L_%g.txt'%k)[1] for k in range(2,7)]