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
52
53
54
55
56
|
#include "cctk.h"
c ==================================================================
SUBROUTINE unphysical_to_physical(grr,grt,grp,gtt,gtp,gpp,dgtt,
& dgtp,dgpp,Psis,dPsis,Psi_power)
c ------------------------------------------------------------------
c
c Convert unphysical metric components and radial (eta) derivatives
c on the 2-sphere into physical quantities, using the conformal factor
c on the sphere
c
c ------------------------------------------------------------------
IMPLICIT NONE
c Input variables
INTEGER ::
& Psi_power
CCTK_REAL,INTENT(INOUT),DIMENSION (:,:) ::
& grr,grt,grp,gtt,gtp,gpp,dgtt,dgtp,dgpp
CCTK_REAL,INTENT(IN),DIMENSION (:,:) ::
& Psis,dPsis
c Output variables
c WE ARE CHANGING THE INPUT VARIABLES !!!!
c ------------------------------------------------------------------
dgtt = Psis**4*dgtt + 4.0d0*Psis**3*dPsis*gtt
dgtp = Psis**4*dgtp + 4.0d0*Psis**3*dPsis*gtp
dgpp = Psis**4*dgpp + 4.0d0*Psis**3*dPsis*gpp
grr = Psis**4*grr
grt = Psis**4*grt
grp = Psis**4*grp
gtt = Psis**4*gtt
gtp = Psis**4*gtp
gpp = Psis**4*gpp
END SUBROUTINE unphysical_to_physical
c ==================================================================
|