From af9d4aacef8eb641775018def1ed55e5aee42968 Mon Sep 17 00:00:00 2001
From: ceschoon <ceschoon@ulb.ac.be>
Date: Fri, 3 Nov 2023 14:42:24 +0100
Subject: [PATCH] Fix in get_ave_background_density()

---
 src/Density.cpp | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/src/Density.cpp b/src/Density.cpp
index 848b5cf..115aee4 100644
--- a/src/Density.cpp
+++ b/src/Density.cpp
@@ -455,11 +455,11 @@ double Density::get_ave_background_density() const
       {d += get(ix,iy,0); n++;}
 
   for(int ix=0;ix<Nx_;ix++)
-    for(int iz=0;iz<Nz_;iz++)
+    for(int iz=1;iz<Nz_;iz++) // Cedric: start from Nz=1 so that the Ny=Nz=0 boundary points are counted twice in the average. Same remark below.
       {d += get(ix,0,iz); n++;}
 
-  for(int iy=0;iy<Ny_;iy++)
-    for(int iz=0;iz<Nz_;iz++)    
+  for(int iy=1;iy<Ny_;iy++)
+    for(int iz=1;iz<Nz_;iz++)    
       {d += get(0,iy,iz); n++;}