[cig-commits] r16035 - in seismo/3D/SPECFEM3D_GLOBE: tags/v4.0.4 trunk

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2009-11-25 06:25:45 -0800 (Wed, 25 Nov 2009)
New Revision: 16035

Modified:
   seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.4/auto_ner.f90
   seismo/3D/SPECFEM3D_GLOBE/trunk/auto_ner.f90
Log:
Brian Savage fixed a bug detected by Rosa Martin from University of Granada, Spain:
an array of size 1:9 instead of 1:8 was copied to a subarray of size 1:8


Modified: seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.4/auto_ner.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.4/auto_ner.f90	2009-11-25 01:53:44 UTC (rev 16034)
+++ seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.4/auto_ner.f90	2009-11-25 14:25:45 UTC (rev 16035)
@@ -35,14 +35,20 @@
 !  <savage13 at gps.caltech.edu>
 !  <savage13 at dtm.ciw.edu>
 !
-!   It is based partially upon formulation in the following references:
+!  It is based partially upon formulation in:
 !
-!   Komatitsch and Tromp, 2002 Part I
+! @ARTICLE{KoTr02a,
+! author={D. Komatitsch and J. Tromp},
+! year=2002,
+! title={Spectral-Element Simulations of Global Seismic Wave Propagation{-I. V}alidation},
+! journal={Geophys. J. Int.},
+! volume=149,
+! number=2,
+! pages={390-412},
+! doi={10.1046/j.1365-246X.2002.01653.x}}
 !
-!    and
+!  and the core determination was developed.
 !
-!   The Core determineation was developed
-!
 
   subroutine auto_time_stepping(WIDTH,  NEX_MAX, DT)
     implicit none
@@ -342,7 +348,7 @@
     integer npts, ispec
     integer, parameter :: NBNODE = 8
     double precision pts(NBNODE+1,2), points(npts,2)
-    pts(1:8,:) = points( ( (ispec-1) * NBNODE)+1 : ( (ispec) * NBNODE )+1, : )
+    pts(1:8,:) = points( ( (ispec-1) * NBNODE)+1 : ( (ispec) * NBNODE ), : )
     pts(NBNODE+1,:) = pts(1,:)  ! Use first point as the last point
   end subroutine get_element
 

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/auto_ner.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/auto_ner.f90	2009-11-25 01:53:44 UTC (rev 16034)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/auto_ner.f90	2009-11-25 14:25:45 UTC (rev 16035)
@@ -348,7 +348,7 @@
     integer npts, ispec
     integer, parameter :: NBNODE = 8
     double precision pts(NBNODE+1,2), points(npts,2)
-    pts(1:8,:) = points( ( (ispec-1) * NBNODE)+1 : ( (ispec) * NBNODE )+1, : )
+    pts(1:8,:) = points( ( (ispec-1) * NBNODE)+1 : ( (ispec) * NBNODE ), : )
     pts(NBNODE+1,:) = pts(1,:)  ! Use first point as the last point
   end subroutine get_element