diff --git a/src/BC-Ellipsoid.f90 b/src/BC-Ellipsoid.f90
index a919ae23a8563ec72d57be47cd704dd8aa4d228e..711ebac7c7ba92dadcdf5ac9c837f2fe0daebcd2 100644
--- a/src/BC-Ellipsoid.f90
+++ b/src/BC-Ellipsoid.f90
@@ -45,7 +45,9 @@ subroutine geomcomplex_ellip(epsi,nxi,nxf,ny,nyi,nyf,nzi,nzf,dx,yp,dz,remp)
     zeromach = ten*zeromach
     is_inside=.false.
     !  orientation=[oriw, orii, orij, orik]
-    call NormalizeQuaternion(orientation)
+    do i = 1,nbody 
+        call NormalizeQuaternion(orientation(i,:))
+    enddo
     !  shape=[shx, shy, shz]
     !  write(*,*) shape, 'SHAPE'