[cig-commits] r21624 - seismo/3D/SPECFEM3D/trunk/src/decompose_mesh

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-03-24 15:10:15 -0700 (Sun, 24 Mar 2013)
New Revision: 21624

Modified:
   seismo/3D/SPECFEM3D/trunk/src/decompose_mesh/decompose_mesh.F90
Log:
added a stop statement if an acoustic material has a non-zero shear-wave velocity Vs


Modified: seismo/3D/SPECFEM3D/trunk/src/decompose_mesh/decompose_mesh.F90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/decompose_mesh/decompose_mesh.F90	2013-03-24 19:23:55 UTC (rev 21623)
+++ seismo/3D/SPECFEM3D/trunk/src/decompose_mesh/decompose_mesh.F90	2013-03-24 22:10:15 UTC (rev 21624)
@@ -361,6 +361,10 @@
        if(idomain_id == 1 .or. idomain_id == 2) then
          ! material is elastic or acoustic
 
+         ! check that the S-wave velocity is zero if the material is acoustic
+         if(idomain_id == 1 .and. vs >= 0.0001) &
+                stop 'acoustic material defined with a non-zero shear-wave velocity Vs, exiting...'
+
          mat_prop(1,num_mat) = rho
          mat_prop(2,num_mat) = vp
          mat_prop(3,num_mat) = vs