Error pops up when initializin HeterogeneousBorefield. Solved by changing line: `initial_ΔT::Vector{S} = zeros(S, Nb)` into: `initial_ΔT::Vector{S} = zeros(Nb)` in HeterogeneousBorefield.jl