N AL N
no                          Ill"4).. . ..  ..

PROGRAM RANWLK

74/175  OPT=0

FIN 5.0+508

PAGE      3

SI b
1 I-
1.17
1  1
121
121
122
1234
S24
12
1 4'3
121
1 4 71
1  Pe
12?
1  3
11 4
126
I .
5 8
136
140
140
161
1458
1469
147
10 1
112
123
I 54
155
187
lb 9
160
161
1623
1 65
186
170
171

9      ~         ..,. *~ r4uw'*-~'~   ~
9

100 CONTINUE
00 110 I = INK
00 105 J=INY
HOLOt(I J)  =  I[EAf)IIJ1
1 0   CONTINUE
110  CONTINUE
C***PLACE RESTART PARTICLIS***
C
IF(IHEDIN.EO.0) GO TO 111
nO 112 NP =  ,NTOTAL
CALL FIND( X(NP) Y4NP) ,e)XPARDYPAR,NXPARNYPARNSMS)
NOUT(NSMS)=N0UT(NS,4S)+1
112   CONTINUE
111   CONTINUE
C     INITIAL IZE COUNTER FR OUTPUTTING HEAOS LVERY "OUTFRO" TIMESIEPS
OUTKNT = I
C
C  ****  START  TIME  STEP  LOOP  .**.
C
NT=D
120 CONTINUE
TIME = TIME * TSTEP
NT=NT*1
IF(MTYPE .EQ. 3) GO TO 140
C
C       SOLVE FOR HEADS IN UPPER LAYER
C
ERACOO = .
L = 1
CALL IAOI (L,'TITIM4,NUMIT*EEMAX
C
IF(ERRCO0 .LF. 0) GOTD 130
CALL ERPROC
IF(FRRCO)  .GE.  100)  GOTO  900
130   CONTINUE
140   CONTINUE
IFtMTYPL .E0. 2) COTO 150
C
C       SOLVE FOR LOWER LAYER HEADS
C
ERRCO0 = 0
L = ?
CALL  IAOl (LiNTiTIME  1NUMIToEEMAX)
IF(ERRCOD ILE. 0) GOTO 150
CALL ERPROC
IF(ERRCOO .GE. 130) GOTO 900
150   CONTINUE
C
C       COMPUTE INTER-LAYER LEAKAGE
C
IF(MTYPE *E. 1) GOTO 160
ERRCOO = 0
CALL LEAK
IF(ERRCOU  .LE.  0)  GOTO  160
CALL FRPROC
IFtERACOD *GE. 100) GOTO 900
160   CONTINUE

0

0

40
ag

is