diff options
Diffstat (limited to 'src/ADMConstraints.F')
-rw-r--r-- | src/ADMConstraints.F | 27 |
1 files changed, 20 insertions, 7 deletions
diff --git a/src/ADMConstraints.F b/src/ADMConstraints.F index 83d3aa0..8c7658a 100644 --- a/src/ADMConstraints.F +++ b/src/ADMConstraints.F @@ -157,16 +157,29 @@ c This may be needed for CalcTmunu c Initialize stress-energy tensor components. - Ttt = 0.0D0 - - Ttx = 0.0D0; Tty = 0.0D0; Ttz = 0.0D0 - - Txx = 0.0D0; Tyy = 0.0D0; Tzz = 0.0D0 - Txy = 0.0D0; Txz = 0.0D0; Tyz = 0.0D0 + if (stress_energy_state .ne. 0) then + + Ttt = eTtt(i,j,k) + + Ttx = eTtx(i,j,k); Tty = eTty(i,j,k); Ttz = eTtz(i,j,k) -c Include macro for stress energy tensor. + Txx = eTxx(i,j,k); Tyy = eTyy(i,j,k); Tzz = eTzz(i,j,k) + + + else + + Ttt = 0.0D0 + + Ttx = 0.0D0; Tty = 0.0D0; Ttz = 0.0D0 + + Txx = 0.0D0; Tyy = 0.0D0; Tzz = 0.0D0 + Txy = 0.0D0; Txz = 0.0D0; Tyz = 0.0D0 + +c Include macro for stress energy tensor. #include "CalcTmunu.inc" + + end if c Calculate the hamiltonian constraint c ------------------------------------ |