BEGIN TEST = 0
-
DO 1 I=2, NTOT1
-
DO 1 J=2, NTOT1
-
TEMP=0.25 * (PHIOLD(I, J + 1) + PHIOLD (I, J-1) + PHIOLD (I+1,J) +
-
PHIOLD (I-1,J))
-
TEST = AMAX1(TEST,ABS(TEMP-PHIOLD(I,J)))
-
1 PHINEW (I,J) = TEMP
-
DO 2 I=2, NTOT1
-
DO 2 J=2, NTOT1
-
2 PHIOLD (I,J) = PHINEW (I,J)
-
IF (TEST>CONVG) GO TO BEGIN
-
ELSE STOP
|