final final fix...

src/ellip_utils.f90

@@ -221,8 +221,8 @@ contains
       if (radius /= radius) then
         write(*,*) "Got an error in grid check!"
         write(*,*) "Radius = ", radius
-        write(*,*) "point = ", centre
-        write(*,*) "Body centre = ", position
+        write(*,*) "point = ", point
+        write(*,*) "Body centre = ", centre
         write(*,*) "Translated point = ", trans_point
         write(*,*) "Orientation = ", orientation
         write(*,*) "Rotated point = ", rotated_point