add compuatation and printing of mean radius for each AH found

output format for BH info in stdout log is now

INFO (AHFinderDirect): AH 1/1: r=1.86525 at (-0.000000,0.000000,-0.000000)
INFO (AHFinderDirect): AH 1/1: area=45.23876324 m_irreducible=0.9486815054

and in BH_diagnostics file is now

# apparent horizon 1 of 1
#
# column 1 = cctk_iteration
# column 2 = cctk_time
# column 3 = centroid_x
# column 4 = centroid_y
# column 5 = centroid_z
# column 6 = mean radius
# column 7 = area
# column 8 = m_irreducible


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@879 f88db872-0e4f-0410-b76b-b9085cfa78c5

1 files changed, 4 insertions, 3 deletions

diff --git a/src/patch/patch_system.cc b/src/patch/patch_system.cc
index da4027c..0aefef6 100644
--- a/src/patch/patch_system.cc
+++ b/src/patch/patch_system.cc
@@ -1560,12 +1560,13 @@ norms.reset();
 // overlap, i.e. that patch_overlap_width=1
 //
 // Arguments:
-// src_gfn = (in) The gridfn to be integrated.
+// unknown_src_gfn = (in) The gridfn to be integrated.  This may be
+//			  either nominal-grid or ghosted-grid.
 // ghosted_radius_gfn = (in) The surface radius.
 // g_{xx,xy,xz,yy,yz,zz}_gfn = (in) The xyz 3-metric components.
 // method = (in) Selects the integration scheme.
 //
-fp patch_system::integrate_gridfn(int src_gfn,
+fp patch_system::integrate_gridfn(int unknown_src_gfn,
 				  int ghosted_radius_gfn,
 				  int g_xx_gfn, int g_xy_gfn, int g_xz_gfn,
 						int g_yy_gfn, int g_yz_gfn,
@@ -1577,7 +1578,7 @@ fp integral = 0.0;
 	for (int pn = 0 ; pn < N_patches() ; ++pn)
 	{
 	const patch& p = ith_patch(pn);
-	integral += p.integrate_gridfn(src_gfn,
+	integral += p.integrate_gridfn(unknown_src_gfn,
 				       ghosted_radius_gfn,
 				       g_xx_gfn, g_xy_gfn, g_xz_gfn,
 						 g_yy_gfn, g_yz_gfn,