[cig-commits] r21860 - seismo/3D/SPECFEM3D/trunk/utils

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-04-13 10:10:53 -0700 (Sat, 13 Apr 2013)
New Revision: 21860

Modified:
   seismo/3D/SPECFEM3D/trunk/utils/compute_expression_of_potential_energy.f90
Log:
added integration weight and 1/2 coefficient to utils/compute_expression_of_potential_energy.f90


Modified: seismo/3D/SPECFEM3D/trunk/utils/compute_expression_of_potential_energy.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/utils/compute_expression_of_potential_energy.f90	2013-04-13 16:36:29 UTC (rev 21859)
+++ seismo/3D/SPECFEM3D/trunk/utils/compute_expression_of_potential_energy.f90	2013-04-13 17:10:53 UTC (rev 21860)
@@ -9,13 +9,19 @@
 
   integer :: i,j
 
-  print *,'potential_energy = potential_energy + &'
+! compute potential energy 1/2 sigma_ij epsilon_ij
 
+  print *,'potential_energy = potential_energy + integration_weight * 0.5 * ( &'
+
 ! expression in 3D
 ! 1 = x, 2 = y and 3 = z
   do i = 1,3
     do j = 1,3
-      print *,'sigma_',i,j,' * epsilon_',i,j,' + '
+      if(i == 3 .and. j == 3) then
+        print *,'sigma_',i,j,' * epsilon_',i,j,')'
+      else
+        print *,'sigma_',i,j,' * epsilon_',i,j,' + '
+      endif
     enddo
   enddo