[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