conv

DiffractionGratings/grating3D_data_conv.geo

@@ -18,6 +18,7 @@ DefineConstant[
     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)"]},
+    xsideg        = {0    , Name StrCat[pp2,"/9skew angle [deg]"],Visible 0},
 
     tag_geom      = {  4      , Name StrCat[pp3,"/0Shape"], Choices {1="Pyramid",2="Cylindrical Hole",3="Torus",4="HalfEllipspoid",5="Checkerboard",6="bi-sinusoidal",7="2D lamellar"}},
     rx            = {107  , Name StrCat[pp3,"/1rx"]},
@@ -99,14 +100,24 @@ For k In {1:5}
 EndFor
 PML_bot_hh = hh_L_6-PML_bot;
 PML_top_hh = hh_L_1+thick_L_1;
+hh_L_7     = PML_bot_hh;
+hh_L_0     = PML_top_hh;
+thick_L_7  = PML_bot;
+thick_L_0  = PML_top;
 
 theta0 = thetadeg*Pi/180;
 phi0   = phideg*Pi/180;
 psi0   = psideg*Pi/180;
+xsi    = xsideg*Pi/180;
+
+dyc = period_y*Cos[xsi];
+dys = period_y*Sin[xsi];
 
 DomainZsizeSca  = PML_top_hh+PML_bot-(hh_L_6-PML_bot);
 DomainZsizeTot  = PML_top_hh-hh_L_6;
 npts_interpX    = period_x/InterpSampling;
 npts_interpY    = period_y/InterpSampling;
+npts_checkpoyX  = 50;
+npts_checkpoyY  = 50;
 npts_interpZSca = DomainZsizeSca/InterpSampling;
 npts_interpZTot = DomainZsizeTot/InterpSampling;