[cig-commits] r19153 - seismo/2D/SPECFEM2D/trunk/src/specfem2D

[liuqy at geodynamics.org]

Author: liuqy
Date: 2011-11-06 14:16:51 -0800 (Sun, 06 Nov 2011)
New Revision: 19153

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.f90
Log:
Fix bug related to back-reconstruct plane wave field (SIMULATION_TYPE=2)



Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.f90	2011-11-06 02:02:36 UTC (rev 19152)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.f90	2011-11-06 22:16:51 UTC (rev 19153)
@@ -543,8 +543,8 @@
 
                  if(SAVE_FORWARD .and. SIMULATION_TYPE ==1) then
                     if(p_sv)then !P-SV waves
-                       b_absorb_elastic_left(1,j,ib_left(ispecabs),it) = tx*weight
-                       b_absorb_elastic_left(3,j,ib_left(ispecabs),it) = tz*weight
+                       b_absorb_elastic_left(1,j,ib_left(ispecabs),it) = (tx + traction_x_t0)*weight
+                       b_absorb_elastic_left(3,j,ib_left(ispecabs),it) = (tz + traction_z_t0)*weight
                     else !SH (membrane) waves
                        b_absorb_elastic_left(2,j,ib_left(ispecabs),it) = ty*weight
                     endif
@@ -634,8 +634,8 @@
 
                  if(SAVE_FORWARD .and. SIMULATION_TYPE ==1) then
                     if(p_sv)then !P-SV waves
-                       b_absorb_elastic_right(1,j,ib_right(ispecabs),it) = tx*weight
-                       b_absorb_elastic_right(3,j,ib_right(ispecabs),it) = tz*weight
+                       b_absorb_elastic_right(1,j,ib_right(ispecabs),it) = (tx - traction_x_t0)*weight
+                       b_absorb_elastic_right(3,j,ib_right(ispecabs),it) = (tz - traction_z_t0)*weight
                     else! SH (membrane) waves
                        b_absorb_elastic_right(2,j,ib_right(ispecabs),it) = ty*weight
                     endif
@@ -731,8 +731,8 @@
 
                  if(SAVE_FORWARD .and. SIMULATION_TYPE ==1) then
                     if(p_sv)then !P-SV waves
-                       b_absorb_elastic_bottom(1,i,ib_bottom(ispecabs),it) = tx*weight
-                       b_absorb_elastic_bottom(3,i,ib_bottom(ispecabs),it) = tz*weight
+                       b_absorb_elastic_bottom(1,i,ib_bottom(ispecabs),it) = (tx + traction_x_t0)*weight
+                       b_absorb_elastic_bottom(3,i,ib_bottom(ispecabs),it) = (tz + traction_z_t0)*weight
                     else!SH (membrane) waves
                        b_absorb_elastic_bottom(2,i,ib_bottom(ispecabs),it) = ty*weight
                     endif
@@ -820,8 +820,8 @@
 
                  if(SAVE_FORWARD .and. SIMULATION_TYPE ==1) then
                     if(p_sv)then !P-SV waves
-                       b_absorb_elastic_top(1,i,ib_top(ispecabs),it) = tx*weight
-                       b_absorb_elastic_top(3,i,ib_top(ispecabs),it) = tz*weight
+                       b_absorb_elastic_top(1,i,ib_top(ispecabs),it) = (tx- traction_x_t0)*weight
+                       b_absorb_elastic_top(3,i,ib_top(ispecabs),it) = (tz- traction_z_t0)*weight
                     else!SH (membrane) waves
                        b_absorb_elastic_top(2,i,ib_top(ispecabs),it) = ty*weight
                     endif