[cig-commits] r11919 - seismo/3D/SPECFEM3D_GLOBE/trunk

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2008-05-06 13:53:20 -0700 (Tue, 06 May 2008)
New Revision: 11919

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/compute_forces_crust_mantle.f90
   seismo/3D/SPECFEM3D_GLOBE/trunk/compute_forces_inner_core.f90
Log:
added comment about attenuation for P waves (Q_\alpha)


Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/compute_forces_crust_mantle.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/compute_forces_crust_mantle.f90	2008-05-06 14:40:30 UTC (rev 11918)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/compute_forces_crust_mantle.f90	2008-05-06 20:53:20 UTC (rev 11919)
@@ -803,6 +803,12 @@
 ! term in xz = 4
 ! term in yz = 5
 ! term in zz not computed since zero trace
+! This is because we only implement Q_\mu attenuation and not Q_\Kappa.
+! Note that this does *NOT* imply that there is no attenuation for P waves
+! because for Q_\Kappa = infinity one gets (see for instance Dahlen and Tromp (1998)
+! equation (9.59) page 350): Q_\alpha = Q_\mu * 3 * (V_p/V_s)^2 / 4
+! therefore Q_\alpha is not zero; for instance for V_p / V_s = sqrt(3)
+! we get Q_\alpha = (9 / 4) * Q_\mu = 2.25 * Q_\mu
 
   if(ATTENUATION_VAL) then
 

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/compute_forces_inner_core.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/compute_forces_inner_core.f90	2008-05-06 14:40:30 UTC (rev 11918)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/compute_forces_inner_core.f90	2008-05-06 20:53:20 UTC (rev 11919)
@@ -562,6 +562,12 @@
 ! term in xz = 4
 ! term in yz = 5
 ! term in zz not computed since zero trace
+! This is because we only implement Q_\mu attenuation and not Q_\Kappa.
+! Note that this does *NOT* imply that there is no attenuation for P waves
+! because for Q_\Kappa = infinity one gets (see for instance Dahlen and Tromp (1998)
+! equation (9.59) page 350): Q_\alpha = Q_\mu * 3 * (V_p/V_s)^2 / 4
+! therefore Q_\alpha is not zero; for instance for V_p / V_s = sqrt(3)
+! we get Q_\alpha = (9 / 4) * Q_\mu = 2.25 * Q_\mu
 
     if(ATTENUATION_VAL) then