diff --git a/DiffractionGratings/grating3D.pro b/DiffractionGratings/grating3D.pro
index 14bf789b86c7df77d6f2ddedd1dca9eda4119433..1a40f0dc9ec5e13378530dd21ba9dfe4fe2119ad 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.;
diff --git a/DiffractionGratings/grating3D_runall.sh b/DiffractionGratings/grating3D_runall.sh
index b200ef0a72f4f55c3f35c380dd384208b28f5bcd..67d992f5840df4d856dc6cbf2f532d6256832ab3 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 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