seasonality function

Calculate seasonality of discharge