diff --git a/src/BC-Cylinder.f90 b/src/BC-Cylinder.f90
index 410391a82b727c9642fee9a607ffe858bdbe6fc9..2bcac4c10c4ead5f6e1cb7a1c06ba0f860e335ab 100644
--- a/src/BC-Cylinder.f90
+++ b/src/BC-Cylinder.f90
@@ -82,11 +82,12 @@ 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
                 cycle
              endif
             !  write(*,*) i, j, k
-             epsi(i,j,k)=remp
+              epsi(i,j,k)=remp
           enddo
        enddo
     enddo