[cig-commits] r12912 - seismo/2D/SPECFEM2D/branches/BIOT

[cmorency at geodynamics.org]

Author: cmorency
Date: 2008-09-17 12:17:36 -0700 (Wed, 17 Sep 2008)
New Revision: 12912

Modified:
   seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90
Log:
Bug fixed for acoustic rmass


Modified: seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90	2008-09-17 19:00:29 UTC (rev 12911)
+++ seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90	2008-09-17 19:17:36 UTC (rev 12912)
@@ -1594,9 +1594,8 @@
           cpsquare = vpext(i,j,ispec)**2
         else
           rhol = density(2,kmato(ispec))
-          lambdal_relaxed = poroelastcoef(1,2,kmato(ispec))
-          mul_relaxed = poroelastcoef(2,2,kmato(ispec))
-          cpsquare = (lambdal_relaxed + 2._CUSTOM_REAL*mul_relaxed) / rhol
+          kappal = poroelastcoef(1,2,kmato(ispec))
+          cpsquare = kappal / rhol
         endif
           rmass_inverse_acoustic(iglob) = rmass_inverse_acoustic(iglob) + wxgll(i)*wzgll(j)*jacobian(i,j,ispec) / cpsquare
         endif