IF(HIN.LT.0.0) THEN !为非饱和状态
IF(INDEX.EQ.0) THEN !排水
DO I=1,N
IF(B(I,1).GT.BC) THEN
IA(I,1)=-1
ENDIF
IF(B(I,M).GT.BC) THEN
IA(I,M)=-1
ENDIF
ENDDO
DO J=2,M-1
DO I=1,N
IF(IA(I,J-1).EQ.-1.AND.B(I,J).GT.BC) THEN
IA(I,J)=-1
KI=I
10 KI=KI-1
IF(KI.GE.1) THEN
IF(B(KI,J).GT.BC) THEN
IA(KI,J)=-1
GOTO 10
ENDIF
ENDIF
KI=I
欢迎光临 Fortran Coder (http://bbs.fcode.cn/) | Powered by Discuz! X3.2 |