SoC_m = abs((SoC1 - SoCk1));

    if h0 == h1
        we = we0;
        tau_e = (tau_e1 - tau_e0)*(SoCk1 - SoC_vec(idx0))/(SoC_vec(idx1)-SoC_vec(idx0)) + tau_e0;
        h = h0;

            else
                if SoC_m < 0.05;
                    mod = 1;
                elseif SoC_m > 0.05;
                    mod = 0;
                end

                if mod == 0;
                    we = we0;
                    tau_e = tau_e0;
                    h = h0;
                elseif mod == 1;
                    we = we1;
                    tau_e = tau_e1;
                    h = h1;
                end
    end