Gaussian Processes in modern Fortran