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

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2012-08-08 18:22:39 -0700 (Wed, 08 Aug 2012)
New Revision: 20560

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
Log:
delete some used variables


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90	2012-08-09 01:10:54 UTC (rev 20559)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90	2012-08-09 01:22:39 UTC (rev 20560)
@@ -42,7 +42,7 @@
 !========================================================================
 
   subroutine pml_init(nspec,nglob,anyabs,ibool,nelemabs,codeabs,numabs,&
-                    nspec_PML,is_PML,which_PML_elem,which_PML_poin,spec_to_PML, &
+                    nspec_PML,is_PML,which_PML_elem,spec_to_PML, &
                     icorner_iglob,NELEM_PML_THICKNESS)
 
 
@@ -60,7 +60,6 @@
   integer :: i,j,k,i_coef
 
   logical, dimension(4,nspec) :: which_PML_elem
-  logical, dimension(4,nglob) :: which_PML_poin
   integer, dimension(nglob) ::   icorner_iglob
   integer, dimension(nelemabs) :: numabs
   logical, dimension(4,nelemabs) :: codeabs
@@ -97,9 +96,6 @@
               ncorner=ncorner+1
               icorner_iglob(ncorner) = iglob
 
-              !array to know which PML it is
-              which_PML_poin(ibound,iglob) = codeabs(ibound,ispecabs)
-
            enddo
         enddo
         endif ! we are on the good absorbing boundary
@@ -142,7 +138,6 @@
                     end do
                     ncorner=ncorner+1
                     icorner_iglob(ncorner) = iglob
-                    which_PML_poin(ibound,iglob) = .true.
                  end do
               end do
               nspec_PML=nspec_PML+1

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90	2012-08-09 01:10:54 UTC (rev 20559)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90	2012-08-09 01:22:39 UTC (rev 20560)
@@ -997,7 +997,7 @@
    rmemory_dux_dx,rmemory_duz_dx,rmemory_dux_dz,rmemory_duz_dz
 
   integer, dimension(:), allocatable :: spec_to_PML,icorner_iglob
-  logical, dimension(:,:), allocatable :: which_PML_elem, which_PML_poin
+  logical, dimension(:,:), allocatable :: which_PML_elem
 
   real(kind=CUSTOM_REAL), dimension(:,:,:,:,:), allocatable :: rmemory_displ_elastic
 
@@ -2815,19 +2815,16 @@
       allocate(is_PML(nspec))
       allocate(icorner_iglob(nglob))
       allocate(which_PML_elem(4,nspec))
-      allocate(which_PML_poin(4,nglob))
       allocate(spec_to_PML(nspec))
 
       is_PML(:) = .false.
       which_PML_elem(:,:) = .false.
-      which_PML_poin(:,:) = .false.
 
       call pml_init(nspec,nglob,anyabs,ibool,nelemabs,codeabs,numabs,&
-                  nspec_PML,is_PML,which_PML_elem,which_PML_poin,spec_to_PML, &
+                  nspec_PML,is_PML,which_PML_elem,spec_to_PML, &
                   icorner_iglob,NELEM_PML_THICKNESS)
 
       deallocate(icorner_iglob)
-      deallocate(which_PML_poin)
 
       if (nspec_PML==0) nspec_PML=1
 !!!!!!!!!!!!! DK DK added this