| |
| | [No title] (Site not responding. Last check: 2007-10-21) |
 | | Princeton shall not be liable * C for any direct, consequential, or other damages suffered by any user * C or any others resulting from the use of this research software. |
 | | C SLMIN=0.2 C DO 100 LOOP=1,10 C sweep right DO 3 J=2,JM-1 DO 1 I=2,IM-1 IF(FSM(I,J).EQ.0..OR.FSM(I+1,J).EQ.0.) GOTO 1 SL(I,J)=ABS(H(I+1,J)-H(I,J))/(H(I,J)+H(I+1,J)) IF(SL(I,J).LT.SLMIN) GOTO 1 DH=0.5*(SL(I,J)-SLMIN)*(H(I,J)+H(I+1,J)) SN=-1. |
 | | H(I+1,J)=H(I+1,J)-SN*DH H(I,J)=H(I,J)+SN*DH 1 CONTINUE C sweep left DO 2 I=IM-1,2,-1 IF(FSM(I,J).EQ.0..OR.FSM(I+1,J).EQ.0.) GOTO 2 SL(I,J)=ABS(H(I+1,J)-H(I,J))/(H(I,J)+H(I+1,J)) IF(SL(I,J).LT.SLMIN) GOTO 2 DH=0.5*(SL(I,J)-SLMIN)*(H(I,J)+H(I+1,J)) SN=-1. |
| www.phys.ocean.dal.ca /~kelley/linux/benchmark/pom97.f (355 words) |
|