diff options
Diffstat (limited to 'src/D3_extract.F')
-rw-r--r-- | src/D3_extract.F | 11 |
1 files changed, 6 insertions, 5 deletions
diff --git a/src/D3_extract.F b/src/D3_extract.F index 4fb6899..9f6fcad 100644 --- a/src/D3_extract.F +++ b/src/D3_extract.F @@ -3,7 +3,7 @@ c ================================================================== - SUBROUTINE D3_extract(cctkGH,do_ADMmass,do_momentum,do_spin,igrid, + SUBROUTINE D3_extract(cctkGH,conformal_state,do_ADMmass,do_momentum,do_spin,igrid, & origin,myproc,interpolation_order,Nt,Np,all_modes, & l,m,x,y,z,Dx,Dy,Dz,Psi_power,Psi, & g00,gxx,gxy,gxz,gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz, @@ -71,7 +71,7 @@ c Input variables CCTK_POINTER :: cctkGH INTEGER,INTENT(IN) :: - & igrid,l,m,Psi_power,myproc + & conformal_state,igrid,l,m,Psi_power,myproc CCTK_INT,INTENT(IN) :: & Nt,Np,all_modes,do_momentum,do_spin,interpolation_order INTEGER,INTENT(IN) :: @@ -161,7 +161,7 @@ c 2. Project quantities onto the 2-surface c c ------------------------------------------------------------------ - CALL D3_to_D2(cctkGH,do_ADMmass,do_momentum,do_spin, + CALL D3_to_D2(cctkGH,conformal_state,do_ADMmass,do_momentum,do_spin, & Psi_power,origin,myproc,interpolation_order,Dx,Dy,Dz,Psi, & g00,gxx,gxy,gxz,gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz, & x,y,z,eta,Nt,Np,theta,phi,Psis,g00s,gxxs,gxys,gxzs, @@ -189,9 +189,10 @@ c 4. Calculate physical quantities on 2-surface c c ------------------------------------------------------------------ - CALL unphysical_to_physical(grr,grt,grp,gtt,gtp,gpp,dgtt,dgtp, + if (conformal_state > 0) then + CALL unphysical_to_physical(grr,grt,grp,gtt,gtp,gpp,dgtt,dgtp, & dgpp,Psis,dPsis,Psi_power) - + end if c ------------------------------------------------------------------ c |