diff options
Diffstat (limited to 'src/D3_extract.F')
| -rw-r--r-- | src/D3_extract.F | 29 |
1 files changed, 22 insertions, 7 deletions
diff --git a/src/D3_extract.F b/src/D3_extract.F index 47a7fb3..7ef0f0f 100644 --- a/src/D3_extract.F +++ b/src/D3_extract.F @@ -4,7 +4,8 @@ c ================================================================== SUBROUTINE D3_extract(cctkGH,conformal_state,do_ADMmass,do_momentum,do_spin,igrid, - & origin,myproc,interpolation_operator,interpolation_order,Nt,Np,all_modes, + & origin,myproc,interpolation_operator,interpolation_order,Nt,Np, + & nx,ny,nz,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, & eta,ADMmass,momentum,spin,mass,rsch,Qodd,Qeven,Extract_temp3d,dtaudt) @@ -73,17 +74,18 @@ c Input variables INTEGER,INTENT(IN) :: & conformal_state,igrid,l,m,Psi_power,myproc CCTK_INT,INTENT(IN) :: - & Nt,Np,all_modes,do_momentum,do_spin,interpolation_order + & Nt,Np,nx,ny,nz,all_modes,do_momentum,do_spin,interpolation_order INTEGER,INTENT(IN) :: & do_ADMmass(2) CCTK_REAL,INTENT(IN) :: & origin(3),Dx,Dy,Dz,eta - CCTK_REAL,INTENT(IN),DIMENSION(:,:,:) :: + CCTK_REAL,INTENT(IN),DIMENSION(nx,ny,nz) :: & Psi,g00,gxx,gxy,gxz,gyy,gyz,gzz, & hxx,hxy,hxz,hyy,hyz,hzz - CCTK_REAL,INTENT(IN),DIMENSION(:) :: - & x,y,z - CCTK_REAL,INTENT(INOUT),DIMENSION(:,:,:) :: + CCTK_REAL,INTENT(IN),DIMENSION(nx) :: x + CCTK_REAL,INTENT(IN),DIMENSION(ny) :: y + CCTK_REAL,INTENT(IN),DIMENSION(nz) :: z + CCTK_REAL,INTENT(INOUT),DIMENSION(nx,ny,nz) :: & Extract_temp3d CCTK_STRING,INTENT(IN) :: & interpolation_operator @@ -122,6 +124,11 @@ c 1. Specify polar coordinates on chosen 2-surface c c ------------------------------------------------------------------ +c print*,'Entering D3_extract' +c print*,nx,ny,nz +c print*,'D3_extract Qodd' +c print*,lbound(Qodd,1),lbound(Qodd,2),lbound(Qodd,3) +c print*,size(Qodd,1),size(Qodd,2),size(Qodd,3) Pi = two*ASIN(one) c Set polar gridspacing @@ -163,14 +170,17 @@ c 2. Project quantities onto the 2-surface c c ------------------------------------------------------------------ +c print*,nx,ny,nz,Nt,Np +c print*,'Calling D3_to_D2' CALL D3_to_D2(cctkGH,conformal_state,do_ADMmass,do_momentum,do_spin, & Psi_power,origin,myproc,interpolation_operator,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, + & x,y,z,eta,Nt,Np,nx,ny,nz,theta,phi,Psis,g00s,gxxs,gxys,gxzs, & gyys,gyzs,gzzs,dPsis,dgxxs,dgxys,dgxzs,dgyys,dgyzs,dgzzs, & ADMmass_int1,ADMmass_int2,momentum_int1,momentum_int2, & momentum_int3,spin_int1,spin_int2,spin_int3,Extract_temp3d) +c print*,'Exited D3_to_D2' if (myproc .eq. 0) then @@ -203,12 +213,17 @@ c 5. Extract gauge invariant quantities on 2-surface c c ------------------------------------------------------------------ +c print*,'Calling D2_extract' CALL D2_extract(eta,igrid,Nt,Np,theta,phi,all_modes,l,m,g00s,grr, & grt,grp,gtt,gtp,gpp,dgtt,dgtp,dgpp, & do_ADMmass,ADMmass_int1,ADMmass_int2, & do_momentum,momentum_int1,momentum_int2,momentum_int3, & do_spin,spin_int1,spin_int2,spin_int3, & mass,rsch,Qodd,Qeven,ADMmass,momentum,spin,dtaudt) +c print*,Qodd +c print*,Qeven +c print*,'Exiting D2_extract' + c End of myproc = 0 conditional endif |