Hybrid gaussian process approach to robust economic model predictive control with an application to solar thermal systems