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

[cmorency at geodynamics.org]

Author: cmorency
Date: 2009-04-01 14:24:22 -0700 (Wed, 01 Apr 2009)
New Revision: 14561

Modified:
   seismo/2D/SPECFEM2D/branches/BIOT/compute_forces_acoustic.f90
Log:
BIOT2D: corrected bug on bottom and top forces saved to reconstruct the forward wavefield for adjoint calculations.


Modified: seismo/2D/SPECFEM2D/branches/BIOT/compute_forces_acoustic.f90
===================================================================
--- seismo/2D/SPECFEM2D/branches/BIOT/compute_forces_acoustic.f90	2009-04-01 20:32:06 UTC (rev 14560)
+++ seismo/2D/SPECFEM2D/branches/BIOT/compute_forces_acoustic.f90	2009-04-01 21:24:22 UTC (rev 14561)
@@ -400,10 +400,10 @@
             potential_dot_dot_acoustic(iglob) = potential_dot_dot_acoustic(iglob) - potential_dot_acoustic(iglob)*weight/cpl
 
              if(save_forward .and. isolver ==1) then
-            b_absorb_acoustic_bottom(j,ispecabs,it) = potential_dot_acoustic(iglob)*weight/cpl
+            b_absorb_acoustic_bottom(i,ispecabs,it) = potential_dot_acoustic(iglob)*weight/cpl
              elseif(isolver == 2) then
             b_potential_dot_dot_acoustic(iglob) = b_potential_dot_dot_acoustic(iglob) - &
-                                               b_absorb_acoustic_bottom(j,ispecabs,NSTEP-it+1)
+                                               b_absorb_acoustic_bottom(i,ispecabs,NSTEP-it+1)
              endif
 
           endif
@@ -463,9 +463,9 @@
             potential_dot_dot_acoustic(iglob) = potential_dot_dot_acoustic(iglob) - potential_dot_acoustic(iglob)*weight/cpl
 
              if(save_forward .and. isolver ==1) then
-            b_absorb_acoustic_top(j,ispecabs,it) = potential_dot_acoustic(iglob)*weight/cpl
+            b_absorb_acoustic_top(i,ispecabs,it) = potential_dot_acoustic(iglob)*weight/cpl
              elseif(isolver == 2) then
-            b_potential_dot_dot_acoustic(iglob) = b_potential_dot_dot_acoustic(iglob) - b_absorb_acoustic_top(j,ispecabs,NSTEP-it+1)
+            b_potential_dot_dot_acoustic(iglob) = b_potential_dot_dot_acoustic(iglob) - b_absorb_acoustic_top(i,ispecabs,NSTEP-it+1)
              endif
 
           endif