From 63ff01a2be81c99455e54e37d2920125ad1f1e6a Mon Sep 17 00:00:00 2001 From: Alberto Ramos Date: Tue, 31 Jan 2023 15:50:51 +0100 Subject: [PATCH] Corrected bug for setting boundary fields Thanks to Fernando Perez Panadero --- src/YM/YMsf.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/YM/YMsf.jl b/src/YM/YMsf.jl index f882536..6ecc916 100644 --- a/src/YM/YMsf.jl +++ b/src/YM/YMsf.jl @@ -107,7 +107,7 @@ function krnl_setbnd_it0!(U, phi1, phi2, lp::SpaceParm{N,M,B,D}) where {N,M,B,D SFBC = (B == BC_SF_AFWB) || (B == BC_SF_ORBI) - if (it == 0) && SFBC + if (it == 1) && SFBC for id in 1:N-1 U[b,id,r] = bndfield(phi1,phi2,lp.iL[id]) end