Merge pull request #97 from CaNS-World/fix-hit-oob

Fix out-of-bounds issue in HIT forcing.
CaNS-World · Nov 15, 2024 · 8fa1d31 · 8fa1d31
2 parents 4c0cec4 + f7a1c28
commit 8fa1d31
Show file tree

Hide file tree

Showing 2 changed files with 1 addition and 2 deletions.
diff --git a/src/forcing.f90 b/src/forcing.f90
@@ -20,7 +20,7 @@ subroutine lscale_forcing(ftype,lo,hi,alpha,dt,l,dl,zc,zf,u,v,w)
     integer , intent(in), dimension(3) :: lo,hi
     real(rp), intent(in) :: alpha,dt
     real(rp), intent(in), dimension(3) :: l,dl
-    real(rp), intent(in), dimension(0:) :: zc,zf
+    real(rp), intent(in), dimension(lo(3)-1:) :: zc,zf
     real(rp), intent(inout), dimension(lo(1)-1:,lo(2)-1:,lo(3)-1:) :: u,v,w
     !
     real(rp) :: f0,factor

diff --git a/src/main.f90 b/src/main.f90
@@ -289,7 +289,6 @@ program cans
   if(.not.restart) then
     istep = 0
     time = 0.
-    !$acc update self(zc,dzc,dzf)
     call initflow(inivel,bcvel,ng,lo,l,dl,zc,zf,dzc,dzf,rho12(2),mu12(2),bforce,is_wallturb,time,u,v,w,p)
 #if defined(_SCALAR)
     call initscal(inisca,bcsca,ng,lo,l,dl,dzf,zc,s)