From f7a1c28a5dc81549e59915b896272f946a454ad4 Mon Sep 17 00:00:00 2001
From: Pedro Costa
Date: Fri, 15 Nov 2024 13:19:27 +0100
Subject: [PATCH] Fix oob issue.
And removed superfluous `update self` directive.
---
src/forcing.f90 | 2 +-
src/main.f90 | 1 -
2 files changed, 1 insertion(+), 2 deletions(-)
diff --git a/src/forcing.f90 b/src/forcing.f90
index b4b26c6..3cda469 100644
--- 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
index 04c41fa..197b9c7 100644
--- 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)