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

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2013-03-13 11:07:16 -0700 (Wed, 13 Mar 2013)
New Revision: 21523

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90
Log:
add PML support for viscoelastic simulation. Here the PML is not derived based on viscoelastic wave equation, but on elastic wave equation. However,when approximately Qp is bigger than 70 and Qs is bigger than 50, based on the experiment observation,we found that the PML can still keep good absorbing efficiency when use it directly in viscoelastic simulation


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90	2013-03-13 18:03:53 UTC (rev 21522)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90	2013-03-13 18:07:16 UTC (rev 21523)
@@ -260,6 +260,9 @@
      ! loop over spectral elements
      do ispec = 1,nspec
 
+        if((.not. PML_BOUNDARY_CONDITIONS) .or. &
+          (PML_BOUNDARY_CONDITIONS .and. (.not. is_PML(ispec))))then
+
         do j=1,NGLLZ
            do i=1,NGLLX
 
@@ -394,6 +397,7 @@
 
            enddo
         enddo
+       endif
      enddo
 
   endif ! end of test on attenuation