From 5c6ac062cc1b9a3916a8f1c5f69850c96842a08d Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Max=20L=C3=BCbke?= Date: Mon, 11 Mar 2024 15:47:29 +0100 Subject: [PATCH] fix(btcs): alpha intercell on explicit constant boundaries See MR @ https://git.gfz-potsdam.de/naaice/tug/-/merge_requests/25 [skip ci] --- julia/TUG/src/Core/BTCS.jl | 5 +---- 1 file changed, 1 insertion(+), 4 deletions(-) diff --git a/julia/TUG/src/Core/BTCS.jl b/julia/TUG/src/Core/BTCS.jl index dd9b8af..60320ef 100644 --- a/julia/TUG/src/Core/BTCS.jl +++ b/julia/TUG/src/Core/BTCS.jl @@ -68,12 +68,9 @@ function calcExplicitConcentrationsBoundaryConstant( sy::T, )::T where {T} alpha_center_neighbor = calcAlphaIntercell(alpha_center, alpha_neighbor) - alpha_center_center = - alpha_center == alpha_neighbor ? alpha_center_neighbor : - calcAlphaIntercell(alpha_center, alpha_center) return sy * alpha_center_neighbor * conc_center + - (1 - sy * (alpha_center_center + alpha_center)) * conc_center + + (1 - sy * (alpha_center_neighbor + alpha_center)) * conc_center + sy * alpha_center * conc_bc end