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

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2013-04-03 11:42:56 -0700 (Wed, 03 Apr 2013)
New Revision: 21705

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.F90
Log:
fix small bug in defining ib_bottom,ib_top,ib_left,ib_right using in backward simulation when using Stacey boundary condition


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.F90	2013-04-03 09:11:22 UTC (rev 21704)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.F90	2013-04-03 18:42:56 UTC (rev 21705)
@@ -784,19 +784,20 @@
 
     ! boundary element numbering
     do inum = 1,nelemabs
-      if (typeabs(inum) == IBOTTOM) then
+
+      if (codeabs(IEDGE1,inum)) then
         nspec_bottom = nspec_bottom + 1
         ib_bottom(inum) =  nspec_bottom
 
-      else if (typeabs(inum) == IRIGHT) then
+      else if (codeabs(IEDGE2,inum)) then
         nspec_right = nspec_right + 1
         ib_right(inum) =  nspec_right
 
-      else if (typeabs(inum) == ITOP) then
+      else if (codeabs(IEDGE3,inum)) then
         nspec_top = nspec_top + 1
         ib_top(inum) = nspec_top
 
-      else if (typeabs(inum) == ILEFT) then
+      else if (codeabs(IEDGE4,inum)) then
         nspec_left = nspec_left + 1
         ib_left(inum) =  nspec_left