postpro python

DiffractionGratings/grating3D.pro

@@ -392,10 +392,10 @@ PostOperation {
         EndIf
       EndIf
       // // For DEBUG
-      // Print [ epsr_xx , OnElementsOf Omega, 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"];
+      Print [ epsr_xx , OnElementsOf Omega, 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"]];
       // For i In {0:nb_slice-1}
       //   Print [ Etot , OnPlane { {-period_x/2+small_delta,-period_y/2+small_delta,zcut_sub~{i}} {-period_x/2+small_delta,period_y/2-small_delta,zcut_sub~{i}} {period_x/2-small_delta,-period_y/2+small_delta,zcut_sub~{i}} } {ninterv_integ,ninterv_integ} , File StrCat[myDir,> "./Views/Etot_XYcut.out",Format Table];
       //   Print [ Edif , OnPlane { {-period_x/2+small_delta,-period_y/2+small_delta,zcut_sup~{i}} {-period_x/2+small_delta,period_y/2-small_delta,zcut_sup~{i}} {period_x/2-small_delta,-period_y/2+small_delta,zcut_sup~{i}} } {ninterv_integ,ninterv_integ} , File StrCat[myDir,> "./Views/Edif_XYcut.out",Format Table];

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])
+Q=np.array(Q)
+TOT = Rnm.sum()+Tnm.sum()+Q.sum()
+print('Rtot',Rnm.sum())
+print('Ttot',Tnm.sum())
+print('Atot',Q.sum())
+print('TOT ',TOT)
\ No newline at end of file