accrued_dose = 0
    total_dose = []
    for incr_dose in tid_arrays:
        accumulated_dose = incr_dose + accrued_dose
        total_dose.append(accumulated_dose)
        accrued_dose = accumulated_dose[-1]
    dose_array = np.hstack(total_dose)