-------
e
"
o
10
15
*
1
1
1
1
1
9
'
»
»
'
2
S
E
9
o
1"
IP
n
p
t
i
i
i
i
2
?
2
t
2
>
2
E
8
1C
10
10
3
a
2
2
2
2
2
2
2
2
2
F
9
11
Q
O
9
3
1
2
-
2
2
2
2
2
?
*f 5 S 6
9 8 r
3 it) 9 7
1? 3 7 E
9 3 7 P
o 8 7 S
9875
3 « 1 1
7 p t t
» c 1 I
112?
2272
277'
2 ' 2 3
? ? J 1)
? ' 7 0
2 >'>
2 2 ? 2
2 7 ? 2
7222
e e e
<; s s
c c c
*»*!>
<; s s
c S 5
r * r
s 5 5
1 < 1
11
i i *i
111
Til
0 ? 1
T n i
211
221
722
2 ? »
5 4 tt t*
t» T -» 2
li T ** 7
** * 7
» 7 » 2
2111
- 1 1 1
S
a
2
3
I
"
2
2
7
3
?
2
2
2
2
2
1
1
4
2
3
5
3
2
3
v
5
3
3
3
2
2
2
2
1
1
"-
*
2
3
T
3
2
5
7
3
3
t
2
2
2
2
2
,
2
1
2
2
3
3
2
3
3
3
3
2
2
2
2
2
2
2^
7 i
2*
1 "
lit
1 1 *
1* «
i 1
211
999
7 £. C.
9 9 ->
3 ? 7
7 T **
Z9 *>
/ '
29 *^
£
2*9 «
-r
29 T
' ^
3 » '
3-» V
J *
3 ' '
3 ? »
3 * '
7 ? 7
2 . ,
22?
2 " "*
2 2 ?
27'
f « n
^ 1
-in
c i n
i i '
SIC
nr
n"
Tl"
« n
V 1 '
VI?
i/ « ^
V I '
I/ 1 O
V1L1
V 1 ^
V10
VI"
V10
VI-
VI 0
V!"
vir
VI"
VIC
VI"
V10
CA)
X)
I
oo
-------
» ?
7 ?
7 ?
.
T T
«
en >"
tsj
3 3
9 1
» 1
1 2
» ?
2 2
2
2
T
3
»
3
8
1
2
2
2
2 ?
2 2
J 2
4 4
12 12
9 3
'-»2'2'r??222?''
??22?'''2222233-'-'
?1333T''33333'?3T
97t>HI:C666?r;5"|it«|iI
B7R->2TI4t5ti4(|i«4«i)
8<»ll!lllll'?2''>J»
poiiiiiiii792-»-»-r
22'7111'2T333333
I33'll??T33T33-rT
?-'3?l??''3?3^3J'T3
33333'''333333''-'
<-l'
'll
ril
* 11
Vll
Vll
tfl t
vt
Vll
VT
Vll
00
OO
-------
en
? 2
- ?
2 2
> 7
2 7
- 7
321"".
' 1
? 2
> 7
r c
v b
E 7
? 1C
17 13
f 11
13 13
1' I'
1" I"
11 11
11 10
f 5
3 3
8 8
3 3
2
2
2
2
2
2
1
2
S
7
B
* r*
« u
11
11
11
12
11
10
IS
e
a
8
8
"
2
'
7
'
2
1
2
1
S
fl
* n
1"
11
11
U
1"
11
IP
-
3
8
9
3
"
2
2
7
"
2
1
2
T
,
7
n
1 U
11
11
13
11
1"
o
1
8
3
f
3
2
'
2
?
2
7
1
5
B
e
4 pi
1U
11
13
1?
1"
3
P
3
3
1
e
4
2 3
7 *
2 2
7 r
? >
" ?
c c
7 *
1C 3
o jp
12 3
11 P
IP 3
8 7
7 S
6 ?
5 «=
8 8
3 3
1 7
1 1
3 «
T 3
7 2
' ?
2 2
7 2
j j
S 0
5 "5
7 C
3 e
7 f,
7 F
7 '
r c
S 5
«; «;
r ^
1 1
1 1
" 1
1 1
3
2
2
2
2
«.
3
1
1
5
r
7
5
1
6
E
S
5
5
7
7
2
1
T 7 ^
2 7 "
** » ^>
777
" " "
22"
13'
U 1 ft
II1;
677
'88
778
<; s <;
(I Cl c
c 5 S
? 5 5
S e 5
r 5 c
S 5 r
222
'77
727
1 1 7
7
f
2
7
2
2
T
4
c
T
3
7
c;
G
5
!
"
5
"
2
2
2
2
2
7
2
2
2
7
T
u
T
7
s
5
5
«
5
c
S
7
2
7
2
2
2
2
2
2
2
3
4
1
S
e
K
5
5
S
5
S
c
5
2
7
2
2
2
;
2
7
t
2
3
4
4
S
5
S
"
5
S
5
5
5
?
2
7
2
2
2
2
2
2
2
3
3
4
4
5
5
5
5
5
5
5
5
5
2
2
2
2
2 » "
7 2 r
, , ,
- 7 ?
?»».
77"
3 ' »
J -r t
3 » '
i, T T
* *i J
(4 T- T
413
544
555
5 c r
555
"j * c
5*5
5 e e
72?
2 " "
7 7 "
2 » "
Vll
Vll
Vll
VI!
Vll
Vll
'1-
ri'
SI"
r 1 **
v 1 i
< 12
S« <9
1 '
«!"
^1?
?1?
r!2
SI'
ri2
'f
VI"
VI"
VI"
VI"
CO
»
i
00
45-
-------
tn
1
2
-
?
1
-
2
-
"
'
3V
r
5
7
7
8
Q
11
I'-
ll
11
11
11
1 8
2 1
' 2
2 2
' 2
2 2
2 2
2 2
? 2
2 2
? 2
(P ft
^ *»
5 1
7 S
e s
* 8
9 9
n 10
1' 12
11 11
11 11
11 11
11 11
11 11
1
9
2
t
2
2
'
2
2
2
1
5
f
0
8
1
13
11
11
11
11
11
2 1
2 2
? 2
? 2
- 2
2 7
2 2
2 ?
' ?
2 2
1 1
1 1
; s
S 6
c E
9 10
13 13
11 11
11 11
11 13
11 11
i' in
3
»
2
2
?
2
-
2
'
1
1
E
9
ft
10
13
11
11
12
10
.,
7*»
o
» -
2 3
3 3
T T
2 2
- ,
2 2
- ?
.
4 1
5 5
r ir
S S
! 1-
11 12
12 11
I- in
12 11
11 11
1" 9
q *
14
1
1 7
2 2
3 3
3 3
3 3
3 3
' ?
2 2
2 2
2 2
5 5
5 5
R £
7 8
1C 9
in s
« e
9 8
9 8
in q
9 A
9 7
9 *
* 1
2 ?
> T
» T
< 3
-, 2
2 2
f t
3 3
5 6
e 7
7 8
1" 10
11 1?
11 12
9 11
1 1"
f 9
9 8
7 7
7 7
' Z 2 7
^9^7
t £ £.
2222
'-22
'222
« 2 2 2
- 3 3 3
2333
? 3 T T
'333
3 * * 3
Be c f
Do
C C C C
*> *> *- fa
5779
7 a >
9 *C 10 12
11 12 13 11
1' 13 1? 11
12 11 IP "3
11 10 9 8
3999
» 9 8 7
3 S ' 7
7 7 P e
7 R « s
~
'
2
'
2
2
3
3
3
3
3
9
11
I'
1"
8
7
7
e
E
S
E
S
2.
2
2
2
2
2
3
3
3
3
3
8
10
10
8
E
6
e
S
e
6
6
E
7 ' '
"9 f f
i. .
2 ' -
? ? ?
2 ' '
222
2 ' '
3^7
3 ' '
T T T
33"
3 ' T
6f ^
b I
7^ ^
*
S S S
» i a
"88
77'
S R r
S P «
E e E
S^e
SEE
R C R
E f e
q c *
114^
V I
V12
vi-
vi-
V1 '
VI-
VI?
VI'
VI?
vi-
r i 7
3 1^
ri3
si*
?13
S13
Sl»
S13
ri'
SI*
<13
;i3
El»
«1*
00
pa
OO
-fi.
-p..
-------
C/l
CO
14 14
14 14
4 4
P 8
1 8
6 8
" 7
9 2
* 2
? 2
' 7
? 2
9 '
? 2
' 2
7??0'7.
4 4
* 4
s E
7 7
8 P
14
7
8
3
1
2
2
2
2
2
2
2
2
2
5
E
7
14 1?
13 1'
T B
B 8
o, <»
1 1
2 2
- "
2 9
7 2
2 3
2 »
2 2
2 9
2 7
7 '
4 2
4 7
<; 4
10 9 a
8 R
8 t
119
9 2 2
222
333
733
3 7 '
T 3 4
5 5 *
2 ? 2
7 ^ >
222
27?
723
9 ? o
274
5 T 7
1 1
1 2
7 ?
7 ?
T 7
? 7
3 3
4 4
" 4
2 9
? 2
2 2
2 2
4 4
4 5
S c
R P
7
1
2
2
2
2
7
3
4
4
7
2
7
2
4
4
S
S
8
7
2
7
7
1
1
9
3
a
4
-
2
-
7
4
4
5
7
P
E F 6
"i 1 9
t
7 ** ?
** "5 9
* £
192
172
122
222
"72
373
T * 3
477
922
2 ' 2
99j
2 T 3
4 S 5
5 5 K
B E 7
7 3 9
9 1C 11
e
r.
2
2
2
2
2
2
5
3
2
2
T
3
T
6
7
9
11
S
R
2
2
2
?
7
2
2
2
7
2
3
3
3
3
E
6
8
S
11
S 6
E S
21
t,
2 2
9 9
* £.
7 2
9 2
7 2
2 2
2 2
2 2
7 2
9 2
2 2
2 3
3 3
3 3
' 3
s e
7 7
« 8
8
5
F
'
2
2
?
3
2
2
2
2
t
T
7
3
B
S
e
7
-
r c,
E E
2 7
2 9
9 9
? 9
3 ?
» »
? 7
- -
3 3
T ^
T T
T T
4 4
C E
c ^
? <;
7 7
7 7
V1T
VI'
V1T
VI'
VI'
VI'
VI?
V1T
S14
il° OJ
' w
f
o
-------
f IT 14 14
I' 14 11 11
13 11 11 14
1* 13 13 13
1"* 13 13 13
1T 1* 13 13
I* 1? 13 1'
J* I' 12 1"
4 4 4 P
tj C} n t^
n c s 7
1 9 8 *
111"
19 7 2 2
7 3 3 T
'332
7722
A < w n
* £ £ £
? 2 2
f > 2 ?
332?
» T 2 2
''22
14
14
13
1'
1'
12
17
1?
e
1
1
2
3
2
>
2
?
2
1
1«
14
13
13
12
12
17
1?
S
R
1
T
7
7
?
?
7
?
7
1
11
17
12
1"
12
12
12
17
ft
*1
B
8
1
7
T
T
T
2
7
2
7
1
12
10
12
12
1?
17
1?
17
1
1
7
2
T
7
T
3
2
7
?
^
tt
in
12
17
17
1?
12
12
1
1
i
^
T
5
,
3
3
3
2
2
2
1
11
n
12
1'
1'
1'
11
11
1
1
T
2
^
7
3
«
4
T
2
2
2
1
1C 1*^
1C 11
12 11
11 1*
12 1T
11 11
11 11
11 11
1 2
3 2
7 i
2 1
2 2
4 ?
»
3 '
3 '
2 2
2 '
2 7
2 2
12
11
1C
ir
11
ID
10
10
2
2
'
1
1
9
f
T
^
T
-
y
2
1'
ID
8
f
Q
ir
ir
1C
2
7
1
7
T
T
»
T
3
T
T
7
10
P
B
E
7
8
9
1C
2
1
2
2
3
3
3
3
T
3
T
J
8
r
5
F
*
7
R
o
>
2
2
3
2
T
3
^
3
3
3
T
4
9
E
E
?
B
S
7
8
2
2
2
2
3
3
3
3
3
3
3
3
,.
5
5
5
g
5
«
E
7
2
2
?
3
*
1
3
T
3
3
3
3
3
3
5
6
F
e
E
E
B
E
7
2
2
3
3
3
3
3
3
3
3
3
3
3
5
5
c
S
B
E
f
5
E
2
2
3
3
3
3
3
3
3
3
T
3
T
c
c
F
S
c
c
r
c
r
2
3
V
3
f
j
'
3
3
'
3
E
*
S
«;
B
-
F
c
S
2
3
»
T
T
3
r
T
3
*
3
E
S14
S14
14
-14
S14
114
?!«
W 1 b
v m
V14
V14
VI"
V14
V14
V14
V 14
V14
V14
V14
V14
V14
V14
V14
VI"
«!"
en
X3
i
OO
SI"
-------
1
T
i
T
S
P
11
1'
14
1 4
14
14
^ T
C
4
a
p
1
2
»
2
»
2
,
1
1
0
T
G
R
11
14
14
1 1
14
4
4
4
0
1
2
?
2
2
2
7
2
1
1
0
4
G
9
12
14
14
13
« -T
1 J
13
4
4
4
Q
1
2
2
2
2
2
2
1
C
3
5
7
1"
13
14
4 (l
1 4
1 *
13
13
If
c
1'
4
4
-
T
1
t
2
j
2
>
2
1
0
2
3
S
7
ir
it
14
1 (l
14
1^
*
i*
ir
* ^
1 £
4
4
P
1
1
1
2
2
»
1
1
1
1
2
4
S
c
10
14
14
1"*
1^
12
12
4 «
1 1
11
4
"
0
1
1
1
7
>
2
'
7
0
?
T
T
7
«
1C
14
1?
f A
1Z
1C
1C
11
4
"
8
8
1
'
2
?
2
7
2
0
T
4
7
4 n
11
14
14
11
11
9
t
o
1C
'3
0
S
o
2
'
2
»
3
?
3
,.
4
C
1
'0
12
M
12
1C
c
30
1
1
1
1
2
2
3
T
3
T
3
,.
q
7
"
11
1'
1'
11
"
T
7
0
Q
1
1
1
^
^
2
*»
2
*
T
3
4
e c
8 »
in i"
1C 11
11 I"*
1(3 11
1C 11
ir 11
ID 1**
p 1 b
1C 11
15 1?
9 9
9 9
1 -
1 ?
2 **
2 2
*? 1
2 1
2 1
2 7
3 "
3 3
4 4
4 4
4 T
3
P
11
17
12
12
11
10
11
«
1 1
1?
1
9
-
2
1
-
2
1
2
7
4
4
s S
? 0
11 11
1- 12
12 1'
1C 9
IP
o 8
O C
a o
O T
/
3 3
9 9
n o
* 7
' 2
*9 "9
w t
' 2
2 2
2 2
' 2
» 7
- 2
4 3
4 4
4 4
4 4
T
O
11
17
17
S
7
6
8
p
9
2
2
2
2
2
9
7
3
s
4
4
4
7
9
11
12
11
8
«
6
9
8
Q
2
2
2
2
2
2
3
3
3
3
4
3
7
9
1"
11
*
8
G
6
P
9
7
2
2
2
2
2
2
3
3
3
3
3
T
7
9
1 O
9
3
E
O
G
S
8
e
2
2
2
2
2
2
2
3
3
3
3
3
3
7
«
9
9
7
R
5
E
7
8
n
2
2
2
2
7
7
T
3
T
3
T
3
' *
e o
3 9
8 B
T -T
F 6
* r
G G
7 7
P «
Q O
8 3
7 »
2 ?
7 -
, ~
3 3
^ T
» 3
t T
3 3
t T
7 3
3 3
'I'
I'
SI'
?1*
*1«
El-
r:t-
?i-
*" 1 F
_ 1
Slr
-le
SI-
VI*
vir
Vlr
V15
Vlr
Vlr
vn
n-
VI"!
Vl«
VI-
.
CO
-------
'22
* » 2
o t t
J t c
4CT?.
«. 9 s
* B 8
4 c, 5
I> tt c
** M D
T ^ 7
J c
» T 4
455
c » 8
On
e
5
5
E
9
14
13
12
1?
12
1?
11
4
1
*
2
2
2
T
I.
f
5
c.
9
n
11
13
1'
1?
1'
11
10
0
1
T
2
2
C
2
2
S
T
T
11
13
1?
12
12
11
11
in
f
a
i
>
2
0
T
t
7
0
11
11
17
I?
17
12
12
11
in
»
1
2
?
3
o a
7 ?
2<9
f
- 2
3 3
4c
b
T Q
10 in
12 11
11 1">
If
in «
12 «
13 in
17 10
11 in
in in
* «
L 1
1 '
1 2
? '
- -
^ T
J T
2-»
2 "
2 ?
3 4
5U
K
77
/
0 0
10 1?
10 »
9 9
8 8
7
7 E
7 F
9 7
*
q 7
a *r
?*
l
2 -
7 »
7 1
2 1
7 '.
r
3
ia
s
9
A
7
E
e
E
s
e
i
i
i
i
3 7
7 2
T T
5 5
9 1C
1 9
9 8
3 '
7 e
c. c
E E
E e
e E
c. «
E e
** 9
_ t
?ry
£.
> ?
? 7
2 2
7 2
? 2
T
T
C
1C?
a
e
E
s
F
F
e
B
E
2
?
2
2
2
4
3
4
5
1C
9
7
S
R
=-
e
o
5
r.
F
>
2
2
'
2
4
3
4
5
n
7
C
6
6
E
E
S
E
S
R
2
2
2
2
3
4
4
4
S
7
7
E
E
E
E
E
S
E
e
^
2
2
3
3
3
4
4
4
E
E
E
E
f
S
E
E
E
K
E
'
3
3
3
4
4
4
e
f>
6
E
E
E
c.
C
«!
E
2
T
3
T
9
'
C
e
s
5
c
E
e
e
E
E
c
r
3
T
3
T
T
Vlr
(14 p
V 1
si*-
c i r
i A*
Sir
Tic
Sir
Sic
Sic
Sir
Sic
-ir
Sic
\M C
V 1 *t
I/ t C
V 1K
Vic
V1F
Vic
V1F
Vic
OJ
en
X)
I
00
-------
tn
ID
1
»
2
J
2
-
9
5
>
,
,n
11
1'
1?
Jl
11
11
*
22?
22'
77 9
Z £
!? 2 "»
77 7
i *
? 2 ?
"* 2 2
22?
2 7 >
? 2 2
22"*
777
£ £ £.
29 1
£ 1
n * -9
' 1 -
4 S e.
7 81°
0 11 17
11 17 10
1? 13 17
1? 12 12
12 12 1?
12 12 1?
12 12 1"
11 10 1C
444
2
2
2
'
?
2
(
**
7
17
14
17
12
1?
I"
9
o
7
~
2
'
7
2
1
1
9
14
1'
1?
I'
10
o
8
P
7
7
2
3
'
0
ft
1
^
8
17
12
12
ir
9
8
8
8
0
T T 7
334
774
,
337
1 ' 2
T -. ?
22'
«» »
. J
889
ir in ^
11 10 8
11 10 R
11 10 1
1" IP o
878
7 e 7
9 3 8
? 8 »
811
7
li h 4
3 ' T
2 > ;
22?
Z7 "»
« «.
8 A «
*
755
7 F E
7 E E
8 ' E
8*7
8 8 P
8 o a
9 »
1 ' ?
' 2
7 7
4 4
4 3
2 2
» 2
2 2
7 7
E E
E E
6 S
E E
5 =-
J C
» 7
P 8
2 2
2223
3733
T T T T
333
* ^ » 3
4333
' > ' 2
2222
2272
Sc c c
r> b o
7776
«; s s G
£=,=,£
P B S S
GERE
SERE
f K 6 6
S S S E
'SEE
2 ? ' 2
3 7
t »
V ^
T J
3 T
3T
J
T T
7 7
? 2
? 2
c c
*> b
EC
b
s s
6£
b
E «
E E
5 «!
E S
P 6
e, P
e 5
? e
? 2
7
JS
,
'
'2
6
E
c
E
S
s
G
r
6
c
3
Vie
V1F
Vl«
Vie
SI1
SI1
SI'
Tl1
si->
Vf
03
-fc.
-------
4
0
«
P
1
*
t
7
>
*»
1
?
1
I
3
-
7
9
1-
4
4
1
>
2
7
2
*
^
'OP,
4
4
3
T
1
1
4
5
7
9
in
4
4
».
2
2
2
2
2
z
»
4
4
4
4
1
1
3
5
7
9
ID '
T T 1
4 » 1
> 7 Z
Z Z 7
22?
> > 2
2 ? '
*> ^ ^
' _ c
Z ? 7
? ? ?
' £ £.
443
431
4 Z T
311
1 0 3
Z 3 4
435
S F 7
S 9 9
9 11 IP
I" 13 1C
8
8
Z
T
Z
z
*»
3
1
1
?
3
4
*
0
10
10
1"
i)
P
-
T
-
3
T
9
2
r
i
z
3
4
c
9
9
8
9
1
1
Z
>
2
t
3
T
Z
3
Z
" 0
z
3
4
. Z
2T
J
3 1
3 3
3 3
2 3
2 2
3 3
3 3
3 I
3 3
2 Z
? Z
7 Z
2 2
' 3
« 4
5 4
4 4
4 3
? 3
3 3
Z 3 3
T T 7
T J y
** T T
J J
» ? 3
333
7 ? 3
T 3 3
373
3 3 T
T T T
333
333
27'
z z »
Z Z -
Z Z 2
33'
333
333
T T 7
33'
3 3 <
333
VI'
\t 1 T
V"
If 4 T
VI
V1T
Vt'
V17
VI'
V17
VI'
V17
V T ^
VI
V1T
t/ 1 7
V 1 r
V T ^
1/17
SI"
«i-
TIP
S19
SIP
*l-
«1B
«!!
'!'
SI"
on
oo
-------
n n
n 10
i? IP
n in
Zn 11
3 3
' 3
* t
6 »
n «
1 1
1 1
1 1
I 1
1 1
1 1
"'no.
2 2
' 2
? 2
1 1
a 3
1C
10
10
in
10
3
3
0
1
1
1
1
1
1
2
2
2
2
0
n i"
1" 10
IP 10
11 10
n i"
«£,
4
a 2
2 2
' -
2 1
T 1
13
If?
O
in
0
"
n
a
2
2
1
1
1
2
2
1
0
7
1C in 1"
« - e
8 »< 9
° e e
091
(i ft f
P c 5
7 C 0
889
881
7"> 'y
£ t
?> *
'
3 'T 3
2 T I
IOC
172
2 ? 2
2 ? 2
? , ,
0 " 1
» ;» ?
333
8
9
e
8
A
l|
0
8
a
2
2
T
ft
2
2
'
»
T
1
»
e * r i «
7 c » »
7 7 C 5 *
7 ; fi c <;
a « 0 2 2
0 1 ? 2
B > 1 2 2
8111'
IT 9 ^ ^
J. t ' '
11779
1 C ' t
1 2 2 ? 1
2 " » ' 2
333^3
32333
^ T T T T
J 5 ^ 3
2 3 Z 3 J
? ' ' ' 3
>>> i o
2 2 r i i
1 " 1 1 2
2 1 ? 7
»»>-> 3
T
«.
e
S
2
?
2
2
1
1
T
3
3
1
2
?
2
3
3
3
<4
5
2
2
7
2
2
?
3
3
T
5
3
2
2
2
2
?
3 3
? 3
U 4
V 1 ^
V18
If 1 B
VI"
VI"
Vl«
511
'19
?!
''I"
Tl"
C/l
oo
-------
tv>
1
ft
c
3
t"
1*^
'
i;
T
3
T
4
9
9
3
R
e
T
1
1
1
1
9
1C
1 7
J /
17
12
T
3
T
4
0
8
8
3
8
1
1
1
1
1
2
ft
9
11
12
T "9
1 £
4 )
1 £.
12
11
3
3
3
4
a
B
8
-------
" '
2 2
- ?
1 2
p P
« ^
JL £
» 7
S 7
9 9
1" 1"
J" 11
<* 11 12
O4
11 11
n n
11 in
* 1 U
in ID
3 3
T 7
P n
U
7 7
3 3
8 f
1 I
2 2
2 '
4 2
7 1
C 1
4 4
7 7
9 13
11 12
12 12
12 11
12 11
11 1"
9 9
4 4
3 C
Qf
r s
3 3
8 »
1 1
2
2
n
1
2
4
7
11
12
11
11
IP
q
3
4
p
f
8
«
1
1
'
2
"
1
2
5
S
1"
11
in
n
n
S
7
4
n
8
1
1
1
> t >
222
" 1 >
IIP
T 9 »
e r (,
655
n 7 7
937
877
937
8 ' E
7.5 S
777
4 e 5
C 0 *
88"
1 ?
222
122
121
* 2 ?
, , ,
111
O n n
111
* 3 3
5 4 »
5 4
"554
764
-en
c E r
7 6 -
7 E =
5 »5 c
« 5 "
*^ n o
111
111
111
1 3 1
2
1
1
5
1
2
2
?
?
7
T
I)
f
b
r
e
i
i
i
i
T» ?
i i
p p
i i
9 ^
2 2
2 3
- 3
? 2
> 2
7 2
4 3
4 4
E ?
5£
.
c 5
** n
1^
c
* 7
*. £
» 1
1 1
1 2
1
"
1
2
2
4
6
c
4
3
3
3
4
4
4
5
D
1
1
1
2
1
»
1
2
2
4
S
5
4
3
'
3
3
4
4
5
0
2
2
i
2
1
0
1
2
2
2
T
T
3
3
7
7
3
4
P
2
2
3
2
2
0
1
2
2
Z
2
2
i
3
I
3
3
3
3
3
P
2
2
3
3
3
2
"11
1 ' »
? ' '
2 -. »
2^ 4
»;
22'
, , ,
7, 2 '
T ' '
33?
33,'
3 3 3
777
3* f
J ?
3 * '
> ** 9
7 Z Z
2 *> "
' 2 2
2*) *
,
?5" 9
£. e
2*9 »
333
3 » T
333
2 ?
5-0
20
-------
T11II1??212112??22222*
1111112771111222272?'''
1I1111222111122-22333*
lllll'l22??ll-?223333*"
V2P
V'C
V2"
V7?
00
-------
PLOT or NOZ , LEVEL i" 4:05 CYCLE i ie90 PARCELS
o\
I .6
! .6
3.3
1 .4
3.2
.8
.6
i.s
.7
2.2
1 .5
I.S
1 .5
1.5
1.5
H.O
6.S
4. 1
S.7
3.2
5.5
3.7
3.0
3.0
3.7
3.0
3.7
3.0
2.2
S.O
6.9
*.7
5.3
3.3
3.2
3.0
S-.S
2.2
5.5
2.2
2-2
1 .5
5. 1
7.9
6.7
5.9
1.9
5*8
3.9
1.7
3.8
1.0
S.5
2-1
5.1
2.2
3.6
8.S
6.0
6.7
1.9
S.7
3.2
3.2
s.o
3.2
1.0
1.0
2.2
5.2
6. 1
a. 8
8.7
6.9
5.9
5.7
3.9
3.2
1.0
1.0
1.0
1.8
3.2
.7
1.1
8'8
9.7
6. 1
8.7
3.3
3.2
S.7
s.O
1.0
1.0
3.2
1.0
2.1
3.5
8 t 8
7.9
8.7
7.0
1.7
1.0
3.2
1.0
1.0
1.0
1.0
1.0
1.0
2.6
9.7
8.7
7.1
5.1
3.9
1.8
1.0
3.2
1.0
1.0
1.3
1.0
2.1
2.6
7.8
6.9
6.0
6.8
1.0
2.4
1.0
1.8
1.0
1.8
3.2
3.0
1 .5
3.S
6.9
6.9
5.9
3.3
1.8
1.0
3.2
1.0
1.0
2.2
3.7
.7
1.3
S.I
t.3
2.5
5.8
3.2
1.8
1.8
3.2
3.2
3.0
2.2
.7
1.7
5.0
5.8
7.6
1. 1
5.9
1 . 1
1.0
3.2
1.0
1.0
1.8
1.0
1 .6
2.5
S.O
S.O
1« 1
5.8
8.6
S.O
S.O
5.8
1.8
S.O
S.O
1.8
1.6
2.1
S.O
5.1
5.3
S.O
S.9
6.9
S.O
S.O
1.1
S.O
3.3
5.0
S.O
S.O
.8
I.S
1.2
6.8
6.0
S.O
5.0
5.9
5.9
S.O
5.8
S.O
S.O
S.O
6.0
S.I
.7
S.O
7.0
S.I
k.'a
S.O
S.O
6.0
6.0
6.8
1.1
5.0
6.0
6.0
6.0
2.2
1.2
7.0
7.9
6.0
s-o
5.8
6*8
5. 1
6.0
6*0
6.0
6.0
6.0
6.0
1-5
3.3
7.Q
6.1
6.0
5.8
1.1
6.0
6.Q
5.1
6.0
6.0
6. 0
6.0
6.0
2.2
S.O
7,0
7.0
6.0
S.9
1.0
7.7
6.0
o.O
6.0
6.0
6.0
6.0
6.0
2.6
2*2
5.0
7.9
o.O
6.6
1.Q
1.B
S.I
5. !
6.8
6.S
6.0
6.8
6.0
7.7
.9
l-S
.6
S.S
1.3
1 .7
2.1
2.1
1.7
1.2
2.6
1.7
2.6
2.6
1.3
1.3
1.7
OJ
oo
-------
NEXUS CODE LISTING AND FLOWCHARTS
BFOR'FLO.FORFLO.S NEXlNC
FUO.YCHARTED er FORFLO
ON
SEP
li-IH'ZZ
»T U.IH.ZZ
ooooot
OOT002 *
000003.
ornu;iH.
CCU005.
coocois.
connj7.
ccnooa.
OOUOO?.
oooo i o .
OI'JOl 1 «
ocoiiz.
osnoi3.
OCOQ1 t«
CSOOlS.
oo(ini&.
oonoi 7.
GOCO 1 &
ccoa 1 1
ccnozc.
ocnpzi
OCf!022«
G~ P fi Z 3
OC002M.
conn2S.
0001126.
OCf.027.
oonnza.
00002''
GG0033.
C C (' 0 3 1 «
Cf r n 3 2 .
i* i. " «* fc
COCO 33.
OCCC3H.
CC0035.
OCCTJ36»
GCC037.
OCC038.
OCC039.
CCCP^C'
CCCOH1.
OCI OM2.
CCCO'13.
CCLO"*1*'
CCC015.
CCCOM6.
OCCOM7.
OCOu1*"*
OCC019.
CCCPSO"
CCCCS1.
NEX|I
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
''PARAMETER Nx-22. NY.,*, Nz-i. NSP-I. MAX-.OOOO
PARAMETER NXMI-NX-I, NYMI-NT-I, NZ"''"-' A
rnHMC.I, OX DY. I)Z. TDX, TOY, TOZ , PM1N, PMAX, «N , QA(
."MNSTRU-SPK'PMXSTRU.SP), -MSTORINSP., BCKSTR(NSP),
2 NV, FNA, FNY. FNZ, V. RV, »', T, KTC,
3 r'A-r'^MAx^^^Ax : «;"M UMFLU;. ,'u1;. FLUXZ.N,
; Jl'S?MNT» ?IuT2IN?** BXl(i,NX,NSP), B X2 ( 1 , NX , NSP ) ,
Cj pmjTl'NTii rL "*'*'*** o « r.
6 BYJIH.NY.NSP), BT2lt.NY.NSP) .,,,
NX NUMBER OF CELLS IN THE X DIRECTION
I1Y NUMBER OF CELLS IN THE T DIRECTION
IIZ . NUMBER OF CELLS [N THE Z DIRECTION
NSP NUMBER OF SPECIES FOR l>HOT°"E"^"p,c ...
HAX - MAXIMUM NUMBER OF PA«CFLS FOR ANY SPECIE
OX - CELL SIZE IN THE X DIRECTION '"ETERS)
nY CELL SIZE I" THE T 0IRECTlON (MET£RS)
DZ CELL SIZE IN THE Z DIRECTION (METERS)
TDX » 2«OX
TOY « 2»DY
ToZ 2»OZ
Prlltl MINIMUM PAPCEL *EIGMT ALLOWED «KG)
PHAX - MAXIMUM PARCEL -HEIGHT ALLOWED (KG)
HN . AVERAGE PARCEL "EIGHT (KG)
BACK BACKGROUND CONCENTRATION (PPM)
UT " TIME STEP (SECONDS)
PMNSTR- STORAGE, BT SPECIES, FOR PMIN
PHXSTR. STORAGE, BY SPECIES, F°* PMAX
>iNSTOR" STORW5E, BY SPECIES, FOR «N
BCKSTK- STORAGE, BY SPECIES, FOR BACK
NPSTOR- STORAGE, BY SPECIES, F°" Np
NV IIX»NY»UZ
FNX « fLOATING NX
FNY FLOATING N>
vNZ : VOLUME"FACTOR ^OR CONVERTING WEIGHTS TO PPM
HV i«'v
NP " NUMBER OF PARCELS
T ELAPSED TIME (SECONDS)
ir~f m riiKRFMT CYCLE NUMBER
CC ' , CO-.CENTI(ATION AKRAY (KG UNTIL CONVERTED TO
CH ' HORIZONTAL DIFFUSION A««AY 00
INC02600
INC02700
1NC02800
INCOZ'OO
1NC03030
INC03ICO
INC03200
INC03300
INC03HCO
INCO^bCO
INC03630
IMC03700
Ii4C03dOO
INC039CO
INCOHOOO
PPM FOR OUTPUT)INCOH1CO
INCOtZCO
INC01300
ITY IN X DIRECTION
ITY IN Y UIHECTION
UIK.J.I.3I VELOCITY IN Z b 1 »*ECT ION _ ___
X PARCEL X COORDINATE ARKAY (CELL FRACTIONS)
Y . PARCEL Y COORDINATE ARRAY (CELL FRACTIONS)
Z PARCEL Z COORDINATE ARRAY (CELL FRACTIONS)
._ . PARCEL WEIGHT AhRAY UG)
FLUX1 FLUX OUT Y - |. SIDE OF VERTICAL GRID (KG)
IIMtU^'UU
INC01BCO
INCOH9CO
1NCOSOOC
1UCOSIOO
00
50
i
00
-------
OCCOS2.
cocosi.
CCC05H«
CCOOSS.
CCC054.
CCC057.
CCC058.
CCCOS9t
CC0060.
cccn<. i.
CC0063.
I LUX2 FLUX OUT Y » NT SIDE Of VERTICAL GRID (KG)
FLUYl FL^X OUT X.I. SIDE OF VERTICAL GRID (KG)
FLUT2 " FLUX OUT X » NX SIDE OF VERTICAL GRID (KG)
BXI BORDER CELLS THAT STORE POLLUTANTS GOING OUT Y
BX2 " BORDER CELLS THAT STORE POLLUTANTS GOJNS OUT Y
BYI BORDER CELLS THAT STORE POLLUTANTS GOING OUT X
It
NY
BY2 BORDER CELLS THAT STORE POLLUTANTS GOING OUT X NX
5 SOURCE ARRAY (KGI
DIMENSION siNZtNr.Nx)
EQUIVALENCE IS(DH)
END
INCUSiOO
INCOS300
INCOS100
SIDE OF GRID
(KG)
SIDE OF GRID
(KG)
SIDE Op GRID
(KG)
SIDE OF GRID
(KG)
INCU6JOO
INCOAMOO
INCO&SOO
INC06600
C/5
?3
I
CO
-------
eFOB«FLO.FORFLO,S |0
FLOACHARTED BY FORFLO /XBOS/ ON 23 SEP 7i AT 11:11:27
GOUOOI i
C00002.
C00003.
COOOOM.
OCOQUb.
OOCOOA.
OOU007.
conooSt
OOHOD9.
000010.
oonoi i .
000012'
000013.
OCOOIM.
GOCI015.
OOC1016.
oonni?.
COOOIR.
000019.
000020*
10
c
c
c
c
c
c
c
c
c
c
c
c
c
c
E>
PROC
COMMON
1 CONC,
ISP
SP
NCYC
f.OL
TDELT
HTIME
CONC
PART
«-GHT
UPSP
PSTRT
SAVE
. IO.UST
EFFICIE
INTE6CR
LOGICAL
10
ISP(NSP), SPINSP), NCYC.
i CONC, PART, WGHT, fpsp, RSTRT, SAVE,
NOL, TOELT, PT1ME,
IOLISTI10)
10CD0100
10000200
I000030C
NAMES OF POLLUTANT SPECIES TO BE PRINTED 1COOOHOO
NAMES OF POLLUTANT SPECIES 10000500
FREQUENCY OF PRINT CYCLES (DEFAULT I I IOC00600
NUMBER OF LEVELS UN THE z DIRECTION) TO BE PRINTED (DEFAULT-NZ)
TIME INTERVAL FOR PRINTING (DEFAULT « 0«) lOQOOBoO
ACTUAL TIMF FON PRINTING IOOOOVOO
FLAG FOR DUMPING CONCENTRATIONS ON UNIT 2 (DEFAULT .FALSE.
FLAG F°H DUMPING PARTICLE POSITIONS ON UNIT 28 (DEFAULT .FALSE.
FLAG FOR DOMPINC, PARTICLE WEIGHTS ON UNIT 2B (DEFAULT
NUMBER OF SPECIES TO BE PRINTED
RESTART FLAG (DEFAULT « .FALSE.)
CYCLE TO PE DUMPED ON UNIT 2V FOR RESTARTING
ARRAY NECESSARY FOR USE IN S3TRAN, S3'S VERSjON OF AN
EFFICIENTLY BUFFERED INPUT/OUTPUT PACKAGE
INTEGER SP, SAVE
LOGICAL CONC, PART, WGHT, RSTRT
'FALSE.
10001300
I0001SOO
100015QO
10001600
10001700
10001803
10001900
10002000
i
c»
-p.
-t.
-------
9FCR«FLO.FORFLO,S
FtOftCHARTEO BT FflRFLO /XflCB/ ON 23 SEP 71 AT li:ii:30
00
oooeni
OOOQU2.
000003.
INCLUDE NEXINC. LIST
INCLUDE 10, LIST
OOOH07.
oououa.
000009.
OCU010.
oouoi i.
OCDCM2.
COOD13.
OGOCtt.
OCttOIS.
020016.
000017.
OCL'018.
ocnoi**
CCOD20.
OC002 I .
OC0022.
CC0023.
OCCC2S.
CCCC26.
CC0027.
CC002B.
CCU029.
CCU030.
CC003I.
OC0032.
CC0033.
CC0036,
CCCC36.
OOC037.
OCOCi38.
OCC037.
CCOCHO.
OCC'CII
OCOCH2.
OCC.Of 3.
OCtlOHM.
NEXUS, NUMERICAL EXAMINATION OF URBAN SMOGi IS A 3-0 I HENS IONAL,
PARCEL-IN-CELL METHOD CODE
NEX00100
NE*oo20o
NEX00300
.NEXOOHOO
NEXOOSOO
NMU0600
NEX00700
NEJUOdOO
CATA
i, j,
DX,
K
or, o*. FINAL / 2»ttnQ*, ico.t tiiBQ. /
ARE USED AS INDICES IN THE x, r, z DIRECTIONS, RESPECTIVELY
c THROUGHOUT THF CODE, AND N is THE POLLUTANT SPECIE NUMBER NEXUUOO
C iET UP CONSTANTS NECtSSART TO START THE RUN AND REAO INPUT NE*01200
C CARDS DEFINING PRINT AND DUMP OPTIONS NEX01300
CALL INPUT NEXOIIOO
c INPUT INITIAL CELLULAR CONCENTRATIONS AND SET ARRAYS NEXOISOO
C IF RSTKT FLAG IS SET, SUBROUTINE RESTRT IS CALLED FROM INPUT NEXOI400
IF (.HOT. RSTRT) CALL SETUP N£*OI700
DO 5 N * t * NSP NEAOIttOO
bACK BCKSTRtN) NLX019QO
c NEGATIVE ARGUMENT TO OUTPUT INSURES PRINTING INITIAL CONDITIONS HEXCJZOOD
b CALL OUTPUTI-NI N£ KNSTOR(N) NEX02SOO
PMIN « PMNSTRfN) NEX02400
PMAX FMXSTR(N) NEX02700
BACK * PCKSTR(N) NEX02800
c COMPUTE PSEUDO VELOCITIES NE*02?oo
CALL DIFFUSINI NE^OSOOO
IF (N .EQ. IT T . T + DT NEX03100
C MOVE AND REWElGHT PARCELS NEX03200
CALL PARCLL(N) NEX03300
C KEEP TKACK OF ALL POLLUTANTS MOVING OUT OF THE MAIN GRID N£*03HOO
CALL BURDE"(N) NEXOSSOO
c INPUT NE" PARCELS
CALL SOURCEf'l)
CALL OUTPUTIN> NExmaco
NPSTOK(N) » NP NEXQ3900
20 CONTINUE NEXOMQOO
c CHECK TO SEE IF MAXIMUM RUN TIME is APPROACHING. IF so, GO To JO.N£*O'UOO
CALL SSWARHIS, i, $301 NEXUSZOO
If CICYC '£9. SAVE* GO TO 30 «£XDH300
IF (T «LT. FINAL) GO TO 10 NEXOMHDO
30 IF (SAVE.ES.O) GO To tO NEXOtSOO
OCC.OM6. CALL V.MLE NEX01600
COOOH7. PRINT 60SO, ICYC NE»01700
CCOOS8. 60HO FORHAT (OCYCLE', 16, DUMPED ON DRUM FOR RESTARTING*) N£*01BOO
OC0019. to T . -T NEXOSVOO
OC0050. DO 50 N . I, NSp NEXOSOOO
OCOOSI. C THIS WILL FINISH OFF ANT DUMP TAPES . "EXO&ICO
CO
CO
-------
000052. ' SO CALL OUTPUTt-
-------
NEXUS
FLOKCHARTED BY FOHFLO /XBOS/ ON 23 SEP 7i AT n:is:Jo
: IIEGIN :
\ ./
i
i
INCLUDE NEXJNC, LIST : :
»
I
1
INCLUDE. 10. LIST : !
._C .»«.3
(..........,........]
C'NEXOOHCO 3
:
-C NEXUS, NUMERICAL EXAMINATION OF 3
CURBAN SMOG, IS A 3-DIMENSIONAL, 3
CNEXOC60C 3
C PARCEL-IN-CELL METHOD CODE 3
CNE»007GO 3
C ....................3
C 3
C'NEXOOSOO 3
DATA DX.-oT, DZ, FINAL / 2'HM70., 100.. '
61180. / :
-C It J. K ARE USED AS IND]CES IN
LTHE X. T, Z DIRECTIONS.
CRESPECT1VELNEXC1000
c THROUGHOUT THE CODE. AND N is
CTME POLLUTANT SPECIE HUfietR
CNEXOtIOC
C SET UP CONSTANTS NECESSARY TO
CSTAKT THE RUN AND READ INPUT
CNEX01200
c CARDS DEFINING PRINT AND DUMP
COPTIONS NEXOI300
CALL INPUT
C INPUT INITIAL CELLULAR 3
(CONCENTRATIONS AND SET ARRAYS 3
CNEX01SOO 3
c IF RSTRT FLAG is SET, suBKOuTjN3
OJ
en
oo
c*
-pa.
-------
I i <ini i > IALLC.U rnuri inrvi j
I UEXCI400 3
I
I
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
*
A
1 TKUE
I
: : CALL SETUP : :
i
i
i
i
: BACK ECKSTR(N) :
C NE6ATIVE ARGUMENT TO OUTPUT 3
ClNSURES PRINTING 1NITJAL 3
CCONDITIONS NEX02000 3
C S3
1
C 103 1
: CONTINUE ;
i
i
i
i
NP HPSTOR(N) :
an m rtNSTOR(N)
PMIN « PHNSTR(N)
PMAX » PMXSTRIN)
BACK « BCKSTRIM)
1
I C COMPUTE PSEUOO VELOCITIES 3
I CNCXC2703 3 >
I
I
: : CALL oiFFUS(N) : :
i
i
I TKUE ' 1
*
z
2,
Z
j
I
I
Z
z
z
z
Z ^j
Z CO
z '
z <»
z -^
z
z
z
z
z
^
z
-------
T T + or
I --- C HOVE AND REWE1GHT PARCELS
I
I
0<-
1
CNEXOJZOO
CALL PARCEL(N)
-C KEEP TRACK OF ALL POLLUTANTS 3
C MOVING OUT OF THE MAIN GRID 3
CNEX03M03 3
CALL BONDER « NP
1
203 I
CONTINUE :
I CHECK TO SEE IF MAXIMUM HUN 3
CTlHE IS APPROACHING. IF SO, GO 3
CTO 30.NEXOH100 3
I
CALL 53w/»KHibi i, i jo> : :
i
i
I FALSE
I
1
I FALSE
I
Z'
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
Z
2
Z
Z
Z
Z
Z
2 f/>
z co
Z -P»
T rs
X Z
X Z
z
z
-------
C 301
1 FALSE
I
1
: : CALL *MLE ; s
i
i
/ PRINT 60HO, Ic*c /
I
I C *OfO FORMAT COCVCLE'i 16, 3
I CD^MPED ON DRUM FOR RESTARTING1) 3
I
C H03 I
: T -T :
I
I
---C THIS WILL FINISH OFF ANT DUhP 3
CTAPES NEX05IOO 3
A C 503
I
I
*: STOP :
i
z
z
z
I
z
z
z
z
2
z
z
z
z
O-J
'en
jo
i
00
-------
SFCR.FLO.FORFLO.S BOROE"
FLO»CHARTED BT FQRFLO /*aos/ ON 23 SEP
AT ii:it:35
oooooi
C00002.
000003*
GCOOOM.
OCOC05.
coooo*.
OC0007.
CCOC08.
COC009.
ooooio.
oonoi I
00001 2«
000013.
OCOOIt.
CCOOI5.
CCR014.
CC10I7.
CC0018.
OC0019.
CC0020.
000021 .
OC0022.
C00023.
000021.
GC0025.
OC0026.
CCC027.
CCOo2S.
C00029.
G00030.
0=3031 .
003032.
C00033.
CC3031.
CC0035.
OCC.036.
000037.
000038.
OC0039.
CC5Q10.
OCC01 1 .
CCC012.
COOOS3.
np p i ti u
w> i* i' O «
OOGO-iS.
CC 00 H 6
000017.
OC0018.
CS0019.
OCCOSQ.
OCOoSI.
C .
c
c
20
25
30
35
so
15
.
SO
too
220
22S
SUBROUTINE BOROE.RIN)
THIS ROUTINE UPDATES THE BORDE" CELLSi EACH OF
T»ICE AS HIGH AS THE CENTRAL GRID CELLS. MASS
CELLS IS RIMOVED FROM FURTHER CALCULATIONS.
INCLUDE NEMNC
00 ICO I I. NX
VEL - uii,i,i,2>/3.
IF (VEL 'SE. 0-I GO TO 25
FXI » FLUXI(i)
DO 20 L « I, H
FX >BXl(Ltl,N)*VEL
BXIILiI.MI BXI(LiIiN) - FX * FX|
FXI - FX
FLUXII I) 0.
GO TO 35
FXI BXi(1,1,N)«VEL
FLUXI(I) FXI
CO 30 L » 1, 3
FX « 6*1/3.
U (VEL-LE.O.I GO TO <»S
FX1«FLUX2(I)
00 MC L«l it
FX'BX2(L,1,N)*VEL
BX2(L,I,N)«BX2(L,I,N)-FX»FX1
FII-FX
FLUX2(I)«0«
GO TO ICO
Fxl.-3*21 II,Nl«VEL
FLUX2(I)*FX!
00 SO L'li3
FX<-BX2(L*I!I,NI«VEL
8X2(L.I|NI«BX2(L,I,N)-FX1«FX
f»I«FX
BX2(HlI,N)'Bx2<1iliN)*FXI
CONTINUE
DO 200 J « li NY
VEL > uii,j,i,i1/3.
IF (VEL >GE. 0<) GO TO 225
FXI FLUTl(Jl
DO 220 L « I i 1
FX -BY1(L,J,N)»VEL
BYKLtJiN) BYKLiJiU) - FX « FXI
FXI « FX
FLUYI(J) 0.
GO TO 33S
FXI » 8Y1(I,J,N».VEL
BORUOIOO
IS 3 TIMES AS'WIDE,
REMOVED FKOM THE OUTER
BOROCHCO
BOHOOSCO
eo«o0400
BOKUC7CO
BOROC8QO
BOK00900
BOKOIOOO
BOK01100
BOKOI2QO
BOH01300
BOnOiMOO
80ROISOO
BOrtOI600
BOR01700
BOxOlbOO
BOKOI900
BOHU2000
80K02IOO
BOK02200
BOK02300
BOK02HOO
B0h02&00
BCH02600
BOh027oO
B0h02800
bOK02900
B0h03000
BOHOilCO
80X03200
BOK03300
BUK03HOO
BOH03SOO
BOK03600
BOK03700
60X03800
BONU3900
BOK01000
60HOHIOO
BOK01200
. BOH01300
CO
-Ci
BOKOH600
80KOt7oO
80x01900
BOK05COO
BOH05100
-------
000052.
GCOOS3.
000054.
OCC,OS>5.
GCOOS6.
GC0057.
OCOnSS.
CCOOS9.
OC')06C.
000041.
CCOOA2.
OC0063.
OC0065.
000066.
OC0047.
CCOP68.
OC0069.
CCD070.
CC0071.
OOL'072.
OOU073.
OC0075.
000076.
CC0077.
- FX1
FLUYl(J) FX1
00 300 U I i 3
FX . B»HL»1.J.N)«VEU
BYHU.J.N) BTl(LiJ.N)
300 FX1 « FX
' BYIM.J.N) BUdtJiN) - FxJ
335 CONTINUE
VEL « U( 1 , J,NX. I 1/3.
IF (V£L .LE. 0«) GO TO 215
FX
00 2HO L«lt1
FX»BT2(L. JiN) .VEL
8Y2(L.J|N)«BY2(L.JiN>-FX«FXl
2HO Fxl.Fx
FLUY2(J)-0.
GO TO 200
215 FXIe-BTZI I .JiN>»VEL
FUUY21J)»FX1
00 250 L-l .3
Fx»-BY2(L+ I i J,N) »VEL
BY2(L,J|N)«BY2(l-,JiN)-FXJ+FX
250 FX1-FX
BY2IHlJ,N)-BY2l1»JtNl-Fxl
200 CONTINUE
RETURN
END '
00R05200
{,0X05300
60H05HOO
BOKObbOO
BOK05600
BOK05700
BOKObBOO
BOK05VOO
BOR06000
80*06100
BOR06200
80R06300
BOK06400
B0r<06500
BOKU6600
BOKU6700
BOR06900
BOK07000
BOK07 100
BON07200
BOK07300
BOR07HOO
BOK07500
BOK07600
60M07700
I
oo
-------
BORDER
FLOWCHARTED BV FORFLO /XBQB/ ON 23 SEP /j AT 11:11:37
SUBROUTINE BOROCh(M)
C THIS ROUTINE UPDATES THE BORDER]
CcELLSi EACH OF WHICH IS 3 TIMES 3
CASBOROOZOO
c TAICE AS HIGH AS THE CENTRAL
LGRIO CELLS. MASS REMOVED FROM
CTHE OUBOR00300
C CELLS IS REMOVED FROM FUKTHER
CCALCULATlONSi BOROOHUO
A
A
A
A
fl
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
: ; INCLUDE NEXINC : :
i
i
i
i
: VCL « ui 1 1 1 , i ,21 /.! :
i
i
I FALSE
I
I
: FXI » FLUXI 1 1 1 :
i
i
B 1
B I
B : FX « -BX i (LI i tN)*vEL :
B : BXI(LtltN) PXllL.liNI - FX » FXI :
B 1
8 C 203 1
I
I
: FLUXI ( i) o. i
z
z
z
z
2
z
z
Z C/4
2 '
2 Tp
1 J^
Z
Z
I
z
2
2
Z
Z
-------
n
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
B
B
Q
B
B
g
B
e
£
B
B
g
B
B
6
B
B
I
C 253 1
: rxj exi < i ,I,N)»V£L ;
: FLUxm ) - FXI :
I
i
I
i
i FX . Bxl (L.I ,I,M«VEt :
: BXl(LiI»N) - BXlCLiIiM - FXJ * FX I
1
C 303 I
I
I
; BXlCtil,M)»BXlH,l,N).FXl I
I
C 353 1
: CONTINUE ;
: vtu-iM i ,NY , i ,2 1/3. :
i
i
1 FALSE
I
I
: FXJ.FLUXZU > :
1
J
I
I
: FX'BXZIL,! ,N)»VEL :
: SX2
-------
oo
A
A
A
A
A
A
A
A
A
A
A
A
t B
A 6
A B
A 8
A B
A 8
A B
A B
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
1
: FLuxz(i>>o< :
i
i
C 153 1
: FX|»-BX2( 1 i I ,N).VEL :
: FLUXZI i >«FXI :
i
i
i
i
: FX=-BX2IL»I i 1 »N)«VEL :
: BX2(Li I |N>«CX.ML,1 ,N>-FX |»FX !
I
C BO] I
1
I
; BX7O, I .N1-BX2IH, 1 ,N)-FXl :
1
C 1001 I
I
I
I
J
: VEL uu iJ.i ,1 )/3« :
i
i
I FALSE
1
I
A
X
X
X
X
X
X
4
» X
X
X
X
X
X
X
x
X
X
A
X
A
X
CO
?a
i
oo
-p*
-C*
fc
z
I
2
z
2
I
Z
Z
Z
1
Z
Z
Z
Z
ft - r i it v i i
-------
*
A
A
*
A
A
*
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
B
0
B
8
B
B
B
B
B
B
fl
B
r * i t C" < i t w t *
i
i
1 s
I
: rx . -BTI it, J,N)«VEL :
: BY1(L,J|N) flYllL.J.M - FX » FXl :
I
I 2201 I
I
I
: FLUYiu) o. :
I
I
C 2243 I
t FXl BYl I 1 iJ,N)«VEL :
: FLUYI i j) « FXI :
1
I
I
I
: FX - BYI /3. ;
I
I
X
X
X
x
X
X
X
X
X
X
X
X
X
X
X
X
x
iw
:?
A t*r\
X Ji.
X 4^
X
z
1
2
I
z
1
2
Z
Z
Z
1
2
Z
2
Z
-------
OO
o
M
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
t
A
A
A
A
A
B
B
Q
B
B
Q
B
B
B
g
B
B
B
B
B
B
t . * \ i «\wt»
I FALSE
I
I
: rxi>FLurzui :
I
i
i
i
: FX.BYZIL, J,NI«VEL ;
: BTZILi J,N)«BY2(L.J|N)-M*FX1 :
I
C 2H03 1
1
1
: FLUY2( JJ-O. :
I
I
C 2HS3 I
! FXI.-BY21 1 . J.NI'VEL ;
! FLOY2(J)-FX1 :
I
[
1
1
; FX--pY2tL»I , J|N)»VEL :
! BY21L i JiN)»B»2(L t J.N1-FX l«FX !
I
C 2503 I
I
I
! BY2IS, J,N)»DY2IH, J.N1-FX1 :
I
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
x
X
x
X
.-.--._---_... -__.._0
z
2
z
z
£
z
z
2
Z
I
i
i
Z
z
2
z
z
z
I
z
z
z
z
i
OJ
C/}
73
1
OO
£>
^
-------
oo
CONTINUE
I
I
RETURN
oo
-------
eK>R«FLO.FORFLO,S DlFFUS . ,
FLO*ChARTED BY F0RFLO /'BOB/ ON 23 SEP 71 AT 11.11.13
03
030001 .
COOC02.
000003.
ocnnoi.
cccoos.
OCC004.
COOU07.
conoos.
COCC09.
OOOUIO.
coon i .
ocnui 2.
GCOOI3.
CGCOIH.
oooois.
000014.
00^017.
OOU018.
CCC019.
OOC020.
000021.
003022-
00(:o23.
000021.
C00025.
000026'
GCC(127.
OD0028.
CC0029.
COI1030'
031)031 .
OCl'032.
CCU033.
GGU031"
OC003S.
GC0036.
GC0037.
OC0038.
000039.
OuOOIO.
OOOOH"
on rt i u y
V j Li (j 1 *
000013.
conoii.
OOOQ15.
OD0016.
000017.
000018.
00001'.
c
c
c
c
(
THE DIFFUSION PSEODO VELOCITY
INCLUDE NEXINC
l"LST* SZ ARE THE LARGEST FRACTION OF A CELL IN ThE X
RESPECTIVELY, THAT » PARCEL «ILL BE ALLOWED TO HOVE IN
DATA sx, SY, sz /3».s /
READ <8) u, ov. OH, DT
CALL S3TRANI8, ZS, 1, 10LIST. IS)
IF IIS «LT. 0) RETURN 0
PRINT 60CO, DT
DO « DT
CALL INVKSIHDV, DM, T, oo>
PRINT 6COOl 00
6000 FORMAT I'OOT «, Fio.3)
DO 10 I 2, NXMI
00 10 J « 2, MYM1
DO 10 K " ll N/
IF (CIK|J,liN) .EO. BACK) GO TO 10 t,.,,,, j i
OVX - DMU,J,ll.lCtK,J.l-l.NI - CCK,J,l*l,N»/tC(K,J,l
OVY . DM«KiJ,l)»tC(K,J-lll»N) * C«K.J»'tIfl
IF N) -
- , 30
20 DVZ * DVINZ.J,I)«N)
30 CONTINUE
CONVERSION FROU M/SEC To FRACTION OF CELL/DT
,1,1) » DVX)/DX
D1FOOIOO
DIF00200
01FOC300
01F00100
, Y, Z DIRECTIONS.
ONE TIME STEP.
OIF00700
DIFOObOO
D1F00900
UIF01000
DIF01100
U1F0120C
Ulf01333
D1FOHOO
DIfOlbOO
DIF01600
OIFOI700
DIF01UCO
OIFOI900
UIF02000
OIf02lOO
Ulf02200
OIF02300
OIF02100
UlF02bOO
D1FQ24CO
,NI«TQX)
,N).TOT)
J,1,N)»TOZ)
Dj)
- C(NZ,J,1,N))/(C(HZ,J,l,N)»D2)
01F02700
OIF028CO
U(K,J,I,2)
U(K,J,I,3)
(U(K,J
UIK.J.I.3)
NZ
NX
U(H,1 ,1,1 I/OX
UU,I , I ,2l/0t
U(K,1 ,1,3)/OZ
UIK.NY.I,1J/ox
UIK.NY, i ,2)/OY
50 U UlKiJiliII/O*
UIK,Jili2) U'KiJi
OOOOSI.
UIK',J|NX,I) >'UIK, J,NX ,l I/DX
U(K,J,NX,2) U< K, J,:iX,2 I/DY
Olf03000
UIF03100
OIF03200
DIF03300
OIF03100
D1F03SOO
OIF03400
01F03700
DIf03800
DIF03900
OIF01000
OIF01100
01F01200
01F01300
DIF01100
01F01SOO
OIFOH600
OU01700
01F01800
01F01900
DlFOSOOO
oirosioo
CM
00
-------
000352.
03JO-.3.
000051.
000-355.
003056.
003J57.
000058.
60 uU.JiNXi3>
i Ji NX«3)/D7
00 73 I
DO 70 J
00 70 K
70 C(K,J,I,
RETUHN
END
NX
i. NY
i, NZ
o«
OIF05200
01FOS300
DIFObtOO
UIFOb6QO
OU05730
oirossoo
oo
OJ
CO
73
I
CO
-------
OIFFUS
FLOflCHARfEO
FO"FLO /X808/ ON 23 SEP 7| AT
00
\ SUBROUTINE OlFFUStN) /
\ /
I
1 C THIS ROUTINE CALCULATES THE 3
I CoIFFUSION PSEUDO VELOCITT ' 3
I CDIF00200 3
I
: : INCLUDE NEXIUC : :
* i
i
: : INCLUDE 10 : :
C S», 5Y, SZ ARE THE LARGEST 3
LFKACTION OF A CELL IN THE X, 1, 3
CZ OIREC01F00500 3
c RESPECTIVELY, THAT A PARCEL 3
CW1LL BE ALLOWED TO HOVE IN ONE 3
CTIHE STOIFOO&OO 3
: DATA sx, sr, sz /j«.s / :
i
i---c READ <8> u, DV, DH, DT oirooeoos
i
i
: : CALL S3TKANO, 25, Y, lOLISTi IS) : :
:.: ..............:.:
i
i
/.,..,,. ,,,.,,, , ,,,.,.,i,..i,.,...,t,,,,,,x FALSE
< IF (IS .LT. 0) > ...
\ /
I TRUE
I
/,.»,,, i ,t , » .x
: RETURN o
N.................../
0< . ........................
I
/,....,..,..,..,,,.......,.,.,t.,,t .,,,,,,/
/ PRINT 6000, DT /
/ ......../
>
I
cn
to
CO
DD DT
-------
CALL INVRSIUDVi OH, T, DD)
00
A-
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
~
a
B
D
tl
B
B
B
B
B
B
B
D
B
B
B
B
P
B
B
B
B
B
B
B
B
B
B
B
8
B
B
B
B
8
e
e
r.
e
e
B
B
B
B
e
B
B
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
1
I
/ PRINT 6000, DD /
1
1 C AOOQ FORMAT COOT », FI0.3) 3
I
I
u t t u
j
I
I
* * * * * f
I
I
I FALSE
1
I
OVX DHIK.J, I )*(C(K, Jt I-l |NI - Q(K,J,I+ :
i INI ) / (C(»ci j, i ,N)*TD*) :
DVY « OHlKtJt 1 )«(CU, J-l , 1 ,N) - C(K,J*1, !
I,N))/(C(iC,J,l,M>TDY) :
1
I
I FALSE
1
I
1 FALSE
1
I
: ovz » DVIK , j, 1 1 .(CIK-I i j, i ,NI - ctK»i,j, :
* I |N) )/(C(K. J, I ,N) »ToZ) :
i
i
.
u
u
u
(J
u
u
u
\J
s u
s u
X
X
X
X
X
X
x
X
X
X
X
X
X
z
z
z
z
z
z
z
z
z
z
z
z
z
z
2
Z
z
z
z
z
z
z
z
z
7
CD
1
no
-------
oo
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
A
A B
A B
A
A B
A
A B
A 6
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
APC
A p C
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ABC
ADC
ABC
ABC
L I U» J i - »»
s u
: DVZ ovii,j,ii.ic(iiJ.i,Ni - CIZ.J.I.N)) : s u
: /icii,j,i,N>»DZ) : s u
: : s u
I s u
i s u
/.. I .........I...' I i\ SU
: 1,0 10 30 :----....._---..-....___..-.-_.--.---.-..--.-__---..-.--.-....__v u
\ .../ su
s u
s u
0<: ........ . ........... ............. o
C 203 I ' S
S
; DVZ - DV(NZ,Ji I }*(C(NZM| i Ji I ,N) - CINZiJt ' S
: i|N))/ic2> + DVYI/OY !
U11111,t,,,.,,,,,i,i,,,i.i,i..,,.,,,,,N FAUSE
< IF (ABSIUIKiJiIf I )> -GE. Sx ) > . ....
\ ..../
I TRUE
1
: uUfj.i.i) SIGUISX,UIK,J,i,ii> :
, o< ...............
I
/.tt,....,,,,,,...,.,,,., ,,,..,,,,,N FALSE
< IF (ABSIUIKtJfI.2)) -GE. SY) > ..........
\ /
I TRUE
I
: l' » SIGN(SY,U|K,J,I,2)) :
I
0< .. .. .........
I
/ ,......,.. ............. x FALSE
< IF (ABSIUIK,JiI.3)I »G£. SZ) >-----......
\...... ...,./
I TRUE
I
: U SIGN(SZ.U(K,J,I,3) I :
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
Z 0-J
z c/>
Z X)
z '
Z 00
z *
z *
z
z
z
z
z
z
z
z
z
-------
OO
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
B C
B C
F C
n
B
B
B
D
R
B
B
n
B
R
0
H
B
B
B
e
8
B
B
B
1)
R
B
B
B
B
a
8
3
B C
a c
8 C
B C.
C 103 I
1
1
1
1
I
1
U(K,l ,1 ,1 ) - U(K, | , I , 1 I/DX :
u(Kt 1 1 1 i2i « uu, i , i ,2>/or :
UIK, 1 , I , 31 » UtKilil .3)/Dz :
U » u«i j, i i ii/oz
: UlKiJ.NXil) U(Kl JlNX i 1 ) /OX
: U(K,J,tJX,2| » U(N t J (NX ,2 ) /OY
1
C 603 I
1
1
1
I
I
I
1
C 731 1
.
.
O-l
CO
' 1
CO
£»
*
-------
C
-------
OFO*«Ft_O.FORFLO,S INPUT
FLOHCHARTEO BY FQRFLO /XBOS/ ON 23 SEP 71 AT 11:11:19
oo
10
003001 »
033002.
OOUU03.
303D01.
033005.
033006.
000007.
ocoooa.
C33009.
000010.
C3301 1 «
030012.
00:013.
003011.
003015.
C w J 3 I 6 .
03 J017.
003018.
C030I9.
C3C020'
COOQ21.
C00022.
C30023.
CC002H.
000025.
C00026.
C03027.
C33028.
013029.
003030.
003031 .
C03332.
033033-
003031.
003035.
033036.
C00037.
C03038.
000039.
C00010.
SOOOH i .
000012.
003013.
003011.
003015.
003016.
003017.
030318.
033019.
C00050*
030Q51.
C
C
C
C
C
SUBROUTINE INPUT
THIS ROUTINE SETS UP CONSTANTS AND ARRAYS NECESSARY FOK
AND INPUT/OUTPUT, THEN READS FREE-FORM DATA CARDS WITH AN AI
ANALYZES THEM AND SETS THE APPROPRIATE VARIABLES. SEE POP
ELEMENT 10 FOR DEFINITION OF VARIABLES. DEFAULT VALUES ARE
SUPPLIED FOR ALL NECESSARY PARAMETERS.
INCLUDE NEXINC
INCLUDE: 10
DIMENSION INI80), lNP{80)
DATA NCYC. NOLI CONCI PART, RSTRT / i, NZ, 3».FALSE* /
DIMENSION A(2)> IHSP(NSP)
DATA A /2«'NO' /
DATA ISP /!/
N V N X ;l Y » N 2
FLO(2i16 IDLISTl1) 3*NV
FLOI2,16 IOLIST(2) NV
FLD(2,16 lOLlsT(3) Ny
FLDI2,16 lOLIsT(l) > 1
FLOI2,16 IOLIST(5) 2
FLOI19,13,IOLISTI1 ) * LOCIU)
FLOc18,13.IOLIST(2 I) « LOC(DV)
FLDI16,|3,IOLIST(3)) LOCIOri)
FLDI 18, IB, lOLlSTd) ) LOCIOT)
INPOOIOO
BOTH COMPUTATION
FORMAT.
INPOOIOO
FLDIIB i a, IOLISTIs))
FLDIte is, IOLISTI*) »
FLD(18 18,IOLIST(7))
FLDI 18 14, lOLlSTISI)
FLDI18 18,IOLIST(13))
IOLISTI9) " IOLISTI5)
V»I./(1.23£-«
LOClNPI
LOCIX)
LOCIY)
LOCIZI
LOCI.»T>
RV
T0X
TOY
TDZ
FNX
FNY
FNZ
P«INT
i ./v
2.«QX
2.»OY
2.«OZ
NX
NY
NZ
6000
6000 FORMAT IM1NPUT C*"DS:')
10 READ IS.SOOO.END-IOOOI INP
5000 FORMAT (80*11
PRINT 6001 i INP
iOOl FORMAT I'O', 80*1)
N . 0
00 20 I I, 80
IF IINPII 1 «EQ« 'I GO TO 20
N . N 1
INIM INPI I I
20 CONTINUE
IF IN «EQ. 0) GO TO 10
IF IINII) >Ea. '0') GO TO too
1NP00600
INP00700
INPOCSOO
INP00900
INPOIOOO
INP01 100
INl>01200
INP01300
1NP01 100
INP01SOO
INP01600
INPOI /OO
INP01600
IflPO 1900
1NP02000
INP02100
INP02200
INP02300
IN»02100
INPU2600
INP02700
INP02800
1NP02900
INPU3000
INP03100
INP03200
INP03330
INC03HOO
INP03SOO
INPOJ600
IHP03700
1NP03800
INP03900
INPOIOOO
ItlPOIlQO
INP01200
INP01300
INi'OIMOO
INP01SOO
iNKQIoCO
INP01700
lNi»01dOO
INPOIVOO
INPOSOOO
.INPOblOO
OJ
in
CO
-------
000052.
030053.
OOOObH.
000055.
000056.
030057.
oooosa.
03HOS9.
000060.
033361 .
0330*2.
033363.
03H06H.
033065.
003066.
033067.
OOJ3&8.
030069.
00'JT70.
033371 .
OCJG72.
000073.
00007H.
003075.
C03376.
003077.
003073.
003379.
033083.
OGQU8 1 .
000082.
C03383.
C 0 3 0 8 1
030035.
000036.
000037.
000008.
030!B9«
000090.
C00091 .
0300V2.
000093.
OOOQ9H.
000095.
000096.
033097.
00309").
000099.
030] 00.
oooioi.
000102.
COPI03.
COOinl.
OCOlOS.
CCC106.
CC0107.
CCOlOfl.
IF ( INI 1 1 .E«.
IF i i'
CONC « .TRUE.
GO TO 10
110 IF (I .NE. *P« I
PART .TRUE.
GO TO 10
300 DO 310 1 1 i N
IF (i'i( ) .GC.
IF ( i N ( i . F a .
IF (1N( ) -E9.
IF (1N< 1 .FQ.
GO TO 9??
310 CONTINUE
GO TO 999
320 NCYC » 0
I I - 1
DO 330 K 1 1 , I
330 NCYC NCYC«IO »
IF (NCYC .EQ. 0)
GO TO 10
350 TOELT » o.
IFAC 1
IEXP 0
1 " I - 1
DO 370 K II, 1
IF ( INIK) >EQ.
TDELT « TOELT.IO
IF ( IEXP .LT. C)
GO TO 370
360 IFAC 0
S* ) GO TO 500
R<) GO TO 600
P1 1 GO TO 999
E* .ANo. INI71 .EQ. 'V') GO TO 300
L' .AND. 1N(I»7) .E(|. >E'I GO TO HQO
>,< .OR. K«l ,GE. N»l) GO TO 60
FLDI0.6, 1N(K» I ) )
( I ) I GO TO 70
NOT INCLUDED IN SPECIE LIST.')
TO MO
GO TO 110
GO TO ?99
0* .AND. !N(I) .LE. '?*> GO TO 310
.'1 GO TO 310
S«> GO TO 350
C1 ) 50 TO 320
FLD(0>6, INIK) ) 18
GO TO ?9?
. ) GO To 360
.IFAC » IFLD(0(6tIN(K) ) - M8 ) 1 0 . I EXP
IEXP « IEXP - 1
INP05200
INHOS300
|NP05-(00
INr-05500
INPOS600
!NPOb7CO
INPOS900
INP06000
INP06 1 00
1NPQ62QO
1NPOA300
II4P06MOO
INP06SQO
1NP06600
INP06700
INP06800
I N f 0 6 V 0 0
IMP07000
INP07100
IsPO/200
1NP07 JOO
INP07MOO
INH07SOO
INP07600
INP07700
!Nf07bOO
INP07900
iNPObZOO
INP08300
INf'UfiHOQ
INP08500
1NP08600
iNpoe/oo
ii.poaeoo
1NP08VOO
I NP09000
INP09100
INP09200
INPC9300
INP09HQO
INF'09500
INP09600
INP09700
1NP09800
INP09VOO
INPIOOOO
INPlOlOO
INP10200
1NP1C300
iNPIO'iOO
INPlObOO
INPI0600
IKf 10700
INHIOBOO
ye
CO
-------
000109. IEXP . -i
CCOl10. 370 CONTINUE
oooii! NCYC « o
CCOl12. GO TO 10
COU113. tOO K - I » 5
COOl1H. NOL 0
C001 15. DO 110 i 6, K
000116. "*10 NOL " NOL'IO » r LD 10 16 i I N I I ) I - 18
CCIMI7. NOL MINO(NOL.NZ)
CCf113. GO TO 10
OOl'II9. 500 bAVE 0
000120. DO 510 I 5i N
cccizi* sio SAVE lo'SAvr » FLt>io»6iiNm> - IB
CC0122. dO TO 10
OOM23. 600 RbTKT .TRUE.
C00121. CALL PESTRT
00ri25« GO TO lo
OCCJ26. 999 PRINT 6999, INP
OG0127. 6999 FORMAT I'CBAD INPUT CA«D: ', 80A1I
Ocr,12f>. CAI L E» IT
OR012V. 1000 DO 1010 I » I, NPSP
OOC130. J » IbP(II
000131. 1010 IhSP(I) " SP(j>
OC0132. U IHPSf .oT. 0) PRINT 6010, (1HSPI 1).I«1,NPSPI
oaT[33« 60iO FORMAT (//// SPECIES TO BE PRINTED:* / U3Aio»
OQU13M. If (UOL .GT. c) PRINT 60lli NcL
OCni3S. 6011 FORMAT (// (IUHEER OF L£VELS TO BE PRINTED:', 13)
OOU136. IF (NCYC .GT. 0) PRINT 6012, NCYC
000137. 6012 FORMAT (// CYCLE INCREMENT FO" PRINTING;', 151
OOIH38. If MPELT .GT, C.) PRINT 60i3, TDELT
OOHI39. 6QI3 FOPMAll//' TICE STEP FOR PRINTING:', |PElO«t,
OOLMNC. IF ICONC) All) 'YES'
COflll. IF (P*RTJ A(2) 'YES-
OOCUH2. PRINT 6011, A
oooiM3. 60IH FORMAT c//*x* «DUMP«, lox* 'DUMP PARTICLES' /
COOl-n. » 7H, 'AND WEIGHTS' / / 7Xi A6, 13X, A6 )
DOOMS. RETURN
000116. END
SECONDS')
CONCENTRATIONS',
INP1C930
INP1 ICOO
INPI 1 100
INP1 1200
INPI 1300
INPI MOO
II.P1 15CO
INPI 1600
INPI I 700
INPI I3CO
INPI 1900
INP12000
INP12100
INP12200
INPI2300
INPI2HOO
1NP126GO
INPI2700
INFI 2800
INPI2900
INPI3000
INf 1 3100
INPI3200
INP13300
INPI3HOO
IhP 13SOO
INPI3630
INP13703
INPI3800
INPI 39QO
INPI
-------
INPUT
FLOWCHARTED BY FORFLO /xece/ ON 23 SEP ;i AT n:iH:s2
SUBROUTINE INPUT
C THIS ROUTINE SETS UP CONSTANTS
CAN!) ARRAYS NECESSARY FOR BOTH
CCOMPUINP00200
c AND INPUT/OUTPUT. THEN READS
CFKEE-FORM DATA CARDS KITH AN Al
CFORMAINP00300
C ANALYZES THEM AND SETS THE
CAPPROPR1ATE VARIABLES^ SEE POP
CINPCOHCa
C ELEMENT 10 FOR DEFINITION OF
CVARIABLES. DEFAULT VALUES ARE
CINP00500
C SUPPLUO FOR ALL NECESSARY
CPAKAMETERS. 1NP00600
IO
t-o
INCLUDE NEXINC
INCLUDE 10
:: c
DIMENSION IM80), INPI80)
DATA NCYC, NOLt CONCI PARTI RSTRT
NZ« 3«-FALSE. /
OlHENSIPf, A(2)i IHSP(NSP>
DATA A /2»'NO« /
DATA ISP l\>
NV = NX«NY»N2
FLDI2,U,IOLIST(1)) 3>NV
FLCI2.I*,IOLJSTI2)) NV
FLD(2i1fc,IOL1ST(3)) NV
FLD(2,16,10LIST|SII I
PLD12,U.IOLISTIS)I - 2
/ li
FLOII 6,IB,JOLIST ( i
FLDI|fl»|8,IOLIST(2>
FLClIB,16,ICLIST(31
FLOI 16. 18, IOLISTMI
FLDIisiis,|OL|STisi
FLD(16,16,IOL1ST(6))
FLDI|8.18,|OLl&T<7))
FLDI16,18,JOLlSTIBII
FLDI16,16,ICUI5TI 10))
LOCIUI
LOCIDV)
LOC(OH)
LOC(OT)
LOCINPI
LOCIX)
LOC
> LOCU)
LOCIWTI
10LISTO) - IOLISTI5)
V«1./(1.23E-6»OX«DT«OZ)
RV J./V
TDK 2**DX
TOY 2.«DY
lOi. 2>*02
FNX NX
r w _ n> V
CM
CO
Z3
i
CO
-------
FNZ ' NZ
VD
A
A
A
A
A
A
A
A
A
I
I
/ PRINT 6000 /
I
I C 6000 FORMAT (MINPUT CARDS:*) 3
I
C 103 I
/ WEAD
-------
vo
A
A
A
A
A
A
A
I FALSE
1
1
I FALSE
1
1
I FALSE
1
I
< GO TO 3CO >
J FALSE
1
I
I
I
< EM GO TO **CO >
I FALSE
1
C 303 1
1
I
: NPSP o :
K » 6 :
i
c HO: i
: NAPE :
I
i
i
< GO TO
1 FALSE
1
r r n » «
.
J
J
J
J
J
J
J
H J
n J
N J
L
L
L
L
L
L
L
L
L
L
L
L
I
L
L
L
L
L
i
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
W
Q
a
Q
8
Q
e
Q
g
g
e
Q
6
Q
G
g
Q
Q
Q
Q
Q
Q
Q
0
Q
0
Q
e
g
e
Q
a
n
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
5
s
s
s
s
s
s
s
s
s
s
s
u
u
u
u
u
u
u
u
u
u
u
u
u
u
0
u
\j
u
u
u
u
u
u
u
u
u
u
u
u
u
u
(J
u
u
u
u
u
u
u
u
u
II
*
X
A
X
X
X
X
X
X
X
X
X
X
X
X
X
A
A
X
X
X
X
X
X
X
X
X
X
X
x
X
x
X
X
X
X
X
X
z
z
z
z
z
z
J
z
z
z
z
z
z
z
z
2
2
2
z
z
z
z
z
z
2
2
Z
Z
2
Z
Z
2
Z
1 W
Z *~~*
1 I
z oo
z -^
z ^
z
z
z
z
z
f
-------
10
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
i
i
: i . 6 :
i
C 603 I
: K > K > i » i :
i
i
i
i
1 FALSE
1
I
: NPSP » HPSP » i :
: ISP(NPSP) « i t
i
i
L 703 I
I
1
/ PHJNT 6070, NAME /
C 4070 FORMAT <«3»i A6,-« NOT 3
C1NCLU0ED IN SPECIE LIST.i| 3
: : CALL EXIT : :
I
C 831 1
I FALSE
1
I
F
F
F
t
f
f
F
F
F
F
F
F
F
H
H
M
M
H
H
H
H
H
H
H
H
H
J
J
J
J
J
J
J
J
J
J
J
J
J
J
j
J
J
J
J
J
J
J
J
J
J
J
J
J
J
J
J
J
J
J
J
J
J
J
I
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
t
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
|>l
N
It
N
V
«
Q
g
a
Q
a
Q
Q
e
Q
Q
Q
e
a
Q
a
Q
6
a
a
w
q
w
a
a
a
a
Q
Q
a
9
Q
a
a
a
Q
a
Q
d
g
a
Q
Q
Q
a
e
3
s
s
s
s
s
s
s
s
s
s
s
&
s
s
s
§
s
s
s
s
s
b
5
b
b
s
s
s
s
s
s
s
s
5
5
s
s
s
s
5
s
5
s
s
i
5
S
s
s
0
U
u
U
u
U
U
U
U
U
U
U
U
U
U
u
U
U
U
U
U
U
U
U
U
U
U
14
U
U
U
U
U
U
U
U
U
U
U
U
u
U
U
U
*
X
A
K
X
I
a
*
X
X
X
X
X
X
X
X
X
X
X
n
X
X
X
X
X
X
A
X
X
X
X
X
X
X
X
*
X
X
X
I
2
z
z
2
z
2
z
I
I
I
z
z
z
z
z
z
z
J
z
z
z
z
z
z
z
z
I
I
I
I
I
I
I
I
OJ
en
X)
1
00
fc.
-pi
-------
to
A
A
A
A
A
A
A
A
A
C 100.1 I
: i IN isi :
i
t
I FALSE
I
I
: CONG -TRUE. :
i
i
C 1101 I
1 FALSE
I
I
: PART .TRUE. :
i
i
C 190] 1
I
I
< GO TO 310 >
I FALSE
1
I FALSE
1
1
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
L
N
N
N
h
N
N
H
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
Q
H
H
Q
d
0
u
Q
Q
8
Q
Q
U
U
S
U
a
a
Q
a
a
a
Q
e
Q
a
e
e
e
«
0
«
«
Q
e
(i
(
0
Q
s
5
5
s
b
S
S
s
s
s
s
5
S
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
S
s
s
s
s
b
s
b
s
s
s
s
s
s
u
u
u
u
u
u
u
u
u
u
u
u
u
(J
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
X
X
X
A
I
X
X
X
X
X
X
X
X
X
X
X
X
X
z
2
2
z
2
I
I
I
I
I
I
I
I
z
z
z
z
z
I
1
z
2
z
z
z
z
z
z
z
2
I
Z
z
z
z
I
I
OJ
CO
73
1
oo
-p>.
r*>
-------
A
A
A
A
A
A
A
A
A
A
A
A
A
I FALSE
I
I
I FALSE
I
1
C 3103 I
I
C 3203 I
: (4CTC o :
: i i - i :
i
i
i
C 330] I
I
1
I FALSE
i
i
I 3503 I
TDLLT o.
IFAC 1
IEXP 0
i i - i :
J
J
J
L
L
L
L
L
L
L
L
L
L
U
L
L
L
L
L
L
L
L
L
L
L
L
L
L
t
L
L
I
L
L
L
L
L
L
N
H
N
N
H
N
N
N
N
N
N
ft
H
N
N
N
N
N
H
H
N
N
N
N
N
N
N
N
N
N
U
w
e
Q
e
u
u
Q
Q
g
Q
U
Q
e
u
Q
g
Q
Q
9
Q
Q
Q
Q
Q
Q
g
Q
Q
Q
0
e
q
Q
Q
Q
Q
Q '
5
5
S
S
s
5
S
S
S
s
s
s
5
s
s
s
s
5
s
s
s
S
s
s
s
s
s
s
s
5
S
s
s
s
s
s
5
S
s
s
s
s
u
(J
U
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
g
U
u
\J
u
u
u
u
u
u
u
u
u
u
u
i
2
2
z
z
2
z
z
z
z
x z
z
J
2
I
I
z
1
z
z
z
2
z
z
z
z
z
z
z
z
en
^~t
i
00
-------
CO
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
1
1
1
I
I FALSE
I
I
; TDLLT TOELT«IO.»«IFAC * ( FLD 10 >6 t 1 N (K » :
: - H8i«io.**|cxp :
1
1 TRUE
I
JEXP » IEXP - i :
*
I
l
C 3603 I
IFAC o :
IEXP « -i :
i
C 3703 I
I
I
: HCYC o :
i
i
t ICO 3 I
: K . i » s :
: NOL o :
I
u
L
L
4.
L
t
L
L
L
L
L
L
U
L
L
L
U
L
L
L
L
L
L
L
L
L
L
U
L
IF
U
u
L
C
V
. a
* e
u
g
Q
<]
g
g
e
a
g
a
g
Q
g
a
g
e
Q
Q
N a
N a
N a
N a
N g
N a
N g
N g
g
g
a
g
g
g
g
g
g
g
e
g
g
g
g
g
a
o
s
s
s
s
- s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
i
s
s
s
s
5
s
s
s
s
s
s
s
s
s
s
s
5
i
s
s
s
s
s
s
s
V
u
u
u
u
u
0
u
u
u
u
u
u
u
g
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
0
u
u
u
u
u
u
u
II
X
X
X
A
X
X
X
X
X
X
X
x
X
X
X
A
z
7
7
7
z
z
7
z
z
z
7
z
z
z
z
7
z
z
7
Z
7
Z
z
2
z
Z
z
z
z
2
Z
Z
j
Z
z
7
Z
7
Z
Z
j
Z
z
z
2
Z
z
2
Z
T
-------
vo
A
A
A
A
I
C -MOD 1
I
I
: NOL HIIIOINOL.NZI :
1
C SOtn I
: SAVE o :
i
i
i
C 5103 I
I
I
C 600J 1
: RSTRT » .TRUE. :
1
I
i : CAUU KEiTur : :
i
i
C 9993 I
g
g
a
e
g
g
g
Q
Q
U
Q
a
Q
Q
a
a
Q
a
Q
Q
Q
9
0
g
e
g
Q
Q
Q
4
Q
Q
Q
Q
Q
Q
g
Q
U
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
s
5
s
s
s
5
s
s
s
5
s
s
u
u
o
V
u
u
u
u
u
it
u
u
u
u
u
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
Z CO
Z ff^~^
'z
I J^
PRINT A199, JNP
-------
o
o
i c,
i
i
CALL E»IT
A
A
A
A
A
1
C 1COCJ I
1
: j « ISP< i >
i
C 10101 I
.
,
1
IF (NPSP ,GT. 01
\ FALSE
>.--.-..
,1
I TRUE
1
/ PRINT 6010f I IHSP1 I) ,I«1 ,NPSP) /
I .................. , ...... . ........... ,.../
I
1 --- C 6010 FORMAT I////1 SPECIES TO
1 CBE PRINTED!' / 113/UO) I
I
o< --------------- ........ ---- .......
I
/ t t t . i i i t f i i i i < . i . 1 1 , < i t i i i < .x
IF (NOL -GT. 0) >-
V ......................................... /
I TRUE
I
/< i i > t i i i i i < . > i i . . i i > 1 1 < i i i > i t < i i > < i < t < t/
/ PRINT 401 1 i NOL /
1
I --- C 6011 FORMAT (// NUMBER OF
I LLtVELS TO BE PRINTED!'. 13)
I
0< .... ------ . --- ... ...... ..
I
i i . t i 1 1 . . i , i i . i i . t . . >x FALSE
IF JNCTC .GT. 0) >....-....
I TKUE
1
/ PKINT 6012i NCYC
1
I---C 6012 FORMAT (// CYCLE
I C INCREMENT FOR PRINTING:1, IS)
I
cn
?a
i
00
I
-------
< IF (TDELT »GT. 0.) >
\ /
I TRUE ,-
l
/>t.,...,,. .,.,,,,t, it ... ,,,,,,, it ,,,,,/
/ PRINT 6013, TDELT /
I
I---C 6013 FORMAT)//' TIME STEP FOR
i CHRINTING:', IPEio.t, SECONDS')
i
o< . .-..
i
< IF (CONC) > -
\ . t/
I TRUE
I
: AIM- 'YES* i
I
o< . .............
I
/i.i,,.i,.,,,., ,.. , ,, ,,<. , ,,,. , .,.,..,,,..v FALSE
< IF (PARTI >_-. .__.._...
\ '/
I TRUE
I
! A(2) » 'YEi' '
I
0< ... ................
I
/ PRINT 60M, A /
---C 60H FORMAT (//6X, 'DUMP', |OXi3
C'OUMP PARTICLES' / CONCENTRATI3
CNS', I 7X, 'AND "EIGHTS' / / 7X3
C A6. 13X> A6 ) 3
/"""""""""\
: RETURN :
\ /
en
?3
i
oo
-------
8FOR»FLO*FORFLO,S OUTPUT ,,.e,..t
FLOWCHARTED BY FORFLO /xeoa/ ON 30 SEP rj AT 17:52:37
o
K)
oooooi
000002*
000003.
000001.
000005.
000006*
000007*
000008.
000009*
000010*
C0001 1
000012*
000013.
000011.
oooois.
oooou.
000017.
ooooi a*
0000 1 9.
000020.
OOOQ21 .
000022.
C00023.
000021.
000025.
000026.
000027.
OC0028.
030027.
000030.
000031 .
000032*
0011033*
OOUU31.
00003S.
000036.
00(1037.
000038.
000037.
ooonio.
OOOOH1 .
000012.
000013.
OCOQ11.
000015.
OC0016.
000017.
000018.
C
c
C
c
10
20
30
HO
MB
6000
6010
SO
6050
60
100
110
C
,EQ. 0) NOPR .FALSE.
SUBROUTINE OUTPUTIM)
THIS ROUTINE PRINTS CELLULAR CONCENTRATIONS, IN PPH, *T REflUESTEO
INTERVALS AMD DOES ALL NECESSARY DUMPING OF DATA FOR USE IN OTHER
PROBLEMS, PLOTTING, ETC.
INCLUDE NEXINC
INCLUDE 10
LOGICAL NOPR
DATA SP /*co« /
N ' ABS(M)
IF (M «LT. -100) GO TO loO
IF (M *LT. 0) GO TO tB
oo 10 I > i, NX
DO 10 J « I, NY
DO 10 K » 1, NZ
CONVERT KG TO PPM
C
NP
A6, ', LEVEL*, 12, 10X, A6, 10X,
PARCELS* //)
OUT00100
OUT00200
OUT00300
OUTOOMOO
OUT00500
OUT00600
OUT00700
OUTOOBOO
OUT00900
OUT01000
OUTOI100
OUTOI200
OUT01300
OUTCI1HOO
OUT01500
OUT01600
OUTOI700
ouTOiaoo
OUT01900
OUT02000
OUT02100
OUT02200
OUT02300
OUT02900
OUT03000
OUT031QO
OUT03200
OUT03300
OUTC3100
OUT03500
OUT03600
OUT03700
OUT03800
OUT03900
OUT01000
OUT01IOO
OUT01200
OUTP1300
OUT01100
OUT01500
OUT01600
OUT01700
OJ
CO
73
i
OO
-pa.
-------
0000^91 IF (IS «LT. 0' RETURN 0
OCOOSO. 120 IF (.NOT. »GHT) GO TO 130
OOOOSl. FLD(2,UiIOLIST(10)) NP
000052. C *R1TE (28) T, NP , (WT(I)iI*ItNP)
000053. CALL SiTRANtZe. 2M, 2, IOLlST(9)i IS)
000054. IF (IS .LT. 0) RETURN 0
OOOOSS. 130 IF (.NOT. COI«C) GO TO ISO
OODOSA. C »RITE (2) Ti ICYC. C
000057. CALL S3TRAN(2, 1, NV*NSP « 2i T, IS)
OOOOSB. IF (IS .LT. 0) RETURN 0
DODOS?. ito RETURN
0000*0" END
OUTOH800
OUTOS900
OUT05000
OUTOSlOO
OUTOS200
OUT05300
OUTOSHOO
OUTOSSOO
OUTO&600
OUT05700
ouTOSaoo
OUTOSfOQ
O-l
OO
-pi
-------
OUTPUT
FLOYICHARTED BT FORFLO /XB08/ ON 30 SEP 71 AT 17:Sz:3?
\ SUBROUTINE OUTPUT-...-...-...-.-.--.---- .......................................Q
\ .......,/ Z
I FALSE Z
I Z
I
/ , I I . « , . t I . I . , . , « < I < » t t I . « I . , , , , , X TRO£
< IF (M «LT. 0) GO TO SB >..-......_._._.....-..-...._.._.-.-._.._-._._.._..-,....0
\ / X
I FALSE . X
1 X
I . X
' X
A....-......_.....> DO 10 i i, NX ; x
A : .: x '
x
B...............>; 00 10 J « I, NT : x
B : x
B I X
B I X
B X
a c--*------>: oo 10 K i, NZ : x
ABC :.... : x
ABC I -X
ABC IC CONVERT KG TO PPH 3 X
ABC 1 i
-------
o
Cn
ABC L IUJ I
i
i
I TRUE
I
: ICYC ICYC * i :
i
i
; NOPR .TRUE. :
1
1
I FALSE
1
I
< 0) NOPR > .FALSE. >
I TRUE
I
: oi NOPR . .FALSE. :
I
t
I FALSE
1
t
t PTJHE . PTIME » TDELT !
; NOPR .FALSE. :
I
C 203 t
I FALSE
I
1
I
J
s
s
s
s
s
s
s
5
S
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
y
U
u
u
u
u
u
y
u
u
A
x
A
X
X
X
X
X
x
X
x
X
X
X
x
X
X
X
X
X
X
x
X
X
X
x
X
X
x
X
X
X
x
X
X
X
X
x
X
X
X
x
X
X
X
&
I
z
7
2
z
z
z
z
I
I
I
z
I
z
z
I
I
z
z
z
I
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
7
OJ
rn
X>
i
oo
-pi.
-------
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A B
A B
A B
A B
A
A
A
A
1 FALSE
I
C 30] I
I
I
C 103 1
t CONTINUE :
I
C 18] 1
: IT » T » .s
: IT MODI |T/60,60) + 1 000 < J T/3600 » 3)
: ENCODE (6000tCLOCK) IT
: 6000 FORMAT (16)
! FLDI 18,6, CLOCK) - 13
I
I
I
I
/ PRINT 6010, SPtN), K, CLOCK, IC*C, NP /
I
I C 6010 FORMAT<»|OENSITY PLOT OF «3
I C A6, ', LEVEL', 12, IOX, A6, 10X3
1 1 1 'CYCLE*, 11, IOX. 16, 3
I CPARCELS' //) 3
I
I
I
C 503 I
1
I -C 6050 FORMAT (22F6.II 3
1
C 603 I
1
C 1003 I
v » i
(111
U X Z
u x z
U X 2
U X 2
U * 2
U X 2
U X 2
U X Z
u x z
x z
X 2
X 2
X 2
2
z
z
z
z
z
z
z
\ 2
z
z
z
z
z
z
z
z
z
z
z
2
z w
z£
_ po
* '
Z 00
.£»
Z 4^
z
z
z
2
z
-------
I FALSE
i
i
i FALSE
i
i
A r
A C 1 103 I
l-.-C WRITE 126) T, NP> ( X ( I ) , I 1 ,NP 1 3
1 C CYI I ) ,1-1, NPI ,
-------
o
oo
: RETURN o
\
C 1303
0<-
t
"«\ TRUE
>......
.../
rr (.NOT! COftCI GO TO |10
FALSE
t »HITE 12) T, ICYCi C OUTOSSOO 3
: : CALL SSTRANIZ, i, NV*NSP » 2, T, is) ; :
I
1
I TRUE
I
: RETURN o :
i
i
i
i
z
z
z
z
z
z
z
C 1103 I
/tltttt ,..!. I ,,.,t,S
RETURN :
Co
S3
i
00
-p.
-------
OFOn«FLO.FOKFLO,S pAKCEL
FLO.-/CHAKTEO BY FQRFLO /xaoe/ ON 21 SEP
AT u:is:o?
OonnOl
OCOOU2.
000003.
onooot.
ooooos.
000004.
000007.
oonons.
000007.
ooooio.
oauai i .
00 JO 1 2.
onooi i.
OO'JOM.
OO'JfjlS.
000016.
00'J017,
000016.
000017.
OC0020.
OT)02I .
000022.
000fj23.
GOU02M.
0300ZS.
000026.
000027.
000028.
000029.
OC0030.
000031 .
000032.
000033.
000931.
C0003S.
Q0003&.
030oi7.
030038.
000037.
C3fl010«
OCCI011 .
OCUO'12.
CC0013.
ooooi'f.
COOQH5.
oaooii.
000017,
000019.
GQJ017.
caooso.
GQOOSl.
bUBHOUTjNE PARCtLCHI
C THIS ROUTINE CALCULATES AN INTERPOLATED (DIFFUSION
C VELOCITY FOR EACH PARCELi AND THEN MOVES THE PARCEL
C POSITION. IF IT MOVtS OUT Of THE GRID, THE BOKDER
PAKOOIOO
AND AOVECTIONI PA800200
To ITS NLA PAK00300
CELL FLUXES PAKOOMOO
C ARF UPDATED. ARF.A-NE1SHTED CONCENTRATIONS ARE C^t-C^Lf.J ill AND »OOCD
C TO EXISTING CELLULAR CONCENTRATIONS.
INCLUDE NEMNC
DIMENSION VCL(3I
HOP NP
HUP * 0
L - I
1 CONTINUE
C CALCULATE THE FRACTION OF THE PARCEL IN EACH OF THE
1 X IL)
J Y(LI
K Z(L)
ll I * J
J| . J » 1
Kl . K » 1
Fxl » X(Ll - 1
FYl Y(L) - J
FZI « Z(L) - K
IT (K £« 01 K « J
FX«l»-Fxl
FY-I.-FYI
FZ-l.-FZt
Fl«FZ»FX«FY
F2«FZ*FXl«FY
Fj-FZ'fXI «Fri
F M . F i F X F Y I
F5«FZI»FX«FY
FA«FZ1 *FXI «Fr
F7=FZI «FXI «FYI
F8«FZ1«FX»FTI
C CALCULATE THE AVERAGE PARCEL VELOCITY FOR X-l, Y-2
20 00 30 M « It 2
30 V£L(M) UlK.j, I ,H)»FI * U(K, Jf I 1 ,M)«F2 « U(K,Jl,ll
1 U(K, JJ , I ,H)«F1 » U(K| , j, 1 ,M)»F5 * UIK1 , J,IJ ,M)»F»
2 U « VEL( 1 >
TIL) T(Ll » VELC2I
PAKOH800
PAKOt700
PAKOSOOO
PAHOSIOO
I
CO
-------
030052.
OOOIJS3.
00005-4.
030055.
000054.
C30059.
COU04Q.
C30Q6I .
000062.
CSC043.
C0004S.
030066,
C3i>067.
oooniB.
OODQ67.
COU070.
00007 I .
CQ0072.
COP073.
cono7H.
CCU07S.
000076.
C00077.
C00078.
003079.
030080.
GS0081 .
CO&U32.
000083.
000081.
001085.
on>io86.
CO JQ87.
G0008S.
00 >0d9.
OOi'o?C»
000091 .
OC'1092.
000093.
C0009t.
OC0095.
000094.
CC')097.
CC0098.
001)099.
00(1100.
oou I ni .
COUI02.
000)03.
OOUI01.
oomos.
00010*'
OOUJ07.
OOUIOS.
ZIL) ZIL) « VELI3)
IF (XIL) .LT. I.) GO TO 210
IF (XIL) .GT. FHX) GO TO 220
IF |Y|U| .LT. I.) GO TO 230
IF /lT(L>
C(K,J|,I,N) . C(KtJltliN) » F<**HT|L)
C(KliJiII.N) . C(Kl,J,ll,N) « F6*wTIL)
C(K| ,Jl ,1 I ,N) . ClICI ,JI ,11, N) * F7»*T|L)
C(KI,JI,I,N) . C(K1,JI,|,N| « Fe«WT(LI
100 NNP « NNP » I
XINNPI X(L>
YINNP) « Y(L)
ZINNP) ZIL)
GO TO 2Q1
C UPDATE FLUXES
210 IF (ZIL) .LC. 2.51 FLUY|(J) FLUYIIJ) « WTILI
GO TO 200
220 IF IZIL) .LE. 2.5) FLUY2IJI FLUY2IJ) » WTU)
GO TO 200
230 IF IZ(L) ,LE. 2.5) FLUX|(J> FLUXl(J) » »T(L)
GO TO 200
210 IF (ZIL) ,LE. 2.51 FLUX2IJ) « FLUX2IJ) » OTIL)
200 CONTINUE
201 L L * I
IF (L >LE. NP| G0 TO I
IF INP &£. MAX) RETURN Q
HP m NNP
RETURN
PAR05200
PAK05300
PftKOSfOO
PAROS500
PAK05700
PAN05SOO
PAHOS900
PAROiOOO
PASOilOO
PAN06200
PAH063DO
PAK06500
PAKU4400
PAK04700
P»KO(.UOO
PAK04900
PAKU7000
PAKO/ 100
PAH07200
PAH07300
PAK07400
PAH07SOO
PAK07600
PAK07/00
PAK07800
PORU7900
PAKOSIOO
PAW06200
PAh08300
PANUSIOO
PAKOB5CO
PAR00400
PAKOb700
PAKOHdOO
PAKOBVOO
PAHOVOOO
PAKOVIOO
PAKOV200
PAC09JOC
PAN09100
PAH09SOO
PAN09600
PAK09700
PAH09800
PAK09900
PARIOOOO
PAXIOICO
PAKI0200
PARI0300
PArtlOMGO
PAH10500
PAH10400
PANIQ700
PAK10800
CO
CO
-------
000109.
END
PARI0900
I
co
-------
PARCEL
FLOwCHARTED gl FO"FLO /XBOB/ ON 23 SEP 71 AT li:is:i|
\
SUBROUTINE PARCEL(N) /
---C THIS ROUTINE CALCULATES AN
C INTERPOLATED (DIFFUSION AND
CAOVECTION) PAR00200
C VELOCITY FOR EACH PARCEti AND
3
3
3
3
CTHEN HOVES THE PARCEL TO ITS N£»3
CPAR00300
C POSITION. IF IT MOVES OUT OF
3
3
ITHL GRID, THE BORDER CELL FLUXES3
CPAROOSOO
E ARE UPDATED. AREA-WEIGHTED
(.CONCENTRATIONS ARE CALCULATED
CAND AODLPARDCSGO
3
3
3
3
C To EXISTING CELLULAR CflNCENTRATj
' '
;
:
J
c t
:
CONS. PAROCA00
INCLUDL NEXINC : :
I
I
DIMENSION VELO1 :
NOP » NP
NM.P 0
U 1
1
3 I
CONTINUE :
---C CALCULATE THE FRACTION OF THE
3
3
CPARCEL IN FACH OF THE N£ I CiHBQR I N .)
:
{
t
;
:
:
!
C CEPAR013CO
1 X(L)
J YtL)
K Z(L)
11 1*1
J| J * 1
Kl K » 1
FXI A(L> - I
FT] V(L) - J
FZI « 1(1.1 - K
1
I
3
]
I
Z
Z
Z
I
Z
Z
Z
Z
Z
I
Z
Z
Z
Z
Z
Z
Z
CA)
oo
-------
t
o<
I
z
z
z
z
z
z
z
z
FX.I.-FXI Z
FY.I.-FY1 Z
FZ.t.-FZl Z
F|-FZ«FX«FY Z
F2»FZ»FXI«FY ' I
F3«rZ»F«l«FY| Z
F«i.FZ-FX»FYI Z
F5«FZl»FX»FY Z
F6«FZ1»FXI»FY I
F7-FZI»F*I»FYI Z
! F8.FZl«FX>FYI Z
: : z
i z
I C CALCULATE THt AVERAGE PARCEL 3 Z
I CVELOCITY FOR X>1, Y-2 ANo Z-3 3 Z
I CDIKECTIPAK03500 3 ' Z
I Z
I 203 I Z
Z
A-. .-- ..->: 00 30 M I . 2 i Z
A : : z
A I Z
A C 303 1 Z
A Z
................ VEL(M) u....-................. ...............................0 z
\ / X Z
I FALSE X Z
I X Z
I X Z
x z
! VFL(3I U(K,JiI,3)>F1 * UIK,J,I 1,31«F2 » 1 X Z
U(K.Jl.11>3)«F3 « UIK.JI.1IJI»FH » U( : X Z
Kl i Ji I |3>>F& » U(Kt iJi I 1 |3)>F6 » UlKli ' ' X Z
JliIJ,3)«F7 » UIK1rJIt1,3)*FB : X Z w
; ! I i %
I x z i
/.,,,,,........,..,,v , z ^
: GO TO so j..--.-....--.-...-...--.-._....-.-....-..-.«.-........-.-..-.-._...o x z r».
\. / . u X Z 4^>
U X Z
U X Z
Q<_..-......._[[[Q I
C MQ3 I U Z
u z
VELCJ) (Z(L) - .5)«.5* (UI2. J, 1 ,31 (F I* : U Z
FSI * l/(2, J, | I i3l«
-------
c
t
:
i
i c HOVE PARCEL PARcmoo :
i
503 I
III) * X(L| » VEL( 1 > i
TIL) r (Li » vtL(2i :
Z(LI " Z(LI + VELI3I
I
I
I FALSE
I
I
I FALSE
I
I
I FALSE
1
1
I FALSE
1
I
1 FALSE
I
i
I TRUE
I
Z(L) - .501 :
C ADO APPROPRIATE FRACTION OF
C#EIGHI TO CONCENTHATION ARRAY
CPAR05700
1
1 X(L)
J » TtL)
1C ZtL)
I 1 I » 1
Jl J « I
XI < » 1
N
N
N
N
N
N
N
N
N
N
ft
N
N
N
N
N
N
N
g
q
q
u
q
Q
V
a
a
u
Q
Q
a
a
a
e
Q
Q
Q
g
g
g
Q
s
s
s
s
5
S
s
s
5
S
5
S
s
s
s
s
s
s
s
s
s
s
y
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
7
z
z
z
z
7
z
z
z
z
z
z
z
z
2,
z
z
z
7
z
z
z
z
z
z
2
z
z
z
Z
i
2
Z
Z
2
Z
|
Z
z
z
z
z
z
z
-------
en
FTI » T
N
CIK,jli I,N) o ClMJliIl.N) * F3«wT(U N
C(K,J|>,M)»C(ICiJliIiNI»f'(*WT(L) N
C(Kt,Ji >N) C(KliJ)IiN) * F5««TtL| N
C«itJi 1,N| . c < K. I , J , I 1 iN ) * F6»WT(LI N
CIICI , Jl , I 1 ,N> C IKI , Jl , 1 1 ,N) » F7»»T(L) N
: N
ClKl.Jt.t.N) C(KliJliIiN) » FB»WT
z *
z
z
z
z
z
z
z
z
z
-------
C 2203
: FLUY2U)
C 2303
: FlUXMJ)
C 2*403
: FLuxzu)
t 2003
!
C 2013
!
1 . I
I
1
I TRUE
1
« FLUT2U) » WT
-------
IF (L (.£ NP> GO TO
I FALSE
I
I
IF (NP .GE. MAX*
I*. .*.... ...«!
I TKUE
RETURN o
TRUE
\ FALSE
NP
0<--«
I
NNP
t
I
RETURN
CO
oo
-p.
-pi.
-------
8FO.,«FLO.FOKFLO,S
FLOACHARTEO BT FQRFLO /xsoa/ ON 23 SEP
AT nsis:i8
00
OOODQI
G20Q02 .
OOOIJ03.
000005.
030106.
000108.
000009.
ooooio.
OOOol | .
OOOJ12.
OCOUI3.
OCOQM.
OOOdlS.
COU014.
coauir.
OCCOI8.
COODI 9.
CCJC20.
OCC02I .
OC0022.
C00023.
COOQ21.
000025.
CCD026.
CC0027.
CCOT28.
COU029.
CC0030.
C0003I .
OC0032.
CDU033.
000031.
C00035.
OOU03A.
000037.
C00038.
OUOOHO.
COUOHI .
COOOH2.
OQOQH3.
COUOHH.
CC001S.
CCOQH6.
copot'.
CCU018.
ccoosv.
ccuoso.
oocasi.
c
c
c
c
c
c
c
c
c
c
SUBROUTINE
THIS ROUTINE READS ALL
PROBLEM FROM UNIT 2?.
INCLUDE NEXINC
DIMENSION IOLISTI&)
FLDI 18,18, I3LISTI 1)
>18, IOLISM2I
i 18, lOLlST(31
,IB,IOL|ST(4)
i 18, lOLfSTIS)
.16,I3LI5U&)
16, !OLIST(5))
INFORMATION NECESSARY TO RESTART A
FLDI 13
FLDIla
FLDI 18
FLO<16
FLD(18
FLOI2,
10
LOC(X)
LOClYt
LOC1ZI
LOCUT)
LOCIT)
LOCIPMNSTR)
2
NSP»5
READ 1271 T, ICYC, PHNSTR, PMXSTR, *NSTOR, BCKSTR, NPSTOR
CALL S3TRAHI29, 25, 2, IOLIST(5), |Sfr
IF iis .LT. oi RETURN o
10 N » i, NSP
NPSTOR(N)
INP .EO. 0) GO To 10
S I « I , S
FLDI2,16,IOLIST(I))
READ 129) < x< i >,
IWTI I I ,l»l ,NPl
CALL S3tRAN|29, 25, 1, loLJST,
IF IIS .LT. 0) RETURN 0
CONTINUE
READ 1291 c
CALL S3TRANI29, 2, NV«NSP, C, IS)
IF (IS .LT. 0) RETURN 0
FLP12,U,|OLIST(A)I
uo
NP
IF
UO
,NP), IY( I) ,!! ,NP) , IZm.l-l ,NP) ,
IS)
CALL S3rRA.N<27, 10'
>18,tOLISTI1))
,18,IOLISTI2))
>18,IOLISTI3))
»l8,IOL|STCn )
1)8,10LISTI5))
RES00100
RES03200
RE..OG300
KL500MOO
HESOObOO
JENTRY *FILE
THIS ROUTINE DUMPS ALL INFORMATION NECESSARY TO RESTART A PROBLEM
ONTO UNIT 29.
LOCIX)
LOCIY)
LOCIZ)
LOCUT)
LOCIT)
FLD(18
FLDI IB
FLOI I8
FLDI 18
FLDI 18
FLOI ie
18,
FLDI2, I6i IOL1ST IS) I
FLD|2, U, IOLISTI6) )
TRITE 129) T, ICYC.
CALL S3TRANI29, 21,
LOCIPMNSTR)
2
NSP»S
PMNSTR, PMXSTR, WNSTOR,
2, IOLIST(5|, (SI
BCKSTRt NPSTOR
IF
00
NP
IF
00
I IS -LT. 0) RETURN 0
20 N « 1 , NSP
. NPSTOR(N)
IMP .EG« 0) CO TO 20
|5 t I, t
RES00700
KESOOdOO
HtSOO?00
RES01000
KESOl 100
HEiOlZOO
KLSOI300
RESUlbOO
HESOUOO
I /oo
HE301900
KE5U2000
KES02IOO
RES02JQO
RES02bOO
R£S02600
RES02700
RESO
-------
000052. C WRITE <2?) .J«J,NP>, (2,!!,NPI HESOS300
OCUOSH. CALL S3TRANI29, 2Mi M, IOLIST, IS) _ REbObHQO
CCOnSS. IF US .LT. 0) RETURN 0 ' KESObSOO
QQ0054. 20 CONTINUE KCSOb600
000057. C #RIIE 12?) C RESUS>;00
OOU053. CALL S3TRANC29, I. NV«NSP» Cl ISI «L!>Ob800
CC005». IF (IS «LT. 01 RETURN 0 KESOSVQO
ccooio. c END FILE 29 Rtsobooo
CC0041. c REWIND 29 ' KESo6ioo
OC0062. CALL S3TRAN(29t 9, 10) RES06200
CC0063. RETURN RESObioo
0000*1. ENO H£S06-(00
CO
CO
Ji.
-------
RESTRT
FLOaCHARTEO
FflRFLO /XflOS/ ON 23 SEP 71 AT 11:15:19
SUBROUTINE RESTRT
C THIS ROUTINE READS ALL 3
CINFORMATION NECESSARY To RLSTART3
CA RES00200 1
C PROBLEM FROM UNIT 27. KES00300 J
INCLUDE NEXINC
: DIMENSION IOLISTUI
: FLDC18,18.tOLISTf I I I LOC(XI
: FLOI|8,|8>IOL1STI2)' * LOCIT)
; FLOI18,18,lOLlSr(3)I . LOCIZ)
: FLOI18,I8|IOL|ST(H)) . LOC1«T)
! FLDt|3,18,IQLlSTIS)I LOCIT)
t FLD(I8,I 3,IOLISTI6)I LOCIP1NSTRI
: FLDI2,)/,, jOLISnS) ) « 2
: FLDI2,14,|OL1STI6)) « NSP«5
I :
i
i---c READ (27) r, icrct PHNSTR, 3
I CPMXSTR, WNSTOR, 8CKSTR, NPSTOH 3
I CRES01100 3
I
J
! : CALL S3TRANI29, 23. 2i IOL1STI5). IS) t :
: : : :
:.: .;.:
i
I
x.11i..t...i,,.,,,,,,.,..«,.,.,,.,,11 ,,,,tx FALSE
< IF US .LT. 0) >
\« ../
I TKUE
I
/« «t..\
: RETURN o :
\ /
o<-..............
i
>: oo jo N « I. NSP :
4---.-^^_--^-._^0_^-j
i
i
i HP NPSTORIN) >
I I
I
CO
oo
-fcx
-------
8
B
A
A
A
A
A
A
A
A
A
A
A
A
A
I FALSE.
1
1
I
t 53 1
I
i c READ <2vi ixu> ii«i .NP) , «T(i»i3
1 CI-liNP), IZ ( 1 ) , l-l ,NP) . RES022003
r c t c» isi : :
i
i
1 TRUE
I
: RETURN o s
i---c REWIND z? RESOSOOO
i
i
z
i
i
z
i
i
i
i
2
I
I
I
I
Z
z
Z
£
z
z
2
Z
z
z
z
z
rr\
?3
t
OS
-fa.
-P*
CALL &31RANI29, 101
-------
K>
N)
1
I
/,,,,,,,,,, I,, I.,,, ,N
: RETURN :
\ /
\ ENTRY WFILE /
\ /
---C THIS ROUTINE DUMPS ALL 3
CINFORMATION NECESSARY TO RESTART]
LA PROBLEM RESC3SOO ]
C ONTO UNIT 2?» RESOJ500 3
FLOt I 8 , 18 , IOUST I I > I LOCIXI '.
FLCM |8 , 18 , IOL1STI2)) LOCIY) :
FUD< 1 3, IB , 10L1STO) I LOc(Z) '
FLOdB.18, IOLIST(5)I « LOC(TI :
FLD(18.18,IOLISTI6)) LOC(PMNSTR) ;
FLD(2,14,|OL|ST(5)) 2 :
FLD(2,JA,|OLIST16» NSP'S :
I
I C i'RITE (29) T, ICYC, PMNSTR, 3
1 CPTlXSTRi ANSTORi BCKSTRt NPSTflR 3
I CRESOHMOO 3
I
1
: : CALL s3TRAN(2?i ZH, 2, IOLISTIS), Isi :':
: : : :
i
i
/,..,,.,,,,, >.,,,,,, .,>.t i»if t.,..,,,v
< IF (IS .LT. 0> '
N. .../
1 TRUE
1
/ \
: RETURN o :
\ /
o< ' Co
1 . X)
A....--..-..-._-,: oo 20 N . 11 HSP : °°
A :»««««tt*««««*««»«t«***«*«««**t«tt*««*«»»««t2 *t^'
A I
A I
A »...,.,,...,. .,>>,...>,.>
: NP * NPSTORiN) :
: ,...:
i
i
/.,«<>,,,,<,,. ,«,., f ..,.<<«. f ,.,, ,N Tf,u£
< IF (NP .E8« a' OJ fo ZO >- .-...-. . .» -. -.-. -..--. .-....-.....»........0
\ / I
I r ,i *c t
-------
CM
*
A
A
A
A B
A B
A
A
A
A
A
A
A
A
A
* r »t.-»i.
I
I
c is: i
i
i c WHITE (29) ixi 1 1 ,i«i ,NPI , (run
I Cl-l.NP), (Zl I ) ,1-1 ,NP) , KtbOSZOO]
1 C 1 (»T( 1 1 , I.| ,NP) RESOS300 J
I
I
: : CALL SSTRANU?, 2s, i, IOLIST, isi : :
i
i
I TRUE
I
: RETURN o :
C 20] I
t
( -C
-------
RETURN
CO
JO
I
oo
-pa.
-------
SFO*»FLO.FORFLO,S SLTUP/NEXUS
FLOJICHARTCO BY FOKFLO /XBOB/ ON 21 SEP i\ AT 11:15:27
00(1001 .
000002.
001)003.
OOJOOi.
000005.
C00006.
0011007.
000009.
OTUQ09.
CO. 1010.
OOIIOI 1 .
oonni2.
0000 1 3.
00001 H.
oauoi s>
OOflol*.
OOU017.
COOOIB.
OOtlOl'-
O00'j20.
oonozi .
CCIUI22.
C01IU23.
0001)25.
COCI126.
000027.
OOUC28.
000029.
000033.
OCC03I
COL'032.
OOOOJ3.
00003H.
000035.
OCO(i34.
CC(1037.
OOLU38.
000039.
CC0010*
COCO1*!.
COOOH2.
OOCO'O,
CCOOH-4.
ocoots.
0000*46.
COOOH7.
ccooia.
OCOOt9.
ocrnso.
OCOU5I t
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
SUBROUTINE SETUP
THIS ROUTINE READS INFORMATION NECESSARY TO START A PROBLEM
FROM UNIT 8.
INCLUDE NEXINC
READ RV
IF (» -LE. PMINI GO TO 20
SPREAD « AMONG NO PARCELS, EACH WITH WEIXiHT P«l
I
SLTOI600
SETOI700
SETOI600
StlOlVQO
SET02000
THE CONCENTRATIONS
SLTC2200
SET02300
SE102MOO
SE102500
SJ.T02600
SET02700
W/WN
fit W/NO SET02BOO
00 10 L I, NO SETOZ900
NP . NP « I SET03000
WT P" S£io3ioo
STARTING PARCEL POSITIONS ARE DERIVED "ITH THE AID OF URANO, A UNIFORMLY
DISTWIBUTEO RANDOM NUMBER GENERATOR (OVtR THE INTERVAL F«0« 0, TO 1.)
JUMP) I - .5, * URAIlDIOUMMY | Stl03<(00
' Y(NPI J - .6 + URANO(DUMM,Y) SEI03bOO
Z(NP) K - .5 + URAHD(DUMMY) SL103400
IF PARCELS FALL OUTSIDE THE OR ID LIMITS THEY A«E ThROHN A*AY. StT03;00
lF(ZlNP).GE.FNZ.OK.XINP).GE.FNX.OR.Y(NP).GE«FNY.OR,X|NP|.LE f.'P«NP-l SEI03900
10 CONTINUE SEIO'tOOO
20 CONTINUE SEIOHlOO
NPSTOR(N> * Np SE.TOH2CO
30 CONTINUE SETOtJOO
SETUP bORDER ARRAYS
LACH OOROER CELL is 6 TIMES THE NORMAL CELL SIZEI THREE TIMES AS
AND TnlCE AS HIGH.
INITIAL CONCENTRATION ASSUMES EXPONENTIAL DECAY TO BACKGROUND
LEVEL - 1/2 IN 8 CELLS
FUG-FAC
00 9 L*l
LONG
SE10M6CO
SETOS?OO
SETOH800
SLTOH900
SCTOSOOO
StTOSlOO
CO
CO
-pi.
-p..
-------
000052.
000053.
00005H.
coooss.
000054.
OODOS7.
COCOS8'
000059.
000060*
Do 1 I'UNX
BXHLtliNI-«C(liliI,N)-BACK>»RV»FU6»6.
7 BX2(L,iiN)"(CIt,NY,1,N)-OACKI»RV«FUG»6.
DO 8 J"Ii^Y
Btl(L,J,N)-(Cll.Ji«iN»-BACKI»RV»FU6»*.
6 BT2ll.iJiN)-IC(liJiNXiNI-BACK)»RV»6.»FUG
V FUG»FUG»FAC
RETURN
END
SETObZOO
SUOS300
SETOStOO
StTObSOO
SETCS600
5EIO&700
SETObSQO
SET05900
SETC6000
CO
CO
-------
SETUP
FLOWCHARTEO RY FORFLO /xeoe/ ON 23 SEP 71 AT 11:15:28
\ SUBROUTINE SETUP /
\ /
L THIS ROUTINE READS INFORMATION 3
[NECESSARY TO START A PROBLEM 3
CSET002CO . 3
C FROM UNIT 6. 5ETOOJOO 3
INCLUDE NEXINC : :
---c READ (8) PMNSTRI PHXSTR( UNSTOP
C BCKSTO SET00500 3
: : CALL SJTRANIB, 2* H»NSP, PHNSTR, is) : t
i
i
/.t it t ttt.ttt it t , ,,..i i t>t> tit tit t i it it . tx FALSE
< IF (IS .LT. 01
\.....
I TRUE
I
/'« it itt . ,i,N
\ ...../
I
: RN5ToR(i) 2000. :
*
I
1---C READ (B) C SETOO'OO 3
I
I
: : CALL S3T«AN<8, z, Nv«r;sp, c, is) : :
:.: i :.!
i
i
,t ttt t >tt t it. , . t tf i . i >f t. , , , it. i tit i t IIIIN FALSE
< IF (IS «LT. 01 >
\ ..../
I TRUE
I
/tttttttit tit ittttt\
: RETURN o :
cn
X)
t
00
-------
tvJ
co
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
B
8
a
B
B
a
a
B
B
B
e
a
D
a
a
B
a
a
B
a
8
8
B
3
B
a
a
9
B
8
B
B
a
8
a
B
a
8
a
B
B
8
B
B
B
8
B
a
B
C
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
r
0
0
D
D
D
0
D
D
0
0
0
0
0
0
D
D
0
0
0
0
0
D
0
0
0
D
0
0
0
0
0
0
0 E
D E
D E
0 E
a E
D E
D E
I
I
IIP « 0
UN . WNSTOR(U)
BACK « BCKSTRlNl
PMIN - PMNSTRlll)
PMAX PMXSTR(N)
I
I
* * * *
1
:
* ^
!
.
>
i
i -c STARTING WEIGHTS FOR EACH 3
I [PARCEL ARE DETERMINED 8r THE 3
I CCONCENTRATISET02IOO 3
I t ON THE DATA FILE AND THE 3
I CAVERAGE WEIGHT SET02200 3
I C " IS TOTAL *EIGHT OF POLLUTANT 3
I UN CELL SET02300 3
I
I
: x « (C(KiJ,iiN) - BACK>*RV
|
< IF (IT .LEi PMIN) GO TO 20
FALSE
C SPREAD * AMONG NO PARCELS. EACH3
CWITH WEIGHT PW SET02600 }
! NO W/«N * 1
t PA « tV/NO
I
I
I
1
t NP m WC » 1
J HT(NP> P*
:
'
:
z
z
z
z
z
z
Z CO
z ,
Z CO
Z 4^
Z -P>
z
z
z
z
z
z
2
Z
2
-------
H-1
tsj
a
a
9
a
a
a
a
a
A 8
A a
A a
A a
A a
A a
A a
A a
A a
A a
A 8
A a
A a
A 8
A a
A a
A a
A a
A 8
A a
A a
A a
A a
A a
A a
A a
A a
A a
A a
A a
A
A
A
A
A
A
A
A
A
A
n
o
o
D
o
o
D
0
D
D
0
D E
D E
0 E
0 E
0 E
0 E
D E
0 E
0 E
0 E
0 E
0 E
D E
0 E
0
0
0
0
0
0 E
C D E
C 0
C D
C D
C 0
C 0
C D
i c STARTING PARCEL POSITIONS ARE 3
CDEK1VEO WITH THE AID OF U«ANOi A3
CUN1SET03200 J
I DISTRIBUTED RANDOM NUMBEK 3
CGENERATOR (OVER THE INTERVAL 3
CFROH Ot TO bET03JOO 3
! X(NP| . I - .5 » URANO(OUMMY)
! Y(NP) » J - .5 + UHANOIDUMMY)
S Z(NP( K - .5 URANDIOUMMYI
>:
i
ic IF PARCELS FALL OUTSIDE THE
I CC,RID LIMITS THEY ARE THKQ*N
I C SET03700
I
I
/t>,.it.,.,t,i,,.,,.i>t.ti.,ii. .. , ,,,,N FALSE
< IFIZINPI .GE.FNZ.OR.xINPt .IJE.FNX«OR.Y(NP) . > - .»..
< GE«FMT.OR,X(NP).LE.It.OR.Y(NP),LE.I. I NP- >
\ /
I TRUE
I
C 103
I
0<-
I
CONTINUE
C 203
I
0<-
I
CONTINUE
NPSTOR(II)
NP
C 303
CONTINUE
I
--C SETUP 80R0EH ARRAYS SETOtSOO 3
C EACH BORDER CELL IS 6 TIMES THE3
CNORMAL CELL SIZE! THREE TIMES AS3
CLSETOSSOO 3
C AND TWICE AS HIGH. SETOH600 3
C INITIAL CONCENTRATION ASSUMES 3
[EXPONENTIAL OECAT TO BACKGROUND 3
ISET04700 3
C LEVEL - 1/2 IN 8 CELLS SETOM600}
en
I
00
FAC-.6
Fut«FAC
-------
A-.- -->: DO v L«i(H :
A : .<.............:
A I
A 1
A «...
A B"- --->: DO 7 j-i.NX :
A B : .....:
A B I
A a i
A B
A B I BXl(L,liN)-(C(I i1tI ,N>-BACKI»RV«FUG»6. I
A B : :
A B I
A B C 73 I
A B
A .... ... _ BX2(L,IiNI'ICtIiNYi1,N)-&ACK)»RV«FUG«6. :
A :,,..... .......:
A I
A I
A B-->: DO 8 j«iiNY :
A B :.. , , :
A B I
A B I
A B
A B ! BYlILiJiN)-(CtIiJi1iM-BACK)"RV'FUG'6. :
A B : ...:
A B I
A B C 83 I
A B
A ... ...._.. BY2IL,J,N)«(C(1,J,NX,N1-BACK)»RV»6.«FUG !
A : ....>.......:
A I
A C 93 I
A .
.................... FUG.FUG'FAC :
: ..t
i
i
/ I I I , I III I Illl I I.IIX
: RETURN :
\.. /
CD
50
i
oo
-------
BFOR.FLO.FORFLO.S SOURCE/NLXUS
BY FoRFLO /XB08/ ON 23 SEP 71 AT Ii:i5:38
OC0001 .
OC0002.
CC0003
CCOOOM.
GCOOGS.
ccuoo*.
COL007.
ccnouB.
OC0009.
ocoo io«
CCUO 1 1 t
CCG01 2.
0000 1 3*
OCOC 1 1«
OCDOI S.
OC'UO 1 6.
000017.
ccnc i 8
P r n n 1 9
i* U U C 1 T »
000020.
OC002 1
CCU022.
CC002 3 .
OC 002 M
OCU02S.
000026.
OCU077.
OC0028.
OCH029.
CCOD3C,
000031 .
000032.
OTOU33.
OOU03S.
Ocnn3<>.
030U3&.
OOU037.
OOOQ38.
000039.
orooMo,
00004 i
conot2.
OOOOH3.
CO ('OH 4 .
COUOH5.
coon1**'
CQ0017 *
000018.
000049.
OOJOSO.
OOOOSI .
C
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
SUBROUTINE SOuRCE(U)
THIS SUBROUTINE READS THE SOURCES FROM UNIT 8, AND CONVERTS THEM
INTO PARCELS, USING THL AVERAGE *EIGHT *N. PARCEL POSITIONS
ARE ASSIGNED USING THE RANDOM NUMBER GENERATOR, THE CALCULATED
'EIGHTS ARE AREA WEIGHTED AND ADDED TO CELLULAR CONCENTRATIONS.
INCLUDt NEXIN'C
FRAC IS USED TO DETERMINE THE VERTICAL FRACTION OF THE CELL THE
Nt». PAKCEL(S) HILL BE PUT 1NTO«
DATA F«AC /. 33333333 /
KEAD (Bl S
CALL S3TRAN(8, 2, NV, S, IS)
IF (is .LT. o RETURN o
DO 20 JJ - 2, NYMI
00 20 II 2 , NXMI
NO is THE NUMBER OF NE* PARCELS CREATED
NO S< 1 i JJi I I I/AN » .5
IF (NO .EQ< 0) GO TO 20
Pn » si l JJ,1 I )/NO
Do 10 L 1, NO
NP M NP » 1
UTINP) P«
X(HP) « II - .5 » URAND(OUMMY)
Y(IIP) JJ - .5 + URANOtOUMMY)
ZINP) > FRAC'URANDIOUMMY) * «S
ASSUMING EACH PARCEL HAS A VOLUME OF ONE CELL, CALCULATE FRACTION
PARCEL IN EACH OF THE NEIGHBORING CELLS, AND ADD THAT FRACTION OF
THE PARCEL WEIGHT IN EACH CELL TO THE CELLULAR CONCENTRATIONS.
I t. X(NP)
J Y(NP)
1 1 i + i
ji j » i
FXI X(NP) - I
FYI . Y(NP) - J
M » 1 . - FXI
Ft . L - FYI
cn, j, i, N> cUiJii,"! « PW»FX»FY
C(l.J.IliN) Ctl.Jitl'Nl » P««FX1»FY
C( 1 ,J1 , I 1 ,N| . Cd.Jl.ll.N) « P#»FXI»FY|
Cd.JliI.N) ClI.Jl.IiN) » PW'FX.FYI
10 CONTINUE
20 CONTINUE
SOURCES FOR BORDER
Y » I. EdGE OF GRID
IF THE WIND is BLOWING INTO THE GRID, AND THERE ARE POLLUTANTS
IN THE BOKOER CELLS, ME* PARCELS ARE CREATED AND PUT INTO THE 2
LOWErt VEHT1CAL CELLS. PARCEL *EIGHT IS ADDED TO THE EXISTING
CELLULArt CONCENTRATIONS BY AREA WEIGHTING.
. IF THE HEIGHT OF THE NEW PARCEL IS TOO SMALL, IT IS AODEO BACK
INTO THE BORDER CELL A«RAY.
DO 130 II-l.NX
souooioo
SOU00200
SOU00300
SOUOOHOO
souoosoo
SOUOC600
SOU00700
SOU00800
SOU00900
souoicoo
SOUO 1 1 00
SOU01200
SOU01300
souo i MOO
SOUOISOO
50U01600
SOUOI 700
SOUOI 800
SOU01900
SOUOiOOO
souo/ioo
SOU02200
S0002300
SOUO<:HOO
&UU02SOO
5000/600
SOU02700
SOU02830
SOU02900
SOU03000
SOUOJIOO
SUU03^00
SOU03 JOO
SOUU3100
SUUOJbCC
iOUU3600
SOU03700
SOU03800
SOU03VOO
SOUOMOOO
SOUOH loo
SOU01200
SOUOH300
SOUU1400
SOUptbOO
bOUUI&OO
SOUO-1700
SOUOHBOO
souoivoo
iOUOSOOO
SUUOb 1 00
O-l
CO
oo
Ji.
-------
NJ
000052.
C30053.
OOOU5H.
000056.
OOJOS7.
03UOS8.
030Q63.
030061 .
003062.
030063.
030061.
030065.
030066.
030067.
000068.
0330«9.
030070.
C33371.
COJU72.
000073.
C'JJOM.
033074.
030077.
030Q73.
C30J79.
03J390.
ooooai.
C30032.
030383.
030091.
030085.
C30US6.
030067.
03UU38.
000090.
C3039I
003092.
03.0093.
C00391.
003095.
000396.
C30397.
0300*8.
OOJQ99.
033100.
033101.
C30102.
0.33101.
033105.
030106.
IFIFLUX1(1 I).LE>0.) GO TO IIS
XFLX'FLUXl< til
FLUX1(1 I)>0.
IF (N'J «E(J. 0) 60 TO 111
00 110 L»liHO
NP«NP»1
VVT(NP)»P'H
XINPI.I 1-.S»URAND(DUMMY)
1 XINPI
IF(XINP).GE'FNX.OR.X(NP).LE.l.'GO TO 109
J » Y(NP)
K « Z(NP|
11-1*1
Jl » J « I
Kl K + 1
FXI X(NP) - I
FY1 Y(NPI - J
FZ1 ZINPJ - K
FX * I. - FXI
ft m |. - fYI
FZ 1 .
IF u
CIKiJ,1
CIK,J,I
C(K,Jli
CIK.J1 .
CIKIiJ,
cm i j.
FZl
01 K
1
N) CU.J.ItNl » PW»FX»FY«FZ
C(K,J,11,N> » Ptl»FXl«FY»FZ
* P«»FX«FYI«FZ'
« PiV»FX«FY«FZI
.N)
l.NI . CIK1,J.I1|N) * Prt»FXI«FY»FZl
C(K1,Jl , I 1,N) CIKI ,J1,I 1.N) + P*»FX1«Frl»FZl
C(KliJI,!,N) . CIKl.JlittNl + PK»FX»FYI»FZ1
GO TO 113
|0» NP«NP-1
110 CONTINUE
GO TO 1 IS
111 8X1(1,II,Nl « 9X1(1,11,N) * XFLX
C Y . NY EDGE OF GRID - PARCELS A«E TREATED IN A HANNEK SIMILAR TO
c THAT OF THE BOTTOM CELLS.
US 1FIFLUX2II I).LC.Q.>GO To |3o
XFLX'FLUX2(I I)
FLUX2I I I )>3.
IF (NO .EQ. 0) G0 TO |2»
P*«XFLX/NO
00 12Q L'liNO
000108.
I |-.5*UaAr|0(OUMMY)
Y(NPI«NY-URANO Z(NP)
SOUOS200
SOUOS300
SOUOSSOO
SOUObSOO
SOUOS6CO
S000b700
souobaoo
SOU06030
SOU06IOO
50006200
SOU063QO
SOU04MOO
SOU06SOO
SOU06630
SOU06/00
SOU06800
SOUU6900
SOUU7000
SOUU7 100
SOU0720Q
SOUU7300
SOU07MOO
SOUU7bOO
bOU07600
bOU07700
SOU07SQO
SOU07930
soooaooo
SOU08IOO
bOU08200
souoa joo
SOU08MOO
SOUOQSOO
SUUOB600
SUUOS700
SOuObaOO
SOUOB900
SOU09000
SOUU9IOO
SOU09200
SOUU93QO
SOU09100
SOU09500
SOU09600
SOU09700
SOUO'800
SOUOVVOO
souioooo
SOU10IOO
SOU10200
SOUI0300
souioioo
SOU10SOO
SOU10600
SOUIQ700
SOH10800
CO
I
co
-------
CM
OOJI09.
003110.
3031 I I.
000112.
OOJ1I 3.
0001II.
0001 IS.
003116.
00'31 I 7.
0031 IS.
0031 I*'
000120.
000121.
003122.
003|23.
000121.
000125.
000126.
000127.
00312?.
003129.
000130"
COOI3I.
C30I32.
00^133.
000131.
000135.
000136.
003137.
0001 3,1.
003137.
000110.
000111.
003112.
000111.
000115.
OC3I16.
000117.
000118.
000119.
003153.
000151.
000152.
000|53.
OO'USH.
00 HSS.
00 JISA.
OOI)lb7.
OCUISB.
0001S9.
003160.
003161.
030162.
000163.
OOJlol.
000165.
II
Jl
Kl
FX1
FY1
FZI
FX
FY
FZ
IF
.1*1
. J » 1
K + I
» X(NP) -
« Y(NP) -
« ZINPI -
. I. - FXI
1. - FY1
I. - FZI
(K >E4. 0)
CtK,J, I ,N) *
CIK,J, I I tN) .
C(K,Jl .I 1 ,N)
CIKI,JlI,NI :
119
K 1
C I < . J l I i N I
C«,J,I I ,N> *
CIK,Jl ,1 I |N)
C I K 1 J i I i N I *
C I K 1 . J I i I 1 i M )
CIKI.Jl,I|N) .
SO TO 120
NP.NP-1
Pw»FX«FY«FZ
* P*»FXl»FY«FZ
t PW»FX1«FY1»FZ
P««FX«FYI«FZ
PW»FX«FY»FZl
» PW«FXI»FY»FZI
C(KI,Jl,IIiNl * P*»FXI«FYl»FZI
ClKl.JltliN) * PH»FX»FYl»FZl
120 CONTINUE
GO TO 130
129 BX2(l,IIiN) > 8X2(1,[l.NI * XrLX
130 CONTINUE
00 23Q JJ«IiNY
x.i. EDGE OF GRID - PARCELS ARE TREATED IN A MANNER SIMILAR TO
THAT or THE BOTTOM CELLS.
IFIFLUfI(JJI.LE.OIGO TO 215
XFLX-FLUYI
FLUYt(JJ)»0<
IF ("0 .EJ. 0) GO TO 211
PW-XFLX/NO
00 210 L»liNO
NP«NP»l
«TtNPI»P»
Y(NP)«JJ-.5»URAHD(DUMMY)
Z(NPI».5»(URAND(OUMMY))*2«
1 X(NP|
J YINPI
IF«Y(NP|.GE.FNY.OR.Y(NP).LE.I.)GO TO 209
1 1 . I » I
Jl J * 1
itl K I
FXI XINPI - I
FYI « Y(NP) - J
FZI ZINP) - K
FX = I. - FXI
FY . 1. - FYI
FZ . I. - FZI
IF (K .£Q. 0) K . 1
CIKiJitiN) C'KiJ'tiN' * PN*FX*FY*FZ
ClK.J.Il.N) C(KfJiIliN) » PI'»FX|»FT«FZ
SOUI0900
soui 1000
SOU! I 100
SOUI UOO
SOUI 1300
SUUI I 100
bOUl IbOO
SOUI 1600
SOUI 1700
soui leoo
SOU| |VQO
SOUI 2000
SOU12100
SOU12200
SOU12300
SUU12SOO
SOU] 2600
SOU12700
SOU12800
SOUI 2900
SOUI 3000
SOU13IOO
SOU13200
SOUI JJOO
SOUI 3100
5UU13bOO
SOUI3600
SOUI3700
SOUI 3800
SOUI J900
souiiooo
SOU1HIOO
SOUI1200
SOUI 1300
SOUI 1100
soui iboo
SOUI1600
SOUI 1700
SOUISBOO
SOOI1900
SOUISOOO
SOUI 5 100
SOUIb^OO
SOUlb300
SOUIblOO
SbUlbbOO
bOUIbbOO
SOUI5700
souibeoo
SOU15900
SOU16000
SOUI6IOO
SOU16200
SOUI6300
SUUI6100
SOUI6SOO
CO
I
oo
-------
0001 66*
COO 167*
p n M i t* A
U U '' 1 O O
OOUl 69.
000 1 70"
000 171.
O0')|72.
00 U73.
COO 1 71
000175.
OCDI76.
OO'JI 77.
000178.
COIJ 1 79 .
OOJ 1 80«
ooo i n i .
00.1162.
0001 S3.
OOOI8H.
OC'J 1 85.
0001 8*.
00 J 187.
000 1B8.
OCOI89.
0 0 J 1 9 0 .
COT 191.
OOJ192.
0001 93«
00 0 1 9 H .
0001 95.
ooai 96.
OOJ197.
000191.
00.1199.
00'J200.
000201 .
000202.
03T203.
00 J201.
000 ?U 5 .
00') 20 6.
03IlfN> « C ( K 1 , J , 1 1 ,
CIM ,jl , I 1 ,N) CIK1 ,jl ,1
» Pw»FX»FYl»FZ
« pr«FX«FY»FZl
N) » P»»FXI«FY«FZ1
1,N) « P««FX1«FY1»FZI
C(K|(JI(I(Nt . ClKI(JllllN) « PH»FX«FYI»FZ1
liO TO 210
20' MP-NP-I
210 CONTINUE
GO TO 215
2it BYiiifJjiN) Brill, jj,N) * XFLX
c x . NX EDGI OF GRID - PARCELS ARE TREATED IN A MANNER SIHILA" TO
c THAT OF THF BOTTOM CE.LLS.
215 1FIFLUY2UJI .LE.O" )GO To
XFLX=FLUY2( JJ)
FLUY2I JJI-0.
NO-XFLX/. - FZl
IF (K >E<). 01 K 1
C CK , J, I ,N> C *
C 1 < , J , 1 i N ) C « , J , I 1 , N 1
CIK.J1 , 1 ,N) , C
-------
SOURCE
FLO'VCHARTEO
FORFLO /X808/ ON 23 SEP 71 AT I 1 : 1 5 : M I
SUBROUTINE SOURCElN)
C THIS SUBROUTINE READS THE 3
[SOl/RCES FROM UNIT 8, AND 3
CCONVERTS THEM souoozco 3
C INTO PARCELS, USING THE AYER-AGE3
C*EI(,HT tin. PARCEL POSITIONS 3
CSOU00300 3
C ARE ASSIGNED USING THE RANDOM 3
[NUMBER GENERATOR, THE CALCULATCD3
CSOUOOMOO 3
C ''(CIGHTS ARE AREA WEIGHTED AND 3
CAOCED TO CELLULAR CONCENTHATIUHSJ
c SOUQOSOO 3
INCLUDE NEXINC
C FRAC IS USED TO DETERMINE THE 3
CyERTJCAL FRACTION OF THE CELL 3
CTHE SOUQ0700 3
C NE« PARCELIS) WILL BE PUT INTO.3
CbOUOOBOO 3
DATA FRAC /.33333333 /
I
1 C READ (8) S SOU01000
I
I
CALL S3TRAN(8t 2i NV, S, IS)
: :
IF US .LT. 0>
\ FALSE
*
I TRUE
I
RETURN o :
Oo
DO 20 JJ
NTM|
-------
04
CT>
A
A
A
A
A
A
4
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
u
B
8
a
0
8
B
ft
B
B
3
B
B
8
g
B
B
a
R
B
B
ft
B
B
B
B
B
H
B
B
B
B
B
B
B
B
B
B
8
B
n
B
B
8
B
B
B
B
B
B
a
B
a
B
B
B
B
.» -
c
C
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
~Xt wvtu II- **"""»
I
I C NO IS THE NUMBER OF NEn PARCELS]
I CCREATEO SOU01500 3
I
I
NO » SI I iJJt 1 P/rtN * »5 i
I
I
I FALSE
1
I
: P« si i .JJ. i I I/NO :
i
I
I
i
NP NP « i :
WT(NP) P* :
XCNPI « II - .5 » URAtiOIOUHMY )
Y(NP) JJ - ,5 » URANOlOUMMY)
Z(NP) F"AC*U"»NO ( ounfu i » .5
-c ASSUMING EACH PARCEL HAS * 3
L VOLUME: OF ONE CELL, CALCULATE 3
CFRACTION SOU02500 3
c PARCEL IN EACH OF THE NEIGHBOR)]
CG CELLS) AND ADD THAT FRACTION 3
COF SOUQ2600 3
c THE PARCEL WEIGHT IN EACH CELL 3
CTO THE CELLULAR CoNCtNTHAT IoNS« 3
CSOU02700 3
j . YINP)
II I « I
Jl J * 1
FX1 « MNP) - 1
FYl Y(NP) - J
FX I . - FX]
FY . 1. - FYl
C(l.J.I.N) . C(I.J.I.N) » PM*pX*FY
C(I.J,I1,H) ClliJ.IltN) * PWFXl'FY
CIl.Jl.It.N) Cll.JI.II.M) * P*>FX|»FYl
C(l,Jl.I,H) " ClliJl.liNI » Ptl«FX«FY|
I
C 103 I
. . ......
^
Z
z
z
I
z
I
I
z
I
I
z
z
I
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
J
I
z
z
z
z
I
z
z
z
z
z
z
z
z
z
z
j
£
1
O-)
On
1
CO
-------
A B
A 8
A B
A B
203
I
0<--
I
CONTINUE
1
I C SOURCES K)R BORDER SOUO<42UO 3
I
1 C 1 1. FDGE OF GRID SOUOH400 3
I C IF THE WIND IS BLOMIflG INTO THE}
I CGRID, AND THERE ARE POLLUTANTS 3
I CSOUOHSCC 3
I L III THE BORDER CELLS. NE* 3
I CPAKCELS AKT CRFATED AND PUT INT03
I CTHL Z SOUt/lbCO 3
I C LO.iER VERTICAL CELLS. PARCEL 3
1 LttEIGHT IS ADDED TO THk EXISTING 3
I CSOUOH700 3
I C CELLULAR CONCENTRATIONS BY ARLA3
I (.WEIGHTING. SOUOHBCO 3
I I IF THE (.EIGHT OF THE NE« PARCEL3
i its TOO SHALL. IT is ADDED BACK 3
I CSOU01900 3
I C INTO THE BORDER CELL AKRAT. 3
I [SOU05000 3
I
I
A->: DO ijo ii«iiNx :
A : :
A I
A I
A /ttiiffitti..,..,.,....f............ ...ti.x TRUE
A < IF(FLUX1II II.LE.O.) GO TO IIS' >- -. ._-. .-...-.-.. ...................................
A \ / Z
A I FALSE Z
A I Z
A I Z
A Z
A XFLx-FLUXi(I I) : Z
A FLUXUlll-0. ! , Z
A No»xFLx/»(N*.s : z
A : : t
A I Z
A 1 Z
» /.,......... ...... it..................I.,«N TRUE i
A < IF INO »EQ. 0) 60 TO lit >..-.................-...........^....-..-...................0 j;
A X / Z w
A I FALSE Z ^
A I ' Z ja
A I Z i
A ............................................. Z OO
A : : z *"
A I Z
A I Z
A ' Z
A --------.>: oo no L-I.NO : z
A B Z
A B I Z
A B I Z
A B . Z
A B : up«up»i : z
A B : WT int- I-PII : z
- - _....._ »< . . ,
-------
OJ
CO
A B
« 6
A B
A 0
A B
A B
A 0
A B
A B
A R
A B
A B
A B
A B
A B
A B
A B
A B
A B
A 8
A B
A B
A B
A B
A B
A 8
A B
A B
A B
A 0
A B
A B
A B
A B
A S
A B
A B
A B
A D
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A A
A B
A B
A B
A B
A B
A B
A B
A B
A 8
A B
Y(t|P)»l .*UHANI){DUi1MY)..5
ZFZ
CtK.Jill.N) C(K tJ, I 1 iN> » P«I«FX1»FY»FZ
C(K,J1,1|,N) C(KtJlilltN) « PW«FXI«FYI»
FZ
CIK.Jl.ltN) C(K,J1,1.N) » PW»FX»FY|»FZ
C(Kl,J,I,M « CUI . J, 1 ,N) » P*»FX«>FY»FZ1
CIMiJiII.N) . CIKl.J.IliN) * Pn*FXI>FY*
FZl
CIM .Jl .11 >N) . CIKI .JI.IliM » PW»FXI»
FYl»FZI
CIK1.J1.I.N) CtKtiJUIiM * P««FX»FY1» :
FZl :
:
C 1095
. '
I
I
1
I..1..1.P-I
U
U
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
!> U
s u
s u
s
s
s
A
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
A
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X ,
X .
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X 2
< 1
t
z
z
z
I
I
I
I
I
I
I
t
I
I
I
t
I
I
I
I
I
I
I
I
I
I
I
I
I
I
I
I
I
I
I
I
I
I
I
I
I
I
I
I
t
I
zn
i
CO
-pa.
-------
OJ
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
»
B
B
B
B
1
C 1101 1
I
1
: GO TO MS
C 11-43 1
: 8X1 ( I , I I ,Nt BX1 ( 1 , I I ,N|
s
:
/
* XFLX :
1
I 1 1 " NY EDGE OF GRID - PARCELS 3
1 CfRL TREATED IN A MANNER SIMILAR 3
I CTO SOU09100 3
i c THAT OF THE BOTTOM CELLS* 3
1 CSOU07ZOO 3
B
B
B
8
8
8
B
B
S
B
B
C 1153 t
< IFtFLUx2( I I ) .LE«0« )GO T
1 FALSE
I
I
XFLX-FLUX2I I 1 )
FLUX2I I I )"0.
NO=XKLX/*N»«S
I
1
< IF (HO «EQ« 01 GO TO
I FALSE
1
I
: PW.XFLX/NO
1
1
1
I
:
t X(UP)»] l-.b»UKAND(DUMMr I :
: YlNP)"NY-URAf40
-------
A O
A B
A D
A B
A B
A B
A B
A B
A B
A B
A B
A B
A n
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A 0
A B
A B
A B
A 8
A 8
A B
A B
A B
A B
A B
A B
A B
A B
A B
A 8
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A R
A B
A B
A B
A B
A B
A B
I
I FALSE
I
I
j . Y(NP) :
K z
-------
M O
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A B
I
I
: NC-NP :
i
i
C 1273 I
; UX2ll,Il,N) BXZUlII.N) » XFLX :
I
----.-I CONTINUE !
I
I
---C X » I. EDGE OF GRID - PARCELS 3
CA*£ TREATED IN A MANNER SIMILAR 3
CTO SOU13500 3
c THAT OF T.HE BOTTOM CELLSI 3
CSOUI3600 3
I FALSE
I
1
: XFLX'FLUYI ( JJ) ;
I FLUY 1 < JJ)«0. I
: NO«XFLX/WN»i5 I
1
1
1 FALSE
I
1
: PH.XFLX/NO :
I
1
X
X
X
X
X
X
X
X
X
X
X
X
X
X
1
z
z
z
z
z
z
I
z
z
z
z
z
z
z
z
z
z
z
L
1
Z
Z
z
2
Z
z
Crt
CO
-------
K)
o
a
a
B
a
a
B
B
A B
A 8
A B
A a
A B
A B
A B
A B
A B
A B
A D
A D
A B
A P
A B
A B
A B
A B
A B
A B
A B
A B
A 8
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A 8
A B
A B
A B
A B
A B
A 6
A 8
A B
A B
A B
8
B
B
B
B
B
8
B
B
B
1
NP«NP»1
WT HJI'I'PVt
XO|P)-l«*URAUD(DUMMY)«tS
Y ( HP )«JJ-.b»URANOI DUMMY)
Z(I4P)>.5*(URAND(DUMMY) 1*2.
1 « X » C(KiJl,liN) * P»»FX«FY|»FZ
C(KI.JttiN) C(KliJtliN) » P»«FX.FY.FZI
C(Kl,J,lltN| . C(Kl.JilltH) » P(,«FM.FT»
FZI
CIKI.Jl.II.Nl CtKl , Jl , 1 1 iH) * P««FX1«
FYl'FZl
C(K|,JI,I,N| » c(KliJliliN) * P*«FX«FYI«
FZI
I
1
U
U
U
U
U
U
U
U
U
U
U
U
U
U
U
U
u
U
U
U
U
U
U
U
U
U
U
U
U
0
U
U
U
U
U
U
U
U
U
S U
C It
X
X
X
X
X
X
X
X
x
*
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
t
z
2
z
z
z
z
z
z
z
7
z
z
2
£
z
z
z
I
z
z
z
z
z
z
z
z
z
z
I
I
z
z
2
z
z
z
I
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
2
z
z
I
I
7
t/J
CO
?3
i
(30
-P>
4^
-------
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
8
B
B
B
B
B
B
B
B
B
C 2093 I
: kp.NP-l
I
C 2103 1
1
I
: GO TO
C 2IM3 I
: BT j ( i ,JJ,N) m an ( i ,
i
i
i
i
i
i
i
C 2153 1
< IFIFLUY2IJJ),LE.O.
1
I
I
S U
s
: s
s
:
JJ,N| * XFLX :
-C X » NX EDGE OF GRID - PARCELS 3
CAKE TREATED IN A MANNER SIMILAR 3
CTO SOU 17600 3
c THAT OF THE BOTTOM CELLS* 3
CSOUI7700 3
FALSE
XFLX-FLUY21JJ1 :
FLUT2(JJ).0. I
NO«XFLX/nN».S :
a
B
B
B
B
B
1
1
< IF (NO >EQ. G) G
I
1
1
FALSE
_
: PW«XFLX/NO :
i
i
i
i
: NP«NP»I
NO :
w :
X Z
X I
x z
A Z
x Z
x z
x z
x z
X Z
x z
x z
x z
x z
« z
X Z
x z
x z
x z
D Z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
3 Z
( Z
I Z
1 w
z CO
1 JO
Z I
z oo
z 4a.
z -^
z
z
z
z
z
z
z
z
-------
A B
A B
A 8
A B
A B
A 8
A 8
A 8
A B
A B
A B
A B
A B
A B
A 8
A B
A B
A B
A B
A B
A 8
A B
A B
A B
A B
A B
A 8
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A B
A 8
A B
A B
A B
A B
A 3
A B
A B
A B
A B
A B
A B
: Y(MP >.JJ-,5»UKANO(OUMMY 1
: Z IMP > o,5* f UfiANU (DUMMY ) ) .2.
: i . x ( N p )
: j YI NP>
1
I FALSE
1
1
K > Z(NP1 <
11-1*1
Jl J * 1
K| . K » 1
FXI « X(MP) - 1
FYl . YINP) - J
FZI ZINH) - K
FX 1. - FX|
FY . 1 . - FYl
FZ 1 . - FZ|
I
I
1 TRUE
1
: K - i :
* * *
i
i
C(K,J.I,N1 C(K,JiI,NI » P««FX.FY«FZ !
C(K,J,I1,N) . ClKiJ.HtN) * PW»FXl»FY«FZ
FZ
C(K|J1>I.N) » C(K.Jl.IiN) « P«»Fx»FYl»FZ
C(Kt,J.l,NI « C(KliJ,t|NI * PW*FX«FY*FZ1
CIKI.J.Il.N) C(Kl.J.IllN) » PVT»FXl»rY»
FZI
C(K1 , Jl , I 1 ,N1 C(K1 i Jl . I 1 iM * PW*FXI*
FYl»pZJ
CIKI.JI.I.NI « CIKI.JI.I.N) » p«i»Fx»fri» :
FZI t
1
I
C 2193 I
u 0 m U l> - 1
u
u
U
U
U
U
u
u
y
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
u
s u
s u
s u
s
s
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
X
I
z
I
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
7
en
oo
-------
A B
A e
A B
A B
A B
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
J
0< .... .......
C 220] I
...... ,; CONTINUE 5
1
I
C 229] I
: BUIIiJJ.NI BT2<1,JJ,N> » XfLX '
I 230] 1
I
I
s x z
s z
....... ........o I
z
z
z
z
z
z
x z
X Z
X Z
x Z
....... . . z
z
z
z
z
z
RETURN
CO
x>
I
CO
-pi
-pi
-------
PICFIC SAMPLE OUTPUT
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
a
0
a
0
0
0
0
0 '
0
0
0
0
a
u
0
0
0
~ '0
0
0
0
- 0 "
0
0
0
0 ""
0
a
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
D
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
TIME'
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
... Q . .
0
0
0
0
0
- 0
0
0
0
0
0
o o o o o c
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
3 O O O O O (
0
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
0
0
3 O O O O O <
o
a
0
0
0
0
o~
0
0
0
0
0
O OOO 0 t
> o o o o o e
o
0
0
0
0
0
"0
0
0
0
.001
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
> o o o o o c
0
0
0
0
0
3 O 0 O
}
0
0
0
0
0
0
0
0
0
0
0
0
0
0
* 0
0
0
0
0
0
" 0
0
0
0
0
0
o o o o o c
a
0
0
0
0
0
O 0 O 0 0 C
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0"
0
0
0
0
0
o
0*
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- - o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1 O O O O O C
> o o o o o c
u
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
3 O O O O O 1
u
0
0
0
0
0
0
0
0
0
0
0
- o -
0
0
0
0
0
o o o o o c
u
0
0
0
0
0
0
0
0
0
0
0
0
3 O 0 0 O O C
a
0
0
0
0
0
0
0
0
0
0
> o o o o o c
> o o o o o c
0
0
0
0
0
> O 0 O O O C
0
0
0
0
0
0
0
> o o o o o c
0
0
0
0
0
3 00 O O 0 <
a
0
0
0
0
> 0 0 00 0 C
0
0
0
0
0
0
0
0
0
0
OO O C
0
0
0
0
0
0
0
0
0
> oo o o o c
0
0
0
0
0
3 O O O O O (
3 O O O O O <
> o o o o o c
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o o o o o c
3 OO O O O <
o
0
0
0
0
0
3 O 0 O 0 O 1
0
0
0
0
0
0
0
0
0
0
oo o c
0
0
0
0
0
0
0
0
Q
0
0
0
0
0
3 OO O O O 1
w
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
3 O O O 0 O <
3000
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o-
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
> o o oo o c
3 o o o o o <
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
. 0
- o
0
0
0
0
0
- 0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o/
0
0
0
0
0
0
0
0
0
i o o o o o e
0
0
a
0
0
0
0
0
0
0
0
0
0
0
0
3 O 0 O 0 O C
0
0
0
0
0
0
-s
0
0
3 O O O
0
0
0
0
3 O O O O O <
0
0
0
0
0
0
0
0
a
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
o -
0
0
0
0
0
0 -
0
0
0
0
0
3 O O O O O (
0
0
0
0
0
0
0
0
0
a
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
) O O O O O C
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
> o o o o o c
0
0
0
0
0
0
0 '
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
> o o o o o c
u
0
0
0
0
0
3 O O O O O 1
0
0
0
0
0
0
0
0
0
0
0
0
0
en
JO
I
00
-------
- - - 0
0
0
0
0
0
0
~ Q
0
0
0
0
0
0
0
o
0
o
0
0
0
0
- o
0
0
- 0
0
0
0
0
0
0
_ g
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
TIME-
0 0
0 0
0 0
0 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
Q
0
0
0
0
0
0
0
0
0
0
o --
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
250.
0
0
0
0
o- -
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 "
0
0
0
0
0
0
0
0
0
0
0
0
0
0
000
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0 -
0
0
0
0
0
. o---
0
0
0
0
0
0
0
0
0
0
o ---
0
0
0
0
0
o - -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
- o-
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
5
s
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
oooooooooooooooc
0
0
0
0
0
0
si
26
V
2
0
0
0
0
0
0
0
0
0
-0
0
0
0
0
0
0
0
0
0
0
0
1
oooooooooooooooc
0
0
0
0
0
0
1
is
HO
61
65
2S
5
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- - o
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
J
IB
57
11*
7?
10
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0-
0
3
15
12
76
56
31
5
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-- o-
0
0
0
0
0
- o -
0
0
0
0
0
0
0
- 0
0
0
t
IS
2V
17
V
i
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
H
5
1
0
0
0
0
o -
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
- o -
0
0
0
0
0
0
0
0
0
0
0
1
- 1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 -
0
0
0
0
0
o
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-~0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
- -fl
0
0
0
0
0
0
0
0
0
0
Q
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
o-
0
0
0
0
0
- 0
0
0
0
0
0
--0
0
0
0
0
0
0 -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o - -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 - - - -CM. _
0 00
o '
o ...
0 .£,
-------
CO
TIME*
0
0
0
0
0'
0
- - - o
0
- - o
0
o-
0
o
0
0
0
" 0
0
0
0
- 0
0
- 0
0
0
0
o
0
0~
0
- - - 0
0
-- 0
0
0
0
0
0
- - - o
0
0-
0
o
0
- - o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-- -0 -
0
0
0
0
0
-o-
0
0
0
0
0
- 0 -
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0 -
0
0
0
0
0
0
0
0
0
0
0
0
0
--0
0
0
0
0
0
~ 0
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
0
0
o-
0
0
0
0
0
' 0
0
0
0
0
0
0
0
0
0
0
0
o-
0
0
0
0
0
0
0
0-
0
0
0
0
0
0-
0
0
0
0
0
o -
0
0
0
- o
0
- 0
0
0
0
0
0
0-
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0-
0
0
0
5
-------
TIME"
0
0
- - o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
... 0
0
0
- - o
0
. _ . 0
0
0
0
0
-- - -0
0
0
o
0
0
0
0
0
o
0
o
0
- 0 -.
0
--- - o
0
0 -
0
0
0
0
0
0
0
0
0
0
0
0
0
o --
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
o -
0
o -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
a
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
750«000
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- -o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 --
0
0
0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 --0
0 0
0 0
a o
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
O 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0
0
0
0
o
0
0
0
0
0 ~
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
2
2
t
1 .
1
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
a
0
0
0
0
0
0
0
0
3
H
10
e
5
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-- 0
0
0
0
0
a
--0
0
0
0
0
0
0
0
0
0
0
0
0
J
»
20
3H
25
- 15
3
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
0
5
17
35
57
ts
- 27
10
3
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
0
1
21
M7
77
57
-28
I1*
3
1
0
0
0
0
0
0
9
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-- 0
0
0
0
0
0
0
0
0
0
0
0
0
7
17
36
57
HS
0
0
0
0
0
0
0
0
0
3
6
18
2'
21
-3015
10
3
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
_. g -.
0
0
0
0
0
0
0
0
0
H
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
0
0
3
6
10
6
__T_ _ \
\
0
a
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-- 0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
2
2
3
1
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
. o -
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0-
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
- -0"
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0 -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 -
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0'
0
0
0 -
0
0
0
0
0
0
0
0
a
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
. 0 -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
« _.^^. , _..T
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o - - -
0
o -
o
0
0
-------
on
O
TINE* lQOO'000
- o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
~ ~ 0
0
0
0
D
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
a
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
a
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o --
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0-
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
--o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
1
t
1
I
I
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
2
5
10
a
1
I
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
2
6
IS
30
20
10
2
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
}
13
22
15
37
18
7
2
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
1
t
IS
26
60
38
20
?
2
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
s
12
2S
48
32
21
6
2
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
2
(,
11
27
J7
11
2
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 -
0
0
0
2
s
11
s
3
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0-
0
0
0
1
3
3
i
1
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
J
1
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0 -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 -
0
0
0
0
0
0
0
0
Q
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0.
0
0
0
0
0
0
0
0
0
0
0
0
0
Q
0
0
0
0
0
0
0
0
0
0
0
n
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
Q
0
0
0
0
0
0
0
0
0
0
0
0
0.
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 - - .
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
Q __ ___,
0
0
0
o
0
0
0
0
-- co
to
00
-------
TIME- !2SOtOOO
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
Q
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
" 0
0
0
Q
0
0
0
0
0
0
0
- -o
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
0
0
-- o
0
0
0
0
0
0
0
0
0
0
0
0
o- -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0 --
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
Q
0
0
0
0
Q
0 -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-- o
0
0
0
0
0
-o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
o
0
0
0
0
0
- 0
0
0
0
0
0
1
t
0
i
1
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
Q
0
0
- o
Q
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
3
5
1
2
0
0
0
0
0
0
0
0
0
0
-- o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
_ -o
0
0
0
0
0
0
0
0
0
0
0
s
e
12
10
9
2
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
2
11
IB
31
23
IS
2
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
- 0
0
0
0
0
2
17
32
S3
38
24
5
1
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
J
20
38
56
It
25
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
1
2
17
31
HH
36
26
s
2
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
- -o
0
0
0
0
0
- 0
0
0
0
0
1
|2
IB
27
21
15
i
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
0
-- 0 -
0
0
0
0
0
s
e
13
9
6
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
i
i
-------
en
TIME.
0
0
0
0
0
- - 0
0
0
0
u
0
0
0
- o
0
I)
0
- o
0
0
0
o-
0
0
Q
0
0
- o
0
0
0
0
0
0
- o
0
0
0
0 .
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
Q
0
--o
0
0
- o
0
0
0
0
0
0
D
0
0
0
0
0
0
- 0
0
0
0
0'
0
0
0
0
0
0
0
Q
0
0
0
0
0
0
Q
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
1500.000
0
0
0
0
0
0
0
0
0
" 0"
0
0
0
0
0
0
0
0
0
0
0
--o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
" " 0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
"" 0
0
0
0
0
0
" "0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
D
0
0
0
0
--- 0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
n
0
--o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
~"0 -
0
0
0
1
2
1 -
2
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
Q
0
u
0
0
1
1
6
s
t
2
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
3
8
II
II
10
a
I
a
0
0
0
0
0
' 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
1
3
V
|3
25
- 26
|V
10
2
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
1
1
11
21
38
31-
33
IB
3
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
2
5
15
25
15
- 1H
3H
U
1
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-- 0
0
0
0
0
0
0
0
0
0
0
0
1
1
A
12
24
3S
38
28
21
1
1
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
0
. 0
0
0
" 0
0
0
0
0
0
0
1
3
?
IS
21
- 25
1*
11
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-- o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
1
1
6
12
-13
7
6
1
0
0
0
0
0
0
0
0
0
0
0
-- 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
1
3
- 1
3
2
1
0
0
0
0
0
0
0
0
0
0
0
0
--- 0
0
0
0
0
0
- o
0
0
0
0
0
0
0
b
0
0
0
0
0
0
0
0
0
0
0
0
0
1
1
2
1
2
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o-
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o-
0
0
0
0
0
o-
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
- 0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
_. Q -
0
0
0
0
0
-- 0 -
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
- o
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0'
0
0
0
0
0
0
0
0
0
0
0
0
a
0
0
0
0
0
0
0
0
0
0
0
0
0
a
0
-0
0
0
0
0
0
0
0
0
Q
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
- o
0
0
0
0
0
0- -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
o - -
0
0
0
0
0
0
0
0
0
0
0
0 ...
0
0
0
0
0
0
0
0
0
0 - - -
0
o - -
0
0 '
0
0
0
CO
?3
I
oo
-------
en
TIME.
0
0
- - o
0
U
0
0
0
0
0
0
0
0
0
0
- Q
0
" " 0
0
0
0
0
0
0
0
"- ~ 0
0
0
0
0
0
0
0
0
- - o
c
0
D
0
0
" 0
0
0
0
. 0 "'
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
o--
0
0
0
0
0
0
0
0
0
0
'0 "
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
Q
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1750.000
0
0
0
0
U
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
- - o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
o--
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
U
0
0
0
0
0
0
0
0
0
0
0
0
U
0
0
0
0
0
U
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
Q
0
0
0
0
0
o
0
0
0
0
0
1
2
2
3
1
0
_- , -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
U
1
1
6
s
3
I
0
a
c
0
0
0
0
0
0
0
0
U
a
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
Q
0
0
0
o
0
0
0
1
2
8
?
13
It
9
s
1
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
6
I
IS
27
22
!»
7
- 3
1
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
- o
0
0
0
0
0
0
0
-- 0
0
0
0
0
0
0
0
0
0
0
--o
0
a
1
10
2?
33
31
23
It
- 6
2
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
o
0
0
2
i!
32
13
35
23
>*»
e
2
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-- o
0
0
0
0
0
0
0
o--
0
0
0
0
0
0
0
0
0
0
-- o
0
0
1
?
27
35
30
23
IS
c
3
1
0
0
0
- - 0 ~
0
0
0
0
0
o -
0
0
0
0
0
- o~~
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-o
0
0
0
6
1 «
17
23
\1
|3
7
3
I
0
0
0
0
0
0
0
0
0
U
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
1
3
10
13
11
8
s
1
0
0
0
0
0
0
0
0
0
o - -
0
0
0
0
0
- o -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-o
0
0
0
0
1
3
5
1
3
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
2
2
1
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
U
0
0
0
n
0
0
o
0
0
0
0
0
U
0
0
0
0
0
U
0
0
0
0
0
0
0
0
0
0
0
0
0
0 -
0
0
0
0
0
0
0
0
0
0
n - ,..--
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
_ Q . -
0
0 *
a
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
U
0
a
0
n
0
0
0
0
0
0
0
0
0
0
0
0
0
0
U
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0-
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- 0 -
0
0
0
0
0
0
0
0
0
0
-o -
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
U
0
0
0
0
0
0
0
0
0
0
0
0
0
0
a
0
0
0
0
0
0
0
0
0
0
a
0
0
0
0
- 0 --
0
0
0
0
0
0
0
0
0
0
- 0 -
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0 -
0
0
. 0
- o
0
- o
0
0
0
CO
OO
-(i.
-Pi
-------
tn
0
0
0
0
0
0
0
0
0
0
o
0
- - 0
0
0
0
0
0
0
0
0
0
0
0
0
~ " 0
0
0
0
- o
0
0
0
u
0
0
0
0
0
u
0
. - 0
0
0
0
0
0
o
0
0
0
0
0
- "0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
TIME
0
0
0
0
" 0
0
0
0
0
0
o
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
~ 0 '
0
0
0
0
0
0
0
o -
0
0
0
~ " 0
0
0
0
0
0
0
0
u
0
0
0
0
0
- o
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
2000
o-
0
0
0
0"
0
0
0
o
0
0
0
0
0
0
0
o-
0
0
0
0
0
0
0
0
0
0
- 0"
0
0
0
0
0
0
0
0
0
0
0
0
0
" 0
0
0
0
a
0
0
0
00
0
0
0
0
0
0
0
0
0
0
--o
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
0
0
c
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
- o '
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
-o
0
0
0
0
0
.0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
1
I
3
2
'" 3
1
1
0
0
0
o
0
0
0
0
0
'0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
0
0
0
0-
0
0
0
0
0
1
2
2
6
5
3
0
0
0
a
0
0
0
0
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
1
5
10
|2
12
8
3
2
0
0
" ~ 0
0
0
0
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
0
0
0
--o
0
0
0
0
0
2
t
10
12
23
23
18
10
s
2
1
0
" 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
t
13
21
30
27
17
10
1
2
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
ff
0
0
0
0
0
0
0
- o
0
0
0
0
0
o ~
0
0
0
0
1
2
a
It
23
36
35
16
to
5
1
0
0 "
0
0
0
0
0
0
0
0
0
0
0
-- o -
0
0
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
7
12
2)
2?
32
I?
?
s
2
1
u
0
0
0
0
0
a
0
0
0
0
0
0
0
0
0
0
0
0
0
" ~0
0
0
0
0
0
- o
0
0
0
0
0
0
0
Q
0
0
0
5
9
n
20
22
15
10
6
i
0
0
0
0
0
0
0
0
o
0
0
0
0
0
~0
0
0
0
0
0
0
0
- o "
0
0
0
0
0
- 0"
0
0
0
0
0
o
0
0
0
0
0
2
s
8
11
13
7
3
2
1
0
o
0
0
0
0
0
0 " "
0
0
0
0
0
- -o -
0
0
0
0
0
0
0
0
0
a
0
0
0
o -
0
0
0
0
0
o --
0
0
0
0
0
1
2
3
s
6
3
1
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
o -
0
0
0
0
0
0
0
o
0.
0
0
0
0
o
0
0
0
0
0
o -
0
0
0
0
0
0
1
2
3
2
1
1
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
a
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
0
0
o
0
0
0
0
0
- o
0
0
0
0
0
o -
0
0
0
0
a -
0
0
0
0
0
0
0
0
0
0
- - o
0
0
0
0
0
- o
0
0
0
0
0
'" 0
0
- o
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
a
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
0
0
0
-- o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
- "0
0
- o
0
0
0
0
0
0
a
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
a
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0 .
0
0
0
0
-Q--
0
0
0
0
0
0
0
0
0
a
0
0
Q
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
- o -
0
0
0
0
0
0
0
0
a
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
Q
0
o
0
0
0
0
0
o - - - .-.._--
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
Q
0 " "
0
o -
0 IA!
o on
o to
o - - - '
0 2
o - £
0
-------
Cn
On
TIME"
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
0 "
0
" 0 '
0
0
0
u
0
0
0
0
0
D
0
' 0
0
0
0
Q
0
- o
0
0
0
0
0
0
0
0
0
0
0
u
Q
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
a
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0"
0
0
0
0
0
" 0'
0
0
0
0
0
- o-
0
0
0
0
a
0
0
0
0
0
0
" 0
Q
0
0
0
0
0 '
0
0
0
0
0
"0 "
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
0
a
0
0
Q
0
0
0
0
0
0
0
0
0
a
'0
0
0
0
0
0
o
0
a
0
2250-
0
0
0
0
" 0
0
0
0
0
0
- 0" "
Q
0
0
0
0
o
0
0
0
0
0
-----
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
' 0 "
0
0
0
0
0
0 "
0
0
0
000
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o -
0
0
0
0
0
" 0"
0
0
0
0
0
0 "
0
0
0
0
0
0~
0
0
0
0
0
0"
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0 "
0
0
0
0
0
o
0
0
0
0
0
0 '
0
1
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
1
3
3
3
2
I
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
a
0
0
1
2
3
6
6
H
2
0
a
0
0
0
0
0
0
Q
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
a
0
0
0
0
0
0
0
0
0
0
o
0
0
1
3
1
"10
IS
11
9
7
3
1
0
0
0
0
0
' 0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
' 0
0
0
0
0
0
0
0
0
0
0
0
-- o
0
1
3
7
9
18
22
19
15
a
1
i
0
0
0
0
" 0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
z
H
n
14
25
27
26
21
11
6
._ J
1
0
0
0
0
0
0
0
0
0
0
0"
0
0
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
- o
0
2
»
11
ir
i 1
3H
29
21
13
7
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
0
0
0 '
0
0
0
0
0
0 "
0
0
0
0
0
0
0
0
0
0
0
- 0
1
2
s
8
17
21 ~
29
26
21
IS
5
1
1
0
0
0
0
" 0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
-o
0
1
3
6
|2
1 A
21
18
12
B
1
i
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
I
1
3
7
"1 1
11
12
9
7
2
2'
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
Q
0
0
0
0
Q
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
1
2
3
~ 5
6
6
2
2
1
o~
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
2
3
2
2
1
1
" 0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
0
0
1
,0
0
0
0
- o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o - -
0
a
0
0
0
- 0 "
0
0
0
0
0
o~ ~
0
0
0
0
0
0--
0
0
0
0
0
o ~
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
u
0
0
0
0
0
0
0
o
0
0
0
0
0
" 0
0
0
0
0
0
0
a
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
-o~
0
0
0
0
0
0
0
o
0
0
0
0
0
o
0
0
0
0
0
o -
0
0
0
0
0
o -
0
0
Q
0
0
o
0
Q
0
0
0
o
0
0
0
0
0
o.
0
0
0
0
0
- o
0
0
0
0
0
0
0
0
0
0
0
0
0
u
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
o
0
0
0
0
0
u-
0
0
0
0
0
0
0
0
0
0
Q
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
0
a
0
0
0
-- o
0
0
0
0
0
-- Q--
0
0
0
0
0
-- o -
0
0 .
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o - -
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
o - -
0
0
0
OJ
CO
-------
U-l
TIME. 2500.000
0
0
0
0
U
0
0
0
" ' 0
0
0
0
0
0
0
0
U
0
o
0
0
0
0
0
0
0
- o
0
0
0
- o
0
-- - o
0
0
0
" o
0
- o
0
D
0
- o
0
0
0
U
0
0
0
0
0
0
0
' 0
0
0
0
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
" 0
0
0
0
0
0
o
0
-o
0
0
0
-Q-
0
0
0
0
0
" 0"
0
0
0
0
0
-Q-
0
0
0
0
0
0
0
0 "
0
0
0
0
0
" 0
0
0
0
0
0
o
0
0
0
0
0
cr
0
0
0
0
0
- o-
0
~ 0
0
0
0
"0 "
0
0
0
- o
0
0"
0
o
0
0
0
~0~
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 "
0
0
0
0
0
o
0
0
0
0
0
o
0
0
0
0
0
0 '
0
0
0
0
0
0 "
0
~ o -
0
0
0
- o
0
~ 0 ~
0
0
0
0
0
U
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-Q-
0
0
0
0
n
o -
0
0
0
0
0
U
0
0
0
0
0
U
0
0
0
0
0
o -
0
U
0
0
0
0
0
0 "
0
0
0
0
0
0 ""
0
0
0
0
0
"0
0
0
0
0
0
- o
0
0
0
0
0
o- -
0
0
0
0
0
1 0
0
0
0
0
0
" 0
0
- o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
'0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
0 "
0
0
0
0
0
o
0
0
0
0
0
0
0
U
0
0
0
0
0
0
0
0
0
0
0
U
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
I
1
U
0
0
0
D
0
U
0
0
0
0
0
U
0
U
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
0
0
1
1
1
2
3
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
- o
0
0
0
0
I
- I
3
3
6
5
6
1
1
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
Q
0
' 0
0
0
0
0
0
- o
0
0
0
0
1
2
7
?
13
12
10
8
s
2
I
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
o~
0
0
0
0
0
0
0
0
0
1
2
t
9
|3
21
20
17
'"II"
6
3
1
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
1
3
7-
13
20
25
21
25
- J5
11
5
2
1
0
" 0
0
0
0
0
0
" 0"
0
0
0
0
0
0
0
Q
0
0
0
0
0
0
0
0
0
0
0
" 0
0
0
0
0
0
-o
0
0
0
2
3
- 8
13
20
31
28
25
15
10
5
2
0
0
0 ~
0
0
0
0
0
0 *
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
1
3
- 7
13
19
21
27
21
17"
10
1
3
1
0
0
0
0
0
0
0
" "0
0
0
0
0
0
0
0
U
0
0
0
0
0
U
0
0
0
0
0
U
0
0
0
0
0
0
0
0
0
0
i.
- &
9
13
19
19
11
IU
6
3
0
0
0
U
0
0
0
0
0
U
0
U
0
0
0
0
0
' 0
0
0
0
0
0
0 ' -
0
0
0
D
0
0
0
0
0
0
0
- o~
0
0
0
1
1
2-_
6
8
11
13
10
10
1
2
1
0
0
o
0
0
0
0
0
' 0
0
0
0
0
0
0
0
U
0
0
0
0
0
U
0
0
0
0
0
0
0
0
0
0
0
o
0
0
0
0
0
- 2
3
1
7
6
1
1
1
0
0
0
U
0
0
0
0
0
U
0
U
0
0
0
0
0
0
o-
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
- o
0
0
0
0
n
0
1
2
3
3
3
1
0
0
0
0
- o
0
0
0
0
0
- Q
0
0
0
0
0
0
0
" 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o
1
0
1
0
1
0
0
0
0
0
0
0
0
0
0
0
0
Q-
0
0
0
0
0
0
0
" 0
0
0
0
0
0
" " 0
0
0
0
0
0
o
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
o - -
0
0
0
0
0
- 0
0
0
0
0
0
-- 0
0
0
Q
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
U
0
0
0
0
0
U
0
0
0
0
0
U
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
U
0
0
0
0
0
U
0
0
0
0
0
0
0
U
0
0
0
0
0
0
0
0
0
0
0
0
Q
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
- 0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
* 0
0
0
0
0
0
0 "
0
0
0
0
0
o
0
0
0
0
D
o
0
0
0
0
0
- o
0
0
0
0
0
o -
0
0
0
0
0
- o -
0
0
0
0
0
- o
0
0
0
0
0
0
0
U
0
0
0
0
0
U
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0 "
0
0
0
0
0
0
0
0
0
0
0
o -
0
0
0
0
0
0
0
0
0
0
0
o --
0
o - -
0
0
0
o -
0
o - -
0
0
0
o -
0
0 " " - -
0
0
0
o -
0
0
0
OJ
OO
-------
PICFIC CODE LISTING AND FLOWCHARTS
CFOR«FLO.FORFLO,5 PICFJC
FLOftCHAPTED BY FoRFLO /X8CB/ 0(: 23 SfP 71 AT 1|:16:C1
OCC001
CGUGO?.
COOOU3.
cotiooi .
COOOOS.
000006.
OOOOC7.
CCUUOB.
ocoon'.
ocnnic.
CC CO 1 I
COC012.
OCt'013.
OGU01H.
CGC01S.
CCUOI6.
COOOI 7.
oconis.
OCOD19.
CCt'020.
OOUII2 1 t
OC0022.
CC0023.
OCDU2H.
GC0025.
OCOCI26.
OCU027.
OC0028.
00^029.
GGOQ3G t
CGC031
000032.
OOC033.
OC003H.
COC035.
000036.
CCD037.
OC')o3f).
CCOU39.
0300HC.
OCOOVl
OCOOH2.
OCUOM3.
OOOOHi.
ocnoib.
OGOOM6.
ooont7.
OCOOS8.
OCOOH9.
000050.
OOOOSI
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
^
c
&EAL K IBARtK 1
INTEGER TAPE
PARAMETER NX so* NY« 27, NP icoo
CGM^Of DX|DY,NN,iCI(NPI.KlBAR(NX|NY),DT»T|
lUINXiNT |2) |UINX,NY l2) |X(N'P) >Y (NP)
DIMENSION xniiocoi
PICFIC IS A 2D POLLUTION CODE USING THE PARTICLE IN CELL METHOD
AITH FICKIAN DIFFUSION
NX NUMBER OF CELLS JN X DIRECTION
NY « NUMBER OF CELLS IN Y DIRECTION
NP MAXIMUM NUMBER of PARCLLS
NH m NUMBER Of PARCELS
OX » CELL SIZt IN X DIRECTION
OY » CELL SIZE IN Y DIRECTION
KI « PARCEL MASS ARRAY
KIHA" TOTAL CELLULAR MASS ARRAY (I.E. K I BAR/ ( OX «DY ) . CELL
OT « TIME STcP
T ELAPSED TIME
KC « NUMBER of CYCLES PROGRAM WILL RUN
D DIFFUSION ARRAY D(I,J,II « DIFFUSION COEFFICIENT IN X
D(l,J,2) .-DIFFUSION COEFFICIENT IN Y
u VELOCITY ARRAY uii,j,i> VELOCITY IN x DIRECTION
0(1, J, 2) . VELOCITY IN T DIRECTION
x PARCEL x COORDINATE ARKAY
Y » PARCEL Y COORDINATE ARHAY
x* » RANDOM NUMBER ARRAY
DATA DX,DY,DT,T/I.O,I.O,2S.O,0.0/
SOURCE COORDINATES
DATA xi.YS/o.s, is.s/
CONSTANTS FOR RANDOM NUMBER GENERATOR
DATA X»U>XSG,xX||)iH/O.Cil.OiVZIIt i300D/
O«TA i.s ,NL ,KO/O,O (C/
TAPE IS A FLAG TO OUTPUT PARCEL POSITIONS ON TAPE
DATA TAPE/O/
N II o fj P
-------
in
CO
OCOOS2.
CCUOS3.
OOf-Obl.
OC005S.
OCOOS6.
OC0057.
CCOfJSS.
0000b9.
OC0060'
CCOO*I.
OCC1062.
GCUDfcJ.
orno*i«
OC0065.
CCC.046.
COC067.
CCtl068.
1C CONTINUE
000070.
ODH071.
CCl'072.
CCC073.
SET RANDOM NUMBERS - IHE RANDOM NUMBER GENERATOR RANDN GENERAUS
li RANDOM NUMBERS MTH A NORMAL DISTRlUUTlON
ABOJT A MEAN, XHU, AMD WITH A MGMA, ASGi AND
STORES THESE NUMBERS IN ARKAY XX.
CALL RANDN
DO 301 » I ill t 2
lr(XX( I I GTili >GO T030I
N L « N L I
IFUIL.GT.HP1 GO TO 301
301 CONTINUE
J« l»l
IFIXXIJ).GT.S.IGO T0302
NS-MS»I
IrlMS.GT.riPl 50 TO 302
Ycisi-xxi J)«YS
302 CONTINUE
30 CONTINUE
ccdu76.
OOU077.
000078.
OC0077.
OCOOSQ.
C0008 I .
OCOOB2.
OOCOS3.
conoBt
cctuies.
CCOOSA.
CCt'087.
oct'Ueo.
000087.
OOC090.
OOf'O'l
CCLC92.
OOC'093.
OCl'OVH.
OCC.09S.
OCCO'*'
OC0097.
CCC098.
000099.
CC.CIOO.
ocniui .
COCIL'2.
COC103.
OCC10M.
OOOIOS.
c
c
c
c
c
c
c
c
c
c
c
c
KO IS THE NUMBER OF RECORDS URITTEN ON TAPE
KO=KO*1
*MTE THE » AND Y ARRAYS ON TAPE FOR LATER PLOTTING
IFI TAPC.NE'O) t-RITE If) (X(N> .N«l |NN) , (YINI ,N-1 ,NN)
CO 15 K*« I ,KC
(FOH A SINGLE PUFF NM«NP, FOR A CONTINUOUS SOURC NN« NN » NPS WHERE
NPS IS THE NUMOER OF PARCELS INTRODUCED PER CYCLE)
T«T»DT
THF » ANY ARRAYS WILL BE WRITTEN ON TAPE AND THE SUBROUTINE CELL *ILL BE
CALLED EVERT IOTH CYCLE
IFIKT.NE.OIGO T02S
NO»KO* I
1F1 ,NNI , (Y(N) (N»1 iNN)
CtLL PRINTS THE CELLULAR "ASS ARKAYt KJBAR
CALL CELL
25 CONTINUE
KT=KT»1
IFIKT.EQ. 10) KT-0
THF SUBROUTINE CALCP ADVECTS AND DIFFUSES THE PARCEL*
CAIL CALCP
is continue
If IPLOt.EO.O) STOP
PRINT VOO,KO
900 f CRMATI Ihl ,2CX, 'RECORDES COPIED >,U)
STOP
END
OJ
en
CO
-------
Picric
SY FoRFLO /X808/ ON 23 SEP 7l AT Ii:i6:0b
BEGIN
REAL KIBAR.M
INTEGER TAPE
PARAMETER NX » 50. NY« 27, NP 1000
COMMON OX,DY,NN,KI(NP>,KlDAR DIFFUSION COEFFICIENT
[ IN Y DIRECTION
C U > VELOCITY ARRAY U(I,J,|) .
[VELOCITY IN x DIRECTION
c ud,j,21 VELOCITY IN r
[DIRECTION
[ X PAKCEL X COORDINATE ARRAY
[ Y . PARCEL Y COORDINATE ARRAY
L XX . RANG DI1 NUMBER ARKAY
CO
Xi
i
oo
-------
I
o
DATA DXiOriCTtT/li3i!iCilbC/ :
I
I C SOURCE COOKOINAYES 1
I
I
DATA xS(Ys/B.5,is.s/ '
i
1-C CONSTANTS FOR RANDOH NurfBES 3
l CGENERATOR 3
i
XMU.XSG.XX<|) ,N/0'01.0.921 I>3300/ :
DATA fls,NL,Ko/c,a,o/ :
i
I C TAPE IS A FLAG TO OUTPUT PARCEL!
I IPOSITIONS ON TAPE 3
1
1
DATA TAPE/O/
NN = NP
KC-131
I C SET VELOCITY AND DIFFUSION J
I CARRAY 3
I
I
A-- ---->: oo s i«tiN* :
A :......... .....:
A I
A (
A .<«.< .......
A 0---- >: DO 5 J-ltNY ;
A B :,..... :
A B 1
A B I
A B
AS U(I,J.II-.OI
A B U(I,J,2)*0.
A B 01 IiJi2)>.t-3
AS 01 I.Jill-l.L-3
A o :
A B I
AS t 53 I
A R
...... ..... ..,: CONTINUE !
: :
l
ic SET PARCEL MASS ARRAY 3
i
i
A... . >: on 10 I«I»HP :
A : :
A 1
A I
A
A : KI 111*1 :
A : :
OJ
CO
oo
-------
10]
CONTINUE
c »j
c i
c» ]
I SET RANDOM NUMBERS - THE RANDOM]
Cl,UnllEK fcENTRATOR RANON GENERATES]
c N RANDOM NUMBERS WITH A NORMA
LDlSTRlBUTION
C ABOUT A MEAN, XMU, AND «ITH A
CSIGMA, XSS, AND
C STORES THESE NUMBERS IN ARRA*
CXX.
I
CALL RANDNlXX,N,XMU,XSG>
DO 301
L THIS SECTION SET THE INITIAL
CPAKCEL POSITIONS ABOUT THE
CSOURCE LOCATION
C 1 IF THE PARCEL POSITION IS
CGKEATEK THAN * SIGMA FROM THE
CSOURCt THE
c PAHCEL is REASSIGNED A NE*
CVALUE )
IFtXX(I)GT.t.»GO T0301
N TRUE
>-----.
./
I FALSE
1
NL»NL»i
IF(NL.GT.NP) GO TO 301
\ TRUE
>-----.
,/
I FALSE
I
CO
C 301 J
1
x(NLi«xx< I i«xs :
i
i
00
J^
Ji.
I
I
Z
2
I
CONTINUE
j«i*i
-------
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
I
1 FALSE
I
I
I
I
I FALSE
I
I
: Y(N5)>xxi JI»TS :
i
C 3023 1
: CONTINUE :
I
C 303 1
C KO IS THE NUMBER OF HECOHOS 3
CttRlTTEN ON TAPE 3
I
i
Z
Z
I
z
2
z
1
2
Z
Z
z
I
I
i
z
KO-KO»1
C WRITE THE X AND T ARRAYS ON J
CTAPE FOR LATER PLOTTING 3
IF(TAPE.NE.O)
, >N FALSE
>..-....
/
I TRUE
I
/ WRITE (9) I X .N-l ,MNI , IT ,N-I ,NN ) /
I
I---C
I C
I C»
I
o< ---- .
I
0°
-------
A-
A
A
A
A
A
A
A
A
t
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
4
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
00 IS KK'liKC
: T-T
< IF(KT.NE.(
C (FOR A SINGLE PUFF NN'NPi FOR A3
CCONTINUOUS SOURC NN> SN » NPS 3
C NPS \% THE NUMBER OF PARCELS 3
CINTROOUCEO PER CYCLE) 3
or :
---C THE X ANT ARRAYS *ILL BE 3
ORITTEN ON TAPE. AND THE 3
[SUBROUTINE CELL «ILL BE 3
C CALLED EVERY 10TM CYCLE 3
I FALSE I
I Z
1 Z
1 Z
I Z
< IFtTAPf
1 TrfUE Z
I ' Z
/ 1RITC (9) (X(NI ,N«l ,HNI , (Y(N) ,N«1 ,NN) / Z
I
I
C
Z
---C CELL PRINTS THE CELLULAR MASS I
CARKAY. KIBAR Z
t
i
: : CALL CELL : : z
i
0
C 253 I
i
i CONTINUE : ?3
: KT-K.T*! : i
fALSE
IF(KT.EQ.|0)
1 TRUE
I
KT«3
-------
4 1 C THE SU6ROUIINE CALCf AUVtCTS 1
A I CANO DIFFUSES THE PAKCELS i
A i i
A 0< *
A I
A . .. .to t ...»
A : : CALL CALCP : :
A :.: «.,....,.:.:
A
A L ......3
C................................3
A C« 3
A
A C ISJ
A
--.-..-----.-.-,: CONTINUE :
-c .......]
c................................3
l» 3
,....,,...,..,,.,,,.......,...,.......,,,. x FALSE
< IF(PLOT.EO.O) >
\ /
t TRUE
I
/.i t.. ......... . t ...N
: STOP :
\ /
0<
/tttiit,,..,i, /
/ PRINT 900,KO ' /
/ /
---C '00 FORHAT(IH1iZOX.'RECOHOES 3
CCOHIED « *,I6) 3
("" \
: STOP :
\ /
GJ
CO
CO
-------
9FOn«FLO.FORFLO,S rALCH
FLOWCHARTEO ar HOTLO /xeos/ on 23 SEP 7i AT 11:16:13
ooonoi
030002-
000003.
oooout.
oonons.
000007.
coon,13.
OOOO09.
OO'lOtO.
OCtlUI1.
OODOI2.
000013.
ocoon.
000015.
c;uuiA.
000017.
000019.
C00020.
OOD02I.
C0dfl?2«
COiH)23.
000025.
020026.
OOi'C27.
00'J02?.
020030.
OC 1031.
OCMQ32.
000033.
CCOU3S.
00003&.
C00037.
ooJina.
000039.
03 JoHl.
OC IOH2.
OOOOH3.
oo JOSH.
030UH5.
OOJO'16.
000017.
O0'5oia.
OOTOt9.
003050.
003051.
L CA).CP
RcAL rIBARiKI
PAPAIIETER NX - 5,0, MY- 27, NP 1000
COMMON OX,DY,N8,KI (IIP) ,KlBAR(riX,NY) ,DT iT,
|(MNX.NY,2).U(NX,NY,2),Xi
AND Y 12) DIRECTIONS
N X 1 = N X - I
NYI=UY-|
NX2«NX-2
MY2=NY-2
ZERO THE rELLULAR MASS ARRAY KlBAR
00 15 I"liHX
00 Ib J"l iNY
MflA'« I i Jl«0.
|5 CQt4TIIIUE
THIS SCCTIOM COMPUTES THE KIBAR ARRAY
00 5 1"! iH8
CALCULATE THE CELL THE PARCEL is CENTERED
NnX | 1 )/OX
M,Y( I
IN
CALCULATE THE FRACTION OF THE PARCEL IN EACH OF THE NEIGHBORING CELLS
FY«Y i i I/DY-M-.S
FX»X( I I/OX-N-.5
1F'FX.GT>0> IGO TO I
NH'MN-2
FX=-FX
6 1F(FY.GT.O.)GO TO 7
HH = Mtt-2
FY«-FY
AD" THE FRACTIONAL PARCEL ''ASS TO EACH CELL
*CllAH(NNiMMI-KlC.AK(NN,MH)»FX'FY«KIII)
S CONTINUE
OJ
CO
CO
CALCULATE THL DIFFUSION PSEUOO VELOCITY
OoilO I«2.NXI
001 10 J-2.NY1
IF(K IBARl I . Jl .NE.O. I GO TO III
V( 1 , J, i 1-0.
V| I ,J,2)«0.
-------
C7>
030(152.
033JS3.
0 0 0 T 5 M .
03 JObS.
OOOOS6.
00 3 J 5 7 »
0300S?.
001060.
OC3061 .
00')Q62.
OOD063.
ccuu6'i.
C33065.
000066.
030067.
OC0068.
03006?.
001070.
030071 .
003372.
030Q73.
03307M.
030D/5.
000076.
000377.
000073.
00007?.
003340.
C0308I .
000032.
003033.
031034.
A ft >) n o Q .
U UU U 83
003036.
OOUQ87.
GO TOllO
1 1 I CONTINUE
V(I,J,2)»-0(I>J,2)»(KIUAKII|J*I)-KIBAR(I|J-1))/(2>»OY*KIBAR(I|J))
1 10 CONTINUE
c CALCULATt THE NEW PAKCEL POSTIONS
00 10 I l>l >NQ
C CALCULATE THE CELL THE PARCEL IS CENTERED IN
I-XI 1 I I/DX
J»Y( I I I/OY
IF I I .GT.NX2.0R. J.GT.NY2) GO TO SO
IF( r .LE.7.0H, J.LE.2) GO TO SO
C CALCULATE THE FRACTION OF THE PARCEL IN EACH OF THE NElGHBORlNa CELLS
FX«X( I I I/OX-I-.S
^t1( I I I/DY-J-.S
Mm J »1
M.J«1
1FIFX.GT.O. (GO TO 1 I
Hmtl-2
f Xc-FX
1 1 IFIFY.GT.O. )GO TO 12
MmH-2
FY=-FY
c CALCULATE THE CELL AVERAGE PARCEL VELOCITY FOR x AND Y DIRECTIONS
12 VX=( -Fl)*VlI(Jil))»(|.-FX)»FY»(UtIiMi|l*V(IiMl
111)* |.-FY)«rX«(UIN,J,M»VtN,J,l))*FX»FY»(U»VY*DT
50 CONTINUE
10 CONTINUE
RETURN
CND
CO
I
CO
-------
CALCP
FL<>*CHASTED B* Ff"
...................: CONTINUE :
>
I C ............j
I C 3
CO
f*+J
I
oo
-£*
-------
CO
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
*
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
A
1 C THIS SECTION COMPUTES THE KIBARJ
i CARRAY i
i
i
..........>: DO s i*iiN6l
i
i c CALCULATE THE CELL THE PARCEL J
I US CENTERED IN J
1
I
N-Al 1 I/DX
M.YI 1 1 /DY
NN=.fl»l
NM»,1»I
I
I---C CALCULATE THE FRACTION OF THt J
i CPARCEL IN EACH OF THE NEIGHBOR IN]
I C CELLS ]
i
i
: FY*ri i )/oY-n-iS :
: FX«X< i i /DX-N..S :
I
i
I FALSE
1
I
: NN-NN-2 :
: FX.-FX :
i
C 63 1
I FALSE
I
1
: MM»MM-2 :
: FY=-FY :
i
1---C ADO THE FRACTIONAL PARCEL MASS j
1 CTO EACH CELL 1
i
C 7J 1
KI8AK(M.M)=MPARtri.M)»U.-F«>*tl.-Fr>«KI<
II
KIfUK(NM,M)°KIOAR(NNiM)*l 1 .-FYI*FX*KI ( 1 1
MBAR(N,Ml)-l<:ieAK(N,MM>»U.-FX)»FY»KI(I>
K|BAK<'
-------
VD
A C 53 I
A . . . . . ,. . . .
. ....,: CONTINUE :
i
ic .....3
1 I 3
It* 3
i c CALCULATE THE DIFFUSION HSEUUO 3
I tvELOCITY 3
i
i
A ..,->; ooiio 1-2, NXI i
A : :
A I
A I
A , ....
A B >: 00110 j-2,NU :
A B : .......:
A R I
A B I
A B /lint i ttiit t tt i t t t t t t t t t ft tt t tt t t t t iittt t\ TRUE
A B < IFlKlBARt I iJ) -NE.O. ) 60 TO 111 >- . ....... . . _0
A B \ /
A B 1 FALSE
A B I
A B I
A 0
A B : vi i, j,2>«o. :
A B : , : z
A B I Z
A B I Z
A B /fffitfiiffft ftf fttt\ i
AS : GO TO) 10 : .-.-.-.-... ..................... .. ... ..0 i
A B \..... / ' X Z
A B X Z
A B X Z
A B 0<------------------------------------------------------------------------------------0
A B C 1 111 I X
A B X
A 8 ' CONTINUE ! X
A B VI1 ,J,2).-D( I,J,2)»IKIBAH(I , J«l )-KIBAR(I, : A
A 0 J-l)>/l2.»OY«KIDAK
C« 3
t CALCULATE THE NE* PARCEL 3
tPOSTIONS 3
A-------->: no 10 II-I.NQ :
A ;..... .....;
-------
i
1
;
C
i
L
J---C CALCULATE THE CELL THE PAKCEL 3
I CIS CENTERED IN 3
1
J.Y( i i i/or :
i
i
J FAl.St
I
I
FALSE
c CALCULATE THE FPACTION OF THE 3
CPARCEL IN EACH OF THE NEIGHBOKIND
C CELLS 3
FX«X( i i I/DX-I-.S ;
Fr«r< i i )/or-j-.s
MoJ»l !
I
I
1 FALSE
1
F"«-FX I
1
in i
i FALSL
i
i
#:?» !
i
I C CALCULATE THE CELL AVERA&E 3
I CPARCEL VELOCITT FOR X ANQ r 3
I
I2J I
X
X
X
V
*
X
X
X
X
X
X
X
z
i
I
z
z
z
z
z
z
z
z
z
z
z
I
i
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z
z w
if
z
z
z
z
z
z
z
z
-------
A VX«(I.-FX)
A l.-F*)«rY.
A (U»DT
I
I
CONTINUE
I
1
I
1
I
I
RETURN
OS
-p..
-------
OFOR-FLO.FOKFLO,S CFLL
FLOVCHARTEU o* FOKFLO /xeoe/ on 21 SEP 7i AT 11:16:20
to
connoi
co')on?>
ootinus-
cct'uni.
co'-ons.
ociinn?.
CO'lOClB.
C00009.
coouio.
GOOD i i
000012.
001,013.
COCOI6.
OCU017.
(JOC018.
000019.
SUBROUTINE CELL
RtAL KiDAR.rj
PARAMETER Mx«5Ci»Y.27.Np«loOO
COMhOll DX,r>Y,llNlKI
-------
CELL
i O.'CnArTtO Dr FiA.LO /XOCB, ON 23 iEP 71 AT Ii:i6:2D
\ SURKOUTIHt CELL /
\ /
I
I
KLAL KIB*K,K1
PARAMETER NX«&0,NY»Z7,NP-1CGO
COMMON DX,CY,m<,K](NP).KIBARtNX,fiY),DT,T.
2I . UIKX.fir .21 ,X(f,p) ,Y(NPI
DIMENSION IKlB(NXtNY)
..[ ....t................. 3
£**»>.3
0 3
C THE SUBROUTINE CELL PRINTS OUT 3
tTHE CONTENTS OF THE K1BAR ARRAY 3
: TX-T-DT :
I
I
/.....I,...,...,,,. , 1 I I I t > IK. I ./
/ PRINT IOC.TX /
/ /
I
I C 100 FORMATCIHI,10X,SHTJM£«,F10.3
I C3> 3
I
1
A-. ---->: DO 10 1-l.NX !
A : ,:
A I
A I
A > .
A n --->: DO zc J«I,NY :
A B :..........,....................... ....:
A 8 I
A B C 203 I
A a
A ... .........; 1KIBI 1 , JloMUAHl I |J) :
A : :
A . I
A I
H /..i ........t,i.......n.,.,...f.ft, i, i .,/
A / PRINT 200.1 IKIB)1iJ)|J-1.NYI /
A / /
A I
A I C 200 FORMATCIX,3011 ) 3
A I
A C 103 1
A
................. CONTINUE :
: .....:
i
i
/....,.,......,.,.. >N
: RETURN :
\ /
CO
CO
-------