debug materials

DiffractionGratings/grating3D.pro

@@ -70,7 +70,7 @@ Function{
       EndIf
     Else
       For j In {2:nb_available_materials}
-        If(flag_mat~{i}==j-1)
+        If(flag_mat~{i}==j)
           epsr[L~{i}]  = Complex[epsr_re_interp_mat~{j}[lambda0/nm*1e-9] , epsr_im_interp_mat~{j}[lambda0/nm*1e-9]] * TensorDiag[1,1,1];
           If (i==1)
             epsr1[] = Complex[epsr_re_interp_mat~{j}[lambda0/nm*1e-9] , epsr_im_interp_mat~{j}[lambda0/nm*1e-9]];
@@ -367,8 +367,8 @@ PostOperation {
         Print [ Etot , OnPlane { {0,-period_y/2,hh_L_6} {0,period_y/2,hh_L_6} {0,-period_y/2,hh_L_1+thick_L_1} } {npts_interpY,npts_interpZTot} , File StrCat[myDir,"Etot_cut_X=0.pos"], Name "Etot_cut_X=0"];
       EndIf
       If (FlagOutPoyCut==1)
-        Print [ Poy_tot , OnPlane { {-period_x/2,0,hh_L_6} {period_x/2,0,hh_L_6} {-period_x/2,0,hh_L_1+thick_L_1} } {npts_interpX,npts_interpZSca} , File StrCat[myDir,"Poy_tot_cut_Y=0.pos"], Name "Poy_tot_cut_Y=0"];
-        Print [ Poy_tot , OnPlane { {0,-period_y/2,hh_L_6} {0,period_y/2,hh_L_6} {0,-period_y/2,hh_L_1+thick_L_1} } {npts_interpY,npts_interpZSca} , File StrCat[myDir,"Poy_tot_cut_X=0.pos"], Name "Poy_tot_cut_X=0"];
+        Print [ Poy_tot , OnPlane { {-period_x/2,0,hh_L_6} {period_x/2,0,hh_L_6} {-period_x/2,0,hh_L_1+thick_L_1} } {npts_interpX,npts_interpZTot} , File StrCat[myDir,"Poy_tot_cut_Y=0.pos"], Name "Poy_tot_cut_Y=0"];
+        Print [ Poy_tot , OnPlane { {0,-period_y/2,hh_L_6} {0,period_y/2,hh_L_6} {0,-period_y/2,hh_L_1+thick_L_1} } {npts_interpY,npts_interpZTot} , File StrCat[myDir,"Poy_tot_cut_X=0.pos"], Name "Poy_tot_cut_X=0"];
       EndIf
       If (FlagOutEscaFull==1)
         If (Flag_interp_cubic==1)
@@ -393,6 +393,7 @@ PostOperation {
       EndIf
       // // For DEBUG
       Print [ epsr_xx , OnElementsOf Omega, File StrCat[myDir,"epsr_xx.pos"]];
+      // Print [ epsr_xx , OnElementsOf Scat, File StrCat[myDir,"epsr_xx.pos"]];
       // Print [ Ecm , OnElementsOf Omega_plot, File StrCat[myDir,"Ecm.pos"]];
       // Print [ E1, OnElementsOf Omega_plot, File StrCat[myDir,"E1.pos"]];
       // Print [ Damp_pml_top, OnElementsOf PMLtop, File StrCat[myDir,"Damp_pml_top.pos"]];

DiffractionGratings/grating3D_postplot.py

 import numpy as np
-Rnm = np.loadtxt('eff_r.txt')[:,1]
-Tnm = np.loadtxt('eff_t.txt')[:,1]
-Q = [np.loadtxt('temp-Q_L_%g.txt'%k)[1] for k in range(2,7)]
-Q.append(np.loadtxt('temp-Q_scat.txt')[1])
+myDir = 'run_results3D/'
+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)]
+Q.append(np.loadtxt(myDir+'temp-Q_scat.txt')[1])
 Q=np.array(Q)
 TOT = Rnm.sum()+Tnm.sum()+Q.sum()
 print('Rtot',Rnm.sum())