From f1df1240008be53ef7170456ba5f8685105365f0 Mon Sep 17 00:00:00 2001
From: s2006749 <s2006749@ed.ac.uk>
Date: Tue, 22 Oct 2024 15:25:12 +0100
Subject: [PATCH] tidier outputs
---
src/genepsi3d.f90 | 5 ++++-
1 file changed, 4 insertions(+), 1 deletion(-)
diff --git a/src/genepsi3d.f90 b/src/genepsi3d.f90
index 1d925c69..51b0092a 100644
--- a/src/genepsi3d.f90
+++ b/src/genepsi3d.f90
@@ -317,10 +317,13 @@ contains
do j=1,ny-1
do jraf=1,nraf
ypraf(jraf+nraf*(j-1))=yp(j)+real(jraf-1, mytype)*(yp(j+1)-yp(j))/real(nraf, mytype)
- write(*,*) "At j = ", j, ", jraf = ", jraf, "ypraf = ", yp(j)+real(jraf-1, mytype)*(yp(j+1)-yp(j))/real(nraf, mytype)
+ if (ypraf(jraf+nraf*(j-1)) /= ypraf(jraf+nraf*(j-1))) then
+ write(*,*) "At j = ", j, ", jraf = ", jraf, "ypraf = ", yp(j)+real(jraf-1, mytype)*(yp(j+1)-yp(j))/real(nraf, mytype)
+ endif
enddo
enddo
if(.not.ncly)ypraf(nyraf)=yp(ny)
+ if(.not.ncly)write(*,*) "Changed ypraf (", nyraf, "). To ", yp(ny)
yepsi=zero
call geomcomplex(yepsi,ystart(1),yend(1),nyraf,1,nyraf,ystart(3),yend(3),dx,ypraf,dz,one)
! if (nrank==0) print*,' step 3'
--
GitLab