this is what I am doing:
delta_list = [[ abs(daily_2016[y,1] - avg_max_array[y,1])] for y in day_of_years ] delta = np.array(delta_list) delta.shape # (366,)