Iterative MAP estimation with increasing data points

I answer my question.
Writing the model statement basic_model = pm.Model() inside the for loop made it work.

for est_end in range(est_start, len(df_est[‘T_mean’]):
basic_model = pm.Model()
with basic_model: