Flexible Integro-Differential Equations For Bayesian Modeling Of Spatio-Temporal Data