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

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2012-09-27 08:47:52 -0700 (Thu, 27 Sep 2012)
New Revision: 20789

Modified:
   seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90
Log:
when initial_field and add_Bielak_conditions is on, ADD_SPRING_TO_STACEY should be set to false. Because currently the code to implement initial_field and add_Bielak_conditions is not correct for  ADD_SPRING_TO_STACEY is true


Modified: seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90	2012-09-27 15:18:42 UTC (rev 20788)
+++ seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90	2012-09-27 15:47:52 UTC (rev 20789)
@@ -473,6 +473,8 @@
   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'
 
+  if(add_Bielak_conditions .or. initialfield ) ADD_SPRING_TO_STACEY = .false.
+
   call read_value_logical_p(ADD_PERIODIC_CONDITIONS, 'solver.ADD_PERIODIC_CONDITIONS')
   if(err_occurred() /= 0) stop 'error reading parameter 51b in Par_file'