[cig-commits] r21988 - seismo/2D/SPECFEM2D/trunk/src/meshfem2D

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-05-06 10:39:39 -0700 (Mon, 06 May 2013)
New Revision: 21988

Modified:
   seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90
Log:
added a stop statement if Stacey and PML are both set


Modified: seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90	2013-05-06 17:36:58 UTC (rev 21987)
+++ seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90	2013-05-06 17:39:39 UTC (rev 21988)
@@ -482,6 +482,10 @@
   ! solve the conflict in value of PML_BOUNDARY_CONDITIONS and STACEY_ABSORBING_CONDITIONS
   if(PML_BOUNDARY_CONDITIONS) any_abs = .true.
 
+  ! CPML and Stacey are mutually exclusive
+  if(STACEY_ABSORBING_CONDITIONS .and. PML_BOUNDARY_CONDITIONS) &
+    stop 'STACEY_ABSORBING_CONDITIONS and PML_BOUNDARY_CONDITIONS are mutually exclusive but are both set to .true.'
+
   call read_value_logical_p(ADD_SPRING_TO_STACEY, 'solver.ADD_SPRING_TO_STACEY')
   if(err_occurred() /= 0) stop 'error reading parameter 51a in Par_file'