Python GEKKO中二维Intermediate的正确使用方法求助

• 之前用GEKKO构建含大型反应网络的平推流反应器模型时，卡在了中间计算环节。
• 尝试用m.Intermediate处理中间变量，也参考过二维Intermediate数组的实现方法，但写出来的代码不仅冗长，还无法正常运行——核心问题是最终方程里的数据数组长度没法随时间步动态调整。
• 知道过度复杂地用m.Intermediate可能不是它的正确用法，但因为反应网络规模太大，实在找不到其他可行的计算方式，附上代码求帮助或优化提示。

# 请粘贴你的GEKKO建模代码到此处

内容的提问来源于stack exchange，提问作者MoDi