From 0b09803657d99c3d57c0bbfab2a3e949b73a2d73 Mon Sep 17 00:00:00 2001
From: s2006749 <s2006749@ed.ac.uk>
Date: Tue, 22 Oct 2024 15:37:04 +0100
Subject: [PATCH] tidierer inputs
---
src/BC-Ellipsoid.f90 | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/src/BC-Ellipsoid.f90 b/src/BC-Ellipsoid.f90
index d95e22b0..c8d72d10 100644
--- a/src/BC-Ellipsoid.f90
+++ b/src/BC-Ellipsoid.f90
@@ -23,7 +23,7 @@ subroutine geomcomplex_ellip(epsi,nxi,nxf,ny,nyi,nyf,nzi,nzf,dx,yp,dz,remp)
use ibm_param
use dbg_schemes, only: sqrt_prec
use ellipsoid_utils, only: NormalizeQuaternion, EllipsoidalRadius, EllipsoidalRadius_debug
- use complex_geometry, only: nraf
+ use complex_geometry, only: nraf,nyraf
implicit none
@@ -83,7 +83,9 @@ 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, " (or maybe it should be ", ((real(j-1,mytype))*dy)/real(nraf,mytype), ")"
- write(*,*) "yp = ", yp
+ if (j.lt.(nyraf)) then
+ write(*,*) "yp(j-1) = ", yp(j-1), " yp(j+1) = ", yp(j+1)
+ endif
endif
do i=nxi,nxf
--
GitLab