store BH diagnostics (= centroid x,y,z, area, m_irreducible)

in Cactus arrays so they're accessible to other thorns and/or I/O methods


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

1 files changed, 32 insertions, 0 deletions

diff --git a/interface.ccl b/interface.ccl
index 87366ce..df8b800 100644
--- a/interface.ccl
+++ b/interface.ccl
@@ -3,3 +3,35 @@
 
 implements: AHFinderDirect
 inherits: Grid ADMBase StaticConformal
+
+################################################################################
+
+#
+# diagnostics computed by this thorn which other thorns may want to use,
+# and/or which we may want to output with Cactus I/O methods
+#
+protected:
+
+#
+# all the remaining diagnostics are arrays subscripted by the
+# (1-origin) horizon number  hn ; subscript 0 is unused
+#
+int BH_diagnostics__int TYPE=array DIM=1 SIZE=N_horizons+1 DISTRIB=constant
+{
+# this is a Boolean flag: 0=false, 1=true
+AH_found	# was this AH found the last time we searched for it?
+		# (if this thorn is active, we search at each time step)
+} "CCTK_INT diagnostics calculated for each apparent horizon found"
+
+#
+# these diagnostics are only defined for those apparent horizons
+# where  AH_found[hn]  is true
+#
+real BH_diagnostics__real TYPE=array DIM=1 SIZE=N_horizons+1 DISTRIB=constant
+{
+# FIXME: it would be nice to gather these together in some way
+centroid_x, centroid_y, centroid_z	# centroid position
+
+area		# area of apparent horizon
+m_irreducible	# irreducible mass = sqrt(area/(16*pi))
+} "CCTK_REAL diagnostics calculated for each apparent horizon found"