Как рассчитать суточную эвапотранспирацию по уравнению Харгривза-Самани и с помощью Python?

---

У меня есть десятилетние данные о погоде, включая максимальную температуру (Tmax), минимальную температуру (Tmin), количество осадков и солнечную радиацию (Ra) за каждый день.

Сначала я хотел бы рассчитать эвапотранспирацию (ETo) для каждого дня, используя следующее уравнение:

ETo=0.0023*(((Tmax+Tmin)/2)+17.8)*sqrt(Tmax-Tmin)*Ra

Затем производится расчет среднемесячных и годовых значений всех параметров (Tmax,Tmin, осадков, Ra и ETo) и вывод их на печать в формате Excel.

Я написал несколько частей. не могли бы вы помочь мне с его завершением? Я думаю, что ему может понадобиться круг "для".

Что я уже пробовал:

import numpy as np
import pandas as pd
import math as mh
# load the weather data file
data_file = pd.read_excel(r'weather data.xlsx', sheet_name='city_1')
# defining time
year = data_file['Year']
month = data_file['month']
day = data_file['day']
# defining weather parameters
Tmax = data_file.loc[:,'Tmax']
Tmin = data_file.loc[:,'Tmin']
Rainfall = data_file.loc[:,'Rainfall']
Ra = data_file.loc[:,'Ra']
# adjusting time to start at zero
year = year-year[0]
month=month-month[0]
day=day-day[0]
#calculation process for estimation of evapotranspiration
ET0=(0.0023*(((Tmax+Tmin)/2)+17.8)*(mh.sqrt(Tmax-Tmin))*Ra

---

year = year-year[0]
month=month-month[0]
day=day-day[0]