da = da.assign_coords(year_month=da.time.dt.strftime("%Y-%m")) result = da.groupby("year_month") - da.groupby("year_month").mean("time")
da = da.assign_coords(year_month=da.time.dt.strftime("%Y-%m")) result = da.groupby("year_month") - da.groupby("year_month").mean("time")