[cig-commits] r21864 - seismo/2D/SPECFEM2D/trunk/src/specfem2D

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-04-13 15:18:57 -0700 (Sat, 13 Apr 2013)
New Revision: 21864

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_energy.f90
Log:
remove all the "mu" terms from the expression of energy in an acoustic medium (since mu = 0 in such a medium, these terms are useless because they are multiplied by zero)


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_energy.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_energy.f90	2013-04-13 22:14:18 UTC (rev 21863)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_energy.f90	2013-04-13 22:18:57 UTC (rev 21864)
@@ -358,10 +358,9 @@
 
       ! get density of current spectral element
       lambdal_unrelaxed_elastic = poroelastcoef(1,1,kmato(ispec))
-      mul_unrelaxed_elastic = poroelastcoef(2,1,kmato(ispec))
       rhol  = density(1,kmato(ispec))
-      kappal  = lambdal_unrelaxed_elastic + TWO*mul_unrelaxed_elastic/3._CUSTOM_REAL
-      cpl = sqrt((kappal + 4._CUSTOM_REAL*mul_unrelaxed_elastic/3._CUSTOM_REAL)/rhol)
+      kappal  = lambdal_unrelaxed_elastic
+      cpl = sqrt(kappal/rhol)
 
       ! double loop over GLL points
       do j = 1,NGLLZ