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

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2012-09-01 16:34:19 -0700 (Sat, 01 Sep 2012)
New Revision: 20664

Modified:
   seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90
Log:
added a stop statement if add_Bielak_conditions and not STACEY_ABSORBING_CONDITIONS


Modified: seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90	2012-09-01 23:26:18 UTC (rev 20663)
+++ seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90	2012-09-01 23:34:19 UTC (rev 20664)
@@ -463,8 +463,10 @@
   call read_value_logical_p(any_abs, 'solver.STACEY_ABSORBING_CONDITIONS')
   if(err_occurred() /= 0) stop 'error reading parameter 51a in Par_file'
 
- !DK DK solve the conflict in set the value of PML_BOUNDARY_CONDITIONS and STACEY_ABSORBING_CONDITIONS
- !in par_file
+  if(add_Bielak_conditions .and. .not. STACEY_ABSORBING_CONDITIONS) &
+    stop 'need STACEY_ABSORBING_CONDITIONS set to .true. in order to use add_Bielak_conditions'
+
+  ! solve the conflict in value of PML_BOUNDARY_CONDITIONS and STACEY_ABSORBING_CONDITIONS
   if(PML_BOUNDARY_CONDITIONS) any_abs = .true.
 
   call read_value_logical_p(ADD_SPRING_TO_STACEY, 'solver.ADD_SPRING_TO_STACEY')