[cig-commits] r14934 - seismo/2D/SPECFEM2D/trunk

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2009-05-08 16:58:03 -0700 (Fri, 08 May 2009)
New Revision: 14934

Modified:
   seismo/2D/SPECFEM2D/trunk/assemble_MPI.F90
   seismo/2D/SPECFEM2D/trunk/specfem2D.F90
Log:
bug detected: for now the implementation of an initial plane wave does not work when MPI is used; therefore I added an exit statement when both options are used. We will remove it when we find and fix that bug.


Modified: seismo/2D/SPECFEM2D/trunk/assemble_MPI.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/assemble_MPI.F90	2009-05-08 22:31:57 UTC (rev 14933)
+++ seismo/2D/SPECFEM2D/trunk/assemble_MPI.F90	2009-05-08 23:58:03 UTC (rev 14934)
@@ -601,7 +601,7 @@
 
 
 !-----------------------------------------------
-! Dummy subroutine, to be able to stop the code whether sequential or parallel.
+! subroutine to stop the code whether sequential or parallel.
 !-----------------------------------------------
 subroutine exit_MPI(error_msg)
 

Modified: seismo/2D/SPECFEM2D/trunk/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2009-05-08 22:31:57 UTC (rev 14933)
+++ seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2009-05-08 23:58:03 UTC (rev 14934)
@@ -2079,9 +2079,14 @@
 
   if(initialfield) then
      write(IOUT,*)
-     write(IOUT,*) 'Reading initial fields from external file...'
+!! DK DK reading of an initial field from an external file has been suppressed
+!! DK DK and replaced with the implementation of an analytical plane wave
+!! DK DK     write(IOUT,*) 'Reading initial fields from external file...'
+     write(IOUT,*) 'Implementing an analytical initial plane wave...'
      write(IOUT,*)
      if(any_acoustic) call exit_MPI('initial field currently implemented for purely elastic simulation only')
+!! DK DK for now there is a bug in the initial plane wave when MPI is used
+     if(nproc > 1) call exit_MPI('for now there is a bug in the initial plane wave when MPI is used')
 
       !=======================================================================
       !