[cig-commits] r15674 - seismo/2D/SPECFEM2D/trunk

[cmorency at geodynamics.org]

Author: cmorency
Date: 2009-09-18 07:59:21 -0700 (Fri, 18 Sep 2009)
New Revision: 15674

Modified:
   seismo/2D/SPECFEM2D/trunk/compute_forces_elastic.f90
   seismo/2D/SPECFEM2D/trunk/compute_forces_fluid.f90
   seismo/2D/SPECFEM2D/trunk/compute_forces_solid.f90
Log:
Fixed normal signs for absorbing boundaries. Thanks Yang Luo.


Modified: seismo/2D/SPECFEM2D/trunk/compute_forces_elastic.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/compute_forces_elastic.f90	2009-09-14 23:45:35 UTC (rev 15673)
+++ seismo/2D/SPECFEM2D/trunk/compute_forces_elastic.f90	2009-09-18 14:59:21 UTC (rev 15674)
@@ -501,8 +501,8 @@
           xgamma = - xiz(i,j,ispec) * jacobian(i,j,ispec)
           zgamma = + xix(i,j,ispec) * jacobian(i,j,ispec)
           jacobian1D = sqrt(xgamma**2 + zgamma**2)
-          nx = + zgamma / jacobian1D
-          nz = - xgamma / jacobian1D
+          nx = - zgamma / jacobian1D
+          nz = + xgamma / jacobian1D
 
           weight = jacobian1D * wzgll(j)
 
@@ -769,8 +769,8 @@
           xxi = + gammaz(i,j,ispec) * jacobian(i,j,ispec)
           zxi = - gammax(i,j,ispec) * jacobian(i,j,ispec)
           jacobian1D = sqrt(xxi**2 + zxi**2)
-          nx = + zxi / jacobian1D
-          nz = - xxi / jacobian1D
+          nx = - zxi / jacobian1D
+          nz = + xxi / jacobian1D
 
           weight = jacobian1D * wxgll(i)
 

Modified: seismo/2D/SPECFEM2D/trunk/compute_forces_fluid.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/compute_forces_fluid.f90	2009-09-14 23:45:35 UTC (rev 15673)
+++ seismo/2D/SPECFEM2D/trunk/compute_forces_fluid.f90	2009-09-18 14:59:21 UTC (rev 15674)
@@ -9,8 +9,6 @@
 !               Nicolas Le Goff, nicolas DOT legoff aT univ-pau DOT fr
 !               Roland Martin, roland DOT martin aT univ-pau DOT fr
 !               Christina Morency, cmorency aT princeton DOT edu
-!               Christina Morency, cmorency aT gps DOT caltech DOT edu
-!               Jeroen Tromp, jtromp aT gps DOT caltech DOT edu
 !
 ! This software is a computer program whose purpose is to solve
 ! the two-dimensional viscoelastic anisotropic or poroelastic wave equation
@@ -573,8 +571,8 @@
           xgamma = - xiz(i,j,ispec) * jacobian(i,j,ispec)
           zgamma = + xix(i,j,ispec) * jacobian(i,j,ispec)
           jacobian1D = sqrt(xgamma**2 + zgamma**2)
-          nx = + zgamma / jacobian1D
-          nz = - xgamma / jacobian1D
+          nx = - zgamma / jacobian1D
+          nz = + xgamma / jacobian1D
 
           weight = jacobian1D * wzgll(j)
 
@@ -748,8 +746,8 @@
           xxi = + gammaz(i,j,ispec) * jacobian(i,j,ispec)
           zxi = - gammax(i,j,ispec) * jacobian(i,j,ispec)
           jacobian1D = sqrt(xxi**2 + zxi**2)
-          nx = + zxi / jacobian1D
-          nz = - xxi / jacobian1D
+          nx = - zxi / jacobian1D
+          nz = + xxi / jacobian1D
 
           weight = jacobian1D * wxgll(i)
 

Modified: seismo/2D/SPECFEM2D/trunk/compute_forces_solid.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/compute_forces_solid.f90	2009-09-14 23:45:35 UTC (rev 15673)
+++ seismo/2D/SPECFEM2D/trunk/compute_forces_solid.f90	2009-09-18 14:59:21 UTC (rev 15674)
@@ -9,8 +9,6 @@
 !               Nicolas Le Goff, nicolas DOT legoff aT univ-pau DOT fr
 !               Roland Martin, roland DOT martin aT univ-pau DOT fr
 !               Christina Morency, cmorency aT princeton DOT edu
-!               Christina Morency, cmorency aT gps DOT caltech DOT edu
-!               Jeroen Tromp, jtromp aT gps DOT caltech DOT edu
 !
 ! This software is a computer program whose purpose is to solve
 ! the two-dimensional viscoelastic anisotropic or poroelastic wave equation
@@ -593,8 +591,8 @@
           xgamma = - xiz(i,j,ispec) * jacobian(i,j,ispec)
           zgamma = + xix(i,j,ispec) * jacobian(i,j,ispec)
           jacobian1D = sqrt(xgamma**2 + zgamma**2)
-          nx = + zgamma / jacobian1D
-          nz = - xgamma / jacobian1D
+          nx = - zgamma / jacobian1D
+          nz = + xgamma / jacobian1D
 
 
           weight = jacobian1D * wzgll(j)
@@ -770,8 +768,8 @@
           xxi = + gammaz(i,j,ispec) * jacobian(i,j,ispec)
           zxi = - gammax(i,j,ispec) * jacobian(i,j,ispec)
           jacobian1D = sqrt(xxi**2 + zxi**2)
-          nx = + zxi / jacobian1D
-          nz = - xxi / jacobian1D
+          nx = - zxi / jacobian1D
+          nz = + xxi / jacobian1D
 
           weight = jacobian1D * wxgll(i)