From 37ecaa15d029f7ca274b2db95c6aced9e63b3992 Mon Sep 17 00:00:00 2001 From: dorband Date: Thu, 6 Mar 2003 09:19:16 +0000 Subject: Included calculation of a normalized hamiltonian constraint: hamnormalized git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/ADMConstraints/trunk@83 b7a48df3-cbbf-4440-997f-b4b717c9f7fc --- interface.ccl | 1 + src/ADMConstraints.F | 7 +++++++ src/InitSymBound.F | 7 +++++++ 3 files changed, 15 insertions(+) diff --git a/interface.ccl b/interface.ccl index 4d409ac..8516ed7 100644 --- a/interface.ccl +++ b/interface.ccl @@ -13,6 +13,7 @@ private: real hamiltonian type=GF { ham + hamnormalized } "Hamiltonian constraint" real momentum type=GF diff --git a/src/ADMConstraints.F b/src/ADMConstraints.F index a6a7703..9c30b1b 100644 --- a/src/ADMConstraints.F +++ b/src/ADMConstraints.F @@ -94,6 +94,7 @@ c Fill with zeros. do i=1,nx ham(i,j,k) = 0.0D0 + hamnormalized(i,j,k) = 0.0D0 momx(i,j,k) = 0.0D0 momy(i,j,k) = 0.0D0 @@ -183,6 +184,12 @@ c = (T_00 - 2 beta^i T_{i0} + beta^i beta^j T_{ij})/alpha^2 end if ham(i,j,k) = HAMADM_HAMADM - 16.0D0*pi*m_rho + if ((HAMADM_HAMADMABS + abs(16.0D0*pi*m_rho))==0) then + hamnormalized(i,j,k) = abs(ham(i,j,k)) + else + hamnormalized(i,j,k) = abs(ham(i,j,k))/ + & (HAMADM_HAMADMABS + abs(16.0D0*pi*m_rho)) + end if c Calculate the Momentum constraints diff --git a/src/InitSymBound.F b/src/InitSymBound.F index 00cfefa..1069c8e 100644 --- a/src/InitSymBound.F +++ b/src/InitSymBound.F @@ -48,6 +48,13 @@ sym(1) = -one sym(2) = one sym(3) = one + call SetCartSymVN(ierr, cctkGH, sym, 'admconstraints::hamnormalized') + if (ierr .lt. 0) then + call CCTK_WARN(0,"Error setting symmetries for ADMConstraints::hamnormalized") + end if + sym(1) = -one + sym(2) = one + sym(3) = one call SetCartSymVN(ierr, cctkGH, sym, 'admconstraints::momx') if (ierr .lt. 0) then call CCTK_WARN(0,"Error setting symmetries for ADMConstraints::momx") -- cgit v1.2.3