diff --git a/src/BC-Ellipsoid.f90 b/src/BC-Ellipsoid.f90
index 91a7b3b68bb57cd35a9132f11a7254f30a6d34d3..5272de5290eec1e7fca7c63b8b8d42f2c69636a2 100644
--- a/src/BC-Ellipsoid.f90
+++ b/src/BC-Ellipsoid.f90
@@ -82,6 +82,7 @@ subroutine geomcomplex_ellip(epsi,nxi,nxf,ny,nyi,nyf,nzi,nzf,dx,yp,dz,remp)
         ym=yp(j)
         if (ym /= ym) then
             write(*,*) "ym = ", ym, " should be ", ((real(j-1,mytype))*dy), ", as j = ", j
+            write(*,*) "yp = ", yp
         endif
 
         do i=nxi,nxf