mur fix

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.;

DiffractionGratings/grating3D_runall.sh

 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