diff options
author | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2000-05-03 10:55:14 +0000 |
---|---|---|
committer | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2000-05-03 10:55:14 +0000 |
commit | 77a660590f49e7b3693eb797651603f8e4c26a3b (patch) | |
tree | 02828a0c4e9a3a9f5dc0c48f37a86c55018effae /src/AHFinder_exp.F | |
parent | 8ccf1c90bdc5a6a3ef43150d7173a10bee93c3b7 (diff) |
Cleaning up.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@87 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src/AHFinder_exp.F')
-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 |