From ffdb67bdd639cff129f848d703a5687b4d6835ed Mon Sep 17 00:00:00 2001
From: Guillaume Demesy <guillaume.demesy@fresnel.fr>
Date: Tue, 26 Nov 2019 22:37:17 +0100
Subject: [PATCH] conv

---
 DiffractionGratings/grating3D_data_conv.geo | 11 +++++++++++
 1 file changed, 11 insertions(+)

diff --git a/DiffractionGratings/grating3D_data_conv.geo b/DiffractionGratings/grating3D_data_conv.geo
index 2908df2..9dfef1a 100644
--- a/DiffractionGratings/grating3D_data_conv.geo
+++ b/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;
-- 
GitLab