X = pd.DataFrame(index=LLN_weather.index) X["day"] = LLN_weather.index.day_of_year X["hour"] = LLN_weather.index.hour def ConsoHot(df): i = df.day h = df.hour if (i < 165 or i>257): return VHotLab_annee[h] else: return VHotLab_vac[h] X["ConsoHot"] = X.apply(ConsoHot,axis=1)

X["total"] = X.ConsoHot +Vbureau+Vlabo X.total.plot()

Tint= 21 cp = 1000 rho = 1.2 Tru = X.total * cp * rho * (Tint - LLN_weather.tsa) Tru.plot()

LLN_weather.tsa.plot()