有没有冻土水热耦合的大神来指导下

我在用偏微分方程做冻土的相变过程中水热耦合，卡在这里很久了，我检查了报错的方程好像也没啥毛病，但就是一直这样报错。想问下有没有大神遇到过类似定情况能来指导下？程序在零度以下可以运算，但是只要在0度以上就会这样。按道理0度以上应该更简单呀，咋回事呢。