diff options
Diffstat (limited to 'src/qlm_weyl_scalars.F90')
| -rw-r--r-- | src/qlm_weyl_scalars.F90 | 57 |
1 files changed, 57 insertions, 0 deletions
diff --git a/src/qlm_weyl_scalars.F90 b/src/qlm_weyl_scalars.F90 index d79751c..ade2e31 100644 --- a/src/qlm_weyl_scalars.F90 +++ b/src/qlm_weyl_scalars.F90 @@ -24,6 +24,9 @@ subroutine qlm_calc_weyl_scalars (CCTK_ARGUMENTS, hn) CCTK_REAL, parameter :: two=2, four=4 CCTK_REAL :: gg(3,3), dgg(3,3,3), ddgg(3,3,3,3), gg_dot(3,3), gg_dot2(3,3), dgg_dot(3,3,3) + CCTK_REAL :: kk(3,3), dkk(3,3,3), kk_dot(3,3) + CCTK_REAL :: tt(3,3) + CCTK_REAL :: dtg, gu(3,3), dgu(3,3,3), gamma(3,3,3), dgamma(3,3,3,3), ri(3,3) CCTK_REAL :: g4(0:3,0:3), dg4(0:3,0:3,0:3), ddg4(0:3,0:3,0:3,0:3) CCTK_REAL :: gu4(0:3,0:3), dgu4(0:3,0:3,0:3) CCTK_REAL :: gamma4(0:3,0:3,0:3), dgamma4(0:3,0:3,0:3,0:3) @@ -122,6 +125,48 @@ subroutine qlm_calc_weyl_scalars (CCTK_ARGUMENTS, hn) ddgg(:,:,3,1) = ddgg(:,:,1,3) ddgg(:,:,3,2) = ddgg(:,:,2,3) + kk(1,1) = qlm_kxx(i,j) + kk(1,2) = qlm_kxy(i,j) + kk(1,3) = qlm_kxz(i,j) + kk(2,2) = qlm_kyy(i,j) + kk(2,3) = qlm_kyz(i,j) + kk(3,3) = qlm_kzz(i,j) + kk(2,1) = kk(1,2) + kk(3,1) = kk(1,3) + kk(3,2) = kk(2,3) + + dkk(1,1,1) = qlm_dkxxx(i,j) + dkk(1,2,1) = qlm_dkxyx(i,j) + dkk(1,3,1) = qlm_dkxzx(i,j) + dkk(2,2,1) = qlm_dkyyx(i,j) + dkk(2,3,1) = qlm_dkyzx(i,j) + dkk(3,3,1) = qlm_dkzzx(i,j) + dkk(1,1,2) = qlm_dkxxy(i,j) + dkk(1,2,2) = qlm_dkxyy(i,j) + dkk(1,3,2) = qlm_dkxzy(i,j) + dkk(2,2,2) = qlm_dkyyy(i,j) + dkk(2,3,2) = qlm_dkyzy(i,j) + dkk(3,3,2) = qlm_dkzzy(i,j) + dkk(1,1,3) = qlm_dkxxz(i,j) + dkk(1,2,3) = qlm_dkxyz(i,j) + dkk(1,3,3) = qlm_dkxzz(i,j) + dkk(2,2,3) = qlm_dkyyz(i,j) + dkk(2,3,3) = qlm_dkyzz(i,j) + dkk(3,3,3) = qlm_dkzzz(i,j) + dkk(2,1,:) = dkk(1,2,:) + dkk(3,1,:) = dkk(1,3,:) + dkk(3,2,:) = dkk(2,3,:) + + tt(1,1) = qlm_txx(i,j) + tt(1,2) = qlm_txy(i,j) + tt(1,3) = qlm_txz(i,j) + tt(2,2) = qlm_tyy(i,j) + tt(2,3) = qlm_tyz(i,j) + tt(3,3) = qlm_tzz(i,j) + tt(2,1) = tt(1,2) + tt(3,1) = tt(1,3) + tt(3,2) = tt(2,3) + ll(0) = qlm_l0(i,j,hn) ll(1) = qlm_l1(i,j,hn) ll(2) = qlm_l2(i,j,hn) @@ -140,6 +185,18 @@ subroutine qlm_calc_weyl_scalars (CCTK_ARGUMENTS, hn) ! Calculate 4-metric + call calc_det (gg, dtg) + call calc_inv (gg, dtg, gu) + call calc_invderiv (gu, dgg, dgu) + call calc_connections (gu, dgg, gamma) + call calc_connectionderivs (gu, dgg, dgu, ddgg, dgamma) + call calc_ricci (gamma, dgamma, ri) + + call calc_3metricdot_simple (kk, gg_dot) + call calc_3metricderivdot_simple (dkk, dgg_dot) + call calc_extcurvdot_simple (gg,gu,ri, kk, tt, kk_dot) + call calc_3metricdot2_simple (kk_dot, gg_dot2) + call calc_4metricderivs2_simple (gg, dgg, & ddgg, gg_dot, gg_dot2, dgg_dot, g4,dg4,ddg4) call calc_4inv (g4, gu4) |