vvEPA
United States
Environmental Protection
Agency
Region H Office
26 Federal Plaza
New York, N.Y. 10007
EPA 902/4-79-002
April 1979
Air
User's Guide
to the
NYSDOT Air Quality System
-------
GCA-TR-78-67-G
USER'S GUIDE TO THE MODIFIED
NYSDOT AIR QUALITY SYSTEM
Final Report
Prepared for
U.S. ENVIRONMENTAL PROTECTION AGENCY
Research Triangle Park
North Carolina 27711
Contract No. 68-02-2539
Task Order No. 7
Project Officer
William Beach
Prepared by
Dominic Caracciolo, Jr.
April 1979
GCA CORPORATION
GCA/TECHNOLOGY DIVISION
Bedford, Massachusetts
-------
DISCLAIMER
This Final Report was furnished to the U.S. Environmental Protection Agency
by GCA Corporation, GCA/Technology Division, Burlington Road, Bedford,
Massachusetts 01730, in fulfillment of Contract No. 68-02-2539, Task Order No. 7.
The opinions, findings, and conclusions expressed are those of the authors and
not necessarily those of the Environmental Protection Agency. Mention of company
or product names is not to be considered as an endorsement by the Environmental
Protection Agency.
-------
ABSTRACT
This manual documents the modifications to the New York State Department
of Transportation Air Quality System made by GCA/Technology Division. These
modifications principally involved the inclusion of the most recent motor
vehicle emission factors, as represented by the MOBILE1 program.
111
-------
IV
-------
CONTENTS
Abstract iii
Figures vi
Tables „ vi
1. Program Description 1
2. Input/Output Requirements and Options 4
Input Flags 4
3. Software Documentation 22
GCAINP 22
GCAOUT 31
Appendix
FORTRAN Source Code 32
-------
FIGURES
Number Pa.g£
70
1 Sample input data ........ •
2 Sample input data description ......... 21
3 Job file structure. • • 23
4 Flow-chart for main program ....... 24
5 Flow-chart for subroutine GCAINP. .... • 25
6 Flow-chart for subroutine GCAOUT. .............. 29
TABLES
1 Factored
2 Input Fl;
3 Emission File
4 Land Use File
5 SYMAP File
6 Required 0
7 Program De
ik File
Record Description ..........
Le
Le .
2-Time Data Description. .....
lults
...... 2
...... 5
8
11
1 3
... 15
17
VI
-------
SECTION 1
PROGRAM DESCRIPTION
The purpose of this manual is to familiarize the user with the modifica-
tions made to the NYSDOT Air Quality System and describe input requirements
and options and their effect on program execution and output.
The modifications implemented involve changes to MOBILE1 so it will
accept input data from the NYSDOT Air Quality System and provide output in a
form similar to the NYSDOT Air Quality System's ALGOL-based emissions program
MOBILE1 is a sot of FORTRAN-based routines designed to aid in the analy-
sis of emissions from motor vehicles. It computes emission factors for
hydrocarbons (HC), Carbon Monoxide (CO), and Nitrogen Oxides (NOX) based on
vehicle type, region, calendar year, vehicle inspection and maintenance
parameters, speed, temperature, and various other parameters. For a complete
list consult the MOBILE1 User's Guide."' This manual assumes the reader is
familiar with the MOBILE1 program. Moreover, users of the NYSDOT Air Quality
System must be well acquainted with motor vehicle emission factor concepts to
ensure intelligent use of the System.
In order to modify MOBILE1 to accept input data from the NYSDOT Air
Quality System (a Factored Link File - for a description see Table 1), it
was necessary to construct a new input routine (subroutine GCAINP) which
would replace the standard one (subroutine INPUT). The standard MOBILEl out-
put routine (subroutine OUTPUT) has also been replaced by subroutine GCAOUT.
The MOBILE 1 main program has been replaced by a program which first calls
CCAINP and then calls GCAOUT. A new subroutine has been added (EFPMPB) to
compute emission factors for lead (Pb) and particulate matter (PM).
GCAINP is responsible for reading the input flags and one-time data,
initializing appropriate arrays, reading the Factored Link File, calculating
emissions, printing the emissions for each link and optionally writing an
Emissions file.
GCAOUT controls the printing of emissions by zone and/or district, and
the optional input of a Land Use file. If the Land Use file is requested,
GCAOUT calculates zonal emission densities, and at the user's option, prints
them. The user also has the option of requesting a SYMAP file to be output
Lo tape .
Gutliman, Lewis E., User's Guide to MOBILEl: Mobile Source Emissions Model,
(U.S. Environmental Protection Agency, Office of Air, Noise and Radiation.
Washington, D.C. August 1978.)
-------
TABLE 1. FACTORED LINK FILE
Klo-
mpnt
Nn.
1
2
3
4
5
(,
7
8
9
10
1 1
12
13
I'l
15
16
1 7
18
19
20 '
2 1
Ti
'! 3
?4
25
26
27
28
29
JO
31
12
1 )
»/,
15
16
17
IH
iy
Dose r i p t ion
Kimc t ion rode
Jur i scl LC L j on code
DIsLrlct
A Node
II Node
Federal Aid Code
Link length
l.lnk typt-
W.iy code
Area type
I'ol 1 L leal unit code
K Hi- II)
I'Lanning link code
l-ll li-r
/oiled speed
llnll'widtli (Tc)
Source code
Crnde code
Capac 1 Ly
Average speed
Vul nine/capacity rut Jo
'Co t .'i ] vo 1 unto
Vo 1 IIITIC , vrh Jc 1 e type ]
Volnini1, vt'hiclo type 2
Volume, vehicle type '3
Volume, vehicle type A
Vo 1 uitii- , vehicle type 5
Volume, vehicle type 6
McGinn tnj-, hour
Length ol time period
Avernpe speed
Vo 1 iinic/c<'ip:tc 1 ty ratio
Tol.il volume
Volume, vehicle type 1
Vo 1 ume , veh i c 1 e type 2
Vo 1 ume , veh Lc 1 e type 3
Volume, vehicle type 4
Volume, vehicle type 5
Volume, vehicle type 6
Loca t i on
1
2
•J - It
5-9
10 - H
15
Ih - 18
19
20
21
22 - 25
26 - 29
30 - 33
34 - W
45 - 46
47 - 48
49
50
51 - 56
57 - 58
59 - 62
63 - 68
69 - 7 3
74 - 78
79 - 83
84 - 87
88 - 91
92 - 95
96 - 97
98 - 99
100-101
102-107
108-113
114-1 IH
119-123
124-128
129-132
1 33-136
1 J7-I40
Size
1
1
2
5
5
1
3
1
1
1
4
4
4
11
2
2
1
1
6
2
4
6
5
5
5
4
4
4
2
2
2
6
6
5
5
5
4
4
4
Type
N
N
N
N
N
N
N
N
N
N
N
A
N
A
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
Picture
9
9
99
9(5)
9(5)
9
9v99
9
9
9
9(4)
X(4)
9(4)
X(ll)
99
99
9
9
9(6)
99
9(3)v9
9(6)
9(5)
9(5)
9(5)
9(4)
9(4)
9(4)
99
99
99
9(5)v9
9(6)
9(5)
9(5)
9(5)
9(4)
9(4)
9(4)
Comment fl
Zone stored in
locations 5-7
24 hour
24 hour
24 hour
24 hour
24 hour
24 hour
24 hour
24 hour
24 hour
Time period No.
Time period No.
Time period No.
Time period No.
Time period No.
Time period No.
Time period No.
Time period No.
Time period No.
Time period No.
Time period No.
1
1
1
1
1
1
1
1
1
1
L
(contInued)
-------
TABLE 1 (continued)
II
1)
1
:>
i
'i
K
•1
1)
1
•t
1
he
I!'1)', i Illl i
I.ellC.I ll
Avel ;l)'r
Volume/
Tnl.il v
Vnluine,
V..|iime,
Volume,
Vo 1 time ,
Vol nine.
MT i pt ion
it}1, hour
i>l 1 i me per i od
speed
r.ip.ie i 1 v r.,1 io
t • 1 nine
vehicle t ypo 1
veil i f 1 e 1 ype J
velili le type 4
Veil i r 1 e type ')
veil i r 1 e type 6
He}- i mi i ii}', hour
I.eli)..! ll
Ave l .Ij-.i •
V..I e/
ol I hue per iod
speed
, .ipr i t v r.'il i o
'» 1 < '1 ;i 1 vo 1 ume
h
1
M
'.'
II
1
.'
Vol uiiii'.
Vol •,
Volume,
Volume,
Volume,
Volume,
l-'i 1 lei
veil i e 1 e L ype 1
veil i r 1 1 type 2
veil 1 e 1 e t y pe 3
veli 1 c 1 e type 4
veh i r 1 1 Ly pe 5
veh i c 1 < t ype f>
Location
141-142
143-144
I4'>-14(.
14 /-I 52
IVI-I r>8
I1)')- If) 3
K./. 1 (,U
J't ~~ 1 \>O
174-177
178-181
182-185
186-187
188-189
1 90- 1 9 1
142-197
198-203
L'04-208
209-21 3
:'l/,-2IH
1' 19-222
2 2 3- 2 2 f)
227-2311
2 11-240
Size
2
2
2
6
fa
5
c;
J
5
4
4
4
2
2
2
6
6
5
5
5
4
4
4
10
Typo
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
A
Picture
99
99
99
9(5)v9
9(6)
9(5)
9(5)
9(4)
9(4)
9(4)
99
99
99
9(5)v9
9(6)
9(5)
9(5)
9(5)
9(4)
9(4)
9(4)
X(10)
Comments
T i me
Time
T i me
Time
Time
Time
T 1 mil
l ime
Time
Time
Time
Time
Time
Time
Time
Time
Time
Time
Time
Time
Time
Time
Time
period
period
period
per Iod
period
period
pe r 1 od
period
period
period
period
period
period
period
period
period
period
period
period
period
period
period
No.
No.
No.
No.
No.
No.
No .
No.
No.
No.
No.
No.
No.
No.
No.
No.
No.
No.
No.
No.
No.
No.
2
2
2
2
2
2
2
2
2
2
2
2
3
3
3
3
3
3
3
3
3
3
-------
SECTION 2
INPUT/OUTPUT REQUIREMENTS AND OPTIONS
The input data requested by subroutines GCAINP and GCAOUT is either re-
quired or optional. The required input files consist of: (a) a card deck
containing the Input Flags and Required One-Time Data and (b) a tape or disk
file which contains the Factored Link data. The card deck is assigned to
logical unit 5, the Factored Link file is assigned to unit 4. (See Table 1
for a description of the Factored Link file.)
The optional input files consist of the following: (a) a card deck con-
taining Optional One-Time Data (placed between the input flag record and
Required One-Time Data - consult Section 2, Optional One-Time Data, for a
complete description of required order), and (b) a tape, disk or card file
containing Land Use data. The unit number assigned to the Land Use file is
determined by the value for input variable OUTFL5.
The only required output file is for printed output. It is permanently
assigned to logical unit number 6.
Optional output files consist of the following: (a) a tape or disk file
containing Emissions data and (b) a tape or disk file containing SYMAP data.
The unit number assigned to the Emissions file is determined by the value of
input variable OUTFL2; that for the SYMAP file from the value of input
variable OUTFL9.
INPUT FLAGS
The Input Flags are coded on a single card which has a format of 1711.
This card is the very first record to be read by subroutine GCAINP. A tabu-
lar description of the Input Flag record can be found in Table 2.
® Column No. 1: MYMRFG
MYMRFG is a MOBILE1 flag which indicates whether the
user is entering vehicle registration and/or annual
mileage accrual distributions or is using the default
national average distributions. If the user wishes to
supply either or both distributions, additional input
is required. The Optional One-Time Data section should
be consulted for the placement of the additional data.
-------
TABLE 2. INPUT FLAG RECORD DESCRIPTION
Col.
No.
I
ft
B
9
10
12
U
15
Variable
Name
•
ALTPLt;
AUIFLC
TPJCFU.
OUTFL1
OUTFI.2
OUTFI.'J
OUT Ft 5
OIITFI.7
Description
registration and/or mileage
Flag Tor optional modi feat Ion of model
mnlnLrnnncp program In cnlculnt iona
linn Tor optional use of four Hpi-clal l.l>V
(h) vehicle loading (c) trailer towing
(il) humidity (NOX only)
Flag for optional i|se of HDV power /weight
i-orracl Ion fnrtor
Klag to control pollutants to be printed
KmlaelmiH rile
printed fur each zone
rl
lirJnlt-d tt>r fitch district
Flag to i<>ntrol optional Input of Land
HBP flU-
Zoniil EmlHflion Dcnsltiea are to be
printed
X.onnl I'.mlHHfnn DcnulticB
SYHAI' flic
Hiibroni Im i.CAINP opcratOB
Value
1
3
Ottier
0
Other
1
Other
0
1
Other
0
1
Oilier
]
1
Other
0
1
2
3
4
5
6
7
Other
0
Non-0
0
1
2
Other
1
Other
0
Non-0
1
Other
0
1
2
3
it
5
6
7
Other
0
I
2
3
Other
Non-0
Action
Uaer auppllca mlluaitc dint.
User supplies both diet.
Error
Use default factnrn
Error
.
Include I/M
Error
Four factor, not Included
Error
Do, not uae 1IDV fnrtor
Uaa HDV factor
Rrror
Nonme thane I1C
Crana
Error
CO
NOX
HC
Pb
PM
CO, NOX, HC
All
Error
Emiaalone file not output
Emissions file written to logical
unit number OUTFL2
Emlaelono not printed by zone
period only
Emissions printed for all tlma
periods
Error
Emissions printed for 2^-hour time
period only
periods
Error
Land Use file not Input
Land Use file read from logical unit
number OUTFL5
Zonal Emission Densities printed for
24-hour tfime period only
all time periods
Error
All pollutants (default)
CO
NOX
HC
Pb
PM
CO, NOX, HC
All pollutants
Error
Square feet (default)
Square miles
Acres
Square feet
Error
number OUTFL9
Npt'ctl InirumuntH 1 rom *i to 60 mpli;
Tnblt' look-up scheme used aa each
(No. of calls to EFCALX-152).
Emission factors calculated as each
link la read from Factored Link file.
(No. of calla to EFCALX-A*No. of
-------
Column No. 2: ALTFLG
ALTFLG is a MOBILE1 flag which indicates whether the user
desires to modify any of the emission factors used in the
program. At the current time MOBILE1 does not support a
value of ALTFLG=1. Appendix C of the MOBILEl User's Guide
should be consulted to modify MOBILEl and GCAINP if user-
supplied emission factors are to be used.
Column No. 3: IMFLAG
IMFLAG is a MOBILEl flag which indicates whether the
computed emission factors should include credits for
an Inspection/Maintenance program. If I/M credits are
to be applied, additional input is required. Consult
the Optional One-Time Data section for the placement
of the additional data.
Column No. 4: ALHFLG
ALHFLG is a MOBILEl flag which indicates whether the
composite emission factor computation should include
the effects of four light-duty vehicle correction fac-
tors: (a) air conditioning, (b) vehicle loading,
(c) trailer towing and (d) humidity (NOX only). If
these correction factors are to be included, additional
input data is required. Consult the Optional One-Time
Data section for the placement of the additional data.
Column No. 5: TRKFLG
TRKFLG is a MOBILEl flag which indicates whether the
composite emission factor computation for heavy-duty
gasoline and diesel-powered vehicles should include the
effect of the "additional" correction factor based on
the power/weight ratios of these vehicles. If these
correction factors are to be included, additional input
data is required. Consult the Optional One-Time Data
section for the placement of the additional data.
Column No. 6: NMHFLG
NMHFLG is a MOBILEl flag which indicates whether the
user desires total or nonmethane HC emissions to be
computed.
Column No. 7: OUTFLO
OUTFLO is used to select whether the emissions are to
be computed in pounds (OUTFLO=0 or blank) or grams
(OUTFLO=1).. Any other value results in a fatal error
aborting the run.
-------
Column No. 8: OUTFL1
OUTFL1 controls which pollutants emissions are to be printed
as each link is read from the Factored Link file. If OUTFL1
equals zero or is left blank, emissions will not be printed
for each link; if OUTFL1=1, CO emissions will be printed; if
OUTFL1=2, NOX emissions will be printed; if OUTFL1=3, HC
emissions will be printed; of OUTFL1=4, Pb emissions will
be printed; if OUTFL1=5, PM emissions will be printed, if
OUTFL1=6, CO, NOX and HC emissions will be printed; if OUTFL1=7,
emissions for all five pollutants will be printed. Any other
value results in a fatal error, aborting the run.
Column No. 9: OUTFL 2
OUTFL2 controls the optional output of an Emissions file.
If OUTFL2 equals zero or is left blank, no Emissions file
will be written; if OUTFL2 is nonzero, an Emission file will
be written to the logical unit number determined by the value
of OUTFL2. Note that a check is made to insure that the unit
number for the Emissions file differs from that assigned to the
Factored Link File (Unit 4). A further check is made to insure
that the unit number assigned to the Emissions file also differs
from that assigned to the optionally input Land Use file (which
is defined by the value of OUTFL5). These tests are performed
to protect the user from possibly writing upon an input file.
See Table 3 for a description of the record format of the
Emissions file.
Column No. 10: OUTFL3
OUTFL3 controls which time periods are to be selected when
emissions are printed for each zone. If OUTFL3 is zero or
blank, emissions for each zone will not be printed. If
OUTFL3=1, emissions will be printed for each zone for the
24-hour time period only. If OUTFL3=2, emissions will be
printed for each zone for all four time periods. Any other
value results in a fatal error, aborting the run. Note that
emissions for all five pollutants will be printed; that is,
the user has no control over the selection of pollutants.
Column No. 11: OUTFL4
OUTFL4 controls which time periods are to be selected when
emissions are printed for each district. If OUTFL4 is zero
or blank, emissions for each district will not be printed.
If OUTFL4=1, emissions will be printed for each district for
the 24-hour time period only. If OUTFL4=2, emissions will be
printed for each district for all four time periods. Any other
value results in a fatal error, aborting the run. As with
OUTFL3, emissions for all five pollutants will be printed.
-------
TABLE 3. EMISSIONS FILE
KI. Mill Ml I
lle.Me r l pt I on
1 I'unc t ion code
7 Jurisdiction code
) District
/i A node
5 H node
6 l-'ederal aid code
7 Link length
8 Link type
4 Way code
10 ' Al ea type
11 Political unit code
12 Kile II)
13 1' 1 .inn i ii}.', 1 i nk code
l/i Killer
1 ') /.oned speed
16 Hal 1 width ( ft)
1 / Sou rce code
IK Crade code
1 ') Capac i l y
70 Average' speed
21 Volume/capacity ratio
22 Volume
2 3 lln i t code
74 2/i hour
'.' 'i T i me per iod No . 1
7I> Time period No. 2
7 7 T i me per iod No . 3
78 7/i hour
7l) T i me per i od No . 1
10 T i me pe r i od No . 2
t 1 T l me pe r 1 od No . i
17 '.'-'t hour
1 1 T Hue per iod No . 1
l/i T i me per iod No . 2
1 '> 1 i me per i od No . I
Locat ion
1
2
3-4
5-9
10-14
15
16-18
19
20
21
22-25
26-29
30-33
34-44
45-46
47-48
49
50
51-56
57-58
59-62
63-68
69
70-77
78-85
86-93
94-101
102-104
110-117
118-125
126-13'i
134-141
142-149
150-157
158-165
S i ze
1
1
2
5
5
1
3
1
1
1
4
4
4
11
2
2
1
1
6
2
4
6
1
8
8
8
8
8
8
8
8
8
8
8
8
Type
N
N
N
N
N
N
N
N
N
N
N
A
N
A
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
Picture Comments
9
9
99
9(5)
9(5)
9
9v99
9
9
9
9(4)
X(4)
9(4)
X(ll)
99
99
9
9
9(6)
99
99v99
9(6)
9 ^
9(8) I
/ NO emissions
9(8) 1
9(8)^1
9(8) \
9(8) 1
\ HC emissions
9(8) t
9(8) 1
9(8) \
9(8) I
\ CO emissions
9(8) /
9(8) j
. — _
(continued)
-------
TABLE 3 (continued).
KlemeiH
1 )i ' s c r i
Nn.
. lf> ??\ hour
i / T i mi* per i oil
IK Ti me per i od
I1' TiiiH' per iod
'i(d !'A lion r
A 1 T i in*' per iod
4;! Time per i IK!
'i 1 T iiin1 per ind
•'.'* Killer
ir ion
No. 1
No . :>.
No. 3
No. 1
No. 2
No. 3
Location
1G6-173
174-1H1
182-189
190-197
198-205
206-213
214-221
222-229
230-240
Si ze
8
8
8
8
8
8
8
8
11
Type
N
N
N
N
N
N
N
N
A
Pic turd
9(8)
9(8)
9(8)
9(8)
9(8)
9(8) 1
9(8) |
9(8) t
X(ll)
Commcn ts
PM emissions
Pb emissions
-------
Column No. 12: OUTFL5
OUTFL5 controls the optional input of a Land Use Flle- ^
OUTFL5 is zero or blank, no Land Use file will be read. It
OUTFL5 is nonzero, the Land Use file will be read from the
logical unit number determined by the value of OUTFL5. A
check is made to insure that the unit number for the Land
Use file differs from those assigned to the optionally output
Emissions file (OUTFL2) and the optionally output SYMAP file.
The tests are performed to protect the user from possibly
writing upon an input file. For a description of the record
format of the Land Use file, see Table 4.
Column 13: OUTFL6
OUTFL6 controls the time periods for which Zonal Emission Den-
sities are to be printed. The densities are computed by
dividing the zonal emissions by the total land area for that
zone. The total land area is computed from data input in the
Land Use file. If OUTFL6 is zero or blank, Zonal Emission
Densities will not be printed. If OUTFL6=1, Zonal Emission
Densities will be printed for the 24-hour time period only.
If OUTFL6=2, densities for all four time periods are printed.
Any other value results in a fatal error, aborting the run.
Column 14: OUTFL7
OUTFL7 controls which pollutants are to be selected when Zonal
Emission Densities are printed. If OUTFL7 is left blank or
set to zero, the default of all pollutants is assumed. If
OUTFL7=1, CO is printed; if OUTFL7=2, NOX is printed; if
OUTFL7=3, HC is printed; if OUTFL7=4, Pb is printed; if
OUTFL7=5, PM is printed; if OUTFL7=6, CO, NOX and HC densi-
ties are printed; if OUTFL7=7, desnities for all five pollu-
tants are printed. Any other value results in a fatal error,
aborting the run.
Column 15: OUTFL8
OUTFL8 controls the units of area in which Zonal Emission
Densities are to be computed. If OUTFL8 is left blank or
set to zero, the default of square feet is assumed. If
OUTFL8=1, densities are in grams or pounds per square mile.
If OUTFL8=2, densities are in grams or pounds per acre. If
OUTFL8=3, densities are in grams or pounds per square foot.
Any other value results in a fatal error, aborting the run.
Note that the choice of grams or pounds is determined by the
value of OUTFLO.
10
-------
TABLE 4. LAND USE FILE
Element
„ Description
1 District
2 Zone
3 Vacant land
4 Residential land
5 Nonresidential land
6 Alternate passive land
7 Filler
Location
1-2
3-5
6-13
14-21
22-29
30-37
38-80
Size
2
3
8
8
8
8
43
Type
N
N
N
N
N
N
A
Picture
99
999
99v9(6)
99v9(6)
99v9(6)
99v9(6)
X(43)
Comments
1000 's
1000 's
1000 's
1000 's
of ft2
of ft2
of ft2
of ft2
11
-------
9 Column 16: OUTFL9
OUTFL9 controls the optional output of a SYMAP file. If OUTFL9
is zero or blank, no SYMAP file will be written; if OUTFL9 is
nonzero, a SYMAP file will be written to the logical unit number
determined by the value of OUTFL9. As with the other output files,
a check is made to be sure that the unit number assigned to the
SYMAP file differs from those assigned to the Factored Link file
(Unit 4) and the optional Land Use file (OUTFL5). For a descrip-
tion of the record format of the SYMAP file, see Table 5.
* Column 17: MODE
The MODE flag determines the mode of execution in which GCAINP
operates. If MODE is blank or zero, GCAINP will initially set
up an array of emission factors based on 3 mph speed increments
between 5 and 60 mph (i.e., 5, 8, 11, ... 59 mph). As each
link is processed, the emission factor for the average speed
for that link will be 'looked up' in this array and emissions
computed. Note that average speeds of 7, 8 and 9 mph will all
use the emission factor for 8 mph; 10, 11 and 12 mph use that
for 11, etc. This mode of execution results in some inaccuracy
but faster execution (particularly for large Factored Link
files) since MOBILE1 subroutine EFCALX is called 152 times
during initialization. EFCALX is used to compute emission
factors for CO, NOX, and HC.
If HOLE is nonzero, EFCALX is called four times per link (once
for each of the four time periods). The computed emissions are
more accurate than those calculated when MODE=0, but execution
time will be greatly increased. This mode of execution is
suggested only if accuracy in the calculated emissions is
required or if a small Factored Link file is being read.
Required One-Time Data
The Required One-Time Data consists of a card deck containing six records. A
tabular description can be found in Table 6. The description of each card follows,
© Card No. 1: Calendar Year
Card 1 contains the last two digits of the calendar year for
which emissions are to be calculated. Valid calendar years
are 70 through 99. A value less than 70 will be flagged as
a fatal error causing the run to be terminated. The calendar
year should be coded in columns 1 and 2.
® Cards 2 through 6
Cards 2 through 6 have the same format of 4F5.0. Columns 1
through 5 contain data for the 24-hour time period; columns
6 through 10 contain data for time period 1; columns 11
12
-------
TABLE 5. SYMAP FILE
dill Ml I
NO.
I
:>
3
It
')
()
/
H
')
10
1 1
17
i i
14
IS
id
i /
IM
pi
70
71
77
7 1
7/1
:>':
'.'it
•} 1
7H
7')
10
H
'17
)3
'i/i
15
if)
Dc-scr i
pcion
'/.one
NO
HC
CO
I'M
I'll
NO
IK:
CO
I'M
I'h
NO
uc
CO
I'M
I'h
NO
IIC
CO
I'M
I'h
NO
IIC
CO
I'M
I'll
NO
IIC
CO
I'M
I'l)
. NO
IIC
CO
I'M
I'l)
24
24
74
7.4
7.4
1
I
1
1
•1
v
V
2
2
?
'\
1
!
•j
3
74
24
24 .
24
74
1
1
1
1
1
. 2
2
2
2
2
Location Size Type Picture Comments
1-S 8 N 9(8)
9-19 11
7j-:i0
11-41
42-S2
V.l-6'3
l>4-74
7ri-8r>
Hf)-')(>
•)'/-l07
108-118
1 1()-12')
1 10-140
141-lr,l
lr)2-ir,2
H) '1-173
1 74-184
IK5-195
196-206
207-217
218-228
2 '10-240
241-251
252-262
2(i "5-273
774-284
7 H 5-2 9 5
7')6-'106
'107-'il7
")18-"528
179-T19
'140-'! 50
)r)l-'!61
')62-'372
'17 3-383
384-394
f
9(10)v9 Emissions in grams
/
i
i
Emissions in pounds
I \
(
-
(continued).
13
-------
TABLE 5 (continued)
•', lemf nt
No .
37
18
"19
40
VI
4?
V!
44
45
46
4 7
4H
49
50
51
5;'
53
54
5 5
56'
5/
58
V.i
60
dl
i>2
Description
NO 3
MC 3
CO 3
I'M 3
Pb 3
NO 74
11C 24
CO
PM
Pb
NO
IK:
CO
I'M
Pb
NO
nc;
CO
I'M
Pb
NO
IK:
CO
I'M
Pb
!•• i 1 ) f. r
Location Size Type Picture Comments
395-405 11 N 9(10)v9 Emissions in pounds
406-426
41.7-427
428-438
439-449
450-460
461-471
472-482
483-493
494-504
505-515
516-526
527-537
538-548
549-559
560-570
571-581
582-592
593-603
604-614
615-625
626-f,36
637-647
648-658
659-669 \
\
' \
V
Emissions in pounds or grams/Ml2
(at user's option)
\
670 1 A X
14
-------
TABLE 6, REQUIRED ONE-TIME DATA DESCRIPTION
o
a
• «! Columns Inscription
1 1. - !! Calendar year
1 - !> Freeway hot-starts
24-hour time period
6 - in Freeway hot-starts
2 time period .:$1
11 -15 Freeway hot-starts
time period #2
16 - 20 Freeway hot starts
time period #3
1 - 5 Nonfreeway hot-starts
24-hour time period
6 - .10 Non freeway hot-starts
3 time period //I
II - 15 Nonfreeway hot-starts
time period 112
16 - 20 Nonfreeway hot-starts
time period //3
1 - 'j Freeway cold-starts
24-hour time period
6 - 10 Freeway cold-starts
time period #1
LI - 15 Freeway cold-starts
time period #2
16 - 20 Freeway cold-starts
time period #3
1-5 Nonfreeway cold-starts
24-hour time period
6 - 10 Nonfreeway cold-starts
5 time period #1
1.1 - 15 Nonfreeway cold-starts
time period #2
If/ - 20 Nonfreeway cold-starts
time period //3
r ,
24-hour time period
6 - 10 Ambient temperature
(> time period //I.
II, - 1 'i Ambient temperature
time period //2
!(>.-• 20 Ambient temperature
time period //3
I'oniuit
12
F5.0
F5.0
F5.0
F5.0
F5.0
F5.0
F5.Q.
F5.0
F5.0
F5.0
F5.0
F5.0
F5.Q
F5.0
F5.0
F5.0
F5.0
F5.0
F5.0
('.(11111110111 S
70-CY-91)
Percent
Percent
Percent.
Percent
Percent
Percent
Percent
Percent
Percent
Percent
Percent
Percent
Percent
Percent
Percent
Percent
O-Tp-llO
Degrees F
O^T^llO
Degrees F
Degrees F
0-Tj^llO
15
-------
through 15 contain data for time period 2, and columns 16
through 20 contain data for time period 3,
• Card No. 2: Freeway Hot-Start Data
Card 2 contains the percentage of freeway hot-starts for
each of the four time periods.
e Card No. 3: Nonfreeway Hot-Start Data
Card 3 contains the percentage of nonfreeway hot-starts
for each of the four time periods.
e Card No. 4: Freeway Cold-Start Data
Card 4 contains the percentage of freeway cold-starts
for each of the four time periods.
® Card No. 5: Nonfreeway Cold-Start Data
Card 5 contains the percentage of nonfreeway cold-starts
for each of the four time periods.
« Card No. 6: Temperature
Card 6 contains the ambient temperature (°F) for each of the
four time periods. Valid temperatures range from 0 to 110°F.
Any temperature not in this range will be flagged as invalid
and the run will be terminated.
Optional One-Time Data
The optional One-Time Data should be included in the input deck if any
of the six MOBILE1 flags are turned on. This optional data should be placed
between the Input Flag record and the Required One-Time Data records. The
order of the cards within the Optional One-Time Data deck is the same order
in which the corresponding MOBILE1 flags appear on the Input Flag record.
The flag names used are the same as those found in the MOBILE1 program so
a description of the input data format will not be given in this guide. The
MOBILE1 User's Guide should be consulted for the correct format.
Program Defaults
The original MOBILE1 input subroutine allowed the user more flexibility
in setting input flags than does the current input routine GCAINP. The func-
tions performed by these flags deal "with input formats, output formats and a
few" internal modes of operation and are not required in subroutine GCAINP.
They have been assigned default values in BLOCK DATA to insure that they have
a defined value during execution of the program as they may be encountered in
any of the subroutines containing block /FLGCOM/. A description of the
default values may be found in Table 7.
16
-------
TABLE 7. PROGRAM DEFAULTS
Variable
Original purpose
Default
value
Reason for default
INFLAG Select list-directed or formatted 0
read
SP3FLG Select 1 or 3 input speeds 0
VMFLAG Select optional use of user VMT 1
mix
IFORM Select output format 1
PRTFLG Select which factors are included 0
in output
ICEVFG Controls listing of evaporative 1
HC
1DLFLG Controls output of idle factors 1
QNFFLC Select optional unformatted output 0
IREJN Select low altitude, California 1
or high altitude
Not used
One speed read per time
period per link
Supplied in factored
link file
Not usedv
HC, CO, NOX included
Evaporative HC computed
No idle output
"ff
Not used
Low altitude
The original purpose of the flag was to control input or output formats.
It is not used in any of the other subroutines, but a default value has been
assigned for consistency.
17
-------
Summary of Fatal Errors
Most of the input data required by the program is edit-checked for valid
values before any processing of the Factored Link File is attempted. The
MOBILE1 flags are first checked; if any of the flags are in error, a message
will be printed indicating the total number in error. Note that the program
does not indicate which of the flags are in error. The user should consult
Table 2 for valid values.
Whether or not errors were detected in the MOBILE1 flags, the program
will check the values of OUTFL0 through OUTFL9. As with the MOBILE1 flags,
if ,-my is in error a message will be printed stating how many are in error.
Table 2 should again be consulted for valid values.
One-Time data will be read next. Optional One-Time data is read first
if any is supplied; these values are not edit-checked. The Required One-
Time data is read next. The calendar year is checked (Card 1) to be sure
the value is between 70 and 99, inclusive. The ambient temperatures for
each of the time periods (Card 6) is checked to be sure they are in the range
of 0 to 110°F.
Upon reading Card 6 of the Required One-Time data the program will
test to see if there were any errors in the input data. If no errors were
detected, processing will continue; otherwise the run is terminated.
The following conditions will also cause abnormal termination during the
input process:
1. Premature end of file.
2. Emissions file and/or SYMAP file having the same logical
unit number as the Factored Link file.
3. Emissions file and/or SYMAP file having the same logical
unit number as the Land Use file.
Once the program begins reading the Factored Link file, it will continue
to edit-check the data in the file. If an error exists in the file, a warn-
ing message will be printed and the record (link) in error will be skipped;
but processing of the file will continue. Warning messages will be printed
upon detection of the following errors:
1. Negative district or zone.
2. Link type of 0 or 7.
Within each link, data for each of the four time periods is also checked.
A negative or zero average speed for the time period results in the emissions
for that time period being set to 0. An average speed less than 5 mph will
be set to 5 mph; an average speed greater than 60 mph will be set to 60 mph.
18
-------
If the total vehicle volume for a given time period is negative or zero, the
emissions for that time period will be set to 0. Messages will be provided
i-f any of the above conditions exists.
Sample Data
A sample card deck may be found in Figure 1. Card 1 contains the MOBILE1
flags, OUTFL0 through OUTFL9 and the MODE flag. Card 2 is the calendar year.
Cards 3 through 6 contain the freeway and nonfreeway hot and cold start per-
centages (see Table 6 for description) and Card 7 contains the ambient tem-
peratures by time period. For clarity, the resultant output from the program
is presented in Figure 2.
19
-------
T3
O
-H
t-i
0)
P-i
OJ
-H
H
O
•r-l
13
o
}_J
}-j QJ
3 P-i
O
ffi 0)
1 g
8 59 60 61 S2 63 W6S 66 6788 69 75 71 T2 73 74 75 76'? 78 fflU
Note: Optional One-Time data should be placed between Cards 1 and 2,
Figure 1. Sample input data.
1 1
1 1
??
33
44
55
6b
! i
n
9fJ
1 1 M 1 0
, 1 ll 1 I':
1 i 1 1 1 1
?22???
133333
144444
155555
i (j h G f; f>
n n 1 1
4 B 0 a 8 8
J99999
•1 'i K / K
iO
1 111
1 1
2 -
33
44
55
66
n
88
99
i 11
000
1 1 1
'?-
3 (3
444
555
6 6fi
f i 1
888
999
II \/ 1.1
un
0 ' -
M I'l 16
', 1 1
222
333
444
555
666
11 7
888
999
14 l!i lit
lOP ' '5
00
U IH
:, 1
22
33
44
55
66
1 I
88
99
n if
0000
la fa n n
1111
2222
3333
4444
5555
6666
7777
8888
9999
ia n n n
00000
n /i ?wfi /i
11111
22
33
44
55
66
7 /
88
99
,'.! ?*
2?
33
44
55
66
77
88
99
/'i 76
2
3
1
5
6
7
8
9
"1
0000000000
?»?930 3! J? 33 3-1 153k .17
1111111111
2?
33
44
55
66
77
88
99
?fl t"l
22222222
33333333
44444444
55555555
66666666
77777777
88868888
99999999
,!0 11 3? 33 34 35 3f 37
20
-------
INPUT DATA VALUES
MYMRFG
ALTFLG
IMFLAG
ALHFLG
TRKFLG
NMHFUG
OUTFLO
OUTFL1
OUTFL2
OUTFL3
OUTFL
-------
SECTION 3
SOFTWARE DOCUMENTATION
This section of the manual will give a brief description of the flow
of subroutines GCAINP and GCAOUT. See Figures 3 through 6 for flow-charts
of the routines. Appendix A contains the GCA version of the source code.
CCAINP
Subroutine GCAINP first reads in the Input Flag record, prints the input
values of all the flags as well as a descriptive statement about the course of
action to be followed for the 'OUT' flags only. Any values in error will be
recorded by the routine and a message will be printed stating how many flags
are in error. Note that particular flags will not be referenced; it is the
responsibility of the user to locate any flags in error. Even if errors
exists, GCAINP will continue and read any Optional One-Time data (the values
are not printed by GCAINP) followed by all Required One-Time data (six records).
The Required One-Time data will also be checked for errors. If errors exist
in either the MOBILE1 or OUT flags or in the Required One-Time data, process-
ing will stop. No attempt to return to the main program is made.
When GCAINP is satisfied that all card input data is correct, it begins
processing by calling MOBILE1 subroutines TFCALX, INITEX and LDVIMX (if
requested). A -test is then made to determine which mode of operation was
requested. If the averaging mode was requested, GCAINP will initialize
array EMX by calling subroutine EFCALX 152 times. The subscripts of EMX are
based on the following: time period (4), type of link (2), speed (19), modal
split (6), and pollutant (3). EMX stores emission factors based on 3 mph speed
increments ranging from 5 to 60 mph. Note that there is space for only three
pollutants (CO, NOX, HC) since EFCALX does not compute emission factors for
PM and Pb.
GCAINP then reads a record from the Link file and checks for valid district
zone and link type. Valid link type codes are 1, 2 and 8 (nonfreeway) and
3 through 6, and 9 (freeway). Codes 0 and 7 are invalid. If any of the
above is invalid, the link is completely skipped, and a message is printed
stating why the link was rejected. GCAINP then sets up necessary data for
EFCALX (arrays LNKDTA and MS) and again checks the mode of execution. If
the averaging scheme is utilized, a table look-up into array EMX is performed.
If not, EFCALX will be called to compute emission factors. The program will
then call subroutine EFPMPB to calculate emission factors for PM and Pb. At
the present time this call has been disabled by "commenting it out" since the
emission factor calculation is not being implemented. Emissions for each of
the five pollutants is then calculated (gins) by multiplying the emission factor
22
-------
(REQUIRED)
(REQUIRED)
INPUT FLAGS
AND
ONE-TIME DATA
(OPTIONAL)
(OPTIONAL)
SYMAP AND EMISSIONS
FILE MAY ONLY BE RE-
QUESTED IF LAND USE
FILE IS INPUT
Figure 3. Job file structure,
23
-------
(START)
CALL
GCAINP
CALL
GCAOUT
STOP
Figure 4. Flow-chart for main program,
24
-------
WRITE FLAGS TO PRINTER
WRITE MESSAGE
TO PRINTER
' KKAI) AITKOI'KIATK
OI'TIONAI, ONK-
TIMK DATA
Figure 5. Flow-chart for subroutine GCAINP.
25
-------
INITIALIZE ARRAY EMX BY
CALLING EFCALX 152 TIMES
READ A RECORD
FROM FACTORED
LINK FILE
0
Figure 5 (continued). Flow-chart for subroutine GCAINP
26
-------
no 1=1,4
o
NO
YES
STORE CORRECT DATA IN
ARRAY LNKDTA
COMPUTE MODAL SPLIT
YES
CO, NOX, IIC
CALCULATE EMISSION
FACTORS BASED ON TABLE
LOOK-UP INTO ARRAY EMX
CO, NOX, HC
CALL EFCALX TO CALCULATE
EMISSION FACTORS BASED
ON AVERAGE SPEED FOR LINK
CALL EFPMPB TO CALCULATE
EMISSION FACTORS FOR PB
AND PM
CALCULATE EMISSIONS FOR
ALL FIVE POLLUTANTS
(CONVERT TO LBS IF REQUESTED)
SUM EMISSIONS INTO ARRAY EMIS
BASED ON TIME PERIOD, DISTRICT
AND ZONE AND POLLUTANT
Figure 5 (continued)- Flow-chart for subroutine GCAINP
27
-------
PRINT EMISSIONS
BY LINK
/ WRITE RECORD TO
EMISSIONS FILE
(UNIT NO. OUTFL2)
Figure 5 (continued). Flow-chart for subsroutine GCAINP
28
-------
START
PRINT EMISSIONS
BY ZONE (IF OUTFL3=1,
24 IIR IF=2, ALL TIME PERIODS)
YES
NO
PRINT EMISSIONS
UY DISTRICT (IF OUTFIT
= 1, 2'i IIR TF=2, MI, TIME PERIODS)/
PRINT TOTAL EMISSIONS
FOR EACH POLLUTANT FOR
EACH TIME PERIOD
READ RECORD FROM
LAND HSF, FILE
(UNIT NO. OUTFIT)
Figure 6. Flow-chart for subroutine GCAOUT
29
-------
RETURN
AREA •= VACANT + RESIDENTIAL
+ NONRESIDENTIAL
+ ALTERNATE PASSIVE
YES
COMPUTE ZONAL EMISSION
DENSITY BY POLL AND TIME.
DENS(POLL, TIME)
EMIS(ZONE,POLL,TIME)*1000/AREA
(IN SQMI, SQFT AND ACRES)
YES
/ PRINT ZONAL EMISSION DENSITIES
OUTFL6 TIME PERIOD, OUTFL7 POLL.
OUTFL8 UNITS (SQMI, SQFT, ACRES)
WRITE RECORD TO SYMAP
FILE
(UNIT NO. OUTFL9)
Figure 6 (conrinued). Flow-chart for subroutine GCAOUT
30
-------
times the link length times the total volume. They are stored in array
CALCEM whose subscripts depend on the number of pollutants (5) and the
number of time periods (4). The computed emissions are then summed by zone
and district into array EMIS. The subscripts are (1100, 5, 4). Since valid
zones range from 0 to 999, emissions for each of the zones is stored in
EMIS (1, j, k) through EMIS (1000, j, k); 1 being zone 0, 1000 being zone 999.
Districts range from 0 to 99 and the emissions for each of the districts
occupies locations EMIS (1001, j, k) through EMIS (1100, j, k). The j-subscript
is £or pollutants, the k-subscript for time periods. Note that the calculation
of emissions takes place for each of the four time periods.
The remainder of GCAINP is concerned with printing the emissions for
each link and writing records to the optional Emissions file. A return is
then made to the main program.
GCAOUT
Subroutine GCAOUT is mainly concerned with printing requested output. It
begins by optionally printing emissions by zone and/or district. The user
selects all time periods or the 24-hour time period only. Note that only
zones and districts will be printed which have nonzero emissions for any of
the five pollutants. Total emissions are then calculated and printed for all
four time periods by looping through the district portion of array EMIS
(EMIS (1001, m, k) through EMIS (1100, j, k)). The user has no control over
this printout.
GCAOUT then reads the optional Land Use file. If none is to be read, a
return is made to the main program. The Land Use file is used to compute
Zonal Emission Densities. If a SYMAP file is being written, then GCAOUT will
compute Zonal Emission Densities for each of the four time periods. If no
SYMAP file is requested, GCAOUT bases the number of time periods for which
densities are to be calculated on whether there is a request for printout for
all time periods or the 24-hour period only.
Densities are computed only for zones-and districts which are read from
the Land Use file. They are computed in three units: grams or pounds per
square mile, per square foot and per acre. The option of grams or pounds
depends on the value of OUTFLO.
Zonal Emission Densities may then be printed for one, some or all of the
pollutants and for the 24-hour or all the time periods. The last function of
GCAOUT is to write individual records to the SYMAP file. Records are written
only for the zones which are read in from the Land Use File.
31
-------
APPENDIX
FORTRAN SOURCE CODE
32
-------
iv GI BH.EASE 2.0
DATE : 7833?
7/«5
PAGE 0001
0001
OOOi
on 03
OOOa
DATA SET PTNirOT AT LEVEL 007 AS OF 11/28/78
PGM TO TFST MODIFICATIONS TO
CALL GCAINP
CALL GCAOU1
STQP
00000000
00000010
00000020
00000030
00000040
00000050
ooooo r, 60
-------
FORTKAN IV Gl
0001
0002
000}
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
001U
0015
OOlb
0017
0018
0019
RELEASE 2.0
BI_K DATA
DATE
7833?
12/17/U5
PAGE 0001
BLOCK DATA
»«»»**THIS BLOCK DATA TO BE USED "ITH GCA MODIFICATIONS
COMMON /FLGCOM/ ALTFLG,ALHFLG,TRKFLG, I RDFLG . SP JFLG, NHHFLG , IDLFLG,
1 INFLAG,VMFLAG,UNFFLG,MYMRFG.IMFLG,ICEVFG,PRTFLG,IFORw
INTEGER ALTFLG,sLHFLG,fokFi r,, SPJFLG/VMFLAG.UNFFLG^PRTFLG
COMMON /REtCOM/ IREJN
CUMWON /GCAQTA/ EMIS(1100,5,U)/COLOST(2,U),HOTST(2,^),TEMPTPCa)
COMMON /GCALNK/ FUNCCD,JURIS,DIST.ZONE/NUOEA,NODES,FEDAID,LMKLEN,
1 LNKTYP,*ArcoD,ARE ATY,POLUNT,FILEID,PLANLK,ZONSPD,HALFWD,SRCCOD,
^ GRADCD,CAP4C,AVGSPO('J),VCBfa),T(lTvnL(U),vnLVTl(«),VOLVT2(t(),
3 VCJLVT3(«) ,VOLVT«f«), VOL VT5(a),vOLVTb(U),8EGHB(«), LENT IMta)
INTEGER FUNCCO.DIST,ZO^E,FEOAID,nAYCnD,AREATY,POLUNT,FILEID,
1 PLAN'LK,ZONSPD,HALFI«D,SRCCOD,GRAOCD,CAPAC,AVGSPD,VCR,TOTVOL,
2 VOLVT1,VOLVT2,VOLVTJ,VOLVT«,VOLVT5.VOLVT6,BEGHR
COMMON /GCAIO/ LINKFL,NCRDRD,NPRNTR
COMMON /GCAwRK/ IEMISSC5,«),CALCEM(5,«),1UNITS(2),IDPOLL(5),
1 JHOUR(U),IOENS(3)
DATA INFLAG,SP3FLG»VMFLAG,IFORM,PRTFLG.ICEVFG,IDLFLG,
1 UNFFLG /2*0»2*lf0,2*1,0/
DATA IREJN /!/
DATA CAUCEM,EMIS /20«0„0,11000*0.0,11000*0.O/
DATA BEGHRd ) ,LENT!M{ 1 ) /a»0/
DATA IUNITS /2HL8,2HGM/
DATA IDENS /aHSQMJ,«HACRE,1HSOFT/
DATA LINKt-L,NCkQRD,NPRNTR /«,5,b/
DATA JHOUR X2«,1,2,3'
DATA IDPOLL /1HCO ,UHNOX ,flHHC ,UHPB ,UHPM /
&«»Ba«aa9«aVCia**™«««o»»vwa**rawaw0«v«v*wVa*™vav*>a>a*av«wa
-------
IV .
0001 c. . .SUBK°UT1NE EFPMPB(EMF) 00005230
"""
OOOJ
0001 RETURN 00005280
0005 END 0000529Q
00005300
-------
F 0 P. T « 4 •-
0001
GI
0002
0003
oooa
1005
0006
0007
0008
0009
0010
001 1
0012
0015
001«
0015
0016
0017
0018
0019
0020
0021
0022
0023
0025
:i_e A
c
c
:
c
c
c
c
c
c
c
S£ 2,0 GCAINP 0«TE = 78358 15/21/3"
SUBROUTINE G C 4 I N P
D. CARACC IOLO JR.
GCA/TECHNCLUGY DIVISION
"?ijWLINGT<"'N BO .
RE OF OR H) M 4 , 01730
517-275-9000
HijRRCUGHS 85700 * I8W 370/158
COMMIT-, /SPicO"/ SPl.SPAv
C 0 " M fj •, / I M c c M / ICYIM,iST^l\,IMTFL.G,MODYRl,MCCYR2
C 0 M M (_• \ / M Y H C C M / "•' 1 M , M Y R , T F
COMMPJ, /JU\KI/ AC , XLOAD, TR A IUR , 4HSHIJM
COMMON /JUNK2/ MGnGT, HDlGT , MGC ID, HDC ID
REAL
-------
00
,FO«1R*N_ IV Gl
0001
0002
0003
000"
0005
OOUb
0007
Ou08
000'
0010
001 1
0012
0013
00 1"
0015
OOlo
0017
00 16
0019
0020
0021
0022
0023
002"
002b
RELEASE
C
C
C
C
C
C
c
c
c
c
2,0 GCAINP. DATE i 79089 . _09/2J/50._
SUBROUTINE GCAINP
D. CARACCIULO JR.
GCA/TECHNOLOGY DIVISION
BURLINGTON RO,
BEDFORD MA. 01730
617-275-9000
BUHRQdGnS B570il • IB« 370/158
C.T-IUON /SPICOM/ SPI.SPAV
COMMON /IMCOM/ 1CYIM, ISTRIN, I*TF Lt , MOO YK) ,MCK'YK2
CQMM,jiv /MYMCOM/ MYM,MYR,TF
COMMON /JuNKl/ AC , XLOAO, TRA ILR. AbShUH
CO"HOS /JU.NK2/ HGrsbT , HDrtGT , H(iL ID, HDCIO
"fcAL xLOAD{3).'-YW(Jo,6),tfYR(20,b),TF(20,h)
CQMVU'V /FLGCOM/ ALTFLG,ALHFLG,TR I F L 3 , OuT F L U , OUT F L5 , OUT FLb , 'JUTFL7,
\ DIJTF Lr , vu 1 F L4
CO" '10'. /(,CALNR/ FbNCCD,JuRIS,UIST,ZONE,NODFA,.1-,OOEd,FEOAIU,LNKLEiv,
1 L?jMYP,i>,AYCUO,AREATY,POLU'JT,FILEID,PLANLK,ZOi>tSPD,HALFwD>SRCCOD,
2 GRADCO,CAPAC, AVGSPO(U), VCRl") , TOT VOL ( U ) , VOL V T 1 ( u ) , VOL V T2 ( u ) ,
3 VULVT3C4),VOLVT«(«),VOLVT5CU),VOLV1t>(«),BEGHR(4),LENTI«(u)
1'. TEGE" Fu'.CC3,niST,iu'.L,FtOAID,ftAYCOO,AHEATY,FOLU'.T,FlLEID,
1 PLiVL^£0\SP|;,nALF«C,SHCCL'D,GWAtCD.CAPAC,AV,&bPi;,vCK,IC'rvOL,
CuMM0". /GCA10/ LIN".FL, NCRDRO, NPRNTR
Cn«"0v /C,CA.-»«n/ IE»ISSC5, u) ,CALLE"(5, " ) , lu^'I Tb(2) , H'POLL (5) ,
1 J-!'j.jMfcP 7h
vt*
-------
FORTRAN IV Gl RELEASE 2,0
GCAINP
DATE = 79069
09/21/50
PAGE 000?
u>
CO
0026
0027
0026
0029
OOiO
0031
0032
0033
003«
0035
0036
U037
0038
0039
OOUO
OOul
00^2
UOJ5
00 uo
00«7
0016
0049
0050
0051
0052
U053
uu5
-------
_F_G_RTRAN
.MkEA_3£_!j.O __. G,CAJNP__
i 79089
0066
0067
0066
OJ3e>9
0070
0071
0073
007u
007S
007o
0077
00 !B
0079
0080
0081
0082
0083
008U
0065
0086
0087
0088
0089
0090
J091
0092
0093
009u
009S
0097
009f
0099
0100
0101
01 o2
olo3
01 Ou
*RITE(NPRNTR,797=99o)WY^
WYMRFG .EG. 3)REAO(NCRURD,7010,ENO=990)"YR
i ,
1 ) d t A 0 ( \ C * J « D , 7 f ?.. , t M> = 9 9 C ) I C Y I >•' , I S T R 1 V , I v T F L G ,
1 JREAIMNCRDRO, 7030,
00001310
00001340
00001 J50
00001360
00001370
00001380
00001390
ooooiaoo
00001010
00001U20
00001U30
00001440
00001050
OOOOlUeO
0 0 0 0 1 u 7 0
00001180
00001U90
00001500
00001510
00001520
00001530
00001500
00001550
00001560
00001570
00001580
00001590
00001600
00001610
00001620
00001630
00001600
00001650
000016eO
C000167C
00001680
00001690
00001700
00001710
C0001720
00001730
00001700
00001750
00001760
00001770
00001780
00001790
00001800
-------
FORTH** IV 61 RELEASE 2.0
DATE i 79089
09/21/50
PAGE 0004
-P-
o
0105
JM06
OlUT
" 0108
0109
0110
01 11
0112
0113
01 IK
0115
0116
01 1 7
0116
0119
0120
0121
0122
0123
0121
0125
0126
0127
0128
0129
01 3 (HOTST(lflO,K«l,«)
50_fENps995) (HOTST(2,K).K«l,(i)
.•TRITE I NPRNTR, 7932) (HOTS T(2,K) , K = l ,«)
REAOlNCROR0.70SO,ENO=995)(COLOST(2,K),K=i,a)
(COUOST(2,K) ,K = 1,«)
ARITE(NPKMR,7935)TEHPTP
00 70 n = l , u
IF(TE^PTPtM ,GE. 0.0 .AnD.
TEMPTP(K) .Lt, 110.0)GO TO 70
nRITE(NPWNTR,7905)JnOUR(K)
70 CONTINUE
IK'»EWR1 .GT. 0
NERR2 .GT. 0 .OR. NERR3 .GT. 0)GO TO 9999
C4LL TFCALX
CALL INITEXtNMYALT, ALTKOD)
IFCIMFLAG .EG. DC&LL
IF(LIN.' 77 KP=1 , 3
77
00001810
00001820
00001830
000018UO
00001850
'0"0001860
00001870
00001880
00001890
00001900
00001910
00001920
00001930
ooooi9ao
00001950
000019bO
00001970
00001980
00001990
00002000
00002010
00002020
00002030
000020UO
00002050
00002060
00002070
00002080
00002090
00002100
00002110
00002120
00002130
000021UO
00002150
00002160
00002170
00002180
00002190
00002200
00002210
00002220
00002230
000022UO
00002250
00002260
00002270
00002280
-------
GCAJNP
DATE i 7.9089
0145
0140
01_«7_
0 I ifl
01U9
0150
0151
Ji52
0153
015U
0155
0150
0157
Ol5b
0159
oioo
Olbl
Olb2
0163
0 Ibu
0 lob
0 loo
Olo7
o 1 00
U109
0170
0171
01/2
j 1 ' 3
V 1 I U
01 75
01 7o
017/
ul M
or/ 9
j i au
0181
90 CONTINUE
C
NLINKS=0
NBAD=0
500 READ(LINKFL, 7 1 00 , ENDs2000 ) FUNCCD, JURI 3 , DI ST , ZONE , NODE A, NOOEB ,
i FEPAID,LNKLEN,LNKTYP,WAYCOO,AREATY,POLUNT,FILEIO,PLANU«,ZONSPD/
2 HAl_F«D,SRCCOD.GRADCD,CAPACi_AVGSPC)(J),VCR(l)>TOtVOLU)fVOl.VTl(l)/
3 VOLVT2U ) , VOLVT~3( I) ,VOLvTu(l"f, VOLVT5( 1 ), VOLVTbd ) , (SEGHR(K) ,
J LE'' x F K = 1
IF(LNKlYP ,EU. 1 .OR. LNKTYP ,EU. 2 .OR. LNKTrP .EQ. BHNDXFR:2
I 'jOXZ'< = ZO^Et 1
INDXDS=01ST+ 1 001
C
•jj uco *uvt = i , "
IF ( avGSPO CKT IUE) .GT. 0)GU TO BOO
nRlTf(JPK(vTh,7950)NLlNKS,''lUOEA,NCOE6,AvGSPDtnTlME),JHOUR(KTIMt)
00 790 KP=! , 5
C'LCE^«P,'sTIvt)=0.0
79j CJ'iTI'iUE
GO T j 1 OuG
600 IFCAVGSPD(KTIMt) . GE . 5)GO TO 810
«RITE(NPRNTR, 7951 )M INKS, :jQDEA,NODEb, A VGSPD ( K T I ME ) , JHOOR ( K T I «E )
AvGSPD(KT IMt)«5
GH ro &o j
"l.i IF CtvGSP'J(nTl-t) ,LE. oO)L.U Tu frou
A^ITt(",p^-jlR,7°52).
-------
FORTRAN IV &1 RELEASE 2.0
CCAINP
DATE « 79089
09/31/50
PA6E 0006
0182
6TB3
018U
018b
0186
0(67
0186
"0189
0190
0191
0192
0193
U]9u
0195
0196
0197
0198
0199
0201
0202
020J
020u
0205
0 2 Go
0207
u20o
0209
0210
021 1
0212
0213
021
0215
890 CONTINUE
C COMPIj.E_ APPROPRIATE INPUT DATA FOR EFCALX
"L.NKDTAC I JBAVGSPO(KTIME) "~
LNKOTAl2)sAVGSPO(KTlME)
LNKDTA( 3
LNKOTA(5>HCOLDSTUNDXFR,KTIME)
LNKDTAlb)3MOTST(INDXFR,KTIME)
Lf-iHDTAC7)=COLDSTCTMDXFR,KTlME)
COMPUTE VMT MIX (WQOAL SPLIT) FOR TIME PERIOD
IF(TOTVOLCKTIME) ,GT. 0)GO TO 895
RRITEtNPRNTR, 798 3) BLINKS, NODE A, NODES, JHOURCM IMF.)
DO 892 KP=1,5
CONTINUE
GO TO 1000
CONTINUE
MS(l)=FLOAT(VOLVTl(KTIME))/FLOAT(TOTVOL(KTIMEn
MS(«)=FLOAT(VOLVT2(KT1ME))/FLOAT(TOTVOL(KTIME))
^S(5)=FLOA!(vOLvr3(K]lME))/FLO£.T(TOTVOL(KTI"E))
MS(2)=FLOAT(VOLVTU(KT1W£))/FLOAT(TOTVOL(KTIME))
892
895
MS(6)sFLOAT(VOLvT6(KTIME))/FLOAT(TOTVOL(KTI ME))
C
C«*»**«»*tiDETERMINE MODE OF EXECUTION
IF (LI^KS-H)GO TO 898
C
C
C
JSE E^lSblG^ FACTOR F'OR \EAREST SPEED
-------
_f_0*URAN IV 61
0210
02 1 7
0218
0220
0221
0223
022u
0225
022°
0227
0228
0229
0230
0231
U232
0233
023«
0235
0236
0237
0230
0239
0200
GCAIMP
VEHLNK«FLO»T(TOTVOL(KTIM£))«RLENTH
C -CQ-
CALCEM(1,M1ME)=VEHLNK»"TEDEF(2)
IFCQUTFLO .EQ. 0 ) C ALCEH I 1 , «• T IME ) *C ALCEM ( 1 , K I IME ) /053 . 59
CALCEMT2/KTJ«E)aVEHLNK»*TEDEF(3)
IF(OUTFl_0 ,EO. 0)CALCEN(2,"TIME)«CALCEMtJ,KTIM£)/u53.59
C -HC-
.JLOQJ.-
00003220
IFCOUTFLO ,tU. 0)CALCt"(3,l
-------
FO«TRAM
61 RELEASE 2.0
GCAIMP
-p-
-p-
02«8
02<49
0250
025!
0252
0253
025u
0255
0256
0256
U259
0260
02ol
02e>2
0263
026"
0265
0260
02oo
IF(OUTFL1 .GT. 53
°*Tt_«_7?i89_ ___ 09/21/SO
,NODES,HLENTH,LNK1YP,ABEATY,FILE 10,
PAGE 0008
00003700
OOOQJ710
2 7.0NSPO,HALFrtD,GRADCO,CAPAC,"" 00003720
3 CTOTVOL(*T ).KT=l/1>,((IDPOLL(KP),CALCE*(KP,KT),KT=l,u),KPEl,KPOLLOOOOJ7JO
<• ) 00003740
C
C
0*******»»NRITE R_ECOHD TO, EMISSIONS FILE IF REQUESTED
1500 IFCOUTFL2 . E~Q. 0)GO TO 500
IvCR=VCR(l)*10
IFflvCN ,LE. 9999JGO TO 1550
?,793B)NR
I VCH«9999
1559 CONTINUE
00 |bOO Jsi,«
30 loOO 1=1,5
IEM1SS(I»J)=CALCEM(I,J)+0.5
1600 CONTINUE
C
'
-------
FORTRAN IV Gl RELEASE 2.0
GCAJ.NP.
DATE a 79089
09/21/50
___P_AGE
f>
Ul
0272
027J
0274
0275
0276
0277
027B
U279
0280
0281
0282
0283
0265
0287
0288
0269
0290
0292
0293
029u
0295
029b
0297
0298
0299
0300
0301
u302
7040
7050
7052
7100
7tt90
7495
7500
7510
751 1
7520
'521
7530
7531
7532
FORMAT(DFS.O)
FORMAT(4F5.0)
FORM4UI?)
FORMAT (2I1,12,13,T5,2I5,I1,I3,3I1,I4,A4,I4,11X,2I2,2I1,I6,I2
1 I4j Io,3I5,314,3(3L2,2Jb,3IS, 3J 4) )
FORMATUOH MODE « .I 2,IX,79HEHISS10NS KILL
1LIMK USING EMISSION FACTORS FOR GIVES SPEfcO)
FORMAH10H MODE » ,12,4X,95MEHI S3 IONS HILL 6E COMPUTED FOR EACH
1LINK USING EMISSION FACTORS FOR DEAREST SPEtO IN 3 MPH INC.)
TO BE COMPUTED IN ,A2,1HS)
FORk'ATClH*,15x,28HEMIsSIONS
FLW*A1Clri + ,lSx,43HE«ISS10NS t, ILL
FORMST(1H»,15X,44HEM1SS10NS KILL
1 5A4)
FORMAT(1H+,15X,34HEMISSIONS
F OR MAT (iM-t, 15x, 39HE MISSIONS FILE
FORM*T(|n+,15x,37 rit MISSION 5 MLL
FOR''ATllH+,15X,bOnEMISSIONS A ILL
1 PERIOD ONLY)
FORMATUHt,15X,54HEMISSIONS
1ERIODS)
75ui
75U2
7550
7b5l
75oO
75nl
75&2
7570
7571
7580
7S9u
7591
7oOO
77^0
7710
00004160
00004)9?
00004200
00004210
00004220
BE COMPUTED FOR EACH 00004230
00004240
00004250
000042bO
00004270
00004280
, 00004290
00004300
00004310
00004320
00004330
T JME00004340
00004350
iiILL BE PRINTED BY ZONE FOR ALL TIME POOOOl3bO
00001370
nILL NOT bE PRINTED bY DISTRICT) 00004380
MLL bt PRINTED bY DISTRICT FOR 24 HR 00004390
00004400
WILL BE PRINTED 8Y DISTRICT FOR ALL
M'T BE PRINTED
BE PRINTED FOR
FOR fACH LINK)
EACH LINK FOR ,
FILE WILL NOT BE BRITTEN)
rtILL BE BRITTEN TO UNIT ,
NOT taE PRINTED BY ZONE)
bE PRINTED BY ZONE FOR 24
12)
riR
ITIMt PERIOD ONLY)
FORMATC1H+, 1 5X , 58M£»I SS10NS
1ME PERIODS)
FORMAT C1H+, 15X, 30HLAND USE FILE »»ILL SOT BE READ)
FORM4TC H + , 15X,37HLAND USE FILE A I L L BE READ FROM UNIT ,12)
FOR'-ATl >-», 15*, JuHZO'-iflL E«ISSIOi- DENSITIES AlLL NOT BE PRINTED)
FQRMii; M + , ]b*, 70-ZONAL E«!SSIl'N DENSITIES rlLL bE PRINTED FOR 2u
I HOUR T *t PERIOD ONLY)
FORMAIi ht, 15x,blriZONAL EMISSION DENSITIES *ILL BE PRINTED FOR
1 TIME" PERIODS)
t-OR^AKlh*, 15X,82nlF ZONAL EMISSION DENSITIES ARE PRINTED,
1 OF A'_L P'JLL'.il ANTS '.ILL flt iSSu'-'CO)
KlR^AT ( l~t , lb> ,u5«Zu>.AL t"!SSH'. DENSITIES ^ILL BE PRINTED
1 5A4)
ZO'JAL EMISSION DENSITIES ARE PRINTED,
L 6E ASSUMED)
A|. E^ISSlLi. PENSHIES rtlLL 3E PRIMED
, 15<,58nlF
. A4, le« rtlL
, l5",5o"Zu'i
FORMAT C IH+
1 DF UMTS/
FIIKVAT i 1-. +
IS/, An)
FJR«4T ( l-i+, 15».29M\CJ SY-'AP FILE .'ILL b£
FORMAT i IM+, i5x,3bHSYMAp FILE is TO BE BRITTEN TO UNIT ,12)
FOR"4T(2I1,12,215, 11, 13, 31 1 , 14, Ao, 14,1 IX, 212, 21 1,16, 12,14,
1 Ib, I 1 ,2018, 1 IX)
FORMAT (uO«l»****EOf E>COUMTtRE3 IN LINK FILE«*«*» /
1 34>'C - i.u-'-BEx JF L1N.KS ^"OCcSSED = ,]7/
2 3at-o . j"BE« -JF u!N*S 5"IPFE.n = ,17)
FORWAT152HO MU«BfcR OF RECORDS nRITTEN TO EMISSIONS FILE
T 100004410
00004420
00001130
00004440
00004450
000044bO
00004470
ALL00004480
00004490
DEFAULT00004500
00004510
FOR , 00004520
00004530
DEFALLT00004540
00004550
I'-. L N I T C 0 0 0 4 5 o 0
00004570
00004580
00004590
00004bOO
OOOOUblO
00004b20
00004b30
00004b40
= ,I7)00004b50
-------
IV Gl RELEASE
GCA1NP
DATE E 7908~9
09/21/50
PAGE 0010
OJOJ
0104
0105
0106
0107
l'30b
0109
Oilu
Gill
0312
U313
U31M
0315
olio
Oil 7
U Jib
U319
ui20
0121
ui22
0 J2i
O.i2u
0325
0326
0327
0326-
u3?9
Oiiu
uiil
0332
IN , 42, ! HS, 3«, 0000*660
-» 24 HOUR • •••• 00003670
T.P. 1 ---"/1X, 00004680
1 3b* D A;(GDE BNODE LL T T 10 SP Hr, C,7X,4( BH BEGINS ,12, 7H L0000«690
-ENGTM , I2.2X)/45X,4(2|M POLL EMISSIONS )/ ) 00001700
7810 FORMAT (/li, 12,2(1 X, 15), 1X,F4.2,2(1X,I1),1X,A4,2(1X,I2),1X,I1,1X,I600004710
7800 FORMAT(lHl,aOX,ln
!li- «t C .
-NO.: ,14,
IN ERROR
Lt SS 1 Hu
FOh
ANODE: ,15,
TIME PERIOD
,15, 8H
P E K I 0 L1
Rt ADING
0000u910
READIvG nOT/COLD START 00004940
00004950
8" D'MODE: , I5.16H AVERAG00004960
,I2,21H EMISSIONS SET T000004970
00004980
8H H\ODE: ,I5,lbM AVERAG00004990
1
/952 FQRW4T(
IE SPEtD: ,
i T . b .:• )
F JR-fl T ( ) on
F JS'-'AT ( 1 0-
FORMAT( 10H
"A T(1 OH
FORMAT ( U'H
F j •; •'• 4 T ( 1 ..- •-
f (iR^i M H'"
F j-C- AT ( H f
F-1RMA T ( ) (in
j- t >.
7961
79b2
79bi
79b4
I2,33H
'• Y •• ^ F u
ALTFLG
IMF LAG
ALHFLb
T K n F L &
CUTFL2 =
- I u , rtM
.3 T H A f.
ANODE:
bO FOR
-12)
-12)
- 12)
. 12)
.12)
12)
-12)
-12)
-12)
I 2, I 5 H SPtEO SET 10 00005000
00005010
,15, 6h B\OOE! ,I5,lbH AVERAGOO005020
TI-E PERIuO ,12,loH SPEED SET00005030
00005040
00005050
00005000
00005070
00005080
00005090
00005100
00005110
00005120
00005130
-------
JFORTRAN J.y_Gl__ RELEASfc Z_L0 GCAI.NP
DATE
79069
09/2J_/50.
PAGE 001J.
0113
0335
0336
U33?
0336
0339
0340
03«1
7970
7971
7972
7973
797u
7975
7976
7980
'961
/982
7983
/Q9Q
FORMATCIOH OUTFL3 - ,12)
FORMATUOH OUTFL« » ,I_2)
FORMiT ( IOH OUTFL5 = ,12)
FORMAT ( IOH OUTFL6 = ,12)
FO«"AT(10M OUTFL7 i ,12)
FORMATUOH OUTFL8 = , 12)
FOR*AT( IOH
FOHMAT( 1
1,12, 2<4H
Fu«^AT ( 1
1 , I i, 24M
F U H M A T ( 1
1YPE: ,11,
in
I '">
If
li.
IM
OUTFL9 = , 12)
RtC. NO.:
LK«O« . i. I
REC. M'. :
thht.iH - LI
REC. Nti,:
2UH I'N EK^LJW
F ORMAT ( 1 In
1 vOLL'«t IN
FU^'-« r (21
E -D
REC. .^0.:
tPthn FUt* 11
"0
«*.«*E',0 OF
, I U , ttH
ANODE: ,
15,
8H ((NODE:
,15,
7H
•>« SMPPEu)
,14, 6H
;,r\ Sn IH(
, I U , 6n
- LINK
, I U , 8h
A *v
3 f r
Ar,
SK
Ar.
>••£ ^t"l(j[.
J 0 n » • « i
> • i
uuE : ,
j
'.'DE: ,
1PPEO)
DDE: ,
,12)
IS,
15,
15,
HH oNUDE:
8h 6NODE:
6h BSOOE:
,15,
,15,
,15,
7H
12H
39H
oooo5iao
00005150
00005160
00005170
00005180
00005190
00005200
DIST: 00005210
00005220
/!ONE: OC005230
000052UO
LINK T00005250
00005260
TUTAL 00005270
0000528C
00005290
00005300
-------
FORTRAN IV Gl
0001
RELEASE 2.0 GCAOuT
SUBROUTINE GCAOUT
DATE t 7933?
12/1 7/
-------
Iv 1,1 RELf-ASt 2.0
GCAOUT
DATE = 78332
12/17/U5
P4GE 0002
002"
002S
002b
0027
0028
002°
0030
0031
0032
0033
003<4
0035
0036
0037
0031
0039
oouo
00111
00u2
OOUJ
ooau
OOU5
OOiife
00«7
OOU8
OOU9
0050
0051
0052
0053
0051
0055
0056
0057
POO
nRlTE(NPRNTR,7500)JTYPt(l)
1 IUNI TS(OUTFLO« 1 ) , JTYPt ( 1)
D(J 900 KZ«1 , 1UOO
DO 800 KPsl ,5
IF(EMIS(KZ,KP,M) . G I , O.OJbO TO 850
CUNTINUt
GO TO 900
NPAGEi JHOUKUT ),BEGH«(KT)/LtNri«(KT)
IDPPLL
850
LlNES'LINt 3*1
JZONEiKZ-1
.ME. 0)GO TO 87S
875
900
1000
*RITE(NPRNTR,7500)JTYPE(n
IUNITS(OUTFlO+l),JTYPE(l)
woj ye (\IPRNTR, 7'
CONTINUE
CONTINUE
^JPAGE,JHOUPC^T)/BtGHS(KT),Lt^TIM(K^)
IDPOLL
1050
EMISSIONS BY DISTRICT
.EQ. 0)GO TO 2010
.EQ. 2)KTIME=«
***J*PRINT
IF(OUTFL1
KTIMEsl
IFtOUTFta
00 2000 KT«1 ,KTIW£
LINES=0
NPAGE=1
wRITE(NPRNTR,7500)JTYPE(2)
1 IUMTS(OUTFLO»1), JTYPE(2)
. NP*GErJHOUR(KT),BEGHR(KT),LENTIM(KT),
.IOPOUL
1800
1850
00 1900 K0«1001,l
DO 1800 KPil,5
IF(EMIS(KO,KP,KT)
CONTINUE
GO TO 1900
100
.GT. 0.0)GO TO 1850
,NE. 0)GO TO 1875
1875
1900
2000
LINESrUNES+1
JOIST=KO-1001
IF(MOD(LINES,SO)
NPAGE'NPAGE*!
WRITE(NPRNTR,7500)JTYPE(2),NPAGE,JHOUR(KT),8EGHR(KT),LENTIM(KT),
IUMTS(OUTFLO + 1),JT¥PE(2). IDPDLL
WRITE(NPBNTR,7510)JOIST,(EMISCKD.KP,KT),Kpsl,5)
CONTINUE
CONTINUE
•PRINT TOTAl EMISSIONS FOR EACH TIME PERIOD FOR EACH
LOOP THRU DISTRICT PORTION OF EMIS ARRAY
00005790
00005800
000058 1 0
00005P20
00005830
000058UO
OCO"5fl50
00005860
00005870
00005880
00005890
00005900
OOOC5910
00005920
00005930
000059UO
00005950
00005960
00005970
00005960
00005990
00006000
00006010
00006020
00006030
000060UO
00006050
00006060
00006070
00006060
00006090
00006100
00006110
00006120
00006130
000061UO
00006150
00006160
00006170
00006180
00006190
00006200
00006210
00006220
00006230
000062"0
POLLUTANT00006250
00006260
-------
FORTRAN IV Gl RELEASE 2,(
GCAOUT
DATE = 78332
17/U5
PAGF 000}
k-n
O
005H
0059
0060
0061
0062
0063
006U
0065
006b
0067
0068
0069
0070
0071
0072
0073
007U
0075
0076
0077
0078
0079
0080
0081
0082
0083
0084
0085
0086
0087
2010 NPAGt»l
DO 2025 HTsl ,U
DO 2025 KP=l,5
Od 2025 KDslOO 1,1100
SUM(H.P,KT)=SUM(KP,KT)+EMIS(KD,KP,KT)
20JS L
01GO TO JOSO
nRITE CNPRNTR, 7900) lUNITS(Oi.iTFLOtl),NPAGE,IDPOLL»JHOuR(n,
1 (SIJW(KP,1),KP=1,5), (JHOUR(KT),HFGHW(KT),LtNTI"(nT),fSUM(KP,KT)i
2 KP=1,5},KT=2,u)
IF(OUTFL5 .ME
RETURN
2050 NLANDsO
LINES=0
IF(DUTFL7 •ER. b)KPOLL«3
BASE THE MjMBEfi OF TIME PERIODS FOR «HICH DENSITIES ARE
TO Bfc COMPUTED ON REQUEST FOR SYMAp FILE
IFIUUTFL9 ,6T. 0)GO TO 2060
NO SYMAP FILE REQUESTED.
IF(OUTFL6 .EG. 0)GO TO 9999
1FCOUTFL6 ,EQ. 1)KT]ME=1
LAND USE FILE
BASE * OF TIME PERIODS ON OUTFL6
2060 READ(OUTFI,5,7000,END = 5000)DIST,ZONE/ V AC ANT , RESDEM . NONRES , ALT PAS
NLANOiNLAND*1
+ RESOEN + NONRES + ALTPAS
DENSl 1,V.Z) - DENSITY IN UNITS/SO MI
OEMS(2,Y,Z) . DENSITY IN UNITS/ACRE
DENS(3,Y»Z) - DENSITY IN UNITS/SQ FT
DO JOOO KTzl.KTIME
DO 3000 KP=1,5
DENSC1,KP»KT)»-1.0
OENS(?,KPtKT)i-l.0
DENS(3,KP,M)*-I .0
IF(AR£A .tQ. 0.0)GO TO 3000
DENS(3,KP(KT)*EMIS(KZ,KP,^T)«1000.0/AREA
00006270
0000628C
00006290
00006300
00006310
00006320
00006330
00006510
00006350
00006360
'1 U 0 0 6 3 7 u
00006380
00006390
00006UOO
0000fc«10
00006U20
00006U30
00006UWO
00006150
00006a60
00006U70
00006«80
00006U90
00006500
00006510
00006520
00006530
000065UO
00006550
0000656P
00006570
00006580
00006590
00006600
00006610
00006620
000066JO
00006610
00006650
00006660
00006670
00006680
00006690
00006700
00006710
00006720
00006730
000067UO
-------
FC3TH4N IV Gl ReLE«SE
DATE z
12/I7/ub
0088
0080
0090
0 C ^ 3
0091
0095
U097
0098
0099
01 00
0101
0102
0103
010U
0105
0106
0107
0108
0109
01 10
0111
01 1 2
01 13
011U
1000 CONTINUE
IF(OUTFL6.EQ.O) GO TO UOOO
IF(OUTFL7 .GT. 5)GO TO 3050
•ONE POLLUTANT TO BF PRINTED
00006750
00006760
00006770
Ofi 006760
0000b790
0000679!
00006800
00006810
0000bf>20
0000&830
IF (L INES . EQ. 0)HfiITE(NPRNTH, 7700)1 UMTS(»ijTrt,1«M,;TFNSrn'iTrL>4),
1 KiPAGEi IDPOLLCOUTFL7)
IFlOUTFLb .EO. 1 ) LINE S=L I Nf S+2
IFCOUTFLb .EQ. 2 ) L I NE S=L I NES* 5
IFf^CiDtL INES,50) ,NE. 0)GO TO JOJO
wR]TE(fMPRMR,7700)IUf»ITS(OUTFI,0»l),IOENS(OUTFL8).NP4GE.
1 IDPULLtOUTFL?)
1FIOUTFL6 .fcQ. 1)*RITE(NPRNTR, 77101ZUNE, JHOUHf 1 ) .
1 OENS(OUTFL8,OUTFL7, 1 )
IFlOUTFlb ,EQ. 2 ) WR I TE ( MPRNTR , 7720 ) ZONE , JHUUR ( 1 ) ,
1 DENS(OUTFL8,OUTFL7,1),(JHOUR(KT),BEGHR(KT),LENTIM(KT),
2 DENS(OUTFL8,OUTFL7,KT ) ,KT=2,«)
GO TO UOOO
•MULTIPLE POLLUTANTS TO BE PRINTED
0000b850
00006P60
00006870
0000686C
0000b890
00006900
00006910
00006920
00006930
000069UO
00006950
00006960
00006970
00006980
00006990
00007000
00007010
3050 IFtLINES .EQ. 0 ) wR I TE ( NPRN TR , 7700 ) I UN I TS (OIJTFL 0 + 1 ) , I DENS (OUTFL8 ) , 00007020
1 NPiGE,(IDPOLL(KP),KP=|,KPOLL) 00007030
IFCOUTFLb .EQ. ULINESsL1 MES + 2 000070«0
IF(OUTFL6 .EO. 2)LINES=LINES+5 00007050
IF(MOD(LINES,50) .NE. 0)GO TH 3070 00007060
NPAGE=NPAGE+1 00007070
WRITE(NPRNTR,7700)IUN]TS(OUTFLO+1),IDENS(OUTFL8),NPAGE.(IDPOLL(KP)00007080
1 ,KPal,KPOLL) 00007090
C 00007100
3070 WRITE(NPRNTR,77JO)ZONE,JHOURt 1 ), ( DENS (OU TFL8 , KP , 1 ) , KP«1 ,KPOLL) 00007110
1FCOUTFL6 .NE. 2)GO TO
-------
FORTRAN i v GI
0115
0116
0117
0118
0119
0120
n 1 21
0122
0123
0124
0125
0126
0127
0128
0129
0130
01 Jl
0132
0133
0134
0135
0136
0137
0138
0139
0140
0141
01 42
0143
0114
0145
0146
0147
0148
RELEASE 2.0 GCAOUT
IF (OUTFLO ,EQ. 1 )GO TO u] 00
C
C EMISSIONS ARE CURRENTLY IN
00 uo50 K T = l , 4
C .NOX-
JL8SU )=EMIS(KZ,2,KT)*10.0 +0
JGMS1 1 )=EMIS(KZ,2,KT )«U5J5.9 +0
1 = 1 + 1
C -HC-
JLflSI I )=E*IS(HZ, 3,KT )* 10,0 +0
JGWS ( I ) =EMIS ( *Z i 3 , * T ) *4S 35 , 9 +0
1 = 1 + 1
C -CO-
JLBSC I )»EMIS(KZ, 1 .* I )*] 0.0 +0
JG"S( I)sEMIS(KZ, 1 ,KT)*4535.9 +0
1 = 1 + 1
C -pf -
JLBS(I)BEM1S(KZ,5,KT)*10.0 +0
JGMS ( I ) zE^IS (K2i 5. K T ) *4535, 9 +0
1 = 1 + 1
C -PB-
JLbSCI)=EMIS(KZ,4,KT)*10.0 +0
JGMSC I ) =EMIS(KZ( 4 , KT ) »4535.9 +0
1 = 1 + 1
4050 CONTINUE
GO TO 4200
C
C EMISSIONS ARE IN GMS
4100 00 4150 KT=1 ,4
C -NOX-
J&MS(n=EKIS(KZ,2,KT}*10.0
JLBSt I )sEMI3(KZ, Zi KT) 410,0/453.
1 = 1 + 1
C -HO
JGMSII)ȣM1S(KZ,3,KT)*10.0
JLBS(I)=EMIS(KZ,3,KT)*10.0/453.
1 = 1 + 1
C -CO-
JGMS(I)»EMIS(KZ, J ,KT)*10.0
JLBS(n»EMISf"KZ,l,KT)* 10, 0/453.
1 = 1 + 1
C «PM»
JGMS(I)=EH1S(KZ,5,KT)*10.0
JLBS ( I )=E MIS(K I, 5,KT) •10.0/453.
1 = 1 + 1
C -PB-
JGMSf I)=EHIS(KZ,4,KT)«10.0
JLBSU)=EMIS(KZ,4,KT)«10.0/453.
DATE
LBS
.5
.5
.5
.5
,5
.5
^5
.5
.5
,5
+ 0.5
59 +0.5
+ 0.5
59 +0.5
+ 0.5
59 +0.5
+ 0.5
59+0.5
+ 0.5
59 +0.5
12/17/45
PAGE 0005
1)0057220
00007230
00007200
00007250
00007260
00007270
00007280
00007290
0000730C
0000731 0
00007320
00007330
00007340
00007350
00007360
00007370
00007380
00007J90
00007000
00007110
00007420
00007U30
00007
-------
0151
< ' 1 5 5
0155
'-MSS
0156
0157
015«
01S9
0160
0161
0162
0163
016U
0 165
0166
0167
0168
0169
0170
0171
01 12
01 73
01 7U
?ELt-St 2.0
a 150 CONTINUE
= 7H3 Jp
I2/ I 7 /u5
«r- E o o o ft
u 2 u 0 IK1
T * I , «
•
*H ••' t L.S ( 1 / i, • 1 ) • 1 C'. <
1 t^ISSf 1
-I-C-
ltMlsS(2.*T)iDE>jS(l,J.'10.0 < o . "i
-PB-
IE«I SS(5. KT ) JpF« S ( ' . u , " T 1 « i o , n + O.S
U250 CONTINUE
G(1 TO 2060
5000 WR!TE(NPRNTR,7600)NLAMD
IFCOUTFl" .GT. 0)WRITE(NPRMTR, 7650)NSY"AP
RETURN
FORMATU2. IJ.«F8.6)
FHRk-AT ( IH! , aO >• , ^UMSU^MARY OF EMISSIONS BY , AU , 50X , 5HP4GE .!«/
1 JOX,1UHTIME PERIOD - , I2,UX, 17H8EGINNING HOUR • ,lSi
2 MX, 15HCENGTH (MRS) • (12 / U5» , 1 JH£MI SSIONS IN ,A2,1HS //
J IX, Ail, 1X»5(6X, AU, IOX)/)
FORMAT ( I X( IJ,2X,5(5X,1PE15.9))
FORMAT(«3H1*««**EOF ENCOUNTERED ON LAND USE FILE*****/
i JIHO NUMBER OF RECORDS READ = ,17)
FORMATfUflHO NUMBER OF RECORDS WRITTEN TO SYMAP FILE = .17)
FORMATC 1H1 ,UOX, J8HZONAL EMISSION DENSITY - EMISSIONS IN ,42,
7000
7^00
7510
7600
7650
7700
2 21H ZONE —TIME PER IOD-- , 7 X , U ( Afl , 1 6X ) , Aa )
7710 FOPHAT(/2X, 12. 6x, H, 5H HOUR , « X , 5 ( 5X , 1 PE 1 5 . <» ) 1
7720 FORMATC/2X, I2/6X, 12, 5H MOUR,UX, 5X.1PE15.9 /
1 (/5X , 1HO, I 1 ,5H BEG ,12,5H LEN , I 2 , 5X , 1 PE 1 5 . 9 ) )
7770 FQRMAT(5X, 1H», II ,5H BEG ,I2,5H LEN , I 2 , 5 ( 5X / 1 PE 1 5 . 9 ) )
7flOO FORMAT(I8,60I 1 1 ,2x)
7900 FOR»«»T( 1M1 ,50X, 19HTOTAL EMISSIONS IN , A3 , i HS , 38X , 5MPAGE
1 6X.11HTIME PERIOD, 2X,5(7X, A«,9X)///
2 8X,I2,5H hQUR.«x,5(Sx, 1PE15.9)//
K3X, IH«, I 1 ,5H 8EG ,I2,5h LEN , I 2 , 5 (5X , 1 PE 15 . 9 ) // ) )
0000770C
00007710
00007720
u 0 0 0 7 7 3 0
. 7 ': 0 7 7 i c
•'0007760
0 : 0 C 7 7 7 C
00007780
00007790
,j o o o 7 b fj r<
'--00076 l 0
00007(320
nnnn7«30
000078UO
00007850
•i 0 0 0 7 P b c
0000787Q
00007860
00007890
OU007900
00007910
00007920
00007910
000079UC
00007950
00007960
00007970
00007980
00007990
00008000
0000801 0
00008020
00008030
00008010
00008050
00008060
00008070
00008080
00008090
r r r< r\ $ ) fi r-
00008110
00008120
00008130
000081 UO
00008150
00008160
00008170
-------
FORTRAN I, C1 RELfiS£ 2.0 KC.OUT DATE > 7833. 12/l7/«5 PAGE 0007
FORMATC22HO.**..ENO UF LIST*.*.0 00008J80
V 00008190
-------
TECHNICAL REPORT DATA
(Please read Instructions on the reverse before completing)
NO
902/4-79-002
3. RECIPIENT'S ACCESSION-NO.
4- TITLt ANOSUBTITLE
USER'S GUIDE TO THE NYSDOT AIR QUALITY SYSTEM
5. REPORT DATE
April 1979
6. PERFORMING ORGANIZATION CODE
7. AUTHOR(S)
F. Benesh
Dominic Caracciolo, Jr.
8. PERFORMING ORGANIZATION REPORT NO.
GCA-TR-78-67-G
9. PERFORMING ORGANIZATION NAME AND ADDRESS
GCA/Technology Division
Burlington Road
Bedford, Massachusetts 01730
12. SPONSORING AGENCY NAME AND ADDRESS
U.S. ENVIRONMENTAL PROTECTION AGENCY
Region II
26 Federal Plaza
New York, New York 10007
10. PROGRAM ELEMENT NO.
11. CONTRACT/GRANT NO.
13. TYPE OF REPORT AND PERIOD COVERED
14. SPONSORING AGENCY CODE
15. SUPPLEMENTARY NOTES
16. ABSTRACT
This manual documents the modifications to the New York State Department of Trans-
portation Air Quality System made by GCA/Technology Division. These modifications
principally involved the inclusion of the most recent motor vehicle emission factors,
as represented by the MOBILE1 program.
17.
I.
KEY WORDS AND DOCUMENT ANALYSIS
DESCRIPTORS
Mobile Source Emissions
MOBILE1
b.IDENTIFIERS/OPEN ENDED TERMS C. COSATI Field/Group
'>, LJlS rRlBUTION STATEMEN1
UNLIMITED
19. SECURITY CLASS (This Report)
UNCLASSIFIED
21. NO. OF PAGES
61
20. SECURITY CLASS (This page)
UNCLASSIFIED
22. PRICE
EPA Form 2220-1 (9-73)
55
------- |