Как сгруппировать массив по дате для цикла?

Мне нужно найти подходящий период в каждом году, вычислить его продолжительность и сумму температур. У меня есть срез массива за несколько лет, он выглядит так:

Табл.1

Я разбиваю этот массив на года. Каждый год прогоняю через код.

За дату устойчивого перехода весной принимается первый день периода, сумма положительных отклонений которого превышает сумму отрицательных отклонений любого из последующих периодов с отрицательными отклонениями. Осенью наоборот - концом периода будет считаться первый день того периода, сумма отрицательных отклонений которого превышает сумму положительных отклонений любого из последующих периодов с такими отклонениями. Чтобы найти самый большой период с температурой выше 0 я написал такие условия:

# переменные для условий
a_max = 1
a = 0 # сюда записывается искомая сумма из начального периода
b = 0
c = 0
d = 0 # сюда записывается искомая сумма из конечного периода
d_max = 0

# условие для нахождения начала периода
for i in sort_df.loc[sort_df['Полугодие'] == 1]['Сумма температур']:
    if i >= 0:
        if i > a > abs(b):
            a = a
        if i > a and abs(b) > a:
            a = i
        if a > i > abs(b):
            a = a
        if i > a_max:
            a_max = i
        if a < abs(b) < a_max:
            a = a_max
    if i < 0:
        if abs(i) > a:
            b = i
            a = 0
        if abs(i) < a:
            b = i

# условие для нахождения конца периода
for i in sort_df.loc[sort_df['Полугодие'] == 2]['Сумма температур']:
    if i >= 0:
        if i >= abs(d):
            c = i
            d = d_max
        if i < abs(d):
            d = d
            c = i
    if i < 0:
        if i < d_max:
            d_max = i
        if c == 0:
            d = i
        if abs(i) > abs(d) and c > abs(d):
            d = i
        if abs(i) > abs(d) and c < abs(d):
            d = d

Возможно они написаны некрасиво и громоздко, но ничего лучше я пока не придумал и они работают для всех возможных случаев. Далее находим индексы значений в a и d:

idx_a = sort_df[sort_df['Сумма температур'] == a].index[0]  # индекс начальной даты
idx_d = sort_df[sort_df['Сумма температур'] == d].index[0]  # индекс конечной даты

Потом находим начальную и конечную дату, продолжительность и сумму периода для положительных значений:

start = sort_df.at[idx_a, 'Начало']  # начальная дата
end = sort_df.at[idx_d, 'Начало']  # конечная дата
delta = end - start  # продолжительность периода
summa = sort_df[idx_a:(idx_d+1)].loc[sort_df['Сумма температур'] >= 0, 'Сумма температур'].sum()  # сумма температур

После этого записываем эти переменные в уже созданный фрейм с заголовками:

warm_tot.loc[len(warm_tot)] = [start, end, delta, summa]  # записываем строку во врейм

Результатом должно быть:

Табл. 2

Всю эту процедуру неудобно делать отдельно для каждого года. Я пытался проверить код на массиве, содержащий данные за 50 лет, в итоге он записал мне одну строку с периодом с 1971 по 2006 год (хотя что-то подобное стоило ожидать). Как можно вписать этот код в один цикл так, чтобы он отдельно проходил по каждому году массива как показан в Табл. 1, а итог записывал в новый массив, где будет один период для каждого года?