[cig-commits] r21975 - seismo/3D/SPECFEM3D/trunk/src/specfem3D

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2013-05-03 02:13:19 -0700 (Fri, 03 May 2013)
New Revision: 21975

Modified:
   seismo/3D/SPECFEM3D/trunk/src/specfem3D/compute_forces_acoustic_noDev.f90
   seismo/3D/SPECFEM3D/trunk/src/specfem3D/initialize_simulation.f90
Log:
remove unneeded comments and add one more correct condition for judge the place of outer boundary of pml


Modified: seismo/3D/SPECFEM3D/trunk/src/specfem3D/compute_forces_acoustic_noDev.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/specfem3D/compute_forces_acoustic_noDev.f90	2013-05-03 08:50:25 UTC (rev 21974)
+++ seismo/3D/SPECFEM3D/trunk/src/specfem3D/compute_forces_acoustic_noDev.f90	2013-05-03 09:13:19 UTC (rev 21975)
@@ -284,7 +284,7 @@
        do iface=1,num_abs_boundary_faces
            ispec = abs_boundary_ispec(iface)
            if (ispec_is_inner(ispec) .eqv. phase_is_inner) then
-              if( ispec_is_acoustic(ispec) ) then
+              if( ispec_is_acoustic(ispec) .and. is_CPML(ispec) ) then
                  ! reference gll points on boundary face
                  do igll = 1,NGLLSQUARE
 

Modified: seismo/3D/SPECFEM3D/trunk/src/specfem3D/initialize_simulation.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/specfem3D/initialize_simulation.f90	2013-05-03 08:50:25 UTC (rev 21974)
+++ seismo/3D/SPECFEM3D/trunk/src/specfem3D/initialize_simulation.f90	2013-05-03 09:13:19 UTC (rev 21975)
@@ -151,15 +151,15 @@
   if( ATTENUATION ) then
     !pll
     NSPEC_ATTENUATION_AB = NSPEC_AB
-    if(FULL_ATTENUATION_SOLID) then  !ZN
-      NSPEC_ATTENUATION_AB_kappa = NSPEC_AB  !ZN
-    else  !ZN
-      NSPEC_ATTENUATION_AB_kappa = 1  !ZN
-    endif  !ZN
+    if(FULL_ATTENUATION_SOLID) then  
+      NSPEC_ATTENUATION_AB_kappa = NSPEC_AB  
+    else  
+      NSPEC_ATTENUATION_AB_kappa = 1  
+    endif  
   else
     ! if attenuation is off, set dummy size of arrays to one
     NSPEC_ATTENUATION_AB = 1
-    NSPEC_ATTENUATION_AB_kappa = 1  !ZN
+    NSPEC_ATTENUATION_AB_kappa = 1  
   endif
 
   ! needed for attenuation and/or kernel computations