[cig-commits] r21826 - seismo/2D/SPECFEM2D/trunk/src/meshfem2D

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-04-11 15:45:29 -0700 (Thu, 11 Apr 2013)
New Revision: 21826

Modified:
   seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_materials.f90
Log:
put the test on Cs back to zero because acoustic media have Cs == 0


Modified: seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_materials.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_materials.f90	2013-04-11 22:39:22 UTC (rev 21825)
+++ seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_materials.f90	2013-04-11 22:45:29 UTC (rev 21826)
@@ -116,7 +116,8 @@
         QKappa(i) = val5read
         Qmu(i) = val6read
 
-        if(rho_s(i) <= 0.00000001d0 .or. cp(i) <= 0.00000001d0 .or. cs(i) < 0.00000001d0) &
+! for Cs we use a less restrictive test because acoustic media have Cs exactly equal to zero
+        if(rho_s(i) <= 0.00000001d0 .or. cp(i) <= 0.00000001d0 .or. cs(i) < 0.d0) &
             stop 'negative value of velocity or density'
         if(QKappa(i) <= 0.00000001d0 .or. Qmu(i) <= 0.00000001d0) stop 'non-positive value of QKappa or Qmu'
 
@@ -143,7 +144,8 @@
         cp(i) = sqrt(val4read/val0read)
         cs(i) = sqrt(val6read/val0read)
 
-        if(rho_s(i) <= 0.00000001d0 .or. cp(i) <= 0.00000001d0 .or. cs(i) < 0.00000001d0) &
+! for Cs we use a less restrictive test because acoustic media have Cs exactly equal to zero
+        if(rho_s(i) <= 0.00000001d0 .or. cp(i) <= 0.00000001d0 .or. cs(i) < 0.d0) &
             stop 'negative value of velocity or density'
         if(QKappa(i) <= 0.00000001d0 .or. Qmu(i) <= 0.00000001d0) stop 'non-positive value of QKappa or Qmu'