001 R80104
DECEMBER 1980
ADDENDUM/SUPPLEMENTAL
INFORMATION FOR
VALLEY
CRSTER
ISC
The attached information on three models is furnished as part of
the UNAMAP (Version 4) package to enable users to update their user's
guides for these three models. The coding changes mentioned in this
information have already been made to the source codes included as
part of the UNAMAP (Version 4) tape.
MODELING SUPPORT SECTION
MAIL DROP 14 (919) 541-5335
OFFICE OF AIR QUALITY PLANNING AND STANDARDS
ENVIRONMENTAL PROTECTION AGENCY
RESEARCH TRIANGLE PARK, NC 27711
-------
December 1980
ADDENDUM/SUPPLEMENTAL INFORMATION
TO THE
VALLEY MODEL
The computer code modifications addressed herein are the complete set of
changes made to date in the Valley Model code as it appears in the Valley Model
User's Guide, EPA-4/2-77-018, September 1977. All modifications have been in-
corporated in the Valley code as part of the UNAMAP-4 available from the National
Technical Information Service.
Each item in this report is preceded by the date on which it was issued.
MAIN designates the main program.
A Computer Program Modifications
Make the following modifications to Appendix B of the Valley Model User's Guide.
Note that the code line sequence numbers in the Subroutine EBPLT1 code of UNAMAP
are 100 larger than those in the User's Guide; the latter are cited herein as in
the original source code in the User's Guide.
1. 8/78 - Insert the following line of code in MAIN:
IF(TEMP.EQ.O.)TEMP=293. MA02125C
The origrnal code caused all concentrations to be zero when conversion to
standard conditions was specified.
2. 8/78 - Delete the following lines of code from Subroutine EBPLT1:
P=1.E+10 PL02370C
TEMP=1.E+10 PL02380C
The original code caused the ambient temperature and pressure to be incorrectly
assigned very large values for the analyses of the second and successive
STAR decks in a single execution of the program, resulting in concentrations
that are obviously in error.
-------
3. 7/79 - Replace lines MA046200 through MA046600 in MAIN with the following
lines of code:
DUM3=-.5*((-H-DUMl)/SIGZ)**2 MA046150
DUM4=-.5*((H-DUM1)/SIGZ)**2 MA046155
IF(DUM4.LT.-40.)GO TO 68 MA046160
IF(DUM3.LT.-40.)GO TO 64 MA046165
DUM2=EXP(DUM3)+EXP(DUM4) MA046170
GO TO 65 PA046175
64 DUM2=EXP(DUM4) MA046180
This change precludes underflow diagnostics that may occur in some computer
systems, but does not affect resulting concentrations.
4. 10/79 - Insert the following lines of code in MAIN:
NPASS=0 MA008505
SGRID=0 MA008510
IF(GRID.EQ.SGRID.OR.NPASS.EQ.O.OR.SUM.EQ.O.O) GO TO 72 MA013840
WRITE(6,71) MA013845
71 FORMATCl'///////' VALLEY DIAGNOSTIC.1/1 BE AWARE THAT YOU HAVE CHMA013850
AANGED THE VALUE OF GRID BELOW.1/1 A RERUN IS REQUIRED UNLESS ALL RMA013855
BECEPTORS OF CONCERN '/' LIE ON AN ELEVATED, FLAT PLANE.1) MAO!3860
72 NPASS=1 MA013865
SGRID=GRID MA013870
These additions produce a warning of potential errors due to a possible
oversight by the user. See page 3-2 (Section 3.1.4, first paragraph) of the
Valley Model User's Guide.
5. 10/79 - Change the following line of code in MAIN to read:
IF(MWT.NE.O.O) GO TO 3220 MA020700
6. 10/79 - In code line PL002500 of Subroutine EBPLT1 replace the literal data
(.i.e., ' VIA VALLEY VI'l with 'VIA VALLEY VI 10/79'.
7. 10/79 - Insert the following line of code in EBPLT1:
REAL MWT MA000790
8. 10/79 - Change the following code line of MAIN to read:
If(IZBRIG.EQ.O) GO TO 15 MA041200
Precludes division overflow in Subroutine BEH072.
9. 10/79 - Change the following lines of code in Subroutine EBPLT1 to read:
803 FORMAT(T9,'HLIFE=1,F6.2,'HRS. {etc.} PL011500
IF(W.EQ.O.O) W=1.E10 PL013300
-------
10. 10/79 - Insert the following lines of code in Subroutine EBPLT1:
IF(TS(NSOR).EQ.1.E10) TS(NSOR)=0. PL014750
IF(D(NSOR).EQ.1.E10) D(NSOR)=0. PL014751
IF(VS(NSOR).EQ.1.E10) VS(NSOR)=0. PL014752
IF(VF(NSOR).EQ.1.E10) VF(NSOR)=0. PL014753
DO 11540 I=1,ISOR PL016950
IF(GT(I).EQ.1.E10) GT(I)=0. PL016951
IF(TS(I).EQ.1.E10) TS(I)=0. PL016952
IF(D(I).EQ.1.E10) D(I)=0. PL016953
IF(VS(I).EQ.1.E10) VS(I)=0. PL016954
IF(VF(I).EQ.1.E10) VF(I)=0. " PL016955
11540 CONTINUE PLOT6956
These changes preserve the proper values of the variables examined under
two constraints: (1) more than one stability/wind (STAR) array is evaluated
in a single job; and (2) the source-data records are, for convenience,
input only once for the entire job. Without the recommended changes,
concentrations for the second and successive STAR data sets are signifi-
cantly underestimated -- only under the constraints just mentioned. The
revised program code should be used to reexamine the second and successive
STAR data sets evaluated in a single job with the original program code
when DUPSOR, an input control parameter, was assigned a value greater than
zero.
11. The following set of program changes must be made as a group.
12/80 - Delete the following lines of code from MAIN:
SIGZ=SQRT((A(LL)*ABS(XP)**B(LL)+C(LL))**2 + SIGI**2) MA040500
IF(IZBRIG.EQ.O)GO TO 15 MA041200
IF(L.LT.5) GO TO 14 MA041300
H=HS+BRIG MA042000
12/80 - Change the following lines of code in MAIN to read:
C ISOR,GT,HST,WT,VFORM,AXDIR,IUR,ISHORT,NN,MAP,RECHT,HFCTR,K,P,MWT,MA001900
A HLIFE.ISHORT,IBID,126,IFR MA010500
73 FORMAT(F5.1,F7.0,F6.0,513,2F4.0,13,311) MA010600
* TS(IL),VS(IL),D(IL),VF(IL),IT,1.0,XXP,DTDZ(IT)JTEMP,P,I26) MA034000
CC SIGMA IS INCREASED FOR ALL LOW-LEVEL SRCS XCPT FOR RURAL STABLE. MA040300
* TS(IL),VS(IL),D(IL),VF(IL),L,WS,XXP,DTDZ(L),TEMP,P,126) MA041700
12/80 - Replace the subroutine name "BEH072" with "BEHVAL" in lines
MA033900 and MA041600 of MAIN.
-------
12/80 - Insert the following lines of code in MAIN:
CC 1280 - OPTIONS ADDED FOR BUOYANCY INDUCED DISPERSION VIA SUM OF MA000070
CC SQUARES OF PG SIGMA AND (DELTA H/3.5), AND FOR FACTOR 2.6 OK 2.4 IN MA000080
CC STABLE PLUME RISE(FINAL) FORMULA, AND FOR SELECTIVE EXCLUSION OF MA000090
CC TRANSITIONAL PLUME RISE. P AND T CHANGED FOR PPM. MA000091
D IBID,126,IFR MA001910
CC IBID = 0 NO BUOYANCY-INDUCED DISPERSION. MA013710
CC = 1 VALLEY WITH BUOYANCY-INDUCED DISPERSION. MA013720
CC 126 = 0 MODEL USES 2.4 IN STABLE PLUME RISE EQ MA013730
CC 1 MODEL USES 2.6 IN STABLE PLUME RISE EQ MA013740
CC IFR = 0 TRANSITIONAL PLUME RISE AVAILABLE FOR ALL STABILITIES. MA013750
CC = 1 NO TRANSITIONAL PLUME RISE. MA013751
CC =2 TRANSITIONAL RISE AVAILABLE FOR STABLE ONLY. MA013752
IF(IFR.GT.O)XXP=0. MA033711
IF(XXP.EQ.O.)BRIGUN(IL)=DUM5 MA034110
IF(IBID.EQ.O)GO TO 800 MA040410
SIGZ=((A(LL)*ABS(XP)**B(LL)+C(LL))**2 + SIGI**2) MA040420
GO TO 9 MA040430
800 SIGZ=SQRT((A(LL)*ABS(XP)**B(LL)+C(LL))**2 + SIGI**2) MA040440
IF(IZBRIG.EQ.O.AND.IBID.EQ.O) GO TO 15 MA041110
IF(IZBRIG.EQ.O.AND.IBID.NE.O) GO TO 802 MA041120
IF(L.LT.5.AND.IBID.EQ.O) GO TO 14 MA041130
IF(L.LT.5.AND.IBID.NE.O) GO TO 804 MA041140
XXP=XP/1000. MA041535
IF(IFR.EQ.1)XXP=0. MA041540
IF(XXP.EQ.O.)BRIG=DUM5 MA041750
IF(IBID.NE.O)GO TO 802 MA041910
H=HS+BRIG MA041920
GO TO 15 MA041930
804 BRIG=BRIGUN(IL)/WS MA041940
802 IF(G.GT.O.O)BRIG=G MA041950
H=HS+BRIG MA041960
DUM=BRIG/3.5 MA041970
DUM=DUM*DUM MA041980
SIGZ=SQRT(SIGZ+DUM) MA041990
12/80 - Change the following lines of code in Subroutine EBPLT1 to read:
DIMENSION TS(50),VS(50)fD(50),VF(50),RECHT(112),KFIRI(8),NFIRI(8) PL000900
DIMENSION VERBL (12),URBRUR(2),THREED(12),KBID(8),KCON(8),KBLK(8) PL001000
C ISOR,6T,HST,WT,VFORM,AXDIR,IUR,ISHORT,NN,MAP,RECHT,HFCTR,K,P,MWT,PL001600
C T27,F5.1,T93,F5.1,' VALLEY '/T2,'1280', PL002500
L T42,F5.1,T78,F5.1,/T5,8A4,T56,'0 KM',3X,F7.3,5('KM ',F7.3),'KM1/ PL005200
M T5,8A4,T56,6(1]....+....I),I]'/T5,8A4,T60,F5.1, PL005300
G V(68),V(54),V(77),V(49),V(61),V(69),V(55),(KFIRI(I),I=1,8), PL014500
G S(68),S(54),S(77),S(49),S(61),S(69),S(55),(KFIRI(I)=1,8), PL021700
-------
12/80 - Insert the following lines of code in Subroutine EBPLT1 :
D IBID,T26,IFR
DATA KCON/'BRIGVGS S'.'TABL'.'E CO1 , 'NSTA', 'NT 1','S 2.'.
A '6. '/
DATA KBID/'BUOY','ANCYV-INDVUCED',' DIS ',' PERS ' ,' ION. ' ,
A ' '/
DATA KBLK/1
A ' '/
DATA KFIRI/'FINAVL RI'.'SE 0','NLY,',1 FOR1,1 ALL1,1 CAS1
A 'ES. '/
DATA NFIRI/'ONLY',1 FIN'.'AL R'.'ISE ','FOR ','NONS','TABL1
A 'E. '/
IF(IFR.EQ.1)GO TO 509
IF(IFR.EQ.2)GO TO 502
DO 501 1=1,8
KFIRI(I)=KBLK(I)
501 CONTINUE
GO TO 509
502 DO 503 1=1,8
KFIRI(I)=NFIRI(I)
503 CONTINUE
509 CONTINUE
IF(IBID.NE.O)GO TO 500
DO 510 1=1,8
KBID(I)=KBLK(I)
510 CONTINUE
500 IF(I26.NE.O)GO TO 520
DO 525 1=1,8
KCON(I)=KBLK(I)
525 CONTINUE
520 CONTINUE
* (SCALE(N1),N1=1,6),(KCON(I),I=1,8),(KBID(I),I=118),
* (SCALE(N1),N1=1,6),(KCON(I),I=1,8),(KBID(I),I=1,8),
DO 3805 1=1,8
KCON(I)=KBLK(I)
KBID(I)=KBLK(I)
KFIRI(I)=KBLK(I)
3805 CONTINUE
PL001610
PL002150
PL002155
PL002160
PL002165
PL002170
PL002175
PL002176
PL002177
PL002178
PL002179
PL010214
PL010215
PL010216
PL010219
PL010222
PL010225
PL010228
PL010231
PL010234
PL010237
PL010250
PL010255
PL010260
PL010265
PL010270
PL010275
PL010280
PLOT 0285
PL010290
PLOT 4550
PL021710
PL022450
PL022451
PL022452
PL022453
PL022459
12/80 - Change the subroutine name "BEH072" to "BEHVAL" in Subroutine
BEH072 line BE000100. This subroutine will hereafter be referred
to as "BEHVAL."
12/80 - Change the following lines of code in Subroutine BEHVAL to read:
1 KST,U,X,DTHDZ,T,P,I26) BE000200
DHA = PRC*(F/(U*S))**0.333333 BE010300
-------
12/80 - Insert the following lines of code in Subroutine BEHVAL:
C BEHVAL - BEH072 WITH OPTION TO PERMIT 2.4 OR 2.6 AS COEFFICIENT BE000298
C IN STABLE FINAL PLUME RISE EQN. X =0 FORCES FINAL RISE, AS IN BEH072 BE000299
PRC=2.4 BE005150
IF(I26.NE.O)PRC=2.6 BE005151
These changes permit the user to optionally designate treatment of the effects
of buoyancy-induced dispersion (BID), to designate whether the program will
assign the constant coefficient in the final stable plume-rise equation a value
of 2.4 or 2.6, and to designate conditions for the use of transitional plume
rise. See Section B below for discussion of the options. They also allow
for the continued identification of the revisions made to the Valley Model
computer code.
The use of the optional constants for the plume-rise formula are self-explanatory.
Briggs' latest recommendation is to assign the constant a value 2.6.
Default is zero in all cases, and such default relegates the analysis to that
controlled by the original Valley code. Statements are printed in the output
indicating when the new options have been invoked.
12. Make the following modifications to MAIN:
12/80 - In statement MA009000 delete "OR TO PPM."
12/80 - In statement MA009300 change "CHI IS AT AMBIENT COND." to
"CHI (UG/M3 ONLY) IS AT AMBIENT COND."
12/80 - In statement MA012300 change "TEMP/CMWT*P)" to "298./(MWT*1013.2)".
12/80 - Delete statements MA012400 and MA012600.
12/80 - Change the following line of code to read:
ALWT=8.31E04*298./(MWT*1013.2) MA021900
These modifications describe and assure the proper values of temperature and
pressure used in the conversion from yg/m3 to ppm. Concentrations in ppm
generated by the original computer code can be in error by up to about
+ 15 percent in rare cases, but most often the error would be near zero.
13. 12/80 - Change the following lines of code in Main to read:
3420 IF(PLHT.GT.O.O) GO TO 3430 MA043500
IF(PLHT.LE.-400.0) GO TO 99 MA043700
WS=WS*(400./(400.+PLHT) MA043800
3430 IF(PLHT.LT.IO.O) PLHT=10. MA044000
-------
These changes eliminate insignificant positive errors in calculated
concentrations. The largest of these errors occurred at receptors
no more than 10 m below the undisturbed plume center!ine elevation.
B. Text Corrections to the Valley Model User's Guide
1. 8/78 - Page 3-2, line 5: Change 293 to 298.
2. 10/79 - The discussion of lateral plume deflection on Page 2-2, last paragraph,
is incomplete. A statement to further define the treatment of emissions from
very low heights during stable conditions is as follows: "The height-dependent
attenuation factor for concentrations on elevated terrain is not applied to
contributions from any stable plume with an undisturbed height of 10 m or less
above stack base."
3. 10/79 - Page 3-1, Section 3.1.1, insert the following:
No attempt should be made to force alignment of the printed concentrations
exactly on 16-point radials from map center or exactly on circles about map
center. Consider, for example, the first printed concentration from center
about the southwest direction. The decimal point lies 9.8 print spaces from
center at a bearing of 227.5 degrees from center (using graphical techniques).
With the variable GRID equal to 100 m, the distance of this receptor from center
of the plot is 980 m. This receptor does not lie exactly on the southwest
radial (225 degrees), nor does it lie on a circle about center that passes
through, say, the first receptor to the west of center. The latter lies
8 print spaces from center, or, for this example, 800 m from center. Hence, the
printed concentration decimal point represents, without relocation, the geo-
graphical position of a receptor in the real world.
4. 12/80 - Page 1-9, add the following lines to the list in Section 1.4.4:
• BUOYANCY-INDUCED DISPERSION
• COEFFICIENT 2.4 OR 2.6 IN FINAL STABLE PLUME-RISE EQUATION
• TRANSITIONAL AND/OR FINAL PLUME RISE
5. 12/80 - Page 2-4, Equation 2.1, and page 2-7, Equation 2.la:
Change h to h0.
6. 12/80 - Page 2-5, last line: Change
"T/(MP)" to "298./(M-1013.2)".
-------
7. 12/80 - Page 2-6: Delete lines 2 and 4 (i.e., definitions of T and P).
12/80 - Page 2-8, last line of Section 2.1.3: Delete "during S < 4" and insert
"except that SIGI is set to zero during stable conditions in the'rural mode. The
vertical dispersion coefficient can be adjusted internally, at the user's dis-
cretion, for buoyancy-induced dispersion so that the net coefficient becomes
o = {(aX + d)2 + SIGI2 + (Ah/3.5)2)1/2
where Ah is either the calculated or user-assigned plume rise at the receptor of
concern."
9. 12/80 - Page 2-8, next to last line: Insert as a reference "1975." Also, add to
list of references on page 5-1 the following paper by Briggs:
, 1975: Plume Rise Predictions. Lectures on Air Pollution and Environ-
mental Impact Analysis. September 29 - October 3, 1975, Boston, Massachusetts.
D. A. Haugen, Workshop Coordinator. AMS, Boston. pp59-104.
10. 12/80 - Insert the following at bottom of page 2-8:
The user may designate a value of 2.4 or 2.6 for the constant coefficient in the
final plume rise equations for stable conditions. The value 2.6 is consistent
with other EPA models. The user also selects one of three techniques which
cause the program to invoke the use of : (1) the original code, which selects
the lesser plume rise from the transitional and final rise formulae for all
stabilities; or (2) only the final rise formula for all stabilities; or
(3) only the final rise formula for nonstable conditions and the lesser
result from both formulae for stable conditions. The second of these three
techniques affords consistency with other EPA models; the third may be
preferred for sources near elevated receptors.
11. 12/80 - Page 2-12, line 3: Delete "made."
12. 12/80 - Page 2-14, replace lines 8 through 12 with the following:
cloudy conditions, and can exist in both the daytime and nighttime. In
rural areas the afternoon mixing height is halved in the Valley Model,
and in urban areas it is averaged with the nighttime mixing height, for 40
percent of the class S = 4 occurrences. The remaining 60 percent of the
neutral class (urban or rural) utilizes the mean afternoon mixing
13. 12/80 - Page 2-16, seventh line from bottom, delete the phrase "for the long-
term option only,"
14. 12/80 - Page 3-9, last line, add the following:
An input value of zero (or blank) in effect defaults to infinity for the
calculations, and is printed as 999.99 in the output.
15. 12/80 - Page 4-15, in @HDG card, replace the comma between X and M with a period.
8
-------
December 1980
ADDENDUM/SUPPLEMENTAL INFORMATION
TO THE
USER'S MANUAL FOR SINGLE SOURCE (CRSTER) MODEL
The computer code modifications addressed herein are the complete
set of changes made to date to the CRSTER code as it appears in the User's
Manual for Single-Source (CRSTER) Model. EPA-450/2-77-013, July 1977. All
modifications have been incorporated in the CRSTER code as part of the
UNAMAP-4 available from the National Technical Information Service (NTIS).
Each item in this report is preceeded by the date on which it was issued.
A. 9/79 - Additional Options
Several options have been added to the CRSTER computer program.
These options potentially make the model more useful for a variety of
applications. However, please note that the inclusion of the
options does not constitute a change in the recommended procedures.
The options are therefore to be exercised for regulatory applications
only with the concurrence of the appropriate Regional Office.
The following list identifies the options that have been added to
the CRSTER program:
- transitional plume rise
- anemometer height input
- wind profile power-law exponent input
- stack downwash
momentum plume rise
The new namelist variables associated with these options are as follows:
Variable Name
IPLAS
Type
Integer
ANHT
Real
Description
0 = Transitional plume
rise not accounted for.
1 = Transitional plume
rise accounted for.
Default = 0
Current recommended usage:
IPLAS = 0
Anemometer height.
Default = 7. meters.
Current recommended usage:
ANHT = Actual anemometer
height.
-------
IDWN
Integer
0 = Stack downwash not
accounted for.
Real
1 = Stack downwash
accounted for.
Default = 0
Current recommended usage:
IDWN = 0
Array of dimension 6.
Positions 1-6 correspond
to the set of wind profile
law exponents as a function
of stability.
Default
.1, .15, .2,
.25, .3, .3
Current recommended usage:
P = Default values
IMOM
Integer
0 = Momentum plume rise
not accounted for.
1 = Momentum plume rise
accounted for.
Default = 0
Current recommended usage:
IMOM = 0
The momentum plume rise option allows the momentum plume rise
equations to be used for cold plume sources. This option has
no effect for sources with buoyant plumes - the sources for
which CRSTER was originally designed to be used. Hence, the
default for this option is to omit the momentum plume rise
equation. However, since it may be necessary to consider
sources with non-buoyant plumes, the option which includes
the momentum plume rise equations may be used for such cases.
-------
03
O
•H
4J
P.
O
O O
O O
fl r-
t—t i—i
O O
O O
> >
U 0
o o
o o
CJ CJ CJ
O
O
r-l H
•. — 33
O 2
o- <
CM •>
o c/o
CM «
o *
-* a
CM O
^ m
o u
CM c/o
»i M
CM -H
-H U
v_/ W
O «
-a- u
CM CO
O- M
* *
(U
O
CJ
Z
1
rH
H
O
Pi
ca
w
1
M
Q
M
M
hJ
PM
M
* /-^ *
(sj ^O ^O
"CO
vO v^ Pu
s_- Q .
PM Q H
SO 33
- Z
h-3
u-i M >^s 33
^^ 33 O Z
O U CM <
CO ^-N CM /— .
••* rH CO
rH CO O ^
r-« X -3" PM >-3
'^ .rf CM 33 O
S r| ^ ^ f£
o •x-x5r' i
^x O CO ,—,
Z — < CM S >-5
O r~- " ?< "
M rH CM II 33
CO ^ rH s- , M
53 Z N*-' co ^-^
W rH O 1— 1 U ^!
S PU
vO
\— X
PM^
H
•<,
O
CM
CM
O
»^f
CM
H
x_.
O
CM
i— I
O
CM
M
^^
O
CM
CM
rH
s^
0
cr
CO
0)
o
rH
rH
O
U-i
0)
0)
4J
0)
CO
U
w
H
O
• CO
•• M
rH Ld
-^ II
PM ^
CO
^G
r^S
-------
oooooooomcm
OS OS OS OS OS OS OS OS OS Os OS
ooooooooooo
oooooooooco
ooooooooooooccooooooooooooococoooooooooo
OOCDOOOCQOOCDCOOOOOOOO
*•
01
TJ
0
O
14-1
O
tO
£^
i-l
iH
00
c
•H
i— 1
i-l
O
H-l
01
f
4-1
2
Tj M
•a <;
/ PH
W W <£
to H 2
HH h"H Q^
PS pi o
- S fa
1-1
CM
-a-
00
^~*
•.
PS
o
fa
p
w
H
p"l
O
o
o
<
^ EC x-s
co to P PH M
M Z II P W H
fa O P fa fa O
H U P M M U
U
•< OS
* r^-
* t>»
* 1-1
O
Pi
0
fa
g
Z
— i 3
0 O
OO CJ
i— 1 CJ
o
H W
CO
O l-l
O W
• y
o- p
• pXj
S
§1
pu W
M o
*
*
*
a
/— S
CO
M
'*^'
P
•X
y^
CO
l-l
^— '
CM P
O *
^ CO
H
0 ^
H CO
o£
O 00
EH
*
CO
i-H
CM
*— 1
•
CO
II
1
r^
\D
\B
VO
VO
VO
•
O
*
*
/— \
CO
M
N_X
P
CO
CO
CO
CO
CO
CO
co O
00 *
r~ *
i-H /'-v
to
O t— )
H ^
to
0 >
o *
x-\ tO
• M
m co
• ^<
W r^
O Os
• CM
fa O
2 *
P ii
P PP.
fa" S
CO
CO
CO
CO
CO
CO
•
o
*
*
x~v
CO
M
~o
r^
vO
\D
vO
vO
vD
0
*
*
s~^
CO
1— 1
N^^
CO
>
*
to
l-l
CO
V
m
oo m
r^ o W
-< 0 p
• J3
OHM
H pa H
OHO
CO ^~
oo oo
l-l T~t
to
M
to
3
i-H ^ —
O /"•>
00 CO
'-< M
^^
0 P
H *
O to
O H_
« CO
H *
P •
• CO
W + co
O to o
• l-l 00
H PH -i
W II O
P to H
**~s hH
fa W O
l-l BC O
-------
moo—i^DO^cMu-ioO'-i^£>o»i>f^-r*~i— oocooOO^'—icNif">-3'-d'u~iv£>rs-ooo'^^
CNJCNCNJCNCSICSCMCMCSCMCVjr-lCNJCNCMOJ
f—4
O
o
II
H
O
00
r-H
o
H
O
O
pa
1
Ed
O
•
H
.J
Ed
to
CO
I-l
to
CO
I-l
X_ '
p
•K
CO
I-l
N_^
CO
>
•K
•
m
n
w
CO
ro
CO
CO
o
*
*
/- N
CO
H
CO
?
CO
I-l
to
H
*
H
*
/ — *
CO
M
p
*
H
Q"
*
M
>
*
CO
H
>
4» vO
m VD
. vo
r— 1 XO
II --i
to •
Ed *
*
3
P
n
sc
Ed
P
y-^
to
SB
Ed
P
E~*
»J
•
<^
EE
Q
fa
fa
1
p
CO
M
PJ
EC
II
[d
Ed
p
2S
1— I
H
O
m
O
00
m
t>H
II
*-)
•^
O
rH
O
CO
M
Ed
33
II
*— 3
^
C/2
M
fa
00
CM
CM
•
ft
m
O '
^
co
Z PH P S
<; M i— i M
o O
-~^ -^
g§
P
co
H
<
p
H H
•<<
3 p
H H
<<;
P p
-------
B. Computer Program Modifications
Several program modifications are described below. Some changes
are cosmetic or produce a code more compatible with a variety
of computer systems. Others may affect the output of the model
slightly; however, the effect is of no consequence for applications
in which CRSTER is most commonly used.
Make the following program line modifications to Appendix A of the User's
Manual and check your computer programs for conformance.
PREPROCESSOR PROGRAM
1. 8/78 - In line MET03040 change ' Day ',12, to ' DAY_',F4_._0_,.
2. 9/79 - Remove the statement number 220 from line MET02040.
3. 9/79 - Insert the following line of code:
220 CONTINUE MET02041
SINGLE SOURCE CRSTER MODEL PROGRAM
1. 8/78 - Insert the following line of code:
IF(ELEV(I7,I8).LT.PELEV)GADJ(I8,I7)=0.0 CV022910
The addition of this line of code will reposition the receptors that
are below ground elevation of the plant to plant elevation. This
prevents concentration estimates from being unrealistically low
for receptors at elevations lower than the plant. This line of
code was inadvertantly omitted from the original program.
2. 8/78 - Delete lines CV002000 and CV002100.
3. 8/78 - Insert the following lines:
IMET=1 CV002710
IPTZ=0 CV002720
4. 8/78 - Remove statement numbers 152, 816, and 515 respectively from lines
CV013200, CV0164QO and CV02770Q, and change the NS^ in line
CV026400 to 12.
5. 8/78 - Remove IHC(41, from line CR000800.
6. 8/78 - Remove IHC, and 6,13,18,24, from line CR002700.*
7. 8/78 - Delete lines CR0086QO, CR068800, CR068900, CR071300 and
CR071400.
*This modification is not required if the options under Section A herein
are employed.
-------
8. 8/78 - Insert a comma after:
TABLE' in line CR045500,
'/' in line CR061400,
'/' in line CR061800 and
CO' in line CR067000.
9. 8/78 - Remove ,KSTL from lines CV027700 and CR000100.
10. 8/78 - Remove TX3(2),TX4(2), from line BL000200.
11. 9/79 - Insert the following line of code:
IF(KST.GE.S) GO TO 2350 CR023210
4/80 - Change the following line of code to read:
IF(CHEK.GE.1.6.AND.KST.LT.5)GO TO 2400 CR022400
12/80 - Change the following line of code to read:
IF(KST.LT.5.AND.H.GT.XMH)GO TO 2800 CR022000
The addition/modification of these lines of code prevent inappro-
priate limited mixing under stable conditions in the rural mode.
This may have a slight affect on concentrations for stable conditions;
however this has no effect on design concentrations for applications
in which CRSTER is most commonly used.
12. 9/79 - Change the following lines of code to read as follows:
C BEH072 (BRIGGS EFFECTIVE HEIGHT) September 79 BE000300
C THIS DIFFERS FROM THE OCT. 1972 VERSION IN STATEMENT 24+1: BE000400
C THE CONSTANT 2.6 PREVIOUSLY WAS 2.4. BE000500
C BE000600
DHA = 2.6*(F/(U*S))**0.333333 BE010300
The code in line BEO10300 is the equation for plume rise for
the stable case as it is currently recommended by Briggs*. The
constant "2.6" represents a change from Briggs' previously
recommended value of "2.4", which was used in the original
CRSTER code. This change causes the plume rise under stable
conditions to be somewhat greater; however, this has no effect
on design concentrations for applications in which CRSTER is
most commonly used.
*Briggs, Gary A., 1972: Discussion on Chimney Plumes in Neutral and
Stable Surroundings. Atmos. Environ. _6_, 507-510 (July 72)
-------
13. 9/79 - Delete lines CR015700 and CR015800 and replace with the following
five lines :
C***IF STABILITY=7,SET KST AND KSTL TO 6. CRO15700
IF(ISTAB(IHR).LT.7) GO TO 1631 CR015800
KST=6 CRO15810
KSTL=6 CRO15820
1631 CONTINUE CR015830
The replacement of these lines of code will effectively reset all
stability class 7 values to 6 and calculate a concentration for the
time period involved. Since CRSTER was originally designed for
relatively tall stacks with hot plumes and since stability class 6
results in nearly zero ground level concentrations, CRSTER was
originally coded to assign a zero for the occurrence of stability
class 7. However, to maintain consistency with the RAM and MPTER
models which have been designed to accommodate, in addition, low
level sources for which the stable conditions of stability class
6 may produce significant concentrations, the calculation of
concentrations in CRSTER for stability class 7 (extremely stable)
conditions is necessary. This change, though, will have no effect on
design concentrations for applications in which CRSTER is most
commonly used.
14. 9/79 - Insert the following lines:
DATA HV1,HVX,HV3,KV1,KVX,KV3/150*0.0,150*0/ CR002410
IG=51-KIG CRO31610
IG=51-KIG CRO36810
IG=51-KIG CR042210
15. 9/79 - In lines CR031600, CR036800 and CR042200 change IG=50,2,-1
to KIG=1,49.
16. 12/80 - Change the following lines of code to read as follows:
C— CRSTER—UNEVEN TERRAIN MODEL—REVISED-12/80 CV000100
*' EMISSION UNITS: GM/SEC ', ' AIR QUALITY UNITS: GM/M**3 12/80') CV003800
These modifications allow for the continued identification of the revisions
made to the CRSTER computer code.
17. 12/80 - Change the following lines of code to read as follows:
*HLH(2,24),WS(20),HP(20),TS(20),TITLE(20),IDENT(5),VF(20),D(20), CR000600
*TAPIN(123),JMAX3(4),CHI(180),HE(20,5),DDD(366) CR001200
8042 IF(IMET.EQ.O)GO TO 1199 CR012600
WRITE(6,1500)TAPIN,((HLH(I,J),J=1,24),1=1,2) CR012900
XMH =HLH(IUR,IHR) CRO14500
-------
12/80 - Delete lines CR011500 through CR012500
These modifications eliminate the need to restructure the mixing height
arrays as read from the preprocessor file.
18. 12/80 - Delete the following line of code;
IDIR=DIR/10. CR019100
12/80 - Add the following lines of code:
IF(DIR.LT.O.O)IDIR=(DIR/10.)-0.005 CR019110
IF(DIR.GE.O.O)IDIR=(DIR/10.)+0.005 CR019120
These modifications eliminate the possibility of calculating incorrect
receptor directions due to computer hardware limitations which do not
allow for the exact representation of whole numbers in decimal notation.
19. 12/80 - Insert the following lines of code:
1VT3=0 CR005010
IVTV=0 CR005020
12/80 - Change the following lines of code to read as follows:
8007 IVT3=IVT3+1 CR035400
IF(IVT2.NE.3)GO TO 8008 CR035500
IVT3=0 CR035600
IVTl=IHR/3 CR035700
12/80 - Delete lines CR035800 and CR035900
12/80 - Change the following lines of code to read as follows:
CHIH=CHIH/3.0 CR036500
8052 IVTV=IVTV+1 CR040800
IF(IVTV.NE.IVT)GO TO 8054 CR040900
IVTV=0 CR041000
AVT=IVT CR041100
IVT1=IHR/IVT CR041200
12/80 - Delete line CR041300.
These modifications eliminate the possibility of the 3-hour and variable
hour counter being miscalculated due to computer hardware limitations
which do not allow for the exact representation of whole numbers in
decimal notation.
-------
C. 8/78 - Input Form Correction
The format of the PREP INITIALIZATION CARDS form at the top of
Figure C-l, PREPROCESSOR INPUT DATA CODING FORMS, page -C-3, is in
error and should be modified to be compatible with the format
described in TABLE 4-1, PREPROCESSOR INITIALIZATION CARD
FORMAT, on page 4-11 of the User's Manual.
D. 8/78 - Clarification of the Highest 50 Concentration Tables
Users should note that the highest 50 maximum daily concentration
table allows only the highest of the 180 24-hour average receptor
concentrations on each day to be a candidate for entry into the
table. However, the highest 50 concentration tables for 1-hour,
3-hour and variable hour time periods allow all of the 180 receptor
concentrations for the specified time period to be candidates for
entry into their respective tables.
If one desires a high 50 maximum daily concentration table where
all of the 180 receptor concentrations are candidates for entry,
then the variable time option should be used. To accomplish
this, the namelist variable IVT should be set equal to 24 and program
line CV011600 on page A-13 should be modified to read:
(OR.IVT.EQ.12.OR.IVT.EQ.24) CV011600
E. 8/78 - Random Number Alternative
The preprocessor program requires a random number subroutine not
supplied with the program. The subroutine calculates numbers to
be used in the development of the hourly randomized wind flow
vectors. However, the set of numbers in Appendix A of this
Supplement should be used routinely. These numbers were
used to derive the randomized wind flow vectors for the test
case. They are the set used by EPA for regulatory applications.
These numbers were computed by the Sperry Rand Corp. random
number generator subroutine, RANDU, installed on EPA's UNIVAC
1110 computer; the random number seed used is 65549.
Each line of data in Appendix A contains the day number and 24
integers. The positions of the integers(1-24) identify the hours
of that day for which the integers are to be used in the derivation
of the randomized wind flow vectors. Using the numbers will require
the user to effectively modify the preprocessor program such that
the given integers get transferred into the integer variable IRAND
in line MET01870 to properly coincide with the KHR hour of the IDY
day.
10
-------
F. 9/79 - Running Averages
A new subroutine has been developed for the purpose of calculating
running averages. The subroutine, entitled RUN, has been written
for the express use with the Single Source (CRSTER) Model. The use of
the procedure is optional and allows for the selection of running
3-hour, 8-hour and/or 24-hour averaging times. For each averaging
time selected the subroutine will produce a table of highest
concentrations and a table of second highest concentrations for
each of the 180 radial receptors defined. The second highest table
values are determined using the criteria that there is one other
non-overlapping value that is at least as high as the second highest
value. This procedure is consistent with that stated in Issue 7B,
OAQPS Guidelines for The Interpretation of Air Quality Standards,
No. 1.2-008 (Revised February 1977). Also, for each averaging
time selected a high-50 table will be produced. This table allows
concentrations for different starting hours for the same day,
direction and distance to be candidates for inclusion.
A note of caution, though, is in order for prospective users.
Selection of all three averaging time options will cause an
increase in the CRSTER model run time by approximately a factor
of thirteen (13). Selection of one or two running averages will
show proportionate increases.
Instructions for the installation and use of Subroutine RUN are
given in the first 47 lines of code for the subroutine. The
computer code can be found in Appendix B of this Supplement.
11
-------
Appendix A
12
-------
1 5J87/o>:>7150073oObl014o4M
2 000401304595314825698912
3 4084130328o25ol435121351
4 Ou235b020117965473131362
5 98l469o932C23249230l16426
6 0797o0262290596o69256060
7 274970384076501988562659
8 352919-354431386324767240
9 900803784851094758053409
10 574929000657900221739625
11 48707065o914901570380053
12 230012267523225303690194
13 327384175765415010402542
14 872460604421751060970114
15 4624943olL64847071560734
16 031380026156(1459 05340'.')36
17 123071450609028623312253
10 433796327611595933943570
19 1190726nl503445782401520
20 371017924920746753193832
21 343415319469304228030123
22 6548961o2827119830914957
23 414435319308434769851114
24 031230172812317226669973
25 46895905865276o405542414
26 692280707860060495820082
27 565525145674018025002500
28 14186103308029565224152C
29 169528787374759841713682
30 9ol742691bo3842765872b52
31 826571960820183023035207
32 721310304800830026039656
33" 697937009050706428217437
34 59020046879422oll9067279
35 110962138097719606980862
3o 036o30324966908971963601
37 2740218302208396o6253bll
38 298254688725045587294978
39 459546595807866723917306
40 380152014558301261130995
41 93450953145o837928419815
42 070610727955269981575000
43 3.585367->06c'125282b946929
44 09310ol43585314349803320
45 460029487728619770315227
46 7oo4215802092b2080095048
47* 956853942896790463688358
48 Oo585195259820,-',633024670
49 777814482231015037482769
50 332032418547661181642312
51 27852769ol76248021959004
52 4022020000605722.?5683760
53 044084940100388675200391
54 024357956530090100966427
55 76658103564359«768885911
56 892197076070049525214088
57 5246947o3544250618709598
08 217698779353299905230921
59 330143041033496367143969
-------
oO 3.52l82l-)628>',19'->895^3 33374
62 32379832781900 7763706259
o3 4«70201318407273G8723431
c>4 4'.0l4l9b0835984b37380281
6b O4'll45492b24615353091464
60 8o501b0502073624otS532494
o7 622068244854795638845018
56 372646946090604899002112
69 343415319469304228030123
70 654896162B271198309149b7
71 414435319308434769851114
72 031236172812317226669973
73 4689b9058652766405542414
74 692280707860060495820082
75 56b52bl45674018025002500
76 141861033880295652241522
77 109528787374759841713682
78 961742091663842765872552
79 826571960820183023035207
80 721310304800830026039656
81 697937009050706428217437
82 59020046879422611^067279
83 110962138097719606980862
84 036630324966908971963631
8b 27402183822083966o253611
86 298254688725645587294978
87 45954ob95807»66723917306
8U 380152014558301261130995
89 934509531456837928419815
90 0786107279b'5269r>81575008
91 3385367t>0691252826946929
92- 893186143585314349803320
93 460029487728619770315227
94" 766421580209262080095848
95 956853942896790463688358
96 Oo5851952598288633824670
97 777814482231015037482769
98 3320324185476hll81642312
99 278527696176248021959004-
100 482282808068572225683760
101 044084948180388675200391
102 024357956530090100966427
103 76658103564359"768885911
104 892197o76070049525214088
105 524694763544250618789598
10o" 217698779353299^05230921
107 330143041033496367143969
108 332182962881998959333374
109 933449097112500240729843
110 323798327819007763706259
111 497020131840727388723431
112 440141950835984537380281
113 044l454925246153b3091464
114 865015050207362466532494
lib 622068244854795638845018
116 37264094609068489^002112
117 343415319469304228030123
118 654896162827119830914957
11Q 414435319308434769851114
14
-------
120 G 3123o 172012 3172.26609973
121 4oO°>59050fc>527ov-.4 05542414
122 o92200707860Oo0495820082
123 5ob52'H45674010025002500
124 1410610330M029565224152;:
12b Io9520707374759041713602
L2o 961742691oo30427o50725b2
127 026571960820183023035207
128 721310304800830026039656
129 697^37009050706428217437
IJO 590200468794226119067279
131 110962138097719606980862
132 036630324960908971963601
133 274021838220839666253611
134 29025460U725645b07294970
135 459546b958070bf;>72391730t>
13o 3801b2014bb8301261130995
137 934509531456837928419815
138 070610727955269981575000
139 338536760691252026946929
140 09318ol43585314349003320
141 460029487728619770315227
142 766421580209262000095840
143 956053942096790463688350
144 065851952590238633024670
145 777814482231015037402769
146 332032410547661181642312
147 270527696176248021959004
148 482282800060572225603760
149 044004940180388675200391
150 024357956530f)P0100966427
151 76u58103564359l'7600a5911
153 092197676070049525214080
153 5246947o35442506l0789590
154" 217690779353299'?05230921
155 330143041033496367143969
156 332102962031990959333374
157 933449097112500240729843
150 323790327019007763706259
159 497020131040727308723431
160 440141950835904537300281
lt>l 044145492524615353091464
162 865015050207362466532494
163 622068244054795630845018
164 37264694609060489^002112
165 343415319469304220030123
I6o* 654096102027119030914957
167 414435319300434769051114
loO 031236172012317226669973
169 460959058652760405542414
170 692200707860060495020002
171 565525145674010025002500
172 141861033080295652241522
173 169520707374759041713602
174 961742o916o3042765072552
175 026571960020103023035207
176 721310304000030026039656
177 697937009050706420217437
170 590200468794226119067279
179 11096213809771^606980862
15
-------
180 U3o63032496o008971963o81
181 27402l838220H39t:>o<>253611.
182 298254oaJ725oHb;:>u72yHy7a
183 4b9b46595807866723917306
184 380152014558301261130995
185 934509531456837928419815
186 078ol072795b269'>8157500a
187 338536760091252820916929
188 89318ol'13585314349803320
189 4o0029487728619770315227
190 766421580209202080095848
191 956853942896790463688358
192 065851952598280033824670
193 777814482231015037482769
1^4 33203241854766U81642312
195 278527o9617o248021959004
19o 402282808068572225683760
197 044084948180338675200391
198 024357956530090100966427
19« 766581035643590768885911
200 892197676070049525214088
201 524694763544250618789598
202 217698779353299905230921
203 330143041033496367143969
204 332182902881998959333374
205 933449097112500240729843
206 323798327819007703706259
207 497020131840727383723431
208 440141950835984537380281
209 044145492524615353091464
210 805015050207362460532494
211 622068244854795638845018
212- 372640946090684899002112
213 343415319469304228030123
214' 654896162827119830914957
215 414435319308434769851114
216 03123ol72812317226669973
217 468959058652766405542414
218 692280707860060495820082
219 565525145674018025002500
2?0 141861033880295652241522
221 169528787374759841713682
222 961742691663842765872552
223 826571960820183023035207
224 721310304800830026039656
225 697937009050706428217437
226* 590200468794226119007279
227 110962138097719606980862
228 036630324966908971903631
229 274021838220839666253611
230 298254688725645587294978
231 459546595807866723917306
232 380152014558301261130995
233 934509531456837928419815
234 078610727955269981575008
235 338530760691252820946929
236 893186143585314349803320
237 460029487728619770315227
238 760421580209262080095848
239 9bo85394-28967904o3688358
16
-------
C40 Oo50519525982d8b33024b70
J41 777814482231035037482769
242 332032418547601101642312
243 278527696176248021950004
244 4022820000605722J5683760
245 044OQ4940100300o75200391
24o 02435795653009Q10096o427
247 7ob58103564359<->768085911
248 89219767o070049525214088
249 524694763544250618789598
250 21769877935329^905230921
251 3301430410334963o71439b9
252 332182962881998959333374
253 933449097112500240729843
254 323798327819007763706259
255 497020131840727388723431
256 440141950835984537380281
257 044145492524615353091464
258 8b501505020736246o532494
259 622068244854795638845018
260 3726469<4b090684899002112
2bl 343415319469304228030123
2b2 65489bl628271198309l4957
263 414435319308434769851114
2b4 031236172812317226b69973
265 468959058652766405542414
2o6 692280707060060495820082
2o7 5o5^25145674018025002500
2b8 1418ol033880295652241522
2o9 169528787374759841713682
270 961742o916b3842765872552
271 826571960820183023035207
212 721310304800830026039656
273 697937009050706428217437
274' 590200468794226119067279
275 110962138097719606980862
27b 036b3032496b908971963681
277 2740218382208396o6253611
278 298254688725645587294978
279 459546595807866723917306
280 380152014558301261130995
281 934509531456837928419015
282 078610727955269901575008
283 338536760691252826946929
284 09318ol43585314349803320
285 4b0029487720619770315227
280* 7664215802092b2080095048
287 95685394209b790463688358
288 065851952598288633824670
289 777814402231015037482769
290 332032418547661181642312
291 278527b9617b248021959004
292 482282808068572225683760
293 044084948180380675200391
294 02435795653009010096b427
295 76b581035643599768885911
296 892197o76070049525214008
297 524694763544250618789598
298 217690779353299905230921
299 3301430410334963671439b9
17
-------
30 u 3321829b2Bai9°8^b9333374
J01 9334490971l2b00240729843
302 32379832781^0117/03700259
303 497020131840727388723431
304 4HOl41950a3b9-84b373802Bl
30b 044l4b'492524615353091464
306 865015050207362466532494
307 622008244854795638845018
308 37264o946090o84899002112
309 343415319469304228030123
310 6b48961628271198309149b7
311 414'135319308434769851114
312 03123ol72812317226o69973
313 4689590b865276o405b42414
314 692280707860000495820082
31b bob52514b67401802b002500
31o 141861033880295652241522
317 169528787374759841713682
318 961742691663842765072552
319 826571960820183023035207
320 721310304800830026039656
321 697937009050706428217437
322 590200468794226119067279
323 110962138097719606980862
324 03o630324966908971963681
325 274021838220839666253611
326 298254688725045587294978
327 459546595807866723917306
328 380152014558301261130995
329 934509531456837928419815
330 078610727955269981575008
331 338536760691252826946929
332- 893186143585314349303320
333 460029487728619770315227
334' 760421580209262080095848
335 956853942896790463088358
336 065851952598288633824670
3.57 777814482231015037482769
338 3320324185476oll81642312
3J>9 278527696176248021959004
340 482282808068572225683760
341 044084948180388675200391
342 024357956530090100966427
343 76o581035643599768885911
344 892197676070049525214088
345 524694763544250618789598
346% 217698779353299905230921
347 330143041033496367143969
348 332182962881998959333374
349 933449097112500240729843
350 323798327819007763706259
351 497020131840727388723431
352 440141950835984537380281
353 044145492524615353091464
354 865015050207362466532494
355 622068244854795638845018
356 372646946090684899002112
357 343415319469304228030123
358 654896162827119830914957
359 414435319308434769851114
18
-------
3ol 408959058652766405542414
3o2 692280707860060495820082
363 565525145674010025002500
14186103^880295652241522
169528787374759841713682
3oo 961742691663842765872552
19
-------
Appendix B
20
-------
txaeacac.aiatececacacaeacactKccai:ai
o. u o i-i
co cc < o
O Z 3 Z
Z t-i O >• 1-1
OJ 9-
Z
O «C
Z
a.
CO
»
«•
*^
z
3
Ct
< O
O X
I
O v*
If) (X
i-i O
>»
UJ O
z z
f-i «t
3 Ce
O 3
cr O
m x
3 I
to K"
CO »
l-l CC
x 3
K O
X
U- I
O l»-|
z cc
o o
IAJ O Uj
-J Z
LS d I-i
Z UJ t-
1-1 LU 3
to Z O
Cr
UJ UJ
X 2 Or
»- O UJ
K
x *• to
K U, CL
i-i Z <-)
t- uJ
UJ 3 X
10 O t-
3 IX
CD O
cc 3 t-
O f
u.
o i-< z
UJ X U.
o. K z:
O uj
-J X h-
LU to <
UJ «J tO
o cc
< -I
to »- -J
< to «t
>-i o
O O-
a;
Z lu
3 CD
u- «t
o:
LJJ uj
x >
u, o u;
oe >- z
o a _j
ct o o
a. *• u.
X X >•
•9 "> >
O O i-i
< < t-
o o «->
» » UJ
^> *^ Q.
KI ro to
v^ <-< UJ
LLJ uj a
CJ U)
u: £C O
3 3 in
O O f^-
in r*-
» » (\j
^% *-s O
UJ LU
Ct O O
cr or z
3 3 <
O O
10 10 o
» » 1/1
*^ ^s r—
T- r- CC
w ^^ r*
UJ UJ O
u u >
c: 40
=j r> of.
o o u
to to uj
o
K
cu
3
CO
LU uj «:
X to cc
K 3 uj
>
u. to «
O »H
cc I
O X
Z •-!
« • X
to
• k UJ O
to •J) 2
LU < «c
1C CC
ce > to
UJ < UJ
> X
^ to tc
»i < <
^T ^^ ^^
E
U) J> OJ
2 Z >
oe tr »-i
(/> UJ
i-i a.
«: «t or
o a
• » O
rvj f\j ir<
x x c:
H- t- r-
»• •• (v.
r- T- O
X X ft:
K h- W
r »
i-i i-i a
x x z
»- r- O
x x u-t
x x c
H- t- at
»*•(_)
3 3 -/
CC £C
-J —I
I- * Z Z Z
K UJ Z
ri O ^
u. 3
h- cc
LU 3
x a ct
(-23
»-. o
•• X
uj a I
Z 3 Z
Z K O
3 tj
Or Cr UJ
3 to
or o
3 X •
O I t-
LU
O O
LU •-(
Z 3 cc
tc to •
uj a. u. r
t- o o »•>
UJ «t fv.
a o z r-
o o
UJ UJ •-« tr-
ee X K
< H- < >-
K CC
O Z uj <
LU »-! cc 3
t- o. ce
< CD CC JO
_j r- uj uj
3 K u.
o UJ Z
_J 3 •-» O
< t/i UJ
l-l tO LU to
I
au
UJ u>
e- to LU
3 K <_>
O «J 3
a. o
cc uj o
3 X ct
f- a.
x O
u 3 X
3 a
o o <
O CC
cc a. »
CL Z
_l LU
I -J -I tO
_J l-l O
i-i 35 X
rg
CO T- UJ
UJ >• r~
to cc
3 UJ O
X l-l >H >-l
Z i-i 03
UJ CC
12 iu
< O
CC Z
UJ 3
>
< O
UJ
13 Z
Z ft
»-l -J
z >-
Z 3
30
ce
LU
X UJ
to o
>^ O
X CC
a,
z cc
O K UJ
M < K
H- U)
CL to «t
O o cc
ce <
O < X
< -J
a cc cc _i
UJ UJ UJ M
K K > 3
»_>«-)«
«* < LU
CC CC X _J
«t •£ l_J Oj
X X « <
<_> t_> UJ I—
Z UJ
O X
<->
UJ
ct or
O ^
U.
00
to o
UJ O
Z I
•-« f\J
-J •
UJ «-
30
o to
uj o
_J LC
•— <
1-^ o
K- z
z <
LU K
CO
I- Z JV
to U>
Z
X to
t- 3
3 _J
o «t
O 3
to
z
o
to
LU
ce
u.
o
o
o (-
3 3
K <
U.
to uj
K «
to 3
a.
«J Z
O
Z •-!
>-i h-
O.
z o
o
l-l UJ
o cc
Lu
Z >
o <;
3 Z
a, »-•
M z
ac z
K 3
z cc
o
cr. »
3 *>
O
-------
0000300000000000000000000000000000000000
«— f•« h >T i- ) O I- £1 O- O «— ''- ». ^ Lr -<' f"- UJ C.' O T— fO f"1 >J ITi O Is- <• J O «— f\J Wl -O U~v -O r- O.j (? O
^^^^»»^»jN3>»Li^u^i/-inu^u'iu^untriir\<><;^i<)<>HD'O'O^)<5f^-f^-r-r>-f^f^.fs.f^f^.|v.uj
O p C, p O O c: O O O C, O C- O O C' O O O C O C. O O C.' O O C. C C- O C' C.' O O C' O C' C.' O
O .:: o o i_; ^ •. c- _> ^ t • • • _ . . .• _ .j i. ._ _ c.i • • O _ O C -li ..i - j ci '. • >j :_j o ,_; o •-- - • _• _.
2_JU.
>-H uj o
O oo -J O »-«OOOO-OO^OOioO(\4>t»r-t-
Q.S in IM»»-»»
UJ «C u. » O • irt * » IA * * *»»»OO-O
xioo <-»>-i o »• aj r~ • uc. h- r^o^T-roT-rsiio
^ Z ^sO» (\J-*(\JT~«*(\|T~r-ro^Oj»«~r-
O r- • JNJ r- » CM T- r-
oo »oj»»(\j»» »
_J>- «-• iPioOOoO«*
UJ<«LU d > in r- i-» ••.»»^•»•»K•••»•r-»-^•>
aa OCD 1-1 < x • «^ ••••2 f^otn^-oir, -4-»>ot3oo»T-«—
K » » 1-1 <-« z ^x^^rj T- K. f\jT-r-">r\j»-r-'«-Kitr\o'. *-»••
>- O u. O ^^ /> »t~ 3 iAi/>»-i »»»»»»»»»»»»Of<^K>
«£T. Oh- irvO*^»AX rvi (\J » Ol/^•*K1u^«*Mu^K1l^^^>-C>•T-r"IO
5; <^oOr— >^< •> » ^% t— K1(\JTT— tOtrif»-»T-»—
cc O O •— z'wv^rga acuorg »
«t cc i-i _j a. a: f> o *> »-»-••
ZO.UJ r\J<>--(>-( ct wfcj*^ou^ooa>^<^' »«i»»N-'Omr~»j*:^-lv.f\4r-»r\ii^' »f\t
i_> <— i < >H » >— -rf ••*»•» 4/> »• t-i * Kt £D»Ki»fxi^«— f\jr-T-rvj
zz »«t-^>oo»rvJO30'<*tri •.••.»»»»
IM
< x 1-1 I- a: »-i>>-»e»o*»«-r>j
-------
C COOOOOCJOC-OCOC C-GOOOO OOOOOC. OOOOOOCOOOO
o o o o o
o
h -o i • H- r--
<-< uu -^ cu uo
^ c c_ c. o
_.- r -3 ^-. >r.-_; ' . > i
r 2if22'ZZZ?>22
O O O O O
K-> ~j u, <• r-
CJOOOOOOOOOOOOOOOOOOOOiJ ,
OOOOO-O ;> O- &• O- O t>OOOOOOOOOO»-r-*— T-«— T-T-T-T— r-fv,
C O l_> C C. t— O C. C- O C_ »- t- «- t- *- «- «- r- «-»-«- T- r- »- «- T- r- ^ ^- ,- r-
r . •.; ' . • i a L >_. r: i c. "j o -;j D O O o O o O 3 c. _> o o o o o in ij o
2222222222ZZ22Z:22:22Z;222222-
Ql f-' Of LIT Q^ LU U- frj ^v Q£ ry .•• OiO^£EQiQCttC£uCOlft:LC«iQCiJl
1/1
u.
o
o
o
CJ
cc
UJ
>
u
LU
_J
UJ
u:
• o
13
Lu
o
a..
at
UJ
z
o
UJ
U)
CXI
r
c
o
O f
^ r- »— II
2" II II >*
tr O «o co ro
3 t- ^^*ccijCJr/i^ o *•*
»—» (Mf\J o I— LU LU uj II >-i ro ro LU
r~(_> UJC522Z>-i*2<;ZI •
»-tOT-»»»»» O Ct CE CC 2 O
»->-i*Z2<£CtJ<— 1-1 3 i-i 3 II i-i
CJ <*/
in M
«- X
» o
T- II
I! <-.
K ee
Z
Z
Z)
a
a:
3
o
X
i
Kl
ee
O *^ <->
3 i-«
2 II
1-1 cc
t- JE
Z 3
O to
O M
ro
«
or
*
*
k 3 O •• O O
C3 *- is> r- r- r- K
• II ii II II II II
Ortt: »*^iuOCC^->^>
£>jcca:3«4^crcc
3 -J _J 2 _l _J
toOw''-'t-tOiri*'w
>-i«-rorO»-»-ruKirO
»•» _j r- z i- o
u.c
i-io>i-ou3c.i-«i:
o o
«- ro
23
-------
C.OOC-OOOOOOOOC'COC'OCOCOOOOOOOOOOCOOOOOOC.-OC-
OdOOOOOOOpCJOOCOOOOOOOOCDOOOOOOOOOOOOOOCiOO
r- <\i f -J W\ -LJ f^- t: • c? C. T~ f\j r«~) *j u"v x-. P— wj ^- I—' v- rsj I* i *J i/~ »O r*- cc u O «— f^J r*"i «*7 u~ -o r^ ex, O O
ru r-j rg ro rvi fu r\i ru rxi »«~i ro ro KI ro KI to to fi to •* ^ •* -J •*•*•*>» >T»J-UIU-, tnu-\mmi/".minu-\ O
(_
oo
(NJ
cc
s
C
M
CD
^ -J
to
CM
O
h-
O
I I
V) V)
•w w ro
ro KI
CD A
O
ou
K •-
a; »-i
n n
n
<^
O
II
a:
rs
O
(M
I
ac
c
=>
M
n
K
O
«-s t- er. CL.
II
tr
w >-,«-' 1-1 o~
»- to
Z -J
O <
L-ji ^^
O O
o e.
O «•*
•-• K)
co >
«•>
O »fc
o
IU CC
* s
at o
3 *^
O I
»•* «*•
w n
u. CL
IM
II < M Z < rl
IM (M
IA
•»•
O
a
O
U)
O -* III
u^ ts* t/} tsi **
>-i o _i _J -i cc rg
x i- w w *» ^ I
• to to to «— tr
UlO CL O O Cj Kl£l
O _j CD >-i trttnin OS
in • /^ » 1-1 I d Ik «/>
»*sO r» Xi-t<-" < >i
r- cc in II It II II II II
ll_j» «ov>^«« ~* **•
_JtOUJ_l_JJ-l_l_l3-J^
U)»| |WW^ZW»*
«* >H-OK)T~MfOtOt-<»nto
«* «CJin«*-U^OOOK
«^ *^ n n in ir> in z
OU,U-CLO«- I-I
z >»
O U.
CJ »*4
m o
in
UJ
IS
-------
COGOC C_C GOOGOCOCOOOOtOOOOOO OO OOOOOOOOOOC-O
o c; o o o o o cr. 00000000000000000000000000000000
T— <"NJ '• > vj i/. £> r^ 'X a t_: <— rvj f -Tli.olv»Ci.Ji,>-Oi— OjKl»Tinor>-uiur O r- f\J If i -T u'l <3 t»- uj Cf LJ
•o •-fv.(^rN-f^fs-wcc. oco.ioDa>aoaou.oouo-s>oocr o> o o o o
»-»-»-»—«—»— T- ^- t— r- «— »- T— «- T- r- «— «— »— »—«—»-»— r" *— r- «— r- «— r- r- v- r- ^-T-T-T— <^ «— (X
r '- _ '. .' C •_ JJ C- '. • C- • CJ '_• O O O O O O O CJ O O O O O O O O O O O O O O O O O a O
ZZ'ZZZfZZrZZZZZZZZ'ZZZZZZZZZZZZZZT'ZZ'ZZZZZ'ZZZ
1/1
o
(X _J
_» «-»
O »-
u:
o
z
z: o
o u
u ••»
O ~ II
or Cr » ^
£ «J oc cc
2 -J -i
z
o
cc.
_J
s^
CO
-J
•4 •
> 00
I -V
f* >^
ct Or
O
t-<
9-
ac.
oo
-J
«t
IX
Z3
O
»^ r-
OC- tf\
CD
« O
»- >-
oc
• o
I I
M 10
U. N-
'-'•*> 00
Ou OO U)
01
II
uj O a:
=» rv j
z
I-I O CO
CO CO
H- >- N.
O O »
K- K T-
II II II
^ ^H O
IX Ct r-l
CO CO
u. O
i-i O
000
^J ^2 /^'
•o o
»» in
CO EG kit
I- UB
O > O
>- «t o
*. u-
oc. 2
-J Z
s* »-t
00 K
O
II
OO
9
K
_l
oo
o
oo
to
> v>
41-1
II II
r* r*
a: QC
S3
O O
O
on o
II —1
t/> to •
UJ LU
II II
a: H
o
h-
o
OU f^
o *>
v\
•-1 O
Z K
•
LU O
_i to
• •*
f* O
a: in
fM
in
• II <•> •
«- a. DC. o
II a -j LU.
a o »* •
^**10OKO O )GCOoOOO>-*COOOO>-<>-ilA^>
*O W^ *^ ^^ ^^ V% ir™ pp C2 H» rp ^^ ^«~ f««» ^^ *o *^ ^J
w w || II < i-l Z < i-i ZZ **
O O u» ^* &- O V5 f* h— O ^*~ ^ O O O O O
OO ^- O O
II
V) V)
I I
Crt CO
cC ao ao
C O C*
u-> m m
M X A
X •-• M
II II II
I
o
u- u. a.
I
Kl «-
o in
U
O ts>
CO OO OC 1-4
O O O K
u~\ in m z
h- I a O
«0
25
-------
Ot-'C-C CJOC.OOOC OOOC. OCCiCOOOOOOOOOOOOCjOOC-DOOOO
ooooooooooooooooooooooooooooocoooooo
^j LT\ ^TJ ("*- oc t^ c ^i ^: zi 3
r-j r i f\i <\> t\-
II II II
K K H- »n
-J _J -J -O
CO OO OO O
CJ O O K
m m in
l-i I O O
z
z
o
X.
I
o
t-
o
-»
•
Ul
u
•
C£
c
*~ z
O -f
r >—
«- w II
II *^«
cr cr ft
f\J Ct r^
<\J
03
< o
rvj o
<.
a:
O u.
»n o
fVJ
o « o
O > h-
' CC « l-l «t l-l
«^ «*• J Z »O I- K II M
>»»*• < O OmUJO Qt»-i n r+
or\i u cu oj _j <-\ t. u u croc
oo »— »— KI u u »— • • m •-• «-» ^ 33
T- O C> fSJ *> »^ r-^- »N ro » t/> I/I O O
»h-K»OI •* r-lla:» r- _J .J (-!•-•
«-«-nii«-i-i tv n cc _i o a ii ••• »»
oo II « ^ u •• •> =)wuj3t^c^(£a: occe
Luoet£tctcsa:ujo:ujQ:O>» «OKKJ-Jujj_j
ZI^-J-J-J1-!—I3-J3 Jt-i«MOtf«-iJ-J***'3*''-'
Z ww wzwz U33I I«*''**^'^OO' >omoc*O f\i f\j ^H rsj CNJ
Z 1-0 _J Z _J Z w ^ II II < 1-1 Z < >-l
ooo>o<:oo^-«to>u><--'>-i>--iuQ;t-«
o
on
*
*
ZZ
O
OO
(_> W «->
26
-------
OOCGOOirOCOOOOOGC'C-OCOOOOOOOOOOOOOOCOC. O O O O
ooooooooopopooopopooooooooaoopoooooooopo
r-
r\. r\i «\j
O <_> C3
i rv f>j r>j f^J rv fv^ ru cv f\/ r\i t\i (\j f\ t\* fvj r^ r\* f\- cv (\i cu r\j fd (\^ r\j
'-J O O w.1 O O O o O O O CJ O "1 _' _ O O CJ O O wj C3 'O O
ZiSZZZZZZZZZZZZZZZZZZZZZZZ
a
z
C*
UJ
1C
CC
Ul
>
o
»-
o
X.
z
oo
o
ct
o
UP
z
to
UJ
C 00
in
1-1 O
UJ O
_i u>
o • **
lf\ **: O
» at in
I I I
f> V}
CL
l~
o
n w a
K *J LL; K H t- ^
_i rvj • -J -I -i w
Ut K t I «*
co < ^** in oo in O
^f II II in
O U. U, £L O Crt »-•
r- n II II II II II
H r* ^ r\ r*> r\ r*
V) (/> 1S>
*f UJ (\J
r\j 3
•K «» Z CC
ce (\i •-» ui
i-i II K u.
** £ Z 1/1
u. Z O Z
ce
O r- h-
O CJ *
C r- *
O
O
in
O
0
Ui
CJ
O
a
rj
••o
CC O
.J -J r- u;
w w (J. 3
v* Z
Pd O r-i
O K t-
1-1 hi ce
o o
o
»•
o
27
-------
o <--
c. o o c-- o c-
c. c
o o o o o
o c- o o
c: o
c, c
L-.'Oi'-w o c. T- ro i*. ^ m <<" r^ .,_> u (_• :> C OOOOOOOOOOT-»-t-r-
(\i r j r-, f •- (vj i« , r«~. •"•*•' t-.F.pip>r>-f,t- !• , r» . f-
'
o c, c. c
OOOO
r~- ^ o c
tt
o
in
UJ
_J
CC:
00
U;
UJ
3
Z
(X
1. •
a. tfl
UJ Z
CX I-l
o. I-
Z
O 1-1
I- ee
£L
UJ
2" a:
*i O
O UJ
IX £.
O
UJ
O
o
o <
*
*
<-» (J O O
28
-------
ct
o
ct
O
(O
UJ
QU
X
IS
z
2
CC
ct
ct
u.
HH
II
UJ
z
HH
z
o
UJ
I/O
*
vs •
z u
HI z
Ct r-t
< h-
0- 2
UJ I— 1
CC CE
CL CL
cc ee
O 0
LU UJ
Z £.
O LO
ex z
UJ «t
CHI ce
'.*
hi
1
1
^«
CE
>» «t
fM •-«
11 It
/•» <->
cc , u. <-»
HI ce x a
» O fv- X
f\J ••
CX _J
o _i w
ec ^ •* >»
«— •* r\j i\i
» r\j Q *
•- CTJ HH ^>
n < t- ex
Ct I- HH .J
_l CX II S^
II •-* X
ru *> ex :*
c\j ex _J «c
S£ * -» HH
v a: o c
VI _l • •
ru
o
a o
LU IJJ
h- <
II II
CX CX
II
Of
LU
II
z
> T- —I *-» HH V «rf
CX
LU X
Z w
UJ «C
iJ3 LU
*
*
-> O
O- ^ 1£ II
o a <.
3 2 Z
«t < HH
C i
o »-<
II t> i
u. u. O "» O
Z CD
HH «/)
II II
Z *
at **
!-» J*
II «t
Z II
H-(
=>
r-
X HH
« Z
X. n
DC
a.
cr
O
UJ UJ
z x.
HH HH X Z
T- i:
UJ
z
=>
O 1C
CX Z
UJ <
r>i tx
HH UJ
CX
LU X
in
O
«- •*
co u z: •
IU UJ X >•
X t£ <
(J> » O
HH X • -»
x o «- »
<3 LU tn x
I O Z U. vT
Z O H- ^ «-
Ct LJ «t v» •-•
LU » in
» z x
HH » O t^l -4J » INJ (M « «
LU II M
HH CX <-» • CX O
Ct O X X C u. K
3tU,\ r-«-QHHHH
GO
O
=1
r-
X
<
2.
x
<.
X.
o
*
*
LU
O
z
*
*
*
*
X
CO
in
o tf,
Cj *-
oo O
» UJ »
>O 3 O
UJ HH LU
(- I- K
HH Z HH
CX O CC
302
Kl
ru
cc
ro
u
-------
OC. C- CCOOOOOC-OtjOC/OC-OCCC/OOOOOOOtiOOC^OCjOC.OC.CO
oooo-jncjoooooooaoooooooooooooooocoooooooo
«— ro r» «* f> >v.. f— .'. c c, - P- C' c «~ rvj f < >i u. -o (••• si c<- o T- rj ri ^r LT -c r^ n o o
^^^<>^>j>o^H;r-^N-r-r^r^i^-^-r*r-cccoco3jkOc^at«cbAGot>o-;>o*o. LT-O-O-D-O
t-. r r . r- • f-1 r- , t-, i«, r*> M t-. K r-, i-. K; K M fi i- f« t-i f"> f r«i i« »• t-1 f K^ f r- • •" c i t«, F i r*. r-, r KI -j
-' c ._•> o , j ::< :~> _i • . r - ^ ; c; o c; -'-> --" o o o o -D c: o cr- t.j •;) o J '-.--• '-'• '..' '_. '-^ c"« .- o >.. z> o
ws \ x
<: u. «_>
a. a <
l-l UJ
t/5 3 »
l-l (£ « GC
I • O
O »- 00
«t «j <
Oc UJ t-
UJ CL
> l-« *
«t A *
O /^
O in ro
Z It"-"-!
<3^o i-i la:»
tt LU *>. C N »-• ZO U V X
ui3*r-x z> i-i*r^
>'*»irJ4rttrir* 3« X m »
< \ uJ i-i r- « cc •HCCINJ
u.»-iz«t IS *U5«
O X «-i >• CD OC 2 ••>»"
Zt-CcJu-ccz Dl-t O-*«-
•-H LU uj G < i~< OZ Z cr u.
2 o c c. n ct XZ »-H ^ a.
o
s*l
U,
i
v.
h-*
mj
Ul
CO
«£
C
u.
t-
•-I
OJ
r-
<
o
z
<.
vi
>•
_) »
«t •
UJ
_J
CD
<
*
*
C
1
X
r> o t— "- o: Z OJ£C^^ or *^T-
!- tu a < a.ii »-i o «». Q. >»• t-i O a:
«-OOZUJ»-OK »O» r\i»-o O
u: UJ 1-1 _j UJ UJ o: Z * X tt v~ sr O »j «— u-
XZ_4K-CCWO>-irME 13 I O tr\»<>
Oi-lH-i-iu.i-iu.acxZS:^>%O -< O LU
i-iS:3KcE> o-Kt-iO^*! *»o: xooo _j
XCCOZUJLUUJ •• •• Z V OJ X (M ^ UJ »• LU K ffi
LUUJ>-OZft'r-Ei-«VJ»-(»-«»»illl Z roCt. «*
at-uiZ»*''-iOx2XZ»o:rMr-ir»-rg 1-1 otr-o »-
z uj a. i- i-i (- u- H-»-t«tO»«tO*». I oca: ^ -» t- ^»u> o
OOS^ «1 =1 <^>^»:»^ >— IA « r+ ^«<^^- "^ 3« «hOr% (X 4TI
«_> o u-i uj n o ui ri rx \ »— m»»-> >-i ^> *^ i— i ^ *^ » o LU ir>i-iC7> a »
LUZuJZlXO ix X. in(M«>< »-** w OO<^i^'«^ «Z ro»» O «-
«/j LU u ^> t- I >-i o» ^- ••^aicr^t 1-1 »- x a x u
LUOO OJ O »— O t>- O- t— LL(r-i(\J(Mr\J»-aiOrvJ>-i^i-< O t- O O LU I i-<
liJtt-OSCO:* UJ • • >j Z >- «-> O*O*LUUjO«»<»*' tu »»»• OL
x a. a o T- IMCO «o o ^- LM «t 2 O* fifNju->rsj«»LnoDtDCL NLS >o<^ao o
t- LU u. i-i Z «•'«-•!-<-»<=>< I- >XCcXa:NNCl * LU \ • _j f-i LU LU < z * K»-i-i^i-i^:ccLic>-itoQ:to _j 1-1 LU < a: u oo LU «-
• <: •• Z » O «t U3 »- t- H O » CO O- II II II M 1£ X II II II II I- E l-l 00 II U. O
xxxi-^xz »*j«^Ln«4- oc< >-ia:>^«
rvj ^-_jrvj uia: OTPKO «— oooc2crQ:i».u.2eia:a: LLJCI dCiJ-w-Z z C.
rj \ rj \ Co ^ z LU jt^ttk«<«t oi£H-i^:Ki-ii-(»-i^vOic z LU 3u.>-ii-i»-i «t o
****** LU> **« iu> -».^x "^
LS «t o LO < CVQ »^> t-
* t\J *
• * «- * JE"*" «- *
#£?• * ^»Os O- *
<_». v> l_» CJ UOOfcV <_» VJ VJ
30
-------
OOOC OOOOOOCiOCC-OOOCOCJOOCJOOCJC O C O
OOOCOOOOOOOOOOOOCOOOOOOOOOOOOp
T— (\j f*~i NJ if^ »ij f-^ ti £* r t ^— f vj r*~ -sj i/* ^,~ r— io c C T™ oj f -^j u~t *fs r*- & ij- v_j
^D ^D C3 CD Ct C^ f~~i r~i c> ^- T** T~ ^— T* t~ T* ^~ ^* ^^ (\i f\i ru f>J CNJ f\j CNJ rvj (\j f>j r*)
C2 O O O O ~ ^ O 'J3
Z Z z ir z Z
o: cc oe cc ct cc
a:
a.
O
O
fM
DO
o
X
II
ao OQ
O O
IT* (TV
X C.
IT 7
o
eo
cc
a.
<
cc
II
CO
II II -» -t
o
u»
v- o
OuJ
Kl
K) ro
O O O
O
oc
3 »i-iXO » » II II
X II II II II II M •» -» r-
O
in
GC fO
ID r-
zz
< 1-1
f— t-»fi ii
*'-lZ
O»-iX
-------
December 1980
ADDENDUM/SUPPLEMENTAL INFORMATION
TO THE
INDUSTRIAL SOURCE COMPLEX (ISC) DISPERSION MODEL
The computer code modifications addressed herein are the complete set of
changes made to date to the ISC code as it appears in the Industrial Source
Complex (ISC) Dispersion Model User's Guide, Volume II, EPA-450/4-79-031,
December 1979. All modifications have been incorporated in the ISC code as
part of the UNAMAP-4 available from the National Technical Information Service
(NTIS).
Each item in this report is preceded by the date on which it was issued.
A. Computer Program Modifications
Program modifications are described below. The changes correct errors
in the original source code released in December 1979.
INDUSTRIAL SOURCE COMPLEX (ISCST) PROGRAM
Make the following modifications to Appendix A of the User's Guide and
check your program for conformance.
1. 4/80 - Change the following lines of code in subroutine MODEL to read as
follows:
IYR=YR S0304630
YRS=SINNUM(IYR) S0304640
YRC=COSNUM(IYR) S0304650
IYR=YR S0304800
YRS=SINNUM(IYR) S0304810
YRC=COSNUM(IYR) S0304820
2. 4/80 - Delete subroutine UPWIND from the ISC Short Term model program.
Subroutine UPWIND was originally placed in the ISCST to save computer
time by not calculating concentrations at polar coordinate receptors
upwind of centrally located sources. However, due to an error in
the subroutine, it was incorrectly applied to source/receptor combinations
not covered by the original definition and as a result provided sub-
stantial but erroneous savings in computer run time. Correcting the
subroutine though, has resulted in only a negligible savings in run time.
Since subroutine MODEL provides an in-line upwind check for all
source/receptor combinations and the computer time saved using the
corrected UPWIND subroutine is negligible, the subroutine is being
deleted.
-------
3. 4/80 - Change the following line of code in subroutine MAXOT to read as
follows:
IY=(K-1)/NXPNTS+1 S0500370
4. 12/80 - Change the following lines of code in subroutine MODEL to read:
1010 IF(ISW(4).NE.1.0R.HS+ZS-GRIDZ(IJ).GT.O.O.OR.ITYPE.EQ.2)GO TO 1020 S0304930
IF(ITYPE.GT.O)GO TO 1095 S0305200
IF(VS.LE.O.O)GO TO 1095 S0305220
1095 IF(1SW(4).NE.1.0R.ISW(1).NE.1.0R.NVS.NE.O.OR.ITYPE.EQ.2)GO TO 1100 S0305480
These modifications are required to allow for the treatment of area sources as
flat terrain sources as indicated in the documentation.
5. 12/80 - Insert the following line of code in subroutine MODEL:
IF(l.EQ.l) JDY=JDAY S0301795
12/80 - Delete line S0301930 from subroutine MODEL.
The modifications correct on error in the Julian Day for the selection of
the meteorological card input data/no listing option.
6. 12/80 - Delete line S0306540 from subroutine MODEL.
12/80 - Insert the following line of code in subroutine MODEL:
DO 1290 K=1,NVS S0306515
These modifications are necessary to eliminate the possibility of calculating
an incorrect Vertical Term when ground-level concentrations are calculated
from sources with significant gravitational settling.
7. 12/80 - Change the following lines of code in subroutine MODEL to read:
C TEST FOR WAKE EFFECTS - CALCULATE XPLUME. S0304070
780 DHA=3.*FM*GAMJI*UBARI*SSI S0304110
12/80 - Insert the following line of code in subroutine MODEL: „
tf*
IF(1.570796327*UBARS*SSI.GT.HB+HB)DHA=DHA*SIN(SS*(HB+HB)*UBARI) S0304115
12/80 - Delete line S0304170 from subroutine MODEL.
The modifications are required to correct a wake effects test plume rise
calculation error.
-------
INDUSTRIAL SOURCE COMPLEX (ISCLT) PROGRAM
Make the following program line modifications to Appendix B of the
User's Guide and check your program for conformance.
1. 4/80 - Change the following line of code in subroutine MODEL to read:
IF (RD.GE.XMX) GO TO 1930 S0211230
This change allows sources located 100 m from a receptor to be considered
for contribution to the receptor. This is compatible with the ISC Short
Term Model.
2. 12/80 - Change the following line of code in subroutine MODEL to read:
IF(ZP.LT.ZS.OR.TYPE.EQ.2)ZP=ZS S0211150
This modification is required to allow for the treatment of area sources as
flat terrain sources as indicated in the documentation.
3. 12/80 - Change the following line of code in subroutine MODEL to read:
IF(TYPE.NE.O.OR.JTYPE.EQ.O.OR.XB.GE.HB10)GO TO 1970 S0211740
This modification corrects a mis-coded switch used in the calculation of
o for stacks with wake effects and calculation distances less than 10
building heights.
4, 12/80 - Change the following line of code in subroutine OUTPT to read:
IF(NOCOMB(NG).EQ.1.AND.IDSORC(NGT+1).GE.O)GO TO 260 S0301621
This modification is required to correct an error that allowed the print output
for an individual source to be printed several times.
5. 12/80 - Change the following lines of code in subroutines OUTPT and HEADNG
respectively to read:
COMMON /HEAD/ MSG3(4),MSG4(8),IS,N3,I3,M3,ILN.LSTILN S0300120
COMMON /HEAD/ MSG3(4),MSG4(8),IS,N3,I3,M3,ILN,LSTILN S0400110
B. 4/80 - Computer Underflow in the Short Term and Long Term Model Programs
An underflow interrupt condition occurs on computers whenever an arithmetic
operation between two very small numbers results in a number too small to
be contained in the computer's arithmetic results register. When this
condition arises, computers such as the UNIVAC 1100 series and the CDC 6000
series set the underflow result to zero and continue processing. The ISC
model programs were specifically designed to allow the computer to zero an
underflow condition, therefore an underflow condition is not to be considered
an error within the ISC programs.
-------
However, some IBM computers abort the run when an underflow occurs. In
order to circumvent this problem the user must inform the computer
system that underflows are to be set to zero and program execution to
continue-. For example, the following FORTRAN call may be inserted at the
beginning of either ISC Model program and thus preclude termination of
runs due to underflows on an IBM 360 series computer:
CALL ERRSET (208,256,-1,1,0,208)
To find the appropriate solution to this problem on other computer
systems, the user should contact on-site systems personnel.
C. 4/80 - Receptor Heights Versus Source Elevations
The ISC model allows each receptor to be situated at any height (z).
Both the ISCST and ISCLT programs terminate, however, if the receptor
elevation is above the lowest source height. In addition, if any receptor
is located below the source base height, the ISC programs automatically
reset the receptor height to the source base height (thus, flat terrain
is assigned for that source receptor combination). This will cause a
receptor to be located at different vertical points if two or more
sources affect the same receptor, and are at different base elevations.
To avoid this problem all sources should be located at approximately the
same base elevation, thus maintaining a consistent correction among
receptors that are below the plane of the sources.
The rationale for the above suggestion stems from the applicability
of the ISC model primarily to very localized industrial sources and
complexes, i.e., the sources are not located in radically different
locations or terrain. The model can be used for widely dispersed
sources, but the user's decision to assume similar source elevations
should be an explicit one.
D. Text Corrections to the User's Guide
Volume I
1. 4/80 - All discussion and figures pertaining to the subroutine UPWIND are
to be disregarded.
2. 12/80 - Page 3-6, Receptor Grid System Option. Delete sentence 3 and
replace it with:
3. 12/80 - Page 2-24. The exponential term in Equation (2-12) should be:
exp
r. i m2"
2 a
I \ y I J
Additionally, a "3" value will automatically generate a grid system using
the Cartesian coordinate system and a "4" value will automatically generate
the polar coordinate direction radials with user-defined starting locations
and spacing distances.
-------
A. 12/80 - Page 3-42, Table 3-4. Line 10 in the Description which pertains to
ISW(2), should be changed such that the line reads:
4= program generates polar coordinate grid direction radials
5. 4/80 - Page 3-60. Add to equation (3-1) the term "+D".
6. 4/80 - Page 3-61. To the list of variables and their definitions add:
D=NPNTS if ISW(4) equals "I" in the first card of Card Group (2); otherwise
A equals "0".
7. 12/80 - Page 4-10. Equation (4-1) should be changed to read:
J = 300
and
(E-(N + N +• 2-N ) - K - L - M) , ,. ,.
T = x y xy" ( (4-1)
(N -(N -N + N ))
se x y xy
8. 12/80 - Page 4-11, lines 9 and 11. References to "N " should be changed
to "N ". s
se
9. 12/80 - Page 4-11. To the list of variables and their definitions add:
' 0 ; if ISW(4)=1 & ISW(11)=2 or
if ISW(7)=1 or NSEASN=1 or NGROUP=0
M =
N -N +N ; if ISW(4)=0 or ISW(11)^2 &
x y xy
if ISW(7)^1 & NGROUP^O & NSEASN^l
10. 12/80 - Page 4-16. Equation (4-2) should be changed to read the same as
equation (4-1).
11. 12/80 - Page 4-17, equation (4-3). Reference to "N " should be changed to
"N ". S
se
12. 12/80 - Page 4-77. Equation (4-4) should be changed to read:
J = 300
and
(4-4)
E - (NXPNTS+NYPNTS+2*NXWYPT) - K - L - M
NSEASN*(NXPNTS*NYPNTS+NXWYPT)
-------
13. 12/80 - Page 4-78. To the list of variables and their definitions under
"Condition c" add:
f
0 ; if ISW(4)=1 & ISW(11)=2 or
if ISW(7)=1 or NSEASN=1 or NGROUP=0
M =/
NXPNTS*NYPNTS+NXWYPT; if ISW(4)=0 or ISW(11)^2 &
if ISW(7)=1 & NSEASN=1 & NGROUP=0
14. 12/80 - Page 4-79. Equation (4-5) should be changed to read the same as
equation (4-4).
15. 12/80 - Page 4-80. Delete variable "L" and replace it with:
"where L and M are the same as under Equation (4-4)".
Volume II
1. 4/80 - All discussion and figures pertaining to subroutine UPWIND are to be
disregarded.
2. 4/80 - Page C-4, Figure C-l. The zero in card column 59 in the first card
should be changed to a "6" such that the data value reads "860".
3. 4/80 - Page C-8, Figure C-l. The data in card columns 8 and 9 are to be
transferred to columns 9 and 10 respectively with column 8 becoming blank.
The change should be made for each of the 16 cards illustrated.
4. 4/80 - Page D-8, Figure D-l. The zero in card column 25 on the fourth last
card should be changed to a "6" such that the data value reads ".0076794".
5. 4/80 - Page C-19, paragraph 1. Delete sentence 3 and replace it with:
Also, because ISW(15), ISW(17), and ISW(18) equal "1" and ISW(4) equals "0",
variables A,B,C, and D become 425 times 5(or 2125), 4 times 1 times 425
times 5(or 8500), 201 times 1 times 5(or 1005), and 0, respectively, according
to their definitions given in equation (3-1).
6. 4/80 - Page C-19. Add to equation (C-2):
variable version - "+D"
numerical interpretation - "+0"
7. 12/80 - Page D-24. The numerical solution in equation (D-3) should be changed
to read:
40000 >_ 19+19+2-1 + (2-1+0) • (19-19+1)
> 764
------- |