diff options
Diffstat (limited to 'src/surface_integral.c')
| -rw-r--r-- | src/surface_integral.c | 8 |
1 files changed, 7 insertions, 1 deletions
diff --git a/src/surface_integral.c b/src/surface_integral.c index fa14602..6886981 100644 --- a/src/surface_integral.c +++ b/src/surface_integral.c @@ -190,12 +190,18 @@ void ADMMass_Surface(CCTK_ARGUMENTS) z_max_k = find_closest(cctkGH, cctk_lsh, cctk_delta_space, ghost, z, *ADMMass_box_z_max, 2); + for(i=0; i<cctk_lsh[0]; i++) + for(j=0; j<cctk_lsh[1]; j++) + for(k=0; k<cctk_lsh[2]; k++) + { + ijk = CCTK_GFINDEX3D(cctkGH, i, j, k); + ADMMass_SurfaceMass_GF[ijk] = 0.0; + } for(i=ghost; i<cctk_lsh[0]-ghost; i++) for(j=ghost; j<cctk_lsh[1]-ghost; j++) for(k=ghost; k<cctk_lsh[2]-ghost; k++) { ijk = CCTK_GFINDEX3D(cctkGH, i, j, k); - ADMMass_SurfaceMass_GF[ijk] = 0.0; /* delimit the cube on whose surface we want to integrate */ if ((i >= x_min_i) && |