1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
|
#include "cctk.h"
MODULE D3_extract_int
c ------------------------------------------------------------------
INTERFACE
SUBROUTINE D3_extract(cctkGH,conformal_state,do_ADMmass,do_momentum,do_spin,igrid,
& 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)
IMPLICIT NONE
CCTK_POINTER :: cctkGH
CCTK_INT,INTENT(IN) ::
& conformal_state
INTEGER,INTENT(IN) ::
& igrid,l,m,Psi_power,myproc
CCTK_INT,INTENT(IN) ::
& 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(nx,ny,nz) ::
& Psi,g00,gxx,gxy,gxz,gyy,gyz,gzz,
& hxx,hxy,hxz,hyy,hyz,hzz
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
CCTK_REAL,INTENT(OUT) ::
& ADMmass(2),mass,rsch,Qodd(:,2:,0:),Qeven(:,2:,0:),dtaudt,
& momentum(3),spin(3)
END SUBROUTINE
END INTERFACE
END MODULE D3_extract_int
|
