From e943a66a52228669b89c7403ee8c1eb1b56900a0 Mon Sep 17 00:00:00 2001
From: knarf <knarf@54511f98-0e4f-0410-826e-eb8b393f5a1e>
Date: Tue, 23 Aug 2011 18:58:36 +0000
Subject: initialize ADMMass_SurfaceMass_GF completely

git-svn-id: http://svn.aei.mpg.de/numrel/AEIThorns/ADMMass/trunk@41 54511f98-0e4f-0410-826e-eb8b393f5a1e
---
 src/surface_integral.c | 8 +++++++-
 1 file changed, 7 insertions(+), 1 deletion(-)

(limited to 'src')

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) &&
-- 
cgit v1.2.3