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

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2013-04-10 12:39:47 -0700 (Wed, 10 Apr 2013)
New Revision: 21796

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
Log:
fix one bug


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90	2013-04-10 19:37:04 UTC (rev 21795)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90	2013-04-10 19:39:47 UTC (rev 21796)
@@ -494,8 +494,10 @@
         do i = 1, NGLLX
           do j = 1, NGLLZ
              iglob = ibool(i,j,ispec)
-             if(mask_ibool(iglob))nglob_interface = nglob_interface + 1
-             point_interface(nglob_interface)= iglob
+             if(mask_ibool(iglob))then
+              nglob_interface = nglob_interface + 1
+              point_interface(nglob_interface)= iglob
+             endif
           enddo
         enddo
        endif

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90	2013-04-10 19:37:04 UTC (rev 21795)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90	2013-04-10 19:39:47 UTC (rev 21796)
@@ -889,7 +889,7 @@
 ! to dump the wave field
   integer :: icounter,nb_of_values_to_save
   logical :: this_is_the_first_time_we_dump
-  logical, dimension(:), allocatable  :: mask_ibool
+  logical, dimension(:), allocatable  :: mask_ibool,mask_ibool_pml !zhinan
 
   double precision, dimension(:,:,:),allocatable:: rho_local,vp_local,vs_local
 !!!! hessian
@@ -2944,16 +2944,16 @@
       which_PML_elem(:,:) = .false.
 !   DK DK add support for using pml in mpi mode with external mesh
       if(read_external_mesh)then
-            allocate(mask_ibool(nglob))
+            allocate(mask_ibool_pml(nglob))
       else
-            allocate(mask_ibool(1))
+            allocate(mask_ibool_pml(1))
       endif
 
       call pml_init(nspec,nglob,anyabs,ibool,nelemabs,codeabs,numabs,&
                   nspec_PML,is_PML,which_PML_elem,spec_to_PML, &
                   icorner_iglob,NELEM_PML_THICKNESS,&
                   read_external_mesh,region_CPML,&
-                  SIMULATION_TYPE,PML_interior_interface,nglob_interface,SAVE_FORWARD,myrank,mask_ibool)
+                  SIMULATION_TYPE,PML_interior_interface,nglob_interface,SAVE_FORWARD,myrank,mask_ibool_pml)
 
       if((SIMULATION_TYPE == 3 .or. (SIMULATION_TYPE == 1 .and. SAVE_FORWARD)) .and. PML_BOUNDARY_CONDITIONS)then
          allocate(point_interface(nglob_interface))
@@ -2966,9 +2966,9 @@
 
          call determin_interface_pml_interior(nglob_interface,nspec,ibool,PML_interior_interface,&
                                               which_PML_elem,point_interface,read_external_mesh,&
-                                              mask_ibool,region_CPML,nglob)
+                                              mask_ibool_pml,region_CPML,nglob)
          deallocate(PML_interior_interface)
-         deallocate(mask_ibool)
+         deallocate(mask_ibool_pml)
 
          if(any_elastic .and. nglob_interface > 0)then
            write(outputname,'(a,i6.6,a)') 'pml_interface_elastic',myrank,'.bin'