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

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2012-06-28 09:52:31 -0700 (Thu, 28 Jun 2012)
New Revision: 20429

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/invert_mass_matrix.F90
Log:
fixed two bugs introduced by Zhinan Xie in the C-PML code: in some cases, variable anyabs_local was used without being assigned a value first


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90	2012-06-28 09:30:30 UTC (rev 20428)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90	2012-06-28 16:52:31 UTC (rev 20429)
@@ -248,6 +248,8 @@
   sigma_zz = 0
   sigma_zx = 0
 
+  anyabs_local=.false. ! to avoid having an uninitialized variable
+
   if( PML_BOUNDARY_CONDITIONS ) then
     accel_elastic_PML = 0._CUSTOM_REAL
     accel_elastic_PML_corner = 0._CUSTOM_REAL

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/invert_mass_matrix.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/invert_mass_matrix.F90	2012-06-28 09:30:30 UTC (rev 20428)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/invert_mass_matrix.F90	2012-06-28 16:52:31 UTC (rev 20429)
@@ -141,6 +141,8 @@
   if(any_poroelastic) rmass_s_inverse_poroelastic(:) = 0._CUSTOM_REAL
   if(any_poroelastic) rmass_w_inverse_poroelastic(:) = 0._CUSTOM_REAL
   if(any_acoustic) rmass_inverse_acoustic(:) = 0._CUSTOM_REAL
+
+  anyabs_local=.false. ! to avoid having an uninitialized variable
   if(PML_BOUNDARY_CONDITIONS)anyabs_local=.false.
   if(.not.PML_BOUNDARY_CONDITIONS .and. anyabs)anyabs_local=.true.