fix(btcs): alpha intercell on explicit constant boundaries

See MR @ https://git.gfz-potsdam.de/naaice/tug/-/merge_requests/25

[skip ci]
This commit is contained in:
Max Lübke 2024-03-11 15:47:29 +01:00 committed by Tim Werner
parent 8233d868b2
commit 5c6ac062cc

View File

@ -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