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

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-05-06 10:33:32 -0700 (Mon, 06 May 2013)
New Revision: 21986

Modified:
   seismo/3D/SPECFEM3D/trunk/src/specfem3D/initialize_simulation.f90
Log:
added a stop statement if Stacey and PML are both set


Modified: seismo/3D/SPECFEM3D/trunk/src/specfem3D/initialize_simulation.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/specfem3D/initialize_simulation.f90	2013-05-06 17:28:58 UTC (rev 21985)
+++ seismo/3D/SPECFEM3D/trunk/src/specfem3D/initialize_simulation.f90	2013-05-06 17:33:32 UTC (rev 21986)
@@ -277,6 +277,9 @@
   if( .not. GPU_MODE .and. GRAVITY ) &
     stop 'GRAVITY only supported in GPU mode'
 
+  if(STACEY_ABSORBING_CONDITIONS .and. PML_CONDITIONS) &
+    stop 'STACEY_ABSORBING_CONDITIONS and PML_CONDITIONS are mutually exclusive but are both set to .true.'
+
   ! absorbing surfaces
   if( STACEY_ABSORBING_CONDITIONS ) then
      ! for arbitrary orientation of elements, which face belongs to xmin,xmax,etc... -
@@ -285,26 +288,23 @@
      ! just to be sure for now..
      if( NGLLX /= NGLLY .and. NGLLY /= NGLLZ ) &
           stop 'STACEY_ABSORBING_CONDITIONS must have NGLLX = NGLLY = NGLLZ'
-     if( PML_CONDITIONS ) then
-        print*, 'please modify Par_file and recompile solver'
-        stop 'STACEY_ABSORBING_CONDITIONS and PML_CONDITIONS are both set to .true.'
-     else if( PML_INSTEAD_OF_FREE_SURFACE ) then
-        print*, 'please modify Par_file and recompile solver'
+     if( PML_INSTEAD_OF_FREE_SURFACE ) then
+        print*, 'please modify Par_file'
         stop 'PML_INSTEAD_OF_FREE_SURFACE = .true. is incompatible with STACEY_ABSORBING_CONDITIONS = .true.'
      endif
   else
      if( STACEY_INSTEAD_OF_FREE_SURFACE ) then
-        print*, 'please modify Par_file and recompile solver'
+        print*, 'please modify Par_file'
         stop 'STACEY_ABSORBING_CONDITIONS must be activated when STACEY_INSTEAD_OF_FREE_SURFACE is set to .true.'
      endif
   endif
 
   if( PML_CONDITIONS ) then
      if( STACEY_INSTEAD_OF_FREE_SURFACE ) then
-        print*, 'please modify Par_file and recompile solver'
+        print*, 'please modify Par_file'
         stop 'STACEY_INSTEAD_OF_FREE_SURFACE = .true. is incompatible with PML_CONDITIONS = .true.'
      else if( .not. SUPPRESS_UTM_PROJECTION ) then
-        print*, 'please modify Par_file and recompile solver'
+        print*, 'please modify Par_file'
         stop 'SUPPRESS_UTM_PROJECTION must be activated when PML_CONDITIONS is set to .true.'
      endif
   else
@@ -313,13 +313,13 @@
   endif
 
   if( STACEY_INSTEAD_OF_FREE_SURFACE .and. PML_INSTEAD_OF_FREE_SURFACE ) then
-     print*, 'please modify Par_file and recompile solver'
+     print*, 'please modify Par_file'
      stop 'error: STACEY_INSTEAD_OF_FREE_SURFACE and PML_INSTEAD_OF_FREE_SURFACE are both set to .true.'
   endif
 
   ! checks the MOVIE_TYPE parameter
   if( MOVIE_TYPE /= 1 .and. MOVIE_TYPE /= 2 ) then
-     stop 'error: MOVIE_TYPE must be either 1 or 2! Please modify Par_file and recompile solver'
+     stop 'error: MOVIE_TYPE must be either 1 or 2! Please modify Par_file'
   endif
 
   ! check that the code has been compiled with the right values