diff options
Diffstat (limited to 'src/spin_integrand3D.F')
| -rw-r--r-- | src/spin_integrand3D.F | 24 |
1 files changed, 18 insertions, 6 deletions
diff --git a/src/spin_integrand3D.F b/src/spin_integrand3D.F index cf617d8..7f93a0d 100644 --- a/src/spin_integrand3D.F +++ b/src/spin_integrand3D.F @@ -5,7 +5,7 @@ c ======================================================================== SUBROUTINE spin_integrand3D(origin,x,y,z,gxx,gxy,gxz, & gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz, - & spin_int,Psi,Psi_power) + & spin_int,Psi,Psi_power,conformal_state) c ------------------------------------------------------------------------ c @@ -19,7 +19,7 @@ c ----------------------------------------------------------------------- c Input variables INTEGER,INTENT(IN) :: - & Psi_power + & Psi_power,conformal_state CCTK_REAL,INTENT(IN) :: & origin(3) CCTK_REAL,DIMENSION(:),INTENT(IN) :: @@ -81,11 +81,23 @@ c ------------------------------------------------------------------------ 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) + 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 c Determinant of 3-metric c ----------------------- |