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

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2012-03-18 22:39:18 -0700 (Sun, 18 Mar 2012)
New Revision: 19807

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_initialfield.F90
Log:
fixed a small bug in the initial plane wave in the case of MPI


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_initialfield.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_initialfield.F90	2012-03-19 05:32:33 UTC (rev 19806)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_initialfield.F90	2012-03-19 05:39:18 UTC (rev 19807)
@@ -128,13 +128,12 @@
     else
       call exit_MPI('Unrecognized source_type: should be 1 for plane P waves, 2 for plane SV waves, 3 for Rayleigh wave')
     endif
+  endif
 
   ! allow negative angleforce(1): incidence from the right side of the domain
     angleforce_abs=abs(angleforce(1))
-    if (angleforce_abs > pi/2.d0 .and. source_type(1) /= 3) then
+    if (angleforce_abs > pi/2.d0 .and. source_type(1) /= 3) &
       call exit_MPI("incorrect angleforce: must have 0 <= angleforce < 90")
-    endif
-  endif
 
   ! only implemented for homogeneous media therefore only 1 material supported
   numat_local = numat