diff options
author | schnetter <schnetter@ef6f4158-a8ee-47d1-ba14-cb73256398e0> | 2010-04-07 17:00:34 +0000 |
---|---|---|
committer | schnetter <schnetter@ef6f4158-a8ee-47d1-ba14-cb73256398e0> | 2010-04-07 17:00:34 +0000 |
commit | fd5a4312226842434f81a445910e3284b7a77b43 (patch) | |
tree | a8ef55912cb63293d294f9b6df96ac155a28c1a4 /src/qlm_coordinates.F90 | |
parent | b1d7f9d2ff5a68aac30dcb1666121c91f01fc15b (diff) |
Correct nans
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/QuasiLocalMeasures/trunk@5 ef6f4158-a8ee-47d1-ba14-cb73256398e0
Diffstat (limited to 'src/qlm_coordinates.F90')
-rw-r--r-- | src/qlm_coordinates.F90 | 10 |
1 files changed, 7 insertions, 3 deletions
diff --git a/src/qlm_coordinates.F90 b/src/qlm_coordinates.F90 index 6807bcd..17b7b0c 100644 --- a/src/qlm_coordinates.F90 +++ b/src/qlm_coordinates.F90 @@ -9,6 +9,7 @@ subroutine qlm_calc_coordinates (CCTK_ARGUMENTS, hn) use cctk use constants use qlm_boundary + use tensor2 implicit none DECLARE_CCTK_ARGUMENTS DECLARE_CCTK_FUNCTIONS @@ -43,7 +44,8 @@ subroutine qlm_calc_coordinates (CCTK_ARGUMENTS, hn) qq(1,2) = qlm_qtp(i,j,hn) qq(2,2) = qlm_qpp(i,j,hn) qq(2,1) = qq(1,2) - dtq = qq(1,1) * qq(2,2) - qq(1,2) * qq(2,1) + + call calc_2det (qq, dtq) area = area + sqrt(dtq) * qlm_delta_theta(hn) * qlm_delta_phi(hn) @@ -98,7 +100,8 @@ subroutine qlm_calc_coordinates (CCTK_ARGUMENTS, hn) qq(1,2) = qlm_qtp(i,j,hn) qq(2,2) = qlm_qpp(i,j,hn) qq(2,1) = qq(1,2) - dtq = qq(1,1) * qq(2,2) - qq(1,2) * qq(2,1) + + call calc_2det (qq, dtq) integral_z = integral_z + qlm_inv_z(i,j,hn) * sqrt(dtq) * qlm_delta_theta(hn) * qlm_delta_phi(hn) @@ -172,7 +175,8 @@ contains qq(1,2) = qlm_qtp(i,j,hn) qq(2,2) = qlm_qpp(i,j,hn) qq(2,1) = qq(1,2) - dtq = qq(1,1) * qq(2,2) - qq(1,2) * qq(2,1) + + call calc_2det (qq, dtq) zdot = vv(1) * (- (1/radius**2) * qlm_xi_p(i,j,hn) * sqrt(dtq)) & + vv(2) * (+ (1/radius**2) * qlm_xi_t(i,j,hn) * sqrt(dtq)) |