From 1cde12cf5a24103aae83c7e54ca6eb242213e472 Mon Sep 17 00:00:00 2001
From: Andrew Boyd <adfboyd@gmail.com>
Date: Mon, 18 Sep 2023 20:27:17 +0100
Subject: [PATCH] resolved z-coord problem
---
src/BC-Cylinder.f90 | 12 +++++++-----
1 file changed, 7 insertions(+), 5 deletions(-)
diff --git a/src/BC-Cylinder.f90 b/src/BC-Cylinder.f90
index 5778bf74..ce9b0d6c 100644
--- a/src/BC-Cylinder.f90
+++ b/src/BC-Cylinder.f90
@@ -72,22 +72,24 @@ contains
!
! Define adjusted smoothing constant
! kcon = log((one-0.0001)/0.0001)/(smoopar*0.5*dx) ! 0.0001 is the y-value, smoopar: desired number of affected points
-!
+ write(*,*) nzi, nzf
do k=nzi,nzf
- zm=(real(nzi+k-2,mytype))*dz
+ zm=(real(k-1,mytype))*dz
+ ! write(*,*) k, zm
do j=nyi,nyf
ym=yp(j)
do i=nxi,nxf
xm=real(i-1,mytype)*dx
point=[xm, ym, zm]
! call EllipsoidalRadius(point, ce, orientation, shape, r)
- ! r=sqrt_prec((xm-cexx)**two+(ym-ceyy)**two+(zm-cezz)**two)
- r=sqrt_prec((xm-cexx)**two+(ym-ceyy)**two)
+ r=sqrt_prec((xm-cexx)**two+(ym-ceyy)**two+(zm-cezz)**two)
+ ! r=sqrt_prec((xm-cexx)**two+(ym-ceyy)**two)
if (r-ra.gt.zeromach) then
+ ! write(*,*) i, j, k
cycle
endif
- ! write(*,*) i, j, k
+ ! write(*,*) i, j, k, zm
epsi(i,j,k)=remp
! write(*,*) remp
enddo
--
GitLab