From 0900db59d5e72cefb0fc8fa51ff66261ddc4458f Mon Sep 17 00:00:00 2001
From: Guillaume Demesy <guillaume.demesy@fresnel.fr>
Date: Mon, 13 Feb 2023 19:35:27 +0100
Subject: [PATCH] define mur on SurfIntTop
---
DiffractionGratings/grating3D.pro | 22 ++++++++++++----------
DiffractionGratings/grating3D_runall.sh | 8 ++++----
2 files changed, 16 insertions(+), 14 deletions(-)
diff --git a/DiffractionGratings/grating3D.pro b/DiffractionGratings/grating3D.pro
index 14bf789..3163790 100644
--- a/DiffractionGratings/grating3D.pro
+++ b/DiffractionGratings/grating3D.pro
@@ -117,9 +117,10 @@
EndIf
For i In {1:6}
- mur[L~{i}] = TensorDiag[1,1,1];
+ mur[L~{i}] = TensorDiag[1,1,1];
EndFor
- mur[Scat] = TensorDiag[1,1,1];
+ mur[Scat] = TensorDiag[1,1,1];
+ mur[SurfIntTop] = TensorDiag[1,1,1];
////// PMLS
a_pml = 1.;
@@ -463,15 +464,16 @@
PostOperation {
{ Name postop_helmholtz_vector; NameOfPostProcessing postpro_helmholtz_vector ;
Operation {
-
- Print [ uz , OnElementsOf SurfIntTop, File StrCat[myDir,"uz_ZP.pos"]];
- Print [ uz , OnElementsOf SurfIntBot, File StrCat[myDir,"uz_ZM.pos"]];
+ // // check normal comp of u
+ // Print [ uz , OnElementsOf SurfIntTop, File StrCat[myDir,"uz_ZP.pos"]];
+ // Print [ uz , OnElementsOf SurfIntBot, File StrCat[myDir,"uz_ZM.pos"]];
+ // Print [ CompZu , OnPlane { {-period_x/2,-period_y/2,hh_L_6+0.5*nm} { period_x/2,-period_y/2,hh_L_6+0.5*nm} {-period_x/2, period_y/2,hh_L_6+0.5*nm} } {npts_interpX,npts_interpY} , File StrCat[myDir,"u_cut_ZM.pos"]];
+ // Print [ CompZu , OnPlane { {-period_x/2,-period_y/2,hh_L_1+thick_L_1-0.5*nm} { period_x/2,-period_y/2,hh_L_1+thick_L_1-0.5*nm} {-period_x/2, period_y/2,hh_L_1+thick_L_1-0.5*nm} } {npts_interpX,npts_interpY} , File StrCat[myDir,"u_cut_ZP.pos"]];
- Print [ epsr_xx , OnElementsOf Omega, File StrCat[myDir,"epsr_xx.pos"]];
- Print [ u , OnElementsOf Omega, File StrCat[myDir,"Edif.pos"]];
- Print [ Etot , OnElementsOf Omega, File StrCat[myDir,"Etot.pos"]];
- Print [ CompZu , OnPlane { {-period_x/2,-period_y/2,hh_L_6+0.5*nm} { period_x/2,-period_y/2,hh_L_6+0.5*nm} {-period_x/2, period_y/2,hh_L_6+0.5*nm} } {npts_interpX,npts_interpY} , File StrCat[myDir,"u_cut_ZM.pos"]];
- Print [ CompZu , OnPlane { {-period_x/2,-period_y/2,hh_L_1+thick_L_1-0.5*nm} { period_x/2,-period_y/2,hh_L_1+thick_L_1-0.5*nm} {-period_x/2, period_y/2,hh_L_1+thick_L_1-0.5*nm} } {npts_interpX,npts_interpY} , File StrCat[myDir,"u_cut_ZP.pos"]];
+ // // check : print fields everywhere
+ // Print [ epsr_xx , OnElementsOf Omega, File StrCat[myDir,"epsr_xx.pos"]];
+ // Print [ u , OnElementsOf Omega, File StrCat[myDir,"Edif.pos"]];
+ // Print [ Etot , OnElementsOf Omega, File StrCat[myDir,"Etot.pos"]];
If (FlagOutEscaFull==1)
diff --git a/DiffractionGratings/grating3D_runall.sh b/DiffractionGratings/grating3D_runall.sh
index b200ef0..3b3cf42 100644
--- a/DiffractionGratings/grating3D_runall.sh
+++ b/DiffractionGratings/grating3D_runall.sh
@@ -1,7 +1,7 @@
rm -r res3D*
-# for t in bi_sinusoidal checker half_ellipsoid hole pyramid U retrieve_2D_lamellar nanowire_solarcell convergence skewed_lattice
-for t in half_ellipsoid
+# for t in half_ellipsoid
+for t in bi_sinusoidal checker half_ellipsoid hole pyramid U retrieve_2D_lamellar nanowire_solarcell convergence skewed_lattice
do
for f in 0 1
do
@@ -10,8 +10,8 @@ do
done
done
-# for t in bi_sinusoidal checker half_ellipsoid hole pyramid U retrieve_2D_lamellar nanowire_solarcell convergence skewed_lattice
-for t in half_ellipsoid
+for t in bi_sinusoidal checker half_ellipsoid hole pyramid U retrieve_2D_lamellar nanowire_solarcell convergence skewed_lattice
+# for t in half_ellipsoid
do
for f in 0 1
do
--
GitLab