c  Mahc_Evolve: Numerical Evolution of the General Relativistic Hydro Equations
c  Copyright (C) 2001 Mark Miller

      if(mahc_hydro) then
c     set up the Tmunu

c     (assume zero shift for now)
      enthalpyww = (rho(i,j,k) + rho(i,j,k)*eps(i,j,k) + press(i,j,k))*
     .   w_lorentz(i,j,k)**2

c     assume we have the physical metric calculated
      vlowx = gxx(i,j,k)*velx(i,j,k) + gxy(i,j,k)*vely(i,j,k) +
     .        gxz(i,j,k)*velz(i,j,k)
      vlowy = gxy(i,j,k)*velx(i,j,k) + gyy(i,j,k)*vely(i,j,k) +
     .        gyz(i,j,k)*velz(i,j,k)
      vlowz = gxz(i,j,k)*velx(i,j,k) + gyz(i,j,k)*vely(i,j,k) +
     .        gzz(i,j,k)*velz(i,j,k)

      Ttt = (enthalpyww - press(i,j,k)) * alp(i,j,k)**2
      Ttx = - enthalpyww * alp(i,j,k) * vlowx
      Tty = - enthalpyww * alp(i,j,k) * vlowy
      Ttz = - enthalpyww * alp(i,j,k) * vlowz
      Txx = enthalpyww * vlowx * vlowx + gxx(i,j,k) * press(i,j,k)
      Txy = enthalpyww * vlowx * vlowy + gxy(i,j,k) * press(i,j,k)
      Txz = enthalpyww * vlowx * vlowz + gxz(i,j,k) * press(i,j,k)
      Tyy = enthalpyww * vlowy * vlowy + gyy(i,j,k) * press(i,j,k)
      Tyz = enthalpyww * vlowy * vlowz + gyz(i,j,k) * press(i,j,k)
      Tzz = enthalpyww * vlowz * vlowz + gzz(i,j,k) * press(i,j,k)

      endif