Modeling soil moisture

[Amitava-Soil]

In code: instead of 'S' , I believe it should be S_obs
#Implement markov chain
for t in range(1,df.shape[0]):
alpha[t]=alpha_fn(DOY[t],c,phi)
S_pred[t]=S_min + alpha[t]*(S_obs[t-1]-S_min)+P[t]
S_pred[t]=np.min([S_pred[t],S_max])

I am also not getting the figures
I appreciate any help.

Thanks

Amitava Chatterjee