From 0943c31af058628296872608289c83eed508b6cd Mon Sep 17 00:00:00 2001
From: Grzegorz Kowal <grzegorz@gkowal.info>
Date: Mon, 9 May 2011 18:36:20 -0300
Subject: [PATCH] Calculate properly kinetic energy in shape_binaries().

---
 src/problem.F90 | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/problem.F90 b/src/problem.F90
index bca6206..1f1bb9e 100644
--- a/src/problem.F90
+++ b/src/problem.F90
@@ -1532,7 +1532,7 @@ module problem
 #endif /* GLM */
 #endif /* MHD */
 #ifdef ADI
-            ekin = ekstar * r1 * r1
+            ekin = 0.5d0 * sum(pdata%u(imx:imz,i,j,k)**2) / pdata%u(idn,i,j,k)
             pdata%u(ien,i,j,k) = gammam1i * prstar + ekin
 #endif /* ADI */
           end if
@@ -1559,7 +1559,7 @@ module problem
 #endif /* GLM */
 #endif /* MHD */
 #ifdef ADI
-            ekin = eksat * r2 * r2
+            ekin = 0.5d0 * sum(pdata%u(imx:imz,i,j,k)**2) / pdata%u(idn,i,j,k)
             pdata%u(ien,i,j,k) = gammam1i * prsat + ekin
 #endif /* ADI */
           end if