diff options
-rw-r--r-- | src/RadiationBoundary.c | 15 |
1 files changed, 13 insertions, 2 deletions
diff --git a/src/RadiationBoundary.c b/src/RadiationBoundary.c index a076f4e..b754b6f 100644 --- a/src/RadiationBoundary.c +++ b/src/RadiationBoundary.c @@ -80,6 +80,8 @@ #include "Symmetry.h" #include "Boundary.h" +/* #define DEBUG */ + /* the rcs ID and its dummy function to use it */ static const char *rcsid = "$Header$"; CCTK_FILEVERSION(CactusBase_Boundary_RadiationBoundary_c); @@ -154,7 +156,16 @@ int BndRadiative(const cGH *GH, int num_vars, int *vars, int *faces, int prev_time_level; /* variable (index) which holds the previous time level */ #ifdef DEBUG - printf("BndRadiative(): got passed GH=%p, num_vars=%d, var_indices[0]=%d, table_handles[0]=%d\n", (const void *) GH, num_vars, var_indices[0], table_handles[0]); + printf("BndRadiative() got passed: GH=%p, num_vars=%d:\n", (const void *) GH, + num_vars); + printf("var index var name table handle\n"); + for (i=0; i<num_vars; ++i) + { + printf("%d %12s %d\n", vars[i], CCTK_VarName(vars[i]), tables[i]); + } + printf("end of table\n"); + + /* CCTK_WARN(0, "stopping code"); */ #endif retval = 0; stencil_alldirs = NULL; max_gdim = 0; @@ -1640,7 +1651,7 @@ static int ApplyBndRadiative (const cGH *GH, + boundary is a physical boundary + have enough grid points */ - memset (doBC, 1, sizeof (doBC)); + /* memset (doBC, 1, sizeof (doBC)); doBC is initialized below */ if (sGHex) { for (i = 0; i < 2 * MAXDIM; i++) |