Gaussian Processes routines for astronomical time series