[cig-commits] r19168 - seismo/2D/SPECFEM2D/trunk/src/specfem2D
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Wed Nov 9 13:39:45 PST 2011
Author: dkomati1
Date: 2011-11-09 13:39:45 -0800 (Wed, 09 Nov 2011)
New Revision: 19168
Modified:
seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_initialfield.F90
Log:
fixed a bug that I had introduced when I commented out the corner points removed for plane waves
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90 2011-11-09 17:59:27 UTC (rev 19167)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_acoustic.f90 2011-11-09 21:39:45 UTC (rev 19168)
@@ -333,9 +333,9 @@
ibegin = ibegin_bottom(ispecabs)
iend = iend_bottom(ispecabs)
-!! DK DK not needed! exclude corners to make sure there is no contradiction on the normal
-!! DK DK not needed if(codeabs(ILEFT,ispecabs)) ibegin = 2
-!! DK DK not needed if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
+! exclude corners to make sure there is no contradiction on the normal
+ if(codeabs(ILEFT,ispecabs)) ibegin = 2
+ if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
do i = ibegin,iend
@@ -380,9 +380,9 @@
ibegin = ibegin_top(ispecabs)
iend = iend_top(ispecabs)
-!! DK DK not needed! exclude corners to make sure there is no contradiction on the normal
-!! DK DK not needed if(codeabs(ILEFT,ispecabs)) ibegin = 2
-!! DK DK not needed if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
+! exclude corners to make sure there is no contradiction on the normal
+ if(codeabs(ILEFT,ispecabs)) ibegin = 2
+ if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
do i = ibegin,iend
@@ -696,9 +696,9 @@
j = 1
ibegin = ibegin_bottom(ispecabs)
iend = iend_bottom(ispecabs)
-!! DK DK not needed ! exclude corners to make sure there is no contradiction on the normal
-!! DK DK not needed if(codeabs(ILEFT,ispecabs)) ibegin = 2
-!! DK DK not needed if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
+ ! exclude corners to make sure there is no contradiction on the normal
+ if(codeabs(ILEFT,ispecabs)) ibegin = 2
+ if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
do i = ibegin,iend
iglob = ibool(i,j,ispec)
! external velocity model
@@ -743,9 +743,9 @@
j = NGLLZ
ibegin = ibegin_top(ispecabs)
iend = iend_top(ispecabs)
-!! DK DK not needed ! exclude corners to make sure there is no contradiction on the normal
-!! DK DK not needed if(codeabs(ILEFT,ispecabs)) ibegin = 2
-!! DK DK not needed if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
+ ! exclude corners to make sure there is no contradiction on the normal
+ if(codeabs(ILEFT,ispecabs)) ibegin = 2
+ if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
do i = ibegin,iend
iglob = ibool(i,j,ispec)
! external velocity model
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90 2011-11-09 17:59:27 UTC (rev 19167)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90 2011-11-09 21:39:45 UTC (rev 19168)
@@ -677,9 +677,9 @@
ibegin = ibegin_bottom_poro(ispecabs)
iend = iend_bottom_poro(ispecabs)
-!! DK DK not needed! exclude corners to make sure there is no contradiction on the normal
-!! DK DK not needed if(codeabs(ILEFT,ispecabs)) ibegin = 2
-!! DK DK not needed if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
+! exclude corners to make sure there is no contradiction on the normal
+ if(codeabs(ILEFT,ispecabs)) ibegin = 2
+ if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
do i = ibegin,iend
@@ -737,9 +737,9 @@
ibegin = ibegin_top_poro(ispecabs)
iend = iend_top_poro(ispecabs)
-!! DK DK not needed! exclude corners to make sure there is no contradiction on the normal
-!! DK DK not needed if(codeabs(ILEFT,ispecabs)) ibegin = 2
-!! DK DK not needed if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
+! exclude corners to make sure there is no contradiction on the normal
+ if(codeabs(ILEFT,ispecabs)) ibegin = 2
+ if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
do i = ibegin,iend
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.f90 2011-11-09 17:59:27 UTC (rev 19167)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.f90 2011-11-09 21:39:45 UTC (rev 19168)
@@ -662,11 +662,11 @@
j = 1
-!! DK DK not needed ! exclude corners to make sure there is no contradiction on the normal
+ ! exclude corners to make sure there is no contradiction on the normal
ibegin = 1
iend = NGLLX
-!! DK DK not needed if(codeabs(ILEFT,ispecabs)) ibegin = 2
-!! DK DK not needed if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
+ if(codeabs(ILEFT,ispecabs)) ibegin = 2
+ if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
do i = ibegin,iend
@@ -759,11 +759,11 @@
j = NGLLZ
-!! DK DK not needed ! exclude corners to make sure there is no contradiction on the normal
+ ! exclude corners to make sure there is no contradiction on the normal
ibegin = 1
iend = NGLLX
-!! DK DK not needed if(codeabs(ILEFT,ispecabs)) ibegin = 2
-!! DK DK not needed if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
+ if(codeabs(ILEFT,ispecabs)) ibegin = 2
+ if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
do i = ibegin,iend
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_initialfield.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_initialfield.F90 2011-11-09 17:59:27 UTC (rev 19167)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/prepare_initialfield.F90 2011-11-09 21:39:45 UTC (rev 19168)
@@ -402,11 +402,11 @@
endif
if(codeabs(IBOTTOM,ispecabs)) then
j = 1
-!! DK DK not needed ! exclude corners to make sure there is no contradiction regarding the normal
+ ! exclude corners to make sure there is no contradiction regarding the normal
ibegin = 1
iend = NGLLX
-!! DK DK not needed if(codeabs(ILEFT,ispecabs)) ibegin = 2
-!! DK DK not needed if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
+ if(codeabs(ILEFT,ispecabs)) ibegin = 2
+ if(codeabs(IRIGHT,ispecabs)) iend = NGLLX-1
do i = ibegin,iend
count_bottom=count_bottom+1
iglob = ibool(i,j,ispec)
More information about the CIG-COMMITS
mailing list