From 6de9114dc18766e16d5e33aa357160e10715d9ed Mon Sep 17 00:00:00 2001
From: s2006749 <s2006749@ed.ac.uk>
Date: Tue, 22 Oct 2024 14:51:38 +0100
Subject: [PATCH] modified debug
---
src/BC-Ellipsoid.f90 | 11 +++++++++--
1 file changed, 9 insertions(+), 2 deletions(-)
diff --git a/src/BC-Ellipsoid.f90 b/src/BC-Ellipsoid.f90
index a4f616d6..53dce403 100644
--- a/src/BC-Ellipsoid.f90
+++ b/src/BC-Ellipsoid.f90
@@ -85,11 +85,18 @@ subroutine geomcomplex_ellip(epsi,nxi,nxf,ny,nyi,nyf,nzi,nzf,dx,yp,dz,remp)
! call EllipsoidalRadius(point, position, orientation, shape, r)
do i_body = 1,nbody
if (cube_flag.eq.0) then
- if ((nrank.eq.0).and.(torq_debug.eq.1)) then
+ ! if ((nrank.eq.0).and.(torq_debug.eq.1)) then
+ ! write(*,*) "Point, position, orientation, shape for body", i_body, " = "
+ ! write(*,*) point,position(i_body,:),orientation(i_body,:),shape(i_body,:)
+ ! endif
+ call EllipsoidalRadius(point,position(i_body,:),orientation(i_body,:),shape(i_body,:),r)
+ if (r /= r) then
write(*,*) "Point, position, orientation, shape for body", i_body, " = "
write(*,*) point,position(i_body,:),orientation(i_body,:),shape(i_body,:)
+ write(*,*) "R calculated = ", r
+ write(*,*) "Timestep = ", itime
endif
- call EllipsoidalRadius(point,position(i_body,:),orientation(i_body,:),shape(i_body,:),r)
+
is_inside = (r-ra(i_body)).lt.zeromach
! if (is_inside) then
! call EllipsoidalRadius_debug(point,position(i_body,:),orientation(i_body,:),shape(i_body,:),r)
--
GitLab