diff options
Diffstat (limited to 'interface.ccl')
| -rw-r--r-- | interface.ccl | 71 |
1 files changed, 49 insertions, 22 deletions
diff --git a/interface.ccl b/interface.ccl index 798f96b..5ca0d73 100644 --- a/interface.ccl +++ b/interface.ccl @@ -2,17 +2,20 @@ # $Header$ implements: AHFinderDirect -# thorns from: CactusBase CactusEinstein ................. AEIThorns -inherits: Grid IO ADMBase StaticConformal SpaceMask SphericalSurface +inherits: Grid ADMBase StaticConformal SpaceMask SphericalSurface IO -# we use an include file provided by SpaceMask +# include file provided by SpaceMask USES INCLUDE: SpaceMask.h ################################################################################ -# -# ***** aliased functions used from other thorns ***** -# +# we use a mask in order to find out whether we interpolate from the excised +# region. (we need 3 time levels because Carpet might need to interpolate +# in time.) + +REAL ahmask TYPE=gf TIMELEVELS=3 TAGS='tensortypealias="Scalar" Prolongation="None"' + +################################################################################ # we use an aliased function maybe provided by DriftCorrect SUBROUTINE SetDriftCorrectPosition \ @@ -20,20 +23,9 @@ SUBROUTINE SetDriftCorrectPosition \ CCTK_REAL IN x, CCTK_REAL IN y, CCTK_REAL IN z) USES FUNCTION SetDriftCorrectPosition -# -# This IOUtil function is used to determine whether the diagnostics file -# should be truncated or appended (as in a recovery run ) -# -CCTK_INT FUNCTION IO_TruncateOutputFiles(CCTK_POINTER_TO_CONST IN GH) -USES FUNCTION IO_TruncateOutputFiles - ################################################################################ # -# ***** aliased functions provided by this thorn ***** -# - -# # This thorn provides an aliased-function interface which other thorns # can use to find out about the horizons. Each function returns a status # code which is >= 0 for ok, or < 0 for an error. @@ -83,11 +75,6 @@ CCTK_INT FUNCTION HorizonCentroid \ # or -1.0 if this horizon wasn't found the most recent time AHFinderDirect # searched for it. # -# If this function is to be used in a multiprocessor run on a horizon -# which was found on some other processor, the parameter -# AHFinderDirect::always_broadcast_horizon_shape -# must be set to true to get the correct answer. -# CCTK_INT FUNCTION HorizonRadiusInDirection \ (CCTK_INT IN horizon_number, \ CCTK_INT IN N_points, \ @@ -111,3 +98,43 @@ PROVIDES FUNCTION HorizonRadiusInDirection \ LANGUAGE C ################################################################################ + +# Grid arrays to store the horizon shape, for easy output and checkpointing. +# These are similar and equivalent to SphericalSurface. +# Note that the horizons are numbered [0...N_horizon-1]. + +CCTK_REAL ah_radius[N_horizons] TYPE=array DIM=3 SIZE=max_N_zones_per_right_angle+1,max_N_zones_per_right_angle+1,6 DISTRIB=constant + +CCTK_REAL ah_origin[N_horizons] TYPE=scalar +{ + ah_origin_x ah_origin_y ah_origin_z +} + +CCTK_REAL ah_centroid[N_horizons] TYPE=scalar +{ + ah_centroid_x ah_centroid_y ah_centroid_z ah_centroid_t + ah_centroid_x_p ah_centroid_y_p ah_centroid_z_p ah_centroid_t_p +} + +CCTK_INT ah_flags[N_horizons] TYPE=scalar +{ + ah_initial_find_flag + ah_really_initial_find_flag + ah_search_flag + ah_found_flag + ah_centroid_valid + ah_centroid_valid_p + ah_centroid_iteration + ah_centroid_iteration_p +} + +################################################################################ + +# This IOUtil function is used to determine whether the diagnostics file +# should be truncated or appended (as in a recovery run ) + +CCTK_INT FUNCTION IO_TruncateOutputFiles \ + (CCTK_POINTER_TO_CONST IN GH) +USES FUNCTION IO_TruncateOutputFiles + +################################################################################ |