From bf8e6181a0f63313a1ade4591bf07cd7e098b45c Mon Sep 17 00:00:00 2001
From: s2006749 <s2006749@ed.ac.uk>
Date: Wed, 23 Oct 2024 13:45:47 +0100
Subject: [PATCH] added buffer to ypraf
---
src/BC-Ellipsoid.f90 | 14 +++++++-------
src/genepsi3d.f90 | 2 +-
2 files changed, 8 insertions(+), 8 deletions(-)
diff --git a/src/BC-Ellipsoid.f90 b/src/BC-Ellipsoid.f90
index c8d72d10..0abe7501 100644
--- a/src/BC-Ellipsoid.f90
+++ b/src/BC-Ellipsoid.f90
@@ -84,7 +84,7 @@ subroutine geomcomplex_ellip(epsi,nxi,nxf,ny,nyi,nyf,nzi,nzf,dx,yp,dz,remp)
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), ")"
if (j.lt.(nyraf)) then
- write(*,*) "yp(j-1) = ", yp(j-1), " yp(j+1) = ", yp(j+1)
+ write(*,*) "yp(j-1) = ", yp(j-1), " yp(j+1,..) = ", yp(j+1:)
endif
endif
@@ -99,12 +99,12 @@ subroutine geomcomplex_ellip(epsi,nxi,nxf,ny,nyi,nyf,nzi,nzf,dx,yp,dz,remp)
! 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
+ ! 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
is_inside = (r-ra(i_body)).lt.zeromach
! if (is_inside) then
diff --git a/src/genepsi3d.f90 b/src/genepsi3d.f90
index 51b0092a..6bd7f4b8 100644
--- a/src/genepsi3d.f90
+++ b/src/genepsi3d.f90
@@ -276,7 +276,7 @@ contains
real(mytype),dimension(ysize(1),nyraf,ysize(3)) :: yepsi
real(mytype),dimension(zsize(1),zsize(2),nzraf) :: zepsi
real(mytype),dimension(ny) :: yp
- real(mytype),dimension(nyraf) :: ypraf
+ real(mytype),dimension(nyraf+20) :: ypraf
real(mytype) :: dxraf,dyraf,dzraf
integer :: i,j,k
integer :: ii,jj,kk
--
GitLab