Multilag multivariate autoregressive model for the generation of operational hydrology