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

[cmorency at geodynamics.org]

Author: cmorency
Date: 2009-04-02 07:39:01 -0700 (Thu, 02 Apr 2009)
New Revision: 14567

Modified:
   seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90
Log:
BIOT2D: corrected specfem2D.F90 to account for free surface enforcement for the backward acoustic wavefield.


Modified: seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90	2009-04-02 08:35:28 UTC (rev 14566)
+++ seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90	2009-04-02 14:39:01 UTC (rev 14567)
@@ -3646,6 +3646,11 @@
       call enforce_acoustic_free_surface(potential_dot_dot_acoustic,potential_dot_acoustic, &
            potential_acoustic,acoustic_surface, &
            ibool,nelem_acoustic_surface,npoin,nspec)
+   if(isolver == 2) then
+    call enforce_acoustic_free_surface(b_potential_dot_dot_acoustic,b_potential_dot_acoustic, &
+                b_potential_acoustic,acoustic_surface, &
+                ibool,nelem_acoustic_surface,npoin,nspec)
+   endif
       endif
 
 ! *********************************************************
@@ -3964,6 +3969,11 @@
     call enforce_acoustic_free_surface(potential_dot_dot_acoustic,potential_dot_acoustic, &
                 potential_acoustic,acoustic_surface, &
                 ibool,nelem_acoustic_surface,npoin,nspec)
+   if(isolver == 2) then
+    call enforce_acoustic_free_surface(b_potential_dot_dot_acoustic,b_potential_dot_acoustic, &
+                b_potential_acoustic,acoustic_surface, &
+                ibool,nelem_acoustic_surface,npoin,nspec)
+   endif
     endif
   endif