[cig-commits] r22797 - seismo/2D/SPECFEM2D/trunk/src/specfem2D
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Tue Sep 17 15:36:21 PDT 2013
Author: dkomati1
Date: 2013-09-17 15:36:21 -0700 (Tue, 17 Sep 2013)
New Revision: 22797
Modified:
seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_pressure.f90
Log:
removed a useless comment
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_pressure.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_pressure.f90 2013-09-17 16:14:58 UTC (rev 22796)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/compute_pressure.f90 2013-09-17 22:36:21 UTC (rev 22797)
@@ -513,10 +513,6 @@
endif
! store pressure
-!! DK DK the formula on the line below is not correct because it does not take into account
-!! DK DK the fact that sigma_yy is not zero in the plane strain case
-! pressure_element(i,j) = - (sigma_xx + sigma_zz) / 2.d0
-!! DK DK the right formula is
! sigma_yy is not equal to zero in a 2D medium because of the plane strain formulation
pressure_element(i,j) = - (sigma_xx + sigma_yy + sigma_zz) / 3.d0
! pressure_element2(i,j) = - sigmap
More information about the CIG-COMMITS
mailing list