[cig-commits] [commit] pluggable: Created SPECFEM3D GLOBE v4.0.2 tag, which incorporates the attenuation fix described by Vala. (fac8c40)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/specfem3d_globe

On branch  : pluggable
Link       : https://github.com/geodynamics/specfem3d_globe/compare/64e1b38f0c5ebb4056cce0b15d41c0b9f94ab6e5...099a4d330d5b173b21e51ad441f9f429e5d37842

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

commit fac8c40dff2d33bc4ae676e1be49d7f4517e8952
Author: Leif Strand <leif at geodynamics.org>
Date:   Wed Aug 6 19:21:04 2008 +0000

    Created SPECFEM3D GLOBE v4.0.2 tag, which incorporates the attenuation
    fix described by Vala.


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

fac8c40dff2d33bc4ae676e1be49d7f4517e8952
 auto_ner.f90                    | 2 +-
 compute_forces_crust_mantle.f90 | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/auto_ner.f90 b/auto_ner.f90
index 67482a4..9b7c5a8 100644
--- a/auto_ner.f90
+++ b/auto_ner.f90
@@ -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
diff --git a/compute_forces_crust_mantle.f90 b/compute_forces_crust_mantle.f90
index 5040c4e..9764dcf 100644
--- a/compute_forces_crust_mantle.f90
+++ b/compute_forces_crust_mantle.f90
@@ -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