[cig-commits] r15667 - seismo/2D/SPECFEM2D/trunk

[cmorency at geodynamics.org]

Author: cmorency
Date: 2009-09-14 14:59:09 -0700 (Mon, 14 Sep 2009)
New Revision: 15667

Modified:
   seismo/2D/SPECFEM2D/trunk/specfem2D.F90
Log:
Added SH (membrane) waves calculation (only for elastic media), which yields to a new line in Par_file to specify which kind of calculation you are running: P-SV or SH. iRead the REAME file for more informations.
Also fixed some bugs related to acoustic kernels.


Modified: seismo/2D/SPECFEM2D/trunk/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2009-09-14 21:56:41 UTC (rev 15666)
+++ seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2009-09-14 21:59:09 UTC (rev 15667)
@@ -5103,6 +5103,11 @@
     potential_dot_dot_acoustic = potential_dot_dot_acoustic * rmass_inverse_acoustic
     potential_dot_acoustic = potential_dot_acoustic + deltatover2*potential_dot_dot_acoustic
 
+    if(isolver ==2)then
+    b_potential_dot_dot_acoustic = b_potential_dot_dot_acoustic * rmass_inverse_acoustic
+    b_potential_dot_acoustic = b_potential_dot_acoustic + b_deltatover2*b_potential_dot_dot_acoustic
+    endif
+
 ! free surface for an acoustic medium
     if ( nelem_acoustic_surface > 0 ) then
     call enforce_acoustic_free_surface(potential_dot_dot_acoustic,potential_dot_acoustic, &