diff options
| -rw-r--r-- | src/AHFinder_exp.F | 26 |
1 files changed, 8 insertions, 18 deletions
diff --git a/src/AHFinder_exp.F b/src/AHFinder_exp.F index 880fa94..d93abb5 100644 --- a/src/AHFinder_exp.F +++ b/src/AHFinder_exp.F @@ -5,7 +5,7 @@ @desc Here I calculate the (normalized) expansion of the horizon function. The expansion is defined by the - following expression (a=lapse, b=shift): + following expression: __2 2 / ab __a__b \ exp = \/ f / u + d f d f / u | K - \/ \/ f / u | - trK @@ -69,8 +69,8 @@ ! kdij K ! ij ! -! ddijk 1/2 d g -! i jk +! ddijk d g +! i jk ! ! dfi d f ! i @@ -87,15 +87,16 @@ ! aux,T* auxilliary variables -! *************************************** -! *** DEFINE {zero,half,one,two} *** -! *************************************** +! ******************* +! *** NUMBERS *** +! ******************* zero = 0.0d0 half = 0.5d0 one = 1.0d0 two = 2.0d0 + ! ************************** ! *** FIND EXPANSION *** ! ************************** @@ -134,7 +135,7 @@ kdxz = kxz(i,j,k) kdyz = kyz(i,j,k) -! Find Ds using finite differences. +! Find derivatives of metric using finite differences. ddxxx = idx*(gxx(i+1,j,k) - gxx(i-1,j,k)) ddxyy = idx*(gyy(i+1,j,k) - gyy(i-1,j,k)) @@ -306,17 +307,6 @@ ahf_exp(i,j,k) = T1 + T2 - T3 - T4 -! Check for NaNs. - - if (ahf_exp(i,j,k).ge.zero) then - else if (ahf_exp(i,j,k).lt.zero) then - else - if (myproc.eq.0) then - write(*,*) 'NaN in ahf_exp at point',i,j,k - end if - STOP - end if - 10 continue end do |