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

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2013-05-13 00:07:58 -0700 (Mon, 13 May 2013)
New Revision: 22052

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90
Log:
remove unneeded code in pml_init for checking


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90	2013-05-11 22:04:53 UTC (rev 22051)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/pml_init.F90	2013-05-13 07:07:58 UTC (rev 22052)
@@ -910,64 +910,5 @@
     endif
  enddo
 
- do ispec = 1,nspec
-  ispec_PML=spec_to_PML(ispec)
-    if (is_PML(ispec)) then
-       do j=1,NGLLZ
-          do i=1,NGLLX
-               iglob=ibool(i,j,ispec)
-!!!! ---------- bottom edge
-               if (region_CPML(ispec) == CPML_BOTTOM_RIGHT .or. region_CPML(ispec) == CPML_BOTTOM_LEFT .or. &
-                   region_CPML(ispec) == CPML_BOTTOM) then
-!                   write(*,*)K_x_store(i,j,ispec_PML), K_z_store(i,j,ispec_PML),&
-!                             d_x_store(i,j,ispec_PML), d_z_store(i,j,ispec_PML),&
-!                             alpha_x_store(i,j,ispec_PML), alpha_z_store(i,j,ispec_PML),& "zhinan_2"
-
-                   write(*,'(8f15.8,i2,1x,a30)')K_x_store(i,j,ispec_PML), K_z_store(i,j,ispec_PML),&
-                             d_x_store(i,j,ispec_PML), d_z_store(i,j,ispec_PML),&
-                             alpha_x_store(i,j,ispec_PML), alpha_z_store(i,j,ispec_PML),&
-                              coord(1,iglob),coord(2,iglob),region_CPML(ispec), &
-                            "zn_1,BOTTOM_RIGHT,BOTTOM_LEFT,BOTTOM"
-               endif
-
-!!!! ---------- top edge
-               if (region_CPML(ispec) == CPML_TOP_RIGHT .or. region_CPML(ispec) == CPML_TOP_LEFT .or. &
-                   region_CPML(ispec) == CPML_TOP) then
-
-                   write(*,'(8f15.8,i2,1x,a30)')K_x_store(i,j,ispec_PML), K_z_store(i,j,ispec_PML),&
-                             d_x_store(i,j,ispec_PML), d_z_store(i,j,ispec_PML),&
-                             alpha_x_store(i,j,ispec_PML), alpha_z_store(i,j,ispec_PML),&
-                              coord(1,iglob),coord(2,iglob),region_CPML(ispec), &
-                            "zn_2,TOP_RIGHT,TOP_LEFT,TOP"
-               endif
-
-!!!! ---------- right edge
-               if (region_CPML(ispec) == CPML_BOTTOM_RIGHT .or. region_CPML(ispec) == CPML_TOP_RIGHT .or. &
-                   region_CPML(ispec) == CPML_RIGHT) then
-
-                   write(*,'(8f15.8,i2,1x,a30)')K_x_store(i,j,ispec_PML), K_z_store(i,j,ispec_PML),&
-                             d_x_store(i,j,ispec_PML), d_z_store(i,j,ispec_PML),&
-                             alpha_x_store(i,j,ispec_PML), alpha_z_store(i,j,ispec_PML),&
-                              coord(1,iglob),coord(2,iglob),region_CPML(ispec), &
-                            "zn_3,BOTTOM_RIGHT,TOP_RIGHT,RIGHT"
-               endif
-
-!!!! ---------- left edge
-               if (region_CPML(ispec) == CPML_BOTTOM_LEFT .or. region_CPML(ispec) == CPML_TOP_LEFT .or. &
-                   region_CPML(ispec) == CPML_LEFT) then
-
-                   write(*,'(8f15.8,i2,1x,a30)')K_x_store(i,j,ispec_PML), K_z_store(i,j,ispec_PML),&
-                             d_x_store(i,j,ispec_PML), d_z_store(i,j,ispec_PML),&
-                             alpha_x_store(i,j,ispec_PML), alpha_z_store(i,j,ispec_PML),&
-                              coord(1,iglob),coord(2,iglob),region_CPML(ispec), &
-                            "zn_4,BOTTOM_LEFT,TOP_LEFT,LEFT"
-
-               endif
-
-       enddo
-     enddo
-    endif
- enddo
-
   end subroutine define_PML_coefficients