This allows us to model the Prebotzinger complex and run simulations based on different parameters