diff --git a/src/BC-Cylinder.f90 b/src/BC-Cylinder.f90
index 97a7a353a43ae3ad51e31b9c79c13a801a56cde4..e62c1744856895529e41cabaf49413f208318a98 100644
--- a/src/BC-Cylinder.f90
+++ b/src/BC-Cylinder.f90
@@ -84,6 +84,8 @@ contains
              point=[xm, ym, zm]
              call EllipsoidalRadius(point, ce, orientation, shape, r)
             !  r=sqrt_prec((xm-cexx)**two+(ym-ceyy)**two+(zm-cezz)**two)
+            !  r=sqrt_prec((xm-cexx)**two+(ym-ceyy)**two)
+
              if (r-ra.gt.zeromach) then
                !  write(*,*) i, j, k
                 cycle