diff options
| author | rhaas <rhaas@ef6f4158-a8ee-47d1-ba14-cb73256398e0> | 2014-04-18 04:26:51 +0000 |
|---|---|---|
| committer | rhaas <rhaas@ef6f4158-a8ee-47d1-ba14-cb73256398e0> | 2014-04-18 04:26:51 +0000 |
| commit | f1f2b4771337a195de3c23f440fe32bb3b37689b (patch) | |
| tree | a8fc5f27414c12c8f2b2dc44fa7e6f78e04cdb27 | |
| parent | 9488da761ef86f9356be7acdb7f32735834fe217 (diff) | |
remove commented out code, remove unused variable
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/QuasiLocalMeasures/trunk@47 ef6f4158-a8ee-47d1-ba14-cb73256398e0
| -rw-r--r-- | src/qlm_killing_normalise.F90 | 26 |
1 files changed, 1 insertions, 25 deletions
diff --git a/src/qlm_killing_normalise.F90 b/src/qlm_killing_normalise.F90 index 7b0bb66..253530b 100644 --- a/src/qlm_killing_normalise.F90 +++ b/src/qlm_killing_normalise.F90 @@ -23,7 +23,7 @@ subroutine qlm_killing_normalise (CCTK_ARGUMENTS, hn) logical :: found(ngeodesics) integer :: nsteps - integer :: ii, iii + integer :: ii integer :: i, j character :: msg*1000 @@ -41,29 +41,6 @@ subroutine qlm_killing_normalise (CCTK_ARGUMENTS, hn) if (veryverbose/=0) call CCTK_INFO ("Calculating normalisation factor:") -#if 0 - nsteps = -1 - do iii=1,ngeodesics - ii = (ngeodesics/2+1) + (2*mod(iii,2)-1) * (iii/2) - if (ii<1 .or. ii>ngeodesics) call CCTK_WARN (0, "internal error") - i = 1 + qlm_nghoststheta(hn) & - + ii * (qlm_ntheta(hn) - 2*qlm_nghoststheta(hn) - 1) / (ngeodesics+1) - if (qlm_xi_p(i,j,hn) < 0) then - qlm_xi_t(:,:,hn) = -qlm_xi_t(:,:,hn) - qlm_xi_p(:,:,hn) = -qlm_xi_p(:,:,hn) - qlm_chi (:,:,hn) = -qlm_chi (:,:,hn) - end if - theta = qlm_origin_theta(hn) + (i-1) * qlm_delta_theta(hn) - call killing_factor (CCTK_PASS_FTOF, hn, theta, factor, nsteps) - if (nsteps>0) exit - end do - if (nsteps < 0) then - call CCTK_WARN (1, "Did not manage to integrate along a Killing vector field line loop") - qlm_have_killing_vector(hn) = 0 - factor = 1 - goto 9999 - end if -#else do ii=1,ngeodesics i = 1 + qlm_nghoststheta(hn) & + ii * (qlm_ntheta(hn) - 2*qlm_nghoststheta(hn) - 1) / (ngeodesics+1) @@ -109,7 +86,6 @@ subroutine qlm_killing_normalise (CCTK_ARGUMENTS, hn) factor = 1 goto 9999 end if -#endif |