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

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2012-01-12 15:27:50 -0800 (Thu, 12 Jan 2012)
New Revision: 19350

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_solid.f90
Log:
added comments by Christina Morency, following a discussion with Xie Zhinan


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90	2012-01-11 22:32:26 UTC (rev 19349)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90	2012-01-12 23:27:50 UTC (rev 19350)
@@ -344,6 +344,10 @@
       e13_sum = e13_sum + e13(i,j,ispec,i_sls)
     enddo
 
+! mul_G is the relaxed modulus. Note that it is defined as the
+! frame modulus (in compute_forces_poro_solid.f90), which Christina Morency noted
+! mul_fr, which is in her case equivalent to the solid phase shear
+! modulus, and whose value is entered in Par_file for example
     sigma_xx = sigma_xx + TWO * mul_G * e11_sum
     sigma_xz = sigma_xz + mul_G * e13_sum
     sigma_zz = sigma_zz - TWO * mul_G * e11_sum

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_solid.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_solid.f90	2012-01-11 22:32:26 UTC (rev 19349)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_solid.f90	2012-01-12 23:27:50 UTC (rev 19350)
@@ -345,6 +345,10 @@
       e13_sum = e13_sum + e13(i,j,ispec,i_sls)
     enddo
 
+! mul_G is the relaxed modulus. Note that it is defined as the
+! frame modulus (in compute_forces_poro_solid.f90), which Christina Morency noted
+! mul_fr, which is in her case equivalent to the solid phase shear
+! modulus, and whose value is entered in Par_file for example
     sigma_xx = sigma_xx + TWO * mul_G * e11_sum
     sigma_xz = sigma_xz + mul_G * e13_sum
     sigma_zz = sigma_zz - TWO * mul_G * e11_sum