diff --git a/boundaries/boundaryhandling.py b/boundaries/boundaryhandling.py
index 4f2a0c611ba41136f01fb76e924e683d41a9192c..fb9225dfaf79f53fa756044041b4447b8da8b014 100644
--- a/boundaries/boundaryhandling.py
+++ b/boundaries/boundaryhandling.py
@@ -91,7 +91,7 @@ class LbmWeightInfo(CustomCppCode):
     def __init__(self, lb_method):
         weights = [str(w.evalf()) for w in lb_method.weights]
         w_sym = LbmWeightInfo.WEIGHTS_SYMBOL
-        code = "const double __attribute__((unused)) %s [] = { %s };\n" % (w_sym.name, ",".join(weights))
+        code = "const double %s [] = { %s };\n" % (w_sym.name, ",".join(weights))
         super(LbmWeightInfo, self).__init__(code, symbols_read=set(), symbols_defined={w_sym})
 
 
diff --git a/geometry.py b/geometry.py
index 80ae249786b27597be612b9c517e727c74d17033..e0cec5ae6728c2e8c17452d1ae3f76916d547ae6 100644
--- a/geometry.py
+++ b/geometry.py
@@ -253,7 +253,7 @@ def add_black_and_white_image(boundary_handling, image_file, target_slice=None,
             image_target_slice = shift_slice(intersection_slice, [-s.start for s in image_slice])
             mask_target_slice = [mask_target_slice[i] if i in plane else slice(None, None) for i in range(dim)]
             image_target_slice = [image_target_slice[i] if i in plane else np.newaxis for i in range(dim)]
-            result[tuple(mask_target_slice)] = zoomed_image[image_target_slice]
+            result[tuple(mask_target_slice)] = zoomed_image[tuple(image_target_slice)]
             return result
 
     return boundary_handling.set_boundary(boundary, slice_obj=image_slice, mask_callback=callback,