This strategy is sufficiently general to allow for the use of a large variety of objective functions, global optimization techniques, and surrogate modeling strategies. In this particular study we focus on only single combination: a weighted least-squares objective function, the differential evolution global optimization algorithm, and Gaussian Process (GP) surrogate modeling.