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

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2012-09-21 04:08:25 -0700 (Fri, 21 Sep 2012)
New Revision: 20754

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/Makefile.in
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/invert_mass_matrix.F90
Log:
fixed a small bug in the new (improved) Stacey conditions


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/Makefile.in
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/Makefile.in	2012-09-20 21:58:52 UTC (rev 20753)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/Makefile.in	2012-09-21 11:08:25 UTC (rev 20754)
@@ -1,10 +1,11 @@
 
 #========================================================================
 #
-#                   S P E C F E M 2 D  Version 6 . 2
-#                   ------------------------------
 #
-# Copyright Universite de Pau, CNRS and INRIA, France,
+#                   S P E C F E M 2 D  Version 7 . 0
+#                   --------------------------------
+#
+# Copyright CNRS, INRIA and University of Pau, France,
 # and Princeton University / California Institute of Technology, USA.
 # Contributors: Dimitri Komatitsch, dimitri DOT komatitsch aT univ-pau DOT fr
 #               Nicolas Le Goff, nicolas DOT legoff aT univ-pau DOT fr

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/invert_mass_matrix.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/invert_mass_matrix.F90	2012-09-20 21:58:52 UTC (rev 20753)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/invert_mass_matrix.F90	2012-09-21 11:08:25 UTC (rev 20754)
@@ -331,10 +331,8 @@
                  ty = rho_vs*vy
                  tz = rho_vp*vn*nz+rho_vs*(vz-vn*nz)
 
-                rmass_inverse_elastic_one(iglob) = rmass_inverse_elastic_one(iglob)  &
-                    + (tx)*weight
-                rmass_inverse_elastic_three(iglob) = rmass_inverse_elastic_one(iglob)  &
-                    + (tz)*weight
+                rmass_inverse_elastic_one(iglob) = rmass_inverse_elastic_one(iglob) + tx*weight
+                rmass_inverse_elastic_three(iglob) = rmass_inverse_elastic_three(iglob) + tz*weight
 
               endif
            enddo
@@ -384,10 +382,8 @@
                  ty = rho_vs*vy
                  tz = rho_vp*vn*nz+rho_vs*(vz-vn*nz)
 
-                rmass_inverse_elastic_one(iglob) = rmass_inverse_elastic_one(iglob)  &
-                    + (tx)*weight
-                rmass_inverse_elastic_three(iglob) = rmass_inverse_elastic_one(iglob)  &
-                    + (tz)*weight
+                rmass_inverse_elastic_one(iglob) = rmass_inverse_elastic_one(iglob) + tx*weight
+                rmass_inverse_elastic_three(iglob) = rmass_inverse_elastic_three(iglob) + tz*weight
 
               endif
 
@@ -443,13 +439,9 @@
                    tx = 0
                    ty = 0
                    tz = 0
-                rmass_inverse_elastic_one(iglob) = rmass_inverse_elastic_one(iglob)
-                rmass_inverse_elastic_three(iglob) = rmass_inverse_elastic_three(iglob)
                  else
-                rmass_inverse_elastic_one(iglob) = rmass_inverse_elastic_one(iglob)  &
-                    + (tx)*weight
-                rmass_inverse_elastic_three(iglob) = rmass_inverse_elastic_one(iglob)  &
-                    + (tz)*weight
+                   rmass_inverse_elastic_one(iglob) = rmass_inverse_elastic_one(iglob) + tx*weight
+                   rmass_inverse_elastic_three(iglob) = rmass_inverse_elastic_three(iglob) + tz*weight
                  endif
 
              endif
@@ -508,13 +500,9 @@
                    tx = 0
                    ty = 0
                    tz = 0
-                rmass_inverse_elastic_one(iglob) = rmass_inverse_elastic_one(iglob)
-                rmass_inverse_elastic_three(iglob) = rmass_inverse_elastic_three(iglob)
                  else
-                rmass_inverse_elastic_one(iglob) = rmass_inverse_elastic_one(iglob)  &
-                    + (tx)*weight
-                rmass_inverse_elastic_three(iglob) = rmass_inverse_elastic_one(iglob)  &
-                    + (tz)*weight
+                   rmass_inverse_elastic_one(iglob) = rmass_inverse_elastic_one(iglob) + tx*weight
+                   rmass_inverse_elastic_three(iglob) = rmass_inverse_elastic_three(iglob) + tz*weight
                  endif
             endif
 
@@ -582,3 +570,4 @@
     rmass_inverse_acoustic(:) = 1._CUSTOM_REAL / rmass_inverse_acoustic(:)
 
   end subroutine invert_mass_matrix
+