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

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2012-05-19 18:49:26 -0700 (Sat, 19 May 2012)
New Revision: 20192

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90
Log:
improved the way the stress tensor is made symmetric


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90	2012-05-20 01:46:06 UTC (rev 20191)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90	2012-05-20 01:49:26 UTC (rev 20192)
@@ -391,10 +391,6 @@
                  sigma_zy = mul_unrelaxed_elastic*duy_dzl
                  sigma_zz = lambdaplus2mu_unrelaxed_elastic*duz_dzl + lambdal_unrelaxed_elastic*dux_dxl
 
-! the stress tensor is symmetric by default, unless defined otherwise
-! this can be overwritten below if needed
-                 sigma_zx = sigma_xz
-
 !! DK DK added this for Guenneau, March 2012
 #ifdef USE_GUENNEAU
   include "include_stiffness_Guenneau.f90"
@@ -437,10 +433,6 @@
                  sigma_zz = c13*dux_dxl + c35*(duz_dxl + dux_dzl) + c33*duz_dzl
                  sigma_xz = c15*dux_dxl + c55*(duz_dxl + dux_dzl) + c35*duz_dzl
 
-! the stress tensor is symmetric by default, unless defined otherwise
-! this can be overwritten below if needed
-                 sigma_zx = sigma_xz
-
               endif
 
               ! Pre-kernels calculation
@@ -465,6 +457,11 @@
 
               jacobianl = jacobian(i,j,ispec)
 
+! the stress tensor is symmetric by default, unless defined otherwise
+#ifndef USE_GUENNEAU
+              sigma_zx = sigma_xz
+#endif
+
               ! weak formulation term based on stress tensor (non-symmetric form)
               ! also add GLL integration weights
               tempx1(i,j) = wzgll(j)*jacobianl*(sigma_xx*xixl+sigma_zx*xizl) ! this goes to accel_x