More self-checking.

git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/SphericalSurface/trunk@12 40f6ab95-0e4f-0410-8daa-ee8d7420be1d

1 files changed, 24 insertions, 0 deletions

diff --git a/src/setup.c b/src/setup.c
index 0744cf8..982bbe9 100644
--- a/src/setup.c
+++ b/src/setup.c
@@ -1,6 +1,7 @@
 /* $Header$ */
 
 #include <assert.h>
+#include <stdio.h>
 
 #include "cctk.h"
 #include "cctk_Arguments.h"
@@ -14,13 +15,36 @@ void SphericalSurfaceInfo_Setup (CCTK_ARGUMENTS)
   DECLARE_CCTK_PARAMETERS;
   
   CCTK_REAL const pi = 3.1415926535897932384626433832795028841971693993751;
+  
+  int group;
+  cGroup groupinfo;
+  cGroupDynamicData groupdata;
+  
   int n;
+  int ierr;
+  
+  
+  
+  if (nsurfaces == 0) return;
+  
+  
+  
+  group = CCTK_GroupIndex ("SphericalSurface::sf_radius");
+  assert (group>=0);
+  ierr = CCTK_GroupData (group, &groupinfo);
+  assert (!ierr);
+  ierr = CCTK_GroupDynamicData (cctkGH, group, &groupdata);
+  assert (!ierr);
   
   
   
   for (n=0; n<nsurfaces; ++n) {
     
     /* internal consistency checks */
+    assert (groupdata.dim == 2);
+    assert (groupdata.gsh[0] >= ntheta[n]);
+    assert (groupdata.gsh[1] >= nphi[n]);
+    
     assert (ntheta[n] >= 3*nghoststheta[n] && ntheta[n] <= maxntheta);
     assert (nphi[n] >= 3*nghostsphi[n] && nphi[n] <= maxnphi);