[cig-commits] r12549 - in seismo/3D/SPECFEM3D_GLOBE/tags: . v4.0.2
leif at geodynamics.org
leif at geodynamics.org
Wed Aug 6 12:21:04 PDT 2008
Author: leif
Date: 2008-08-06 12:21:04 -0700 (Wed, 06 Aug 2008)
New Revision: 12549
Added:
seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.2/
Modified:
seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.2/auto_ner.f90
seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.2/compute_forces_crust_mantle.f90
Log:
Created SPECFEM3D GLOBE v4.0.2 tag, which incorporates the attenuation
fix described by Vala.
Copied: seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.2 (from rev 12548, seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.1)
Modified: seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.2/auto_ner.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.1/auto_ner.f90 2008-08-06 16:12:54 UTC (rev 12548)
+++ seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.2/auto_ner.f90 2008-08-06 19:21:04 UTC (rev 12549)
@@ -145,7 +145,7 @@
double precision, dimension(NUM_REGIONS-1) :: ratio_top
double precision, dimension(NUM_REGIONS-1) :: ratio_bottom
integer, dimension(NUM_REGIONS-1) :: NER
- double precision NEX_ETA
+ integer NEX_ETA
! This is PREM in Kilometers, well ... kinda, not really ....
radius(1) = 6371.00d0 ! Surface
Modified: seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.2/compute_forces_crust_mantle.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.1/compute_forces_crust_mantle.f90 2008-08-06 16:12:54 UTC (rev 12548)
+++ seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.2/compute_forces_crust_mantle.f90 2008-08-06 19:21:04 UTC (rev 12549)
@@ -288,7 +288,7 @@
if(ATTENUATION_3D_VAL) then
one_minus_sum_beta_use = one_minus_sum_beta(i,j,k,ispec)
else
- radius_cr = xstore(iglob)
+ radius_cr = xstore(ibool(i,j,k,ispec))
call get_attenuation_index(idoubling(ispec), dble(radius_cr), iregion_selected, .FALSE., AM_V)
one_minus_sum_beta_use = one_minus_sum_beta(1,1,1,iregion_selected)
endif
More information about the cig-commits
mailing list