[cig-commits] [commit] devel: updated comments about Carcione 2007 (736465c)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem2d

On branch  : devel
Link       : https://github.com/geodynamics/specfem2d/compare/a7b59ef7f62c94aa042786bbe47e0d5c8c556b87...67d3f2483297e3de8343c24513ff97b57ef12a80

>---------------------------------------------------------------

commit 736465c53dbb4d84386e5eee5332fc3bd71cab39
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Thu May 22 16:45:16 2014 +0200

    updated comments about Carcione 2007


>---------------------------------------------------------------

736465c53dbb4d84386e5eee5332fc3bd71cab39
 src/specfem2D/compute_forces_poro_fluid.f90   | 2 +-
 src/specfem2D/compute_forces_poro_solid.f90   | 2 +-
 src/specfem2D/compute_forces_viscoelastic.F90 | 5 +----
 src/specfem2D/compute_pressure.f90            | 4 ++--
 4 files changed, 5 insertions(+), 8 deletions(-)

diff --git a/src/specfem2D/compute_forces_poro_fluid.f90 b/src/specfem2D/compute_forces_poro_fluid.f90
index cc50e33..a81cc01 100644
--- a/src/specfem2D/compute_forces_poro_fluid.f90
+++ b/src/specfem2D/compute_forces_poro_fluid.f90
@@ -435,7 +435,7 @@
 
 ! When implementing viscoelasticity according to Carcione 1993 paper, the attenuation is
 ! non-causal rather than causal. We fixed the problem by using equations in Carcione's
-! 2004 paper and his 2007 book.
+! 2004 paper and his 2007 book (Carcione 2007 page 125).
 
 !J. M. Carcione, H B. Helle, The physics and simulation of wave propagation at the ocean
 ! bottom, Geophysics, vol. 69(3), p. 825-839, 2004
diff --git a/src/specfem2D/compute_forces_poro_solid.f90 b/src/specfem2D/compute_forces_poro_solid.f90
index 9e92775..97de32b 100644
--- a/src/specfem2D/compute_forces_poro_solid.f90
+++ b/src/specfem2D/compute_forces_poro_solid.f90
@@ -438,7 +438,7 @@
 
 ! When implementing viscoelasticity according to Carcione 1993 paper, the attenuation is
 ! non-causal rather than causal. We fixed the problem by using equations in Carcione's
-! 2004 paper and his 2007 book.
+! 2004 paper and his 2007 book (Carcione 2007 page 125).
 
 !J. M. Carcione, H B. Helle, The physics and simulation of wave propagation at the ocean
 ! bottom, Geophysics, vol. 69(3), p. 825-839, 2004
diff --git a/src/specfem2D/compute_forces_viscoelastic.F90 b/src/specfem2D/compute_forces_viscoelastic.F90
index 823497a..584a7b1 100644
--- a/src/specfem2D/compute_forces_viscoelastic.F90
+++ b/src/specfem2D/compute_forces_viscoelastic.F90
@@ -746,6 +746,7 @@ subroutine compute_forces_viscoelastic(p_sv,nglob,nspec,myrank,nelemabs,numat, &
             ! J. M. Carcione, D. Kosloff and R. Kosloff, Wave propagation simulation in a linear
             ! viscoelastic medium, Geophysical Journal International, vol. 95, p. 597-611 (1988).
 
+            ! Compute unrelaxed elastic coefficients from formulas in Carcione 2007 page 125.
             ! When implementing viscoelasticity according to the Carcione 1993 paper, attenuation is
             ! non-causal rather than causal. We fixed the problem by using equations in Carcione's
             ! 2004 paper and his 2007 book. See also file doc/problem_attenuation_reference_Specfem2D_fixed_by_Xie_Zhinan.pdf
@@ -755,10 +756,6 @@ subroutine compute_forces_viscoelastic(p_sv,nglob,nspec,myrank,nelemabs,numat, &
             ! J. M. Carcione, Wave fields in real media: wave propagation in anisotropic, anelastic
             ! and porous media, Elsevier, p. 124-125, 2007
 
-            ! compute unrelaxed elastic coefficients from formulas in Carcione 2007 page 125.
-            ! When implementing viscoelasticity according to the Carcione 1993 paper, attenuation is
-            ! non-causal rather than causal. We fixed the problem by using equations in Carcione's
-            ! 2004 paper and his 2007 book. See also file doc/problem_attenuation_reference_Specfem2D_fixed_by_Xie_Zhinan.pdf
             lambdal_relaxed_viscoelastic = (lambdal_unrelaxed_elastic + 2._CUSTOM_REAL*mul_unrelaxed_elastic/3._CUSTOM_REAL)&
                                            / Mu_nu1(i,j,ispec) &
                                            - (2._CUSTOM_REAL*mul_unrelaxed_elastic/3._CUSTOM_REAL) / Mu_nu2(i,j,ispec)
diff --git a/src/specfem2D/compute_pressure.f90 b/src/specfem2D/compute_pressure.f90
index b70dfa5..669038e 100644
--- a/src/specfem2D/compute_pressure.f90
+++ b/src/specfem2D/compute_pressure.f90
@@ -331,7 +331,7 @@
 
 ! When implementing viscoelasticity according to Carcione 1993 paper, the attenuation is
 ! non-causal rather than causal. We fixed the problem by using equations in Carcione's
-! 2004 paper and his 2007 book.
+! 2004 paper and his 2007 book (Carcione 2007 page 125).
 
 ! J. M. Carcione, H. B. Helle, The physics and simulation of wave propagation at the ocean
 !  bottom, Geophysics, vol. 69(3), p. 825-839, 2004
@@ -547,7 +547,7 @@
 
 ! When implement viscoelasticity according to Carcione 1993 paper, the attenuation is
 ! non-causal rather than causal. We fixed the problem by using equations in Carcione's
-! 2004 paper and his 2007 book.
+! 2004 paper and his 2007 book (Carcione 2007 page 125).
 
 ! J. M. Carcione, H B. Helle, The physics and simulation of wave propagation at the ocean
 !  bottom, Geophysics, vol. 69(3), p. 825-839, 2004