diff --git a/DiffractionGratings/grating3D.pro b/DiffractionGratings/grating3D.pro
index 14bf789b86c7df77d6f2ddedd1dca9eda4119433..31637905019f034743a1524964fbc73d3d379218 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 b200ef0a72f4f55c3f35c380dd384208b28f5bcd..3b3cf42452e93acfc8de442ddc9bad4e63004e75 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