[cig-commits] r20205 - seismo/2D/SPECFEM2D/trunk/src/specfem2D
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Mon May 21 18:32:09 PDT 2012
Author: dkomati1
Date: 2012-05-21 18:32:09 -0700 (Mon, 21 May 2012)
New Revision: 20205
Modified:
seismo/2D/SPECFEM2D/trunk/src/specfem2D/define_external_model.f90
Log:
fixed the Q = 0 bug in define_external_model.f90
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/define_external_model.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/define_external_model.f90 2012-05-20 23:08:53 UTC (rev 20204)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/define_external_model.f90 2012-05-22 01:32:09 UTC (rev 20205)
@@ -43,8 +43,8 @@
!========================================================================
- subroutine define_external_model(x,y,iflag_element,rho,vp,vs,QKappa_attenuation,&
- Qmu_attenuation,c11,c13,c15,c33,c35,c55 )
+ subroutine define_external_model(x,y,iflag_element,rho,vp,vs,QKappa_attenuation, &
+ Qmu_attenuation,c11,c13,c15,c33,c35,c55)
implicit none
@@ -66,8 +66,8 @@
rho = 2000.d0
vp = 3000.d0
vs = vp / sqrt(3.d0)
- QKappa_attenuation = 0
- Qmu_attenuation = 0
+ QKappa_attenuation = 9999. ! this means no attenuation
+ Qmu_attenuation = 9999. ! this means no attenuation
c11 = 169.d9
c13 = 122.d9
c15 = 0.d0
@@ -78,9 +78,9 @@
rho = 2500.d0
vp = 3600.d0
vs = vp / 2.d0
- QKappa_attenuation = 60
- Qmu_attenuation = 60
- c11 = 0.d0
+ QKappa_attenuation = 60.
+ Qmu_attenuation = 60.
+ c11 = 0.d0 ! this means no anisotropy
c13 = 0.d0
c15 = 0.d0
c33 = 0.d0
More information about the CIG-COMMITS
mailing list