The Oregonator system, consisting in a set of three Ordinary Differential Equations (ODEs), is widely used in the literature to model the oscillatory Belousov-Zhabotinsky (BZ) chemical reaction. Recently, D’ambrosio et al. used adapted numerical methods to follow the apriori known qualitative behaviour of the solution and improve the quantitative matching with experiments. Budroni et al. employed a modified version of the Oregonator to model a network of diffusively coupled inorganic oscillators, confined in micro-compartments by means of a flow-focus microfluidic technique. This class of networks is effective for predicting and understanding the global dynamics of those systems where the diffusion of activatory or inhibitory signals...