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

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2012-07-11 09:38:08 -0700 (Wed, 11 Jul 2012)
New Revision: 20520

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
Log:
an important modification when using PML in MPI version compiled by the all checked ifort flag. The MPI_version of PML for acoustic/elastic simulation is ready


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90	2012-07-11 11:32:47 UTC (rev 20519)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90	2012-07-11 16:38:08 UTC (rev 20520)
@@ -1311,12 +1311,6 @@
     stop 'PML boundary conditions do not ready for poroelastic simulation'
   endif
 
-#ifdef USE_MPI
-  if(PML_BOUNDARY_CONDITIONS .and. (any_acoustic .or. any_poroelastic) )then
-   stop 'PML_BOUNDARY_CONDITIONS do not ready for mpi version of the code'
-  endif
-#endif
-
   ! allocate memory variables for attenuation
   if(ipass == 1) then
     allocate(e1(NGLLX,NGLLZ,nspec_allocate,N_SLS))
@@ -2872,6 +2866,7 @@
         allocate(rmemory_duz_dx(2,NGLLX,NGLLZ,nspec_PML))
         allocate(rmemory_duz_dz(2,NGLLX,NGLLZ,nspec_PML))
 
+        rmemory_displ_elastic(:,:,:,:,:) = ZERO
         rmemory_dux_dx(:,:,:,:) = ZERO
         rmemory_dux_dz(:,:,:,:) = ZERO
         rmemory_duz_dx(:,:,:,:) = ZERO
@@ -2891,10 +2886,13 @@
       if (any_acoustic .and. nspec_PML>0) then
 
         allocate(rmemory_potential_acoustic(2,NGLLX,NGLLZ,nspec_PML))
-
         allocate(rmemory_acoustic_dux_dx(2,NGLLX,NGLLZ,nspec_PML))
         allocate(rmemory_acoustic_dux_dz(2,NGLLX,NGLLZ,nspec_PML))
 
+        rmemory_potential_acoustic = ZERO
+        rmemory_acoustic_dux_dx = ZERO
+        rmemory_acoustic_dux_dz = ZERO
+
       else
 
         allocate(rmemory_potential_acoustic(1,1,1,1))