diff options
author | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2000-03-30 13:18:54 +0000 |
---|---|---|
committer | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2000-03-30 13:18:54 +0000 |
commit | 07a82a6eac4167e16cb04d9cbf99841e2e9a79d5 (patch) | |
tree | fdef569f3cce63d1d49f08d4bfb327a9476eff97 /src/AHFinder_pow.F | |
parent | 1b675444b0f3b1b478db0e4143fdcbbdba96d396 (diff) |
Small chnages to output and some extra comments.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@57 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src/AHFinder_pow.F')
-rw-r--r-- | src/AHFinder_pow.F | 33 |
1 files changed, 28 insertions, 5 deletions
diff --git a/src/AHFinder_pow.F b/src/AHFinder_pow.F index 1e2794a..b1d3455 100644 --- a/src/AHFinder_pow.F +++ b/src/AHFinder_pow.F @@ -881,33 +881,44 @@ call AHFinder_int(CCTK_FARGUMENTS) +! Messages to screen and logfile. + if (myproc.eq.0) then - if (intarea.ne.zero) then +! Find inverse of area. Be carefull not to divide by zero! +! Also, a value of 1.0D10 indicates an error in the integration +! routine, and all integrals end up with the same value. +! So don't divide by it or I will be confused later. + + if ((intarea.ne.zero).and.(intarea.ne.1.0D10)) then aux = one/intarea else aux = one end if +! Screen. + if (veryver) then write(*,*) write(*,"(A21,ES14.6)") ' Surface area =',intarea write(*,"(A21,ES14.6)") ' Mean value of H =',aux*intexp write(*,"(A21,ES14.6)") ' Mean value of H^2 =',aux*intexp2 - write(*,*) 'number of interpolated points: ', + write(*,*) 'Number of interpolated points: ', . inside_min_count - write(*,*) 'number of those that are negative: ', + write(*,*) 'Number of those that are negative: ', . inside_min_neg_count end if +! Logfile. + if (logfile) then write(1,*) write(1,"(A21,ES14.6)") ' Surface area =',intarea write(1,"(A21,ES14.6)") ' Mean value of H =',aux*intexp write(1,"(A21,ES14.6)") ' Mean value of H^2 =',aux*intexp2 - write(1,*) 'number of interpolated points: ', + write(1,*) 'Number of interpolated points: ', . inside_min_count - write(1,*) 'number of those that are negative: ', + write(1,*) 'Number of those that are negative: ', . inside_min_neg_count end if @@ -920,6 +931,8 @@ if (myproc.eq.0) then +! Headers to screen. + if (veryver) then write(*,*) @@ -939,6 +952,8 @@ end if +! Headers to logfile. + if (logfile) then write(1,*) @@ -958,6 +973,8 @@ end if +! c0 coefficients. + do l=1+stepz,lmax,1+stepz if (veryver) then write(*,"(I4,A6,ES14.6)") l,' ',c0(l) @@ -967,8 +984,12 @@ end if end do +! cc and cs coefficients (only for nonaxisymmetric case). + if (nonaxi) then +! Reflection symmetry on y: some coefficients are not there. + if (.not.refy) then if (veryver) then @@ -997,6 +1018,8 @@ end do end do +! No reflection symmetry on y. + else if (veryver) then |