On fractal interpolation functions in hydrological modelling