From 892f6795a04d899c52b67b2829e41fc69fb2d9f9 Mon Sep 17 00:00:00 2001 From: miguel Date: Tue, 20 Mar 2001 10:35:30 +0000 Subject: Cleaning up. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/ADMConstraints/trunk@45 b7a48df3-cbbf-4440-997f-b4b717c9f7fc --- src/ADMConstraints.F | 13 +++++-------- 1 file changed, 5 insertions(+), 8 deletions(-) diff --git a/src/ADMConstraints.F b/src/ADMConstraints.F index c97828d..453a43c 100644 --- a/src/ADMConstraints.F +++ b/src/ADMConstraints.F @@ -113,12 +113,12 @@ c Calculate constraints. c Calculate the stress energy tensor at this point c ------------------------------------------------ -c Fill with zeroes. +c Initialize stress-energy tensor components. Ttt = 0.0D0; Ttx = 0.0D0; Tty = 0.0D0; Ttz = 0.0D0 - Txx = 0.0D0; Txy = 0.0D0; Txz = 0.0D0 - Tyy = 0.0D0; Tyz = 0.0D0; Tzz = 0.0D0 + Txx = 0.0D0; Tyy = 0.0D0; Tzz = 0.0D0 + Txy = 0.0D0; Txz = 0.0D0; Tyz = 0.0D0 c Inverse metric may be needed for CalcTmunu. @@ -127,8 +127,8 @@ c Inverse metric may be needed for CalcTmunu. det = DETG_DETCG - uxx = UPPERMET_UXX; uxy = UPPERMET_UXY; uxz = UPPERMET_UXZ - uyy = UPPERMET_UYY; uyz = UPPERMET_UYZ; uzz = UPPERMET_UZZ + uxx = UPPERMET_UXX; uyy = UPPERMET_UYY; uzz = UPPERMET_UZZ + uxy = UPPERMET_UXY; uxz = UPPERMET_UXZ; uyz = UPPERMET_UYZ c Call macro for stress energy tensor. @@ -150,7 +150,6 @@ c = (T_00 - 2 beta^i T_{i0} + beta^i beta^j T_{ij})/alpha^2 m_rho = ialp2*Ttt - if (shift_state == SHIFT_ACTIVE) then m_rho = m_rho + ialp2 @@ -164,7 +163,6 @@ c = (T_00 - 2 beta^i T_{i0} + beta^i beta^j T_{ij})/alpha^2 & + betay(i,j,k)**2*Tty & + betaz(i,j,k)**2*Ttz)*2.0D0) - end if ham(i,j,k) = HAMADM_HAMADM - 16.0D0*pi*m_rho @@ -197,7 +195,6 @@ c S_i = - g_{ia} n_b T^{ab} = - g_i^a n^b T_{ab} c = - (T_{i0} - beta^j T_{ij})/alpha m_sx = - ialp*Ttx - m_sy = - ialp*Tty m_sz = - ialp*Ttz -- cgit v1.2.3