aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authoreschnett <eschnett@ef6f4158-a8ee-47d1-ba14-cb73256398e0>2010-07-09 21:51:26 +0000
committereschnett <eschnett@ef6f4158-a8ee-47d1-ba14-cb73256398e0>2010-07-09 21:51:26 +0000
commitb4b4c3f252e551adcca11a114a8fa07f778b10b6 (patch)
tree1d2af21e9829c8af8718ab8039b837fb1efc30a8
parentf8e8d9b5834a4b6e84835f86d9ad754780870f69 (diff)
Declare and define "theta"
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/QuasiLocalMeasures/trunk@23 ef6f4158-a8ee-47d1-ba14-cb73256398e0
-rw-r--r--src/qlm_multipoles.F902
1 files changed, 2 insertions, 0 deletions
diff --git a/src/qlm_multipoles.F90 b/src/qlm_multipoles.F90
index 570f216..d8ab9bb 100644
--- a/src/qlm_multipoles.F90
+++ b/src/qlm_multipoles.F90
@@ -35,6 +35,7 @@ subroutine qlm_multipoles (CCTK_ARGUMENTS, hn)
CCTK_REAL :: delta_space(2)
integer :: i, j
+ CCTK_REAL :: theta
!!$ integer :: a, b, c, d
if (veryverbose/=0) then
@@ -65,6 +66,7 @@ subroutine qlm_multipoles (CCTK_ARGUMENTS, hn)
do j = 1+qlm_nghostsphi(hn), qlm_nphi(hn)-qlm_nghostsphi(hn)
do i = 1+qlm_nghoststheta(hn), qlm_ntheta(hn)-qlm_nghoststheta(hn)
+ theta = qlm_origin_theta(hn) + (i-1)*qlm_delta_theta(hn)
! 2-metric on the horizon
qq(1,1) = qlm_qtt(i,j,hn)