Это должно сработать, но я вполне уверен, что есть векторизованный способ сделать то же самое.
df=pd.DataFrame([1,1,2,2,3,3,4,4,1,2],columns=['phaseM'])
df['phaseMday']=np.linspace(1,len(df),len(df))
for i in range(1,len(df1)):
if df['phaseM'].iloc[i]<df['phaseM'].iloc[i-1]:
df['phaseMday'].iloc[i]=1
else:
df['phaseMday'].iloc[i]=df1['phaseMday'].iloc[i-1]+1