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

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-04-13 08:47:19 -0700 (Sat, 13 Apr 2013)
New Revision: 21858

Added:
   seismo/3D/SPECFEM3D/trunk/utils/compute_expression_of_potential_energy.f90
Log:
added utils/compute_expression_of_potential_energy.f90


Added: seismo/3D/SPECFEM3D/trunk/utils/compute_expression_of_potential_energy.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/utils/compute_expression_of_potential_energy.f90	                        (rev 0)
+++ seismo/3D/SPECFEM3D/trunk/utils/compute_expression_of_potential_energy.f90	2013-04-13 15:47:19 UTC (rev 21858)
@@ -0,0 +1,23 @@
+
+  program expression_total_energy
+
+! Dimitri Komatitsch, CNRS Marseille, France, April 2013
+
+! wrote a small program to avoid risking having typos in the expression if I typed it manually
+
+  implicit none
+
+  integer :: i,j
+
+  print *,'potential_energy = potential_energy + &'
+
+! 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,' + '
+    enddo
+  enddo
+
+  end program expression_total_energy
+