diff options
Diffstat (limited to 'src/momentum_integrand3D.F')
-rw-r--r-- | src/momentum_integrand3D.F | 25 |
1 files changed, 19 insertions, 6 deletions
diff --git a/src/momentum_integrand3D.F b/src/momentum_integrand3D.F index 011fa68..cc481f9 100644 --- a/src/momentum_integrand3D.F +++ b/src/momentum_integrand3D.F @@ -5,7 +5,7 @@ c ======================================================================== SUBROUTINE momentum_integrand3D(origin,Dx,Dy,Dz,x,y,z,gxx,gxy,gxz, & gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz, - & momentum_int,Psi,Psi_power) + & momentum_int,Psi,Psi_power,conformal_state) c ------------------------------------------------------------------------ c @@ -18,7 +18,7 @@ c ----------------------------------------------------------------------- c Input variables INTEGER,INTENT(IN) :: - & Psi_power + & Psi_power,conformal_state CCTK_REAL,INTENT(IN) :: & Dx,Dy,Dz,origin(3) CCTK_REAL,DIMENSION(:),INTENT(IN) :: @@ -81,11 +81,24 @@ c Because other codes evolve Psi**4 c Abbreviations for metric coefficients c ------------------------------------- - p = psi(i,j,k)**ip + if (conformal_state > 0) then + + p = psi(i,j,k)**ip + + dxx = p*gxx(i,j,k); dxy = p*gxy(i,j,k) + dxz = p*gxz(i,j,k); dyy = p*gyy(i,j,k) + dyz = p*gyz(i,j,k); dzz = p*gzz(i,j,k) + + else + + p = 1.0d0 + + dxx = gxx(i,j,k); dxy = gxy(i,j,k) + dxz = gxz(i,j,k); dyy = gyy(i,j,k) + dyz = gyz(i,j,k); dzz = gzz(i,j,k) + + end if - dxx = p*gxx(i,j,k); dxy = p*gxy(i,j,k) - dxz = p*gxz(i,j,k); dyy = p*gyy(i,j,k) - dyz = p*gyz(i,j,k); dzz = p*gzz(i,j,k) c Determinant of 3-metric c ----------------------- |