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

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2012-05-19 18:46:06 -0700 (Sat, 19 May 2012)
New Revision: 20191

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_solid.f90
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90
Log:
fixed comments about the non symmetric form
and fixed a small bug of that non symmetric form in the anisotropic case


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90	2012-05-20 00:17:53 UTC (rev 20190)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_fluid.f90	2012-05-20 01:46:06 UTC (rev 20191)
@@ -403,7 +403,7 @@
 
           jacobianl = jacobian(i,j,ispec)
 
-! weak formulation term based on stress tensor (non-symmetric form)
+! weak formulation term based on stress tensor (symmetric form)
 ! also add GLL integration weights
           tempx1(i,j) = wzgll(j)*jacobianl*(sigma_xx*xixl+sigma_xz*xizl)
           tempz1(i,j) = wzgll(j)*jacobianl*(sigma_xz*xixl+sigma_zz*xizl)

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_solid.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_solid.f90	2012-05-20 00:17:53 UTC (rev 20190)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_poro_solid.f90	2012-05-20 01:46:06 UTC (rev 20191)
@@ -420,7 +420,7 @@
 
           jacobianl = jacobian(i,j,ispec)
 
-! weak formulation term based on stress tensor (non-symmetric form)
+! weak formulation term based on stress tensor (symmetric form)
 ! also add GLL integration weights
           tempx1(i,j) = wzgll(j)*jacobianl*(sigma_xx*xixl+sigma_xz*xizl)
           tempz1(i,j) = wzgll(j)*jacobianl*(sigma_xz*xixl+sigma_zz*xizl)

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90	2012-05-20 00:17:53 UTC (rev 20190)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_forces_viscoelastic.F90	2012-05-20 01:46:06 UTC (rev 20191)
@@ -437,6 +437,10 @@
                  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