Can PyGFunction be used to model horizontal arrays? Any recommendations on where to start or what to adapt would be much appreciated!