Fortran 子程序求助！

zhxbkofkyo 发表于 2022-10-21 13:02
大概明白了，就是当子流类型不止MIXED一种时，例如包括MIXED和CISOLID是，RATES数组长度就会根据子流数量 ...

非常感谢您！解决了困扰我一个月的问题！太感谢了！
我在Fortran 子程序上推进的很困难，还请老师指导一二。如何才能像您一样，理解的这么透彻，如何才能知晓Aspen 帮助文档中没提及的这些规定。（例如刚才的数组长度）

zhxbkofkyo 发表于 2022-10-21 05:02
大概明白了，就是当子流类型不止MIXED一种时，例如包括MIXED和CISOLID是，RATES数组长度就会根据子流数量 ...

大佬 您好。又遇到新问题了，按照自带的煤气化动力学子程序进行模拟，一直报错RCSTR 质量不守恒。

  **  ERROR
      BLOCK B7GISF-1 IS NOT IN MASS BALANCE:
      MASS INLET FLOW = 0.63392644E+00, MASS OUTLET FLOW = 0.72997716E+00
      RELATIVE DIFFERENCE = 0.15151713E+00
      CHECK USER KINETICS OR STOICHIOMETRY/MW'S

在物流结果中发现碳C反应前后没有消耗，但是其他物质CO、H2、CO2、CH4反应前后有变化。检查后发现C的反应速率正常。 很奇怪 不知道怎么办

这个看不出来，你看看你的转化反应有没有问题呗

zhxbkofkyo 发表于 2022-11-24 13:52
这个看不出来，你看看你的转化反应有没有问题呗

反应完全按照  自带的‘固定床煤气化’ 中的动力学子程序 编写。一直报错就很离谱

zhxbkofkyo 发表于 2022-11-24 13:52
这个看不出来，你看看你的转化反应有没有问题呗

RATES(NCOMP_NCC*2-1)       ！我怀疑是这个语句 与 固相C的反应速率 不对应

反应前后就是C没变化 很奇怪

RATES(NCOMP_NCC*2-1)和RATES(NCOMP_NCC*2)，分别是C和S吧，先看看C和S是不是定义成CISOLID类型了，打印下NCOMP_NCC，这个数值是否等于常规+CISOLID的数量和？尝试下RATES(NCOMP_NCC*2)改成RATES(NCOMP_NCC)会是什么结果？多试试呗，也没有好的办法

zhxbkofkyo 发表于 2022-11-28 09:00
RATES(NCOMP_NCC*2-1)和RATES(NCOMP_NCC*2)，分别是C和S吧，先看看C和S是不是定义成CISOLID类型了，打印下N ...

好的  感谢 大佬