[cig-commits] r13025 - seismo/2D/SPECFEM2D/branches/BIOT

[cmorency at geodynamics.org]

Author: cmorency
Date: 2008-10-14 11:01:00 -0700 (Tue, 14 Oct 2008)
New Revision: 13025

Modified:
   seismo/2D/SPECFEM2D/branches/BIOT/compute_forces_acoustic.f90
Log:
correction of adjoint source in compute_forces_acoustic.f90 


Modified: seismo/2D/SPECFEM2D/branches/BIOT/compute_forces_acoustic.f90
===================================================================
--- seismo/2D/SPECFEM2D/branches/BIOT/compute_forces_acoustic.f90	2008-10-14 17:47:17 UTC (rev 13024)
+++ seismo/2D/SPECFEM2D/branches/BIOT/compute_forces_acoustic.f90	2008-10-14 18:01:00 UTC (rev 13025)
@@ -514,16 +514,8 @@
       do j=1,NGLLZ
         do i=1,NGLLX
       iglob = ibool(i,j,ispec_selected_rec(irec))
-          xxi = + gammaz(i,j,ispec_selected_rec(irec)) * jacobian(i,j,ispec_selected_rec(irec))
-          zxi = - gammax(i,j,ispec_selected_rec(irec)) * jacobian(i,j,ispec_selected_rec(irec))
-          jacobian1D = sqrt(xxi**2 + zxi**2)
-          nx = - zxi / jacobian1D
-          nz = + xxi / jacobian1D
-
-          weight = jacobian1D * wxgll(i)
-
-      potential_dot_dot_acoustic(iglob) = potential_dot_dot_acoustic(iglob) + weight*&
-          (nx*adj_sourcearrays(irec_local,NSTEP-it+1,1,i,j) + nz*adj_sourcearrays(irec_local,NSTEP-it+1,2,i,j))
+      potential_dot_dot_acoustic(iglob) = potential_dot_dot_acoustic(iglob) + &
+          adj_sourcearrays(irec_local,NSTEP-it+1,1,i,j)
         enddo
       enddo
             endif