diff --git a/src/BC-Ellipsoid.f90 b/src/BC-Ellipsoid.f90
index 53dce403b8f8992f4a489fc95e1eb8ec550c414b..076c3e7699f0bff478ef3b4c5ab60ced2d0b6312 100644
--- a/src/BC-Ellipsoid.f90
+++ b/src/BC-Ellipsoid.f90
@@ -55,6 +55,7 @@ subroutine geomcomplex_ellip(epsi,nxi,nxf,ny,nyi,nyf,nzi,nzf,dx,yp,dz,remp)
     ! Intitialise epsi
     epsi(:,:,:)=zero
 
+    
 
     ! Update center of moving ellipsoid
     ! if (t.ne.0.) then
@@ -79,6 +80,10 @@ subroutine geomcomplex_ellip(epsi,nxi,nxf,ny,nyi,nyf,nzi,nzf,dx,yp,dz,remp)
         do j=nyi,nyf
         ! ym=(real(j-1,mytype))*dy
         ym=yp(j)
+        if (ym /= ym) then
+            write(*,*) "ym = ", ym, " should be ", ((real(j-1,mytype))*dy), ", as j = ", jj
+        endif
+
         do i=nxi,nxf
             xm=real(i-1,mytype)*dx
             point=[xm, ym, zm]