diff --git a/src/ellip_utils.f90 b/src/ellip_utils.f90 index 83fae98747793698e68a1d4005f7c457f02942f7..17a044016ad22bc9229d8441b95c786ff30227e7 100644 --- a/src/ellip_utils.f90 +++ b/src/ellip_utils.f90 @@ -251,7 +251,7 @@ contains real(mytype) :: r call EllipsoidalRadius(point,centre,orientation,shape,r) - + is_inside = ((r-ra).lt.zeromach) end subroutine is_inside_ellipsoid