[cig-commits] r21931 - seismo/3D/SPECFEM3D/trunk/src/specfem3D

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2013-04-25 01:34:10 -0700 (Thu, 25 Apr 2013)
New Revision: 21931

Modified:
   seismo/3D/SPECFEM3D/trunk/src/specfem3D/compute_forces_acoustic_noDev.f90
Log:
fix error in computing temp1l_new,temp2l_new,temp3l_new


Modified: seismo/3D/SPECFEM3D/trunk/src/specfem3D/compute_forces_acoustic_noDev.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/specfem3D/compute_forces_acoustic_noDev.f90	2013-04-25 00:15:29 UTC (rev 21930)
+++ seismo/3D/SPECFEM3D/trunk/src/specfem3D/compute_forces_acoustic_noDev.f90	2013-04-25 08:34:10 UTC (rev 21931)
@@ -151,19 +151,19 @@
                 temp3l_new = temp3l
 
                 do l=1,NGLLX
-                   hp1 = hprime_xx(l,i)
+                   hp1 = hprime_xx(i,l)
                    iglob = ibool(l,j,k,ispec)
                    temp1l_new = temp1l_new + deltat*potential_dot_acoustic(iglob)*hp1
 !!! can merge these loops because NGLLX = NGLLY = NGLLZ          enddo
 
 !!! can merge these loops because NGLLX = NGLLY = NGLLZ          do l=1,NGLLY
-                   hp2 = hprime_yy(l,j)
+                   hp2 = hprime_yy(j,l)
                    iglob = ibool(i,l,k,ispec)
                    temp2l_new = temp2l_new + deltat*potential_dot_acoustic(iglob)*hp2
 !!! can merge these loops because NGLLX = NGLLY = NGLLZ          enddo
 
 !!! can merge these loops because NGLLX = NGLLY = NGLLZ          do l=1,NGLLZ
-                   hp3 = hprime_zz(l,k)
+                   hp3 = hprime_zz(k,l)
                    iglob = ibool(i,j,l,ispec)
                    temp3l_new = temp3l_new + deltat*potential_dot_acoustic(iglob)*hp3
                 enddo