diff options
-rw-r--r-- | interface.ccl | 6 | ||||
-rw-r--r-- | schedule.ccl | 2 | ||||
-rw-r--r-- | src/radius.c | 2 | ||||
-rw-r--r-- | src/setup.c | 1 |
4 files changed, 11 insertions, 0 deletions
diff --git a/interface.ccl b/interface.ccl index 03af59d..833b7e3 100644 --- a/interface.ccl +++ b/interface.ccl @@ -8,6 +8,12 @@ IMPLEMENTS: SphericalSurface PUBLIC: # To be defined by the setting thorn. +# Suggested meanings: +# zero: The surface data are invalid (garbage) +# one: The surface data are valid in principle, although they may be outdated +CCTK_INT sf_active[nsurfaces] TYPE=scalar + +# To be defined by the setting thorn. # Suggested meanings: # zero: No surface is provided at this time step # negative: No surface could be found at this time step diff --git a/schedule.ccl b/schedule.ccl index ad7de7b..a4aac33 100644 --- a/schedule.ccl +++ b/schedule.ccl @@ -1,6 +1,7 @@ # Schedule definitions for thorn SphericalSurface # $Header$ +STORAGE: sf_active STORAGE: sf_valid STORAGE: sf_info STORAGE: sf_radius sf_origin @@ -12,6 +13,7 @@ STORAGE: sf_shape_descriptors SCHEDULE SphericalSurface_Setup AT basegrid { LANG: C + OPTIONS: global } "Calculate surface coordinate descriptors" diff --git a/src/radius.c b/src/radius.c index 0356259..0958b4c 100644 --- a/src/radius.c +++ b/src/radius.c @@ -40,6 +40,7 @@ void SphericalSurface_Set (CCTK_ARGUMENTS) if (set_spherical[n]) { + sf_active[n] = +1; sf_valid[n] = +1; sf_area[n] = 4 * pi * pow (radius[n], 2); @@ -103,6 +104,7 @@ void SphericalSurface_Set (CCTK_ARGUMENTS) CCTK_REAL const ry2 = pow (radius_y[n], 2); CCTK_REAL const rz2 = pow (radius_z[n], 2); + sf_active[n] = +1; sf_valid[n] = +1; sf_area[n] = 0 * 4 * pi * pow (radius[n], 2); diff --git a/src/setup.c b/src/setup.c index 3bc37d9..d4cd599 100644 --- a/src/setup.c +++ b/src/setup.c @@ -115,6 +115,7 @@ void SphericalSurface_Setup (CCTK_ARGUMENTS) /* mark surface as uninitialised */ + sf_active[n] = 0; sf_valid[n] = 0; } /* for n */ |