三维SIMPLE程序应用中出现的问题

开心果

咱网站上贴的三维simple程序中solve程序块儿中为什么要加
IF(DENOM.EQ.0.) THEN
WRITE(*,998) DENOM,PT(K-1),AKM(I,J,K),AP(I,J,K),NF,I,J,K
      WRITE(2,998) DENOM,PT(K-1),AKM(I,J,K),AP(I,J,K),NF,I,J,K
998FORMAT(1X,';DENOM=';,F15.5,4X,'T(K-1)=';,
     & F15.5,4X,';AKM(I,J,K)=';,F15.5,4X,';AP(I,J,K)=';,F15.5,4X,';NF=';,I3,
& ';I=';,I3,2X,';J=';,I3,2X,';K=';,I3)
DENOM=0的情况不可避免嘛？

apfel

我也在想这个问题，你知道了吗？

apfel

我刚刚调通了一个简单的三维算例，这一部分的影响不大。所以，先暂时不考虑这个问题了。