&EPA
United States
Environmental Protection
Agency
Office of Mobile Sources
Emission Control Technology Division
2565 Plymouth Road
Ann Arbor, Michigan 48105
EPA 460/3-85-007
August 1985
Air
A Program to Calculate Size Specific
Paniculate Emissions for Mobile
Sources — A User's Guide
-------
EPA 460/3-85-007
A Program to Calculate Size Specific Participate
Emissions for Mobile Sources — A User's Guide
by
Energy and Environmental Analysis, Inc.
1655 North Fort Myer Drive
Arlington, Virginia 22209
Contract No. 68-03-1865
Work Assignment 1
EPA Project Officer: Mark Wolcott
Technical Representative: Robert I. Bruetsch
Prepared for
U.S. ENVIRONMENTAL PROTECTION AGENCY
Office of Mobile Sources
Emission Control Technology Division
2565 Plymouth Road
Ann Arbor, Michigan 48105
August 1985
-------
This report was furnished to the Environmental
Protection Agency by Energy and Environmental
Analysis, Inc., 1655 North Fort Myer Drive,
Arlington, Virginia, in fulfillment of Work
Assignment 1 of Contract No. 68-03-1865. The
contents of this report are reproduced herein as
received from Energy and Environmental Analysis,
Inc. The opinions, findings, and conclusions
expressed are those of the authors and not
necessarily those of the Environmental Protection
Agency. Mention of company product names is not to
be considered as an endorsement by the Environmental
Protection Agency.
Publication No. 460/3-85-007
-------
TABLE OF CONTENTS
Page
1. INTRODUCTION 1-1
2. OVERVIEW OF METHODOLOGY 2-1
3. INPUTS 3-1
3.1 The Control Section 3-1
3.1.2 The Flags 3-1
3.1.2.1 IPROMT 3-1
3.1.2.2 PRO3ID 3-4
3.1.2.3 TAMFLG 3-4
3.1.2.4 VMFLAG 3-4
3.1.2.5 MYMRFG 3-4
3.1.2.6 MISFLG 3-5
3.1.2.7 SPCFLAG 3-5
3.1.2.8 IMFLAG 3-5
3.2 One Time Data Section 3-6
3.2.1 Alternate Catalyst Tampering Factors 3-6
3.2.2 Alternate VMT Mix 3-6
3.2.3 Alternate Registration and Mileage Accumulation ... 3-6
3.2.3.1 Alternate Registration 3-7
3.2.3.2 Alternate Mileage Accumulation Rate .... 3-7
3.2.4 Alternate Misfueling Rates 3-7
3.2.5 Alternate Lead Content in Gasoline 3-8
3.3 Parameter Section 3-8
3.4 Summary of User Inputs 3-8
3.5 Program Data File 3-8
-------
TABLE OF CONTENTS (cont)
Page
it. OUTPUTS 4-1
4.1 Report 4-1
4.2 Diagnostic and Error Messages 4-2
5. PROGRAM STRUCTURE AND NOTES 5-1
5.1 The Program I/O Unit Assignments 5-1
5.2 Program Structure 5-1
5.3 Common Blocks 5-1
5.4 Global Variables 5-1
5.5 Other Program Notes
6. EXAMPLES 6-1
Appendix A: Source Code Listing .'.. A-l
Appendix B: Program Data File B-l
-------
I. INTRODUCTION
This report is provided as documentation to a program which calculates the
total areawide lead and particulate emissions from gasoline and diesel
fueled on-road vehicles, trucks and motorcycles for particle sizes up to 10
microns. The program provides a large amount of default (normal) data
and also many options to tailor the conditions to each locality. The period
of interest spans 1980 to 1995. (Years later than 1995 can be calculated but
there is no data specific to years later than 1995).
The methodology and sources of data for this program are described in
detail in the report Size Specific Total Particulate Emission Factors for
Mobile Sources prepared August 198^, for EPA by EEA. That report describes
the calculation of particulate emission factors for each of the six vehicle
classes used in this program. The vehicle classes are: 1) light duty
vehicles, 2) light duty trucks 6000 pounds or less, 3) light duty trucks more
than 6000 pounds, <0 heavy duty gasoline trucks, 5) heavy duty diesel trucks,
and 6) motorcycles. It also provides the emission factors for each major
type of particulate (lead, organic, sulfate, and diese; Composite). A brief
summary of the methods used for the calculation is given in Section 2.
The program inputs and format are designed to parallel those found in
MOBILE 3. The normal user inputs include the scenario year, the vehicle
speed, the driving conditions and the particle size cutoff (i.e., the largest
particle size to be included in the total emissions). By setting control
flags, the user can elect to replace the default data used in calculations
with locality specific data. These include tampering rates for catalyst
equipped light duty vehicles and trucks, misfueling rates, the VMT mix,
mileage accumulation rates, registration fraction by age, the lead content
of gasoline and the implementation of an I/M program. A description of
each input is given in Section 3.
The output report gives a breakdown of particulate emission factors type
for each vehicle class as well as totals and fractions of the total that each
1-1
-------
participate contributes. The output report and any error and diagnostic
messages that are generated are described in Section 4.
The program structure is described in Section 5. A dictionary of subroutines
and common blocks of interest to programmers wishing to make changes to
the program is included. Users may skip this section.
Section 6 provides example runs of the program with various options. Use
of the program in both batch and interactive modes are shown.
The appendices include a source code listing and a listing of the data file
containing the default data.
1-2
-------
2. OVERVIEW OF METHODOLOGY
Particulate emissions consists of lead, organic and sulfate emissions. Dif-
ferent types of vehicles emit varying amounts of each participate type
depending upon fuel type, emission control technology and operating mode.
Lead participates consist mainly of lead salts rather than lead alone.
Organic emissions include elemental carbon and soluble organics. Particu-
late emissions from diesel vehicles are a composite particulate having
elements of the three other types but consisting largely of organic
emissions.
The calculations seek to model the particulate emissions from all types of
on-road vehicles. The procedure is to calculate the emissions factor for
each type of particulate from each vehicle class and then sum over all
vehicle classes to get the total emission factor in grams per mile traveled.
For any given year the procedure may be expressed mathematically
6
EFpm,n,s =Z) ti,nEFi,n,s + EFbrakes (MB> + EF tires
i=l
where EFpm)n>s _ sjze specific all-vehicle class total particulate
emission factor on January 1 of calendar year n at
vehicle speed s (g/mile)
i = vehicle class designator; l=light-duty vehicles (LDV),
2=light-duty trucks I (LDT1), 3=light-duty trucks II
(LDT2), 4=heavy-duty gas vehicles (HDGV), 5-heavy-
duty diesel vehicles (HDDV), 6=motorcycles (MC)
s = vehicle speed; avg. Federal Test Procedure (FTP) =
19.6, avg. Sulfate Emissions Test (SET) = 34.8
(miles/hr); (Note: The FTP and SET are driving cycles
used for the determination of emission factors.)
tj n = area travel fraction of vehicle class i in calendar year
n
EFi n s = particulate emission factor for vehicle class i in
calendar year n at vehicle speed s (g/mile)
EFbrakes = airborne brake wear particulate emission factor com-
ponent = 0.002 grams/mile; this emission factor
2-1
-------
component is assumed to be the same for all vehicle
classes, vehicle speeds and calendar years (all i,s, and
n) due to lack of separate information for each i,s, and
n
= airborne tire wear participate emission factor component
= 0.002 grams/mile; this emission factor component is
assumed to be the same for all vehicle classes, vehicle
speeds and calendar years (all i,s, and n) due to lack of
separate information for each i,s, and n
MB = fraction of airborne particles less than a user-specified
size cutoff (0.1=10 m) that are attributable to vehicle
brake wear, from Table 2-20 or Figure 2-4
Calculating the emissions factor for each vehicle class requires knowledge
of emissions factors and numbers of each type of vehicle within the class.
For instance, light-duty vehicles includes leaded, unleaded and diesel vehicles
each with different emission characteristics.
The emissions factors required are given in the report Size Specific Total
Particulate Emission Factors of Mobile Sources prepared for EPA by EEA.
The emissions factors are functions of vehicle class, model year, fuel type,
catalyst type (for sulfates), vehicle speed (for lead and sulfates), misfueling
and catalyst tampering rates (for unleaded vehicles). Changes in emissions
factors from model year differences occur because of changes in fuel economy,
emission control technology and the shift from leaded to unleaded gasoline.
Each vehicle class is modeled as a fleet of vehicles including vehicles up to
20 years old. The fraction of vehicles by age within each class favors newer
vehicles since older vehicles wear out and are taken out of service. Newer
vehicles also have higher mileage accumulation rates than older vehicles.
Misfueling rates increase with age since there is more opportunity to misfuel
the vehicle over a long period of time. For each vehicle class the registra-
tion fraction, mileage accumulation rate and misfueling rate are assumed
to vary year by year up to 20 years of age. Vehicles 20 years old and older
are modeled as a single age.
For a given scenario year, the total emission factor from a vehicle class is
a sum of the emissions factors over the previous 20 model years weighted
by the fleet fraction of vehicle class. The total emission factor from all
2-2
-------
vehicle classes are further weighted by the areawide travel fraction (i.e.
the ratio of vehicle miles traveled by a vehicle class to the total traveled
by all vehicle classes in the area).
The final results of the calculation are given in grams per mile traveled.
The total emissions can be calculated by multiplying the composite
emission factor by the total number of miles traveled by all classes of
vehicles.
2-3
-------
3. INPUTS
Inputs to the program allow the user to specify and change the scenario
conditions to fit the locality being modeled. Like MOBILE 3, the various
input options and scenario options are set by using control flags. The sec-
tions of user provided data are divided into a control section, one time data
section and parameter section just like MOBILE 3. Section 3.1 describes
the control section, Section 3.2 describes the one time data section and
Section 3.3 describes the parameter section.
3.1 THE CONTROL SECTION
The control section contains the control flags necessary to select the
various options available to the user. All of the flags must be present and
have valid values for the program to run. Inputting "1" for each flag will
result in the program using the default values and options built into the
program. Table 3-1 summarizes the control flags.
3.1.2 The Flags
3.1.2.1 IPROMT
This flag works exactly as it does in MOBILE 3. It is used to select whether
prompts are required and whether horizontal or vertical format is used for
further flag inputs. Horizontal format has all the control flags in a single
record. Vertical format has each flag in a separate record. No prompt is
issued for this input since its purpose is to determine if prompts are required.
(Format = II). The options are:
1 - No prompting; use vertical format
2 - Program prompts; use vertical format
3 - No prompting; use horizontal format
b - Program prompts; use horizontal format
3-1
-------
TABLE 3-1
FLAGS CONTROLLING INPUT AND EXECUTION
Record
Number
1
Variable
IPROMT
IOUNEW
Description
Flag for prompting of
input
Values for output units
1) X=report
2) Y=diagnostic error
reports
3) Z=prompt messages
(Default for these units
is 6)
FMT
Value
1XYZ
2XYZ
3XYZ
fXYZ
Action
No prompts
vertical format
prompt
vertical format
no prompts
horizontal format
prompts
horizontal format
PROJID
80 character title used
as output header
20A4
TAMFLG Flag for optional input
of tampering rates
II
l=normal
2=user data in one-
time section
VMFLAG Selects optional user
supplied VMT mix among
vehicle types
II
1=MOBILE 3 VMT
mix
2=user VMT mix for
each scenario in
parameter section
3=user VMT mix for
all scenarios in
one-time section
-------
Record
Number
Variable Description FMT
MYMRFG Selects optional mileage II
accumulation and regis-
tration versus age
Value
Action
l=normal
2=user supplied
mileage accumulation
rates in one-time
section
3=user supplied
registration in
one-time section
4=user supplied
mileage accumula-
tion and registra-
tion rates in
one-time section
MISFLG
VjJ
Flag for optional input
of misfueling rates
II
l=normal
2=user data for LDV,
LDT1, LDT2 in one-
time section
3=use "simple" rates
SPCFLG
Selects optional lead
content in leaded and
unlead gasoline
II
l=normal MOBILE 3 lead
content
2=user supplied lead
content
IMFLAG
Indicates I/M program
II
l=no I/M
2=I/M
-------
IOUNEW
This variable allows the user to change the output units for the results,
messages and prompts. The default value for all of these is 6. IOUNEW
must be part of the same record as IPROMT and immediately follows it.
Value = XYZ Format = 3A1
Options:
X = Results unit
Y = Diagnostic messages
Z = Prompting messages
3.1.2.2 PRO3ID
This is an 80 character title to be used as a header.
3.1.2.3 TAMFLG
This flag allows the user to input alternate catalyst tampering rates if
desired for light duty vehicles and trucks. If user supplied tampering rates
are selected, the user needs to specify tampering rates only for the case
selected by IMFLAG (i.e., I/M program or no I/M program).
Options: 1 - normal tampering rates
2 - user supplied tampering rates in one-time section (3
values)
3.1.2.4 VMFLAG
This flag selects the areawide vehicle miles traveled fraction. The default
fractions come from MOBILE 3. User supplied fractions should have one
value for each vehicle class.
Options: 1 - use MOBILE 3 VMT mix
2 - user supplied VMT mix for each scenario
3 - user supplied VMT mix for all scenarios in this run (one-
time data)
3.1.2.5 MYMRFG
This flag selects the mileage accumulation rate and registration versus
age. The rates and registration fractions cover ages up to 20 years and are
-------
specified for each vehicle class. User supplied data can be used to replace
the default data for any or all vehicle classes.
Options: 1 - normal mileage accumulation rates and registration
fraction
2 - user supplied mileage accumulation rates (one time data)
3 - user supplied registration fraction (one time data)
4 - user supplied mileage accumulation rates and registra-
tion fractions (one time data)
3.1.2.6 MISFLG
This flag selects the misfueling rates desired. The program normally uses a
set of misfueling rates that varies according to vehicle class and age. The
user supplied rates are generated for each vehicle class from a slope and
intercept with mileage as the dependent variable.
The simple rates are a set of misfueling rates that have been averaged over
vehicle age; hence, only the vehicle class matters.
Options: 1 - normal misfueling rates
2 - user supplied slope and zero mile intercept
3 - "simple" misfueling rates, vehicle class dependent only
3.1.2.7 SPCFLAG
This flag selects normal or user supplied lead content of leaded and unleaded
gasoline. The lead content can be changed for any number of years with
the default lead content used for those years previous to the changes.
Options: 1 - normal MOBILE 3 lead content
2 - user supplied lead content of gasoline
3.1.2.8 IMFLAG
This flag indicates whether an inspection maintenance program is in effect
in the area.
Options: 1 - no I/M program
2 - I/M program in effect
3-5
-------
3.2 ONE TIME DATA SECTION
This section contains data only if the flags in the control section have been
set to expect user supplied data here. This data is input only once and is
used for all scenarios of the run. The data in this section concerns locality
dependent variables.
The one time data section is included only in the following cases in the
order listed.
1) The user wishes to input local catalyst tampering rates (TAMFLG = 2)
2) The user wishes to use an alternate VMT mix for all scenarios
(VMFLAG = 3)
3) The user wishes alternate registration and/or mileage accumulation
rates (MYMRFG = 2, 3 or 4)
4) The user wishes to input alternate misfueling rates (MISFLG = 2)
5) The user wishes alternate lead content in gasoline for a range of
years (SPCFLG = 2)
The order of the data in this section is the same as the order of flags in the
control section.
3.2.1 Alternate Catalyst Tampering Factors
The data to be entered are the tampering rates for LDV, LDT1 and LDT2.
This data should be appropriate for the IMFLAG chosen, i.e. if there is no
I/M program in effect, enter the rates for no I/M program. The format is
3(F4.3,1X).
3.2.2 Alternate VMT Mix
The data to be entered is the areawide vehicle miles traveled fraction for
each vehicle class. The sum of the fractions should equal one. If not, a
warning message will be printed and the VMT fractions will be normalized.
There is one record with format 6(F5.3,1X).
3.2.3 Alternate Registration and Mileage Accumulation
This section allows changes to the registration and mileage accumulation
rates for any or all vehicle classes. Since LDV, LDT1 and LDT2 include
diesel vehicles with different registration and mileage accumulation rates,
3-6
-------
the fuel type must also be specified when indicating which vehicle class the
user wishes to change. For HDGV, HDDV and MC the fuel type that is
input is ignored. The registration and mileage accumulation rates for leaded
and unleaded vehicles are automatically set to be equal; when the user
chooses leaded or unleaded as the fuel type, the registration and mileage
accumulation for the type that was not chosen are set to the same value as
the type chosen.
3.2.3.1 Alternate Registration
The data is entered in three records for each vehicle class. The first record
has format (I1,1X,I1) where the data to be entered is the vehicle class (1-6)
and fuel type (1 or 2 = gasoline, 3 = diesel). The second record contains the
registration information for the first ten years of age in format (10 F5.3).
The third record contains the registration information for years 11 to 20
years or more (11-20/) in format 10 F5.3.
There can be as many sets of three records as desired. After the last set,
there should be a card with three zeroes ('000') to signal the end of
registration information.
3.2.3.2 Alternate Mileage Accumulation Rate
The data is entered in three records for each vehicle class in the same
fashion as the registration data. The first record has format (I1,1X,I1) where
the data is be entered is the vehicle class (1-6) and fuel type (1 or 2 =
gasoline, 3 = diesel). The second record contains the mileage accumulation
rate for the first 10 years in format 10(I5,1X). The third record contains the
mileage accumulation rate for the years 11 to 20 or more (11-20/) in format
There can be as many sets of three records as desired. After the last set,
there should be a card with three zeroes ('000') to signal the end of the
mileage accumulation information.
3.2.4 Alternate Misfueling Rates
The data to be entered is the zero mile misfueling rate and the slope for
the four vehicle classes LDV, LDT1, LDT2 and HDGV1. A negative slope or
3-7
-------
intercept will produce an error message. There is one record for each class
with format (F8.6,1X,F8.6).
3.2.5 Alternate Lead Content in Gasoline
The data entered replaces the default lead content information for future
years. For years prior to the starting year, the default information is used.
The first record gives the starting year for the user supplied data and the
last year in format (I4,1X,I4). Subsequently, there is one record for each
year giving the leaded and unleaded lead content in grams per gallon in
format F5.3,1X,F5.3. The data on the last record is used for the last year
and all later years.
3.3 PARAMETER SECTION
This section contains data that is changed for each scenario. There can be
as many scenarios as required in a single run of the program. Each scenario
consists of four records (five, if VMFLAG = 3). The first record has the
year of the scenario in lb format. The second record has the vehicle speed
in F4.1 format. The third record specifies the driving conditions in II format
where a T indicates transient or cyclic driving and a '2' indicates a steady
cruise. The fourth record contains the particle size cutoff in mirons in
F5.2 format. The particle size cutoff is the largest size of particle to be
included in the calculations. None of the four records have default values;
all of them must be correctly specified. If VMFLAG is equal to three there
is a fifth record to each scenario giving the VMT mix. The format is
described in Section 3.2.2.
3.4 SUMMARY OF USER INPUTS
The user inputs are summarized in Table 3-2.
3.5 PROGRAM DATA FILE
The default data file contains all of the data necessary to run the program
under standard conditions except for scenario year, vehicle speed, driving
conditions and particle size cutoff. The subroutine EMFIN and associated
subroutines reads the file in its entirety. User supplied data which supplants
the default data is read in afterwards from the one time data section or
parameter section.
3-8
-------
TABLE 3-2
INPUT RECORD SEQUENCE SUMMARY
VjJ
vD
Input Record Sequence
1 IPROMT Flag
1 PROJID
1 TAMFLG+
1 VMFLAG+
1 MYMRFG+
I MISFLG+
1 SPCFLG+
1 IMFLAG+
(1 Tampering record if TAMFLG=2)
(1 VMT mix record if VMFLAG=3)
(3/class Alternate registration if MYMRFG=3 or 4)
(1 Alternate registration end card if MYMRFG=3 or 4)
(3/class Alternate mileage accumulation if MYMRFG=2 or 4)
(1 Alternate mileage accumulation end card if MYMRFG=2 or 4)
(4 Misfueling records if MISFLG=2)
(1 Starting and ending years for lead content of gasoline if SPCFLG=2)
(- Lead content of gasoline cards if SPCFLG=2)
Section
Control
Control
Control
Control
Control
Control
Control
Control
One-time
One-time
One-time
One-time
One-time
One-time
One-time
One-time
One-time
-------
1 Year of scenario Parameter*"1"
1 Vehicle Speed Parameter
1 Driving conditions Parameter
1 Particle size cutoff Parameter
(1 VMT mix if VMFLAG=2) Parameter
Notes +: May be on single record if horizontal format is chosen
++: Parameter section may be repeated as many times as desired
-------
The default data file structure is shown in Table 3-3. If changes to the
data file are made care should be taken to change the corresponding array
dimensions in the program. The program comments include references to
tables in the original report to aid making changes.
3-11
-------
TABLE 3-3
DEFAULT DATA FILE SEQUENCE
1. LDV LDT1 LDT2 Organic Emission Factors
a) no. of records - Format (12)
b) catalyst type, fuel type, emission factor - Format=(M, 3X, II, IX, II,
IX, F5.3)
see note 1
2. LDV LDT1 LDT2 Sulfate Emission Factors
a) no. of records - Format (12)
b) year, catalyst type, fuel type, vehicle speed, emission factor -
Format=(I
-------
8. LDV LDT1 Lead Emission Factors
a) no. of records - Format (12)
b) year, fuel type, emission factor - Format=(M, IX, II, IX, F5.3)
see note 1
9. LDT2 Lead Emission Factors
a) no. of records - Format (12)
b) year, fuel type, emission factor - Format=(I4, IX, II, IX, F5.3)
10. Lead Content of Gasoline by Year
a) no. of records - Format (12)
b) year, lead content of leaded, lead content of unleaded - Format=(I4,
IX, II, IX, F5.3)
see note 1
11. LDV Fraction of Fleet by Emission Control Type
see note 2
12. LDV Fleet Sales Fractions by Fuel Type
see note 3
13. LDV Travel Weighting (Registration and Mileage Accumulation Rates)
see note ^
1*. City/Highway Fuel Economy
a) no. of records - Format (Iw)
b) year, fuel economy for LDV, LDT1, LDT2, HDGV1, HDGV2, HDGV
15. Fuel Economy Speed Corrections
a) no. of records - Format (12)
b) speed, transient driving speed correction, steady driving speed
correction
16. LDT1 Fraction of Fleet by Emission Control Type
see note 2
17. LDT1 Fleet Sales Fraction by Fuel Type
see note 3
18. LDT1 Travel Weighting (Registration and Mileage Accumulation Rates)
see note 4
19. LDT2 Fraction of Fleet by Emission Control Type
see note 2
20. LDT2 Fleet Sales Fraction by Fuel Type
see note 3
3-13
-------
21. LDT2 Travel Weighting (Registration and Mileage Accumulation Rate)
see note 4
22. LDDT Travel Weighting (Registration and Mileage Accumulation Rates)
see note 4
23. HDGV Fleet Sales Fraction by Fuel Type
see note 3
24. HDGV Travel Weighting (Registration and Mileage Accumulation Rates)
see note 4
25. HDDT Registration
a) no. of records - Format (12)
b) registration fraction of total fleet by age (repeat)
26. MC Travel Weighting (Registration and Mileage Accumulation Rate)
see note 4
27. Simple Misfueling Rates
I/M program in effect, no I/M program
5 records for LDV, LDT1, LDT2, HDGV, MC
Format=(5X, Ft.2, IX, F4.2)
28. Misfueling by Age
a) no. of records - Format (12)
b) each record contains misfueling rates for:
LDV with I/M
LDV w/o I/M
LDT1 with I/M
LDT1 w/o I/M
LDT2 with I/M
LDT2 w/o I/M
HDGV1 with I/M
HDGV1 w/o I/M
Format=3X, 8 F4.2
29. Particle Cumulative Fraction by Size
a) Particle Sizes
b) particle cumulative fractions
report a) and b) for leaded, unleaded
with catalyst, unleaded w/o catalyst,
diesel, brakes
Format - 5(F5.2, IX)
-------
30. Low Altitude HDDV Conversion Factors
a) no. of records
b) model* year, conversion factor report as specified in 30a)
31. Fraction of Lead Exhausted
non-misfueled, 1975-1980 misfueled
1981+ misfueled vehicles
Format = 3(F4.2, IX)
32. LDV, LDT1, LDT2 Tampering Rate (Catalyst Removal)
a) LDV w/o I/M
LDT1 w/o I/M
LDT2 w/o I/M
b) LDV with I/M
LDT1 with I/M
LDT2 with I/M
Format = 3(F5.3, IX)
33. HDDT Mileage Accumulation by Age and Weight Class
a) no. of records - Format (12)
b) each record contains mileage accumulation for class 2B, Light,
Medium and Heavy HDDT - Format=(3X, 3(15, IX), 15)
repeat as specified in a)
34. HDDT Fleet Sales Fraction by Year and Weight Class
a) no. of records - Format (12)
b) each record contains year and sales fraction for 2B, light, medium and
heavy HDDT - Format=(I4, IX, 4(F5.3, IX))
35. VMT Mix
1 record with the fraction of VMT miles traveled by LDV, LDT1, LDT2,
HDGV, HDDT and MC. This data has been taken from MOBILE 3 except
that LDDT VMT have been split evenly between LDT1 and LDT2 since
MOBILE 3 does not distinguish the two classses of LDDT.
3-15
-------
Note J: repeat this record as many times as specified in a). The "year" is the
earliest year for which the emission factor applies.
Note 2: These records are read by the subroutine GTEMCT. Each set of
records consists of:
a) This record gives the number of years each of the following records
contains. The first year is automatically 1974 and earlier, the
second is 1975, the third is 1976 etc. the last year is that year and
later. Format =12
b) Fraction of leaded, no catalyst by year
c) Fraction of unleaded, no catalyst by year
d) Fraction of unleaded, oxid. catalyst/no air by year
e) Fraction of unleaded oxid. catalyst/air by year
f) Fraction of unleaded 3 way catalyst by year
g) Fraction of unleaded 3 way + oxid. catalyst by year
h) Fraction of diesels by year
j) Fraction of unleaded with catalysts by year
k) Fraction of unleaded with catalysts/no air by year
1) Fraction of unleaded with catalysts/air by year
Format = 15(F5.3, IX)
Note 3: Each set of records is read by subroutine FLTSLE and contains:
a) number of records Format = 12
b) model year, unlead sales fraction, leaded sales fraction, diesel
sales fraction
repeat b) as specified in a)
Format = 14, IX 3(F5.3, IX)
Note 4: Each set of records is read by subroutine TRAVWT and contains 20
records. Each record contains age, registration fraction and mileage
accumulation rate by age. Format = (14, IX, F5.3, IX, 15)
3-16
-------
. OUTPUTS
Output from the program consists of prompting messages, error messages
and the report on results. The report includes information on the scenario
and the calculated emissions from each type of vehicle. The diagnostic
messages inform the user about errors and warn about data errors that are
non-fatal. Prompt messages are available if IPROMT is set to 2 or 4 in the
control section. These messages prompt the user to input data in the correct
sequence. In the CONTROL section there are two formats vertical flag
input(2) and horizontal flag input (4). The prompt messages are the same in
either case.
4.1 REPORT
Example reports are shown in Section 6. The report is divided into three
sections. The first section provides scenario information. The second
section gives the emission factors and the last section gives the fraction of
the total particulates from each vehicle class and the fractions within each
class.
The scenario information reflects information entered in the control section
and parameter section. The first item is the project identification (title)
which was entered as part of the control section. Next the scenario year
and vehicle speed is presented. The flag settings and the driving conditions
are shown next. The flag settings indicate whether the user provided data
or whether the normal default data was used. In the case of the I/M flag,
where there is no default value, the flag indicates whether there is an I/M
program in effect.
In the second section the results of the calculations are shown. The results
show the contribution of the four types of particulates to exhaust emissions
from each vehicle class. For light duty vehicles, both gasoline and diesel
vehicles are included in the same classes. However, it is easy to separate
the contributions from diesels, since only diesel vehicles emit diesel par-
ticulates. The totals show the total grams per mile of all particulates
-------
emitted by vehicles of each class. The VMT weighted totals are the same
totals weighted by the areawide VMT travel fraction. These values are
important for gauging the relative contribution from each class since there
are different numbers of vehicles in each class. Particulates from brake
and tire wear are added to the VMT weighted totals from exhaust emissions
to get the total VMT weighted emissions.
The last section shows the fractional contribution of each paniculate to
the total for each vehicle class.
4.2 DIAGNOSTIC AND ERROR MESSAGES
Messages are of three types: error messages, warning messages and comments.
Error messages can be fatal or non-fatal. Fatal errors cause the program
to halt. Fatal errors are caused by user supplied data that the program
does not understand or which is invalid and cannot be corrected by the
program. Non-fatal errors are errors that do not stop execution. There
occur where values are out of bounds. In some cases, the program adjusts
the out of bound values to be in bounds. The program execution will halt
when there are 50 non-fatal errors. Scenarios with non-fatal errors should
be carefully checked for errors. Comments are useful statements about
options taken. They do not affect the calculations and are present for
informational purposes only. The meanings of the messages are presented
in Table 4-1 for fatal errors, Table 4-2 for non-fatal errors and Table 4-3
for warnings and comments. The error number given in the tables is the
code passed to the QUITER subroutine which writes the message.
4-2
-------
TABLE 4-1
FATAL ERROR MESSAGES
Error
Number
1 "Message Code '#m
This is an unknown of out of bounds message it should never occur. // is the
unknown message number
2 "IPROMT Flag Out Of Bounds - Prompt Mode Type (1 to 4) Must be
Corrected Before Your Next Run"
This message occurs if IPROMT is not equal to 1, 2, 3 or k
3-8 "'Value1 Out of Bounds for 'Name of Flag' Flag"
This message indicates that the flag named has a value that is out of
bounds. It should be corrected for the next run.
16 "Range of Years for Lead Content by Year Cannot be Negative or Zero"
This message occurs when the user has set SPCFLG=2 and incorrectly
specifies the range of years on the first card of the lead content of
gasoline data.
17 "'Value1 intercept must be positive"
This error occurs when MISFLG=2 and the user has entered a negative zero
mile misfueling rate
22 "You Must Correct All Flags Before You Run Again"
This message occurs when one or more of the control flags is out of bounds.
It ends the run. This message occurs in conjunction with errors 3-8.
23 "End of File in OUTGEN Data File"
This error occurs when there is no more data in unit IOUGEN when the
program was expecting data. This message will occur only when using a
data file for user supplied data.
31 "End of File in IOUDAT Data File"
This error occurs when the program finds an end of file in the default data
file. This error occurs only when the default data has been changed.
-------
35 "Driving Condition Flag Must Be 1 or 2 Cannot be "Value"1
The user has incorrectly entered the driving conditions. Entering
l=transient or cyclic, a 2=cruise
36 "Year Cannot Be 'Value'—Must Be Later Than 1979"
The scenario year must be later than 1979
37 "No Such Vehicle Class 'Value' For Registration Data"
This error occurs when MYMRFG=3 or 4 and an incorrect vehicle class has
been entered. The vehicle class must be a number between one and six.
38 "No Fuel Type 'Value' for Registration Data"
This error occurs when MYMRFG=3 or 14 and an invalid fuel type has been
entered. The fuel type must be between one and three. (l=leaded,
2=unleaded, 3=diesel).
39 "No Vehicle Class 'Value* for Mileage Accumulation Data"
This error occurs when MYMRFG=2 or 4 and an invalid vehicle class has
been entered. The vehicle class must be between one and six (1=LDV,
3=LDT2, 4=HDGV, 5=HDDT, 6=MC).
40 "No Fuel Type 'Value' for Mileage Accumulation Data"
This error occurs when MYMRFG=2 or 4 and an invalid fuel type has been
entered. The fuel type must be between one and three (1 or 2=gasoline,
3=diesel)
4-4
-------
TABLE it-2
NON-FATAL ERROR MESSAGES
Error
Number
10 " 'Value' Out of Bounds for VMT Mix (0. to 1.)
The user has specified an incorrect VMT fraction. The VMT weighted totals
in the results will be incorrect.
11 " 'Value' Sum of VMT Mix is not Equal to 1."
The user has specified one or more incorrect VMT fractions. The sum of
the fractions does not add up to one. The VMT weighted totals in the
results will be incorrect.
12 " 'Value' Negative Model Year Mileage"
The user has set MYMRFG=2 or 4 and has entered a negative mileage accumu-
lation rate. The results for this vehicle class and fuel type will be
incorrect.
13 " 'Value1 Negative Model Year Registration"
This message occurs when MYMRFG=3 or 4 and a negative registration
fraction has been entered. The results for this vehicle class will be
incorrect.
19 " 'Value' Speed Must be Positive"
The speed for this scenario is less than zero. The program will take the
absoute value and continue execution.
28 "Lead Content of Leaded Gas Cannot Be 'Value' for Year 'Value of Year' "
This message occurs when the user has set SPCFLG=2 and incorrectly speci-
fies a negative lead content for the given year. The program makes the
value positive and continues.
29 "Lead Content of Unleaded Gas Cannot be 'Value' for Year 'Value of Year' "
This message occurs when the user has set SPCFLG=2 and incorrectly spec-
ifies a negative lead content for the given year. The program makes the
value positive and continues.
30 " 'Value' Tampering Rate Must be Positive or Zero"
This message occurs when TAMFLG=2 and the tampering rate is negative.
The program sets the tampering rate to zero and continues execution.
4-5
-------
32 "Argument Out of Bounds in Function 'Name of Function1"
This error indicates a bad argument was passed to the named function.
This error should not occur. Contact technical assistance.
-------
TABLE 4-3
WARNING AND COMMENT MESSAGES
Message
Number
9 " 'Value1 Speed Reduced to 55 M.P.H. Maximum"
The speed for the scenario was reduced to 55 mph.
14 " 'Value1 Registration with Zero Mileage"
This error will occur only when MYMRFG=2, 3 or 4. There is a vehicle
class with a registration fraction but no mileage accumulation.
15 " 'Value' Mileage with Zero Registration"
This error will occur only when MYMRFG-2, 3 or 4. There is a vehicle
class with mileage accumulation but no registration fraction.
18 " 'Value' Negative Slope for Ageing Vehicles"
This error will occur only when MISFLG=2 and the user has entered a nega-
tive slope, (i.e. less misfueling with mileage accumulation). Where the
derived misfueling rate is negative, it is set equal to zero.
20 " 'Value1 MYR Sum not = 1. (Will Normalize)"
This message will occur only if MYMRFG=3 or 4 and the sum of the
registration fractions do not sum to one. The program will renormalize the
fractions so that they do sum to one.
21 " 'Value1 Speed Increased to 5 M.P.H. Minimum"
The scenario speed was less than 5 mph. The program increases it to 5 mph
before proceeding.
26 " 'Value' Particle Size Cutoff Increased to 0.1 Microns"
The scenario particle size cutoff was less than 0.1 microns. The program
set the particle, size cutoff to 0.1 microns before continuing execution.
27 " 'Value' Particle size Cutoff Decreased to 10 Microns"
The scenario particle size cutoff was greater than 10 microns. The
program set the particle size cutoff to 10 microns before continuing
execution.
4-7
-------
33 "Vehicle Speed = 'Value' Less Than 19.6 M.P.H. for Sulfate Cal."
This message simply warns that the sulfate calculations for light duty
vehicles are being extrapolated.
34 "Vehicle Speed = 'Value' Greater than 34.8"
Same as message 33.
-------
5. PROGRAM STRUCTURE AND NOTES
This section describes the program I/O unit assignments. It also documents
the setup of the default data file which is read in prior to each run.
5.1 The Program I/O Unit Assignments
The MAIN program contains the default input and output unit numbers.
The input unit numbers are stored in the variables IOUDAT and IOUGEN.
IOUDAT is the unit number assigned to the deault data file. It is normally
unit 4. IOUGEN is the unit designated for user input. It is normally unit 5.
These assignments cannot be changed unless the MAIN program is altered.
The output unit numbers are stored as IOUREP, IOUERR and IOUASK.
They are the output unit for reports, error and warning messages and user
prompts respectively. They are all normally assigned to unit 6. The user
can change any or all output unit numbers by using the IPROMT flag to
specify new numbers (see control section inputs). All of the I/O unit
assignments are stored in the common block IOUCOM.
5.2 PROGRAM STRUCTURE
Figure 5-1 diagrams the program structure in terms of subroutines. The
subroutines are described in Table 5-1. Each subroutine is called by the
routine to the left of it on the diagram.
5.3 COMMON BLOCKS
The common blocks and the subroutines that use them are listed in Table
5-2.
5.
-------
FIGURE 5-1
PROGRAM STRUCTURE
CONSEC
OUTHDR
EMFIN
GTEMCT
FLTSLE
TRAVWT
ONESEC
GETVMT
GETLCG
MAIN
PARSEC
CALC
OUT
TWFCAL
PARSIZ
— GETVMT
5ULFAT
— DIESEL
— NLCONS
— GETPB
— GETPBCG
— GETCS
^"•CT^fC
GETEFS
— NLCONS
— GETCFE
—
— GETEFO
— NLCONS
GETCFE
OUTHDR
-------
TABLE 5-1
PROGRAM MODULE DICTIONARY
MAIN Controls overall execution of the program
CONSEC Reads and verifies control flags
EMFIN Reads all default data
ONESEC Reads and verifies one time data
PARSEC Reads parameters for each scenario
GTEMCT Reads fractions of given vehicle class with various types of emission
control technology
FLTSLE Reads fleet sales fractions by fuel type for a given vehicle class
TRAVWT Reads and verifies registration and mileage accumulation by age for
a given vehicle class
GETVMT Reads and verifies VMT fractions by vehicle class for area (also
known as even travel fraction)
GETLCG Reads user supplied lead content of gasoline by year
TWFCAL Calculates travel weighting function for all vehicle classes and sets
default fleet sales functions for those not previously read
PARSIZ Calculates fraction of each particulate type which is below particle
size cutoff
LEAD Calculates total lead emission for given vehicle class, year and fuel
type
NLCONS Finds unleaded emissions factors for given year and vehicle class
GETPB Finds leaded emissions factor for given vehicle class and year
GETPBCG Finds amount of lead in gasoline for given year and fuel type
GETCS Finds speed dependent correction factor for given speed
GETCFE Finds combined city/highway fuel economy (miles per gallon) for
given year and vehicle class
ORGANC Calculates total organic emissions for given vehicle class, year and
fuel type
5-3
-------
GETEFO Finds organic emission factor for given year, vehicle class, catalyst
and fuel type
SULFAT Calculates total sulfate emission for given vehicle class, year, speed
and fuel type
GETEFS Finds sulfate emission factor for given year, vehicle class, catalyst
type, fuel and speed
DIESEL Calculates diesel particulates for given vehicle class
OUT Outputs results
OUTHDR Outputs project name
-------
TABLE 5-2
COMMON BLOCKS
Name of
BLOCK
Description
1. ANSWER Contains particulate totals by type and vehicle class
PBSUM (6) = lead totals
ORGSUM (6) = organic particulate totals
SULSUM (6) = sulfate particulate totals
DIESUM (6) = diesel particulate totals
CLSSUM (6) = particulate totals by class
BRKSUM = brake wear particulate
TIRSLJM = tire wear particulate
SUMALL = sum of all particulates
Block Used by CALC and OUT
2. CONST Contains fraction of lead burned that is exhausted
FPB = fraction of lead exhausted by non-catalyst vehicles
FPBE = 1975-1980 misfueled vehicles
FPBL = 1981+ misfueled vehicles
EMCT
Block used by LEAD, EMFIN
Contains fraction of fleet with emission control devices and catalyst
removal rates
FPBEMC (15,6)
FNLEMC (15,6,8)
FDIEMC (15,6)
IEMCYR (15)
IEMCT
TAMPER (3,2)
fraction of leaded vehicles
fraction of non-leaded vehicles with catalyst
type
fraction of diesel vehicles
model year
number of model years
catalyst removal rates for LDV, LDT1 and
LDT2
For first four variables
1st index = year pointer
2nd index = vehicle class
1=LDV, LDDV
2=LDT1, LDDT1
3=LDT2, LDDT2 6=MC
3rd index = catalyst type (non-leaded only)
1=LDV, LDDV 4=HDGV1, HDDT
2=LDT1, LDDT 5=HDGV2
3=LDT2 6=MC
4=HDGV, HDDT
5=HDGV2
5-5
-------
For TAMPER
= vehicle class
= I/M flag
Block used by FREMCT, GTEMCT, ONESEC, NLCONS, SULFAT
1st index
2nd index
EM DIE Contains Diesel Emission Factors
IDIEYR (20)
IDIEVH (20)
EFDP (20)
NDITDT
- year
- vehicle class
- emission factor
- no. of entries
Block used by EMFIN, DIESEL
EMORG Contains organic emission factors
EMPB
8.
FLAGS1
IORGYR (20)
IORGFT (20)
IORGCT (20)
IORGVH (20)
EFO (20)
NORTDT
year
fuel type
catalyst type
vehicle class
emission factor
no. of entries
Block used by EMFIN, GETEFO
Contains lead emission factors
EFLEAD (20) - lead emission factor
IPBYR (20)
IPBVH (20)
IPBFT (20)
NPBTDT
year
- vehicle class
- fuel type
- no. of entries
Block used by EMFIN, GETPB
EMSUL Contains sulfate emission factors
ISULYR (20)
ISULFT (20)
ISULCT (20)
ISULVH (20)
EFS (20,2)
NSUTDT
year
fuel type
catalyst type
vehicle class
emission factor, speed
no. of entries
Block used by EMFIN, GETEFS
Contains control flags (see control section)
Block used by CONSEC, EMFIN, ONESEC, PARSEC, GETLCG,
GETVMT, NLCONS, OUT
5-6
-------
FLEET Contains fleet sales fractions and VMT mix
FSF (26,6,3) - fleet sales fractions
10. GAS
1st index
2nd index
3rd index
IFSFYR (26,6)
1st index
2nd index
IFSF (6)
VMT (6)
index
vehicle class
fuel type
model year
index
model year
highest index by vehicle class
VMT fraction by vehicle class
Block used by TWFCAL, FLTSLE, GETVMT, CALC, FRFLT, OUT,
EMFIN
Contains lead content of gasoline, city/highway, fuel economy and
vehicle speed correction
LCGSIZ - no. of lead content entries
LCGYR (20) - year for lead content
PBCG (20,2) - lead content of gasoline
1st index = index
2nd index = fuel type
l=leaded
2=unleaded
NCFEYR (27)
CFE (27,6)
Year for fuel economy
miles per gallon by vehicle class
1st index = index
2nd index = vehicle class
1=LDV 4=HDGV1
2=LDT1 5=HDGV2
3=LDT3 HDGV
ICFESZ - no. of entries for city/highway fuel economy
IFECSZ - no. of entries for speed corrections
FECUS (13) - vehicle speed
FEC (13,2) - speed correction
1st index = index
2nd index = driving conditions
l=cyclic
2=steady
Block used by EMFIN, GETUCG, GTPBCG, GETCS, GETCFE
11. HDDT Contains heavy duty diesel truck mileage accumulation
IHTRG - no. of entries
IHDDTMA (20,4) - HDDT mileage accumulation rate
1st index = age
2nd index = wt class
IHDDYR (21)
HDDTRG (21,
1=2B 3=med HDDT
2=Lt. HDDT 4=heavy HDDT
- calendar year
registration fraction
1st index = index
2nd index = wt. class
5-7
-------
12. HDDV
Contains HDDT low altitude gm/bhp-hr to gm/mile conversion
factors
IHDDMY (13) -
HDDVCF (13) -
model year
conversion factor
Block used by EMFIN, DIESEL
13. MSFL Contains misfueling rates
RMISFL (20,4,2) - misfueling
1st index
2nd index
3rd index
TMISFL(6,2) =
1st index
2nd index
= age
= vehicle class
= IMFLAG
simple misfueling rates
= vehicle class
= IMFLAG
Block used by EMFIN, ONESEC and NLCONS
PROJCM Contains project title
PROJID - title
Block used by CONSEC, OUT, OUTHDR
15. PSIZE Contains particle size cumulative fraction for chosen particle size
cutoff
PSZCUT = particle size cutoff
PSL = cumulative fraction for leaded vehicles
PSNL = cumulative fraction for non-leaded vehicles w/o
catalyst
PSNLCT = cumulative fraction for non-leaded vehicles with
catalyst
PSDIE = cumulative fraction for diesel vehicles
PSBRK = cumulative fraction for brake wear
PSTIRE = cumulative fraction for tire wear
Block used by ONESEC, PARSEC, PARSIZ, LEAD, ORGANC,
SULFAT, DIESEL, CALC, OUT
5-8
-------
16. PSDIST Contains particle size distribution
PSDIST (5,2,6) - particle size distribution
1st index = index
2nd index = l=size
2=cumulative fraction
3rd index = source
l=leaded 4=diesel
2=unleaded with catalyst 5=brake
3=unleaded w/o catalyst 6=tire
Block used by EMFIN, PARSIZ
17. SCENR Contains scenario variables
IYEAR - scenario year
VEHSPD - vehicle speed
IDRIV - driving conditions
l=cyclic or transient
2=steady cruise
Block used by EMFIN, PAR5EC, CALC, OUT
5-9
-------
TABLE 5-3
GLOBAL VARIABLES
Purpose
Age of vehicle in years
Fuel Type l=leaded gas
2=unleaded gas
3=diesel
IVHCLS 1-6 Vehicle class
1=LDV, LDDV
2=LDT1, LDDT1
3=LDT2, LDDT2
4=HDGV, HDDT
5=HDGV2
6=MC
Note that this scheme is slightly different than the
vehicle classes the user uses to input data and
interpret results.
IYRPTR 1980+ Calendar Year
5-10
-------
5.5 OTHER PROGRAM NOTES
The vehicle classes inside the program are defined differently in comparison
to the vehicle classes as defined in the original report. In the report the
classes were defined as: 1) LDV and LDDV, 2) LDT1 and LDDT1, 3) LDT2
and LDDT2, 4) HDGV, 5) HDDT, and 6) MC. As far as the user is concerned
the program treats vehicle classes in the exact same way as the report.
Inside the program, the vehicle classes are defined in terms of vehicle class
and fuel type. The first three classes are unchanged and motorcycles are
still the sixth class. However, HDDT are now part of the fourth class when
the emissions factors are called and calculations done. HDGV2 are considered
as the fifth class since they have different fuel economy than HDGV. This
split in the definition at vehicle classes lasts until the subroutine CALC is
called. CALC changes the definitions back to those found in the orignal
report.
5-11
-------
6. EXAMPLES
The following examples were run on the Ann Arbor MTS system to
demonstrate both interactive and batch use of the program. The compiled
source code was in the file "-load" and the associated data file was
"partic.ef.da". The examples are:
1. No tampering (catalyst removal) scenario for 1985 and 1990 with
transient driving conditions, vehicle speed 25 mph and including
particles up to 10 microns in size. Program run interactively with
prompts for user input.
2. No misfueling scenario for 1990 and 1995 with transient driving
conditions, vehicle speed 25 mph and including particle size cutoffs of
10 microns for both years and 5 microns for 1990. Program run
interactively with prompts for user input.
3. Alternate lead content of gasoline for!985 to 1990. A 1990 example
with vehicle speed 25 mph transient driving conditions and a 10 micron
size cutoff was run. Program was run interactively.
4. A user specified VMT mix is used for each scenario. Three cases are
run for 1985 with vehicle speed 30 mph, transient driving conditions
and a particle size cutoff of 10 microns for the first two cases and 5
microns for the last case. The program was run in batch mode with no
prompts printed. The user specified data was a a file "example 1"
assigned to unit 5. The normal exit is flagged by the message "END
OF FILE IN IOUGEN DATA FILE".
5. Alternate tampering rates, LDV registration and mileage
accumulation with age and alternate misfueling rates scenario. The
user data is in file "example 2". The tampering rates for LDV, LDT1
and LDT2 are set to .100, .100 and 0.150 respectively in line 4.
Gasoline powered LDV registration rates are altered in lines 5-8 and
gasoline LDV mileage accumulation rates are altered in lines 9-12.
The zero mileage misfueling rates and slopes with mileage for LDV,
LDT1, LDT2 and HDGV1 are given in lines 13-16. The program was run
for 1987 with vehicle speed 25 mph, cruising driving conditions and
particle size cutoff 5 microns.
6-1
-------
Srun -load 4 = par-tir. .ef .da
4000 . EXAMPLE 1
Enter project id:
test exsmole -- no tamper ing
ENTER TAMFLG, VMFLAG. MYMRFG, MISFLG, SPCFLG, IMFLAG
£ .1, 1 I 1 1
i < i
•ENTER FRACTION OF CATALYST EQUIPPED VEHICLES WITH
CATAI YST RC-MOVPTi FOR I.DV. IJVI i ANN LDT2 FORMAT=3< F4 . 3 , 1 X )
NO I/M PROGRAM IN JN FFFF.OT
.000 .000 .000
ENTER YEAR FOR THIS RUN -- (FORMAT = 14)
1985
ENTER VEHICLE SPEED -- (FORMAT = F4.1)
25.00
ENTER DRIVING CONDITIONS (1= TRANSIENT, 2 = CRUISE)
1
ENTER PARTICLE SIZE CUTOFF (0.1-10.) -FORMAT(F5.2>
10.00
TEST EXAMPLE -- NO TAMPERING
CALCULATIONS FOR YEAR 1985 SPEED 25.0 MPH PARTICLE CUTOFF =10.00 MICRONS
FLAG SETT INGS:
TAMPERING RATES = USER
MII.F.AGF. ?, RFC 1ST. - NORM
DRIVING COhmiTJONS -•• 'IRAN
MISFUELING RATES
I/M PROGRAM
= NORM
= NO
1 PARTICULATE EMISSIONS (GRAMS/MILE) BY
ho
VFH. CLASS
I.HV
LI.1T1
1.1X12
I-WGV
HPDT
MO
TOTAL BRfiKf
U:AV.c ORGANIC
0.0 fan 0.0234
0.024? 0.0303
0.044? 0.0286
0.0987 0.1422
0 . 0 0 . 0
0.0341 0.0
." PART HUH. AT F. 8 - 0
SUI FATE
0.009)
0 . 0066
0 . 0059
0 . 002?
0 . 0
0 . 0
.Of 25
DIESEL
0 . 0 i 28
0.0195
0.0)95
0 . 0
1 .9198
0 . 0
VEHICLE
TOTAL
0.0583
0.0808
0 . 0983
0.2432
1 .9198
0.0311
VMT MIXTURE = NORM
SPECIAL CONDITIONS = NORM
VMT FRAC. VMT WT. TOTAL
0.7845
0.0827
0.0470
0.0419
0.0354
0.0085
0.0457
0.0067
0.0046
0.0102
0.0680
0.0003
TOTAL TIRE PARTICIPATES = 0.0020
TOTAL VMT WEIGHTED PARTICULATE EMISSIONS = 0.1500
PARTICULATE FRACTION BY VEHICLE CLASS
VFH. CLASS
1. PV
Lirn
LPT2
HI1GV
HDDT
MO
I.FAO
0 . 2224
0 , ?007
O.'Vf.OJ
0.405'''
0.0
i . 0000
ORGANIC
0.4019
0 . 3753
0 . £9 1 2
0 . 5846
0.0
0.0
SI. It. FATE
0.1557
0.0822
0.0601
0.0095
0.0
0.0
DIESEL
0.2200
0.2417
0. 1985
0 . 0
1 . 0000
0.0
ENTER YEAR FOR THIS RUN -- (FORMAT = 14)
1990
•ENTER VEHICLE SPEED -- (FORMAT = F4.1)
; '• 3Pt flQ
ENTER DRIVING CONDITIONS (1= TRANSIENT, 2= CRUISE)
PARTICLE SIZE CUTOFF (0.1-10.) -FORMAT(F5. 2)
10.00
-------
EXAMPLE 1 (cont)
TEST EXAMPLE -- HO TAMPERING
C A I cut A 'F ? 0 M S F 0 R YEAR 1 9 9 f i S PEED £ 5 . 0 M P H PARTICLE C I J T 0 F F - i 0 . 0 0 M I r: R o N 5
FLAG Sl-TTJMuS:
TAMPERING RATES = USER MISF1JELING RATES = NORM '-'MT MIXTURE = NORM
HII.FAt"!" if, RFCJST. - N'On.'H J /M PF
-------
*run -load 4=partic.of.da
# I n it a 1 i d MT9 c r
0.
0.
0.
1 .
0.
0174
0434
0434
0
73 i 9
0
=
VE
USER
NO
HICLE
TOTAL
0
0
0
0
1
0
.0456
.0741
.0833
. 1674
.7319
.0294
VMT
MIXTURE
= NOR
SPECIAL CONDITIONS = NOR
CLASS
VMT
0
0
0
0
0
0
FRAC . VMT
.7845 0
.0827 0
.0470 0
. 04 1 9 0
.0354 0
.0085 0
WT. TOTAL
.0358
.0061
.0039
. 0070
.0613
.0003
tf PART ( CUt. A't 1: S -- 0 . 0 i 25
PARTICULATES = 0.
WEIGHTED
F'
1. FAD
0.0525
0.0624
0.1538
0 . 3530
0.0
1 . 0000
0020
PARTICULATE EMISSIONS
ARTICULATE
ORGANIC
0.3481
0.2047
0.1820
0.5745
0.0
0 . 0
FRACTION
SULFATE
0.2188
0.1470
0.1428
0.0726
0.0
0.0
BY
= 0.
1 290
VEHICLE
CLASS
DIESEL
0.
0.
0.
0.
1.
0.
3806
5859
52 1 3
0
0000
0
-------
ENTER YEAR FOR THIS RUN -- (FORMAT = 14)
• 1995
ENTER VEHICLE SPEED -- (FORMAT = F4 . 1 )
"""25 '.Off
• ENTER DRIVING CONDITIONS <1 = TRANSIENT, 2= CRUISE)
1
ENTER PARTICLE SIZE CUTOFF (0.1-10.) -FORMAT (F5. 2)
• 10.000
TEST EXAMPLE -- NO MISFUELING
• CALCULATIONS FOR YEAR 1995 SPEED 25.0 MPH PARTICLE
FLAG SETTINGS:
0 TAMPERING RATES = NORM MISFUELING RATES = USER
MIt.FAGF. fc. RFCTST. - NORM J/M PROGRAM = Nfi
DRIVING CONDITION?; = (RAN
PARTICULATE EMISSIONS (GRAMS/MILE) BY VEHICLE
£ VFH. CLASS l.FAD ORGANIC. SUI. FATE DIESEL TOTAL
I.HV 0.0007 0.0 $49 0.0088 0.02)0 0.0454
U.ITJ 0.000? 0.0)24 0.0104 0.0670 0.0906
O LDT2 0.002J 0.0 $24 0.0)26 0.0670 0.0951
HPGV 0.0349 0.0680 0.0180 0.0 0.1209
HIiDT 0.0 0.0 0.0 1.6606 1.6606
9 MC 0.0294 0.0 0.0 0.0 0.0294
TOTAL RRAKF PAFnm.ll.ATES = 0.0$ 25
TOTAL TIRE PARTICULATES = 0.0020
A V TOTA^ VMT WEIGHTED PARTICULATE EMISSIONS = 0.1262
un
PARTICULATE FRACTION BY VEHICLE CLASS
9 VEH. CLASS t.EAU ORGANIC SHI.. FATE DIESEL
LDV 0.0160 0.3278 0.1948 0.4614
l.DT i 0 . 0085 0 . J 370 0.1152 0 . 7394
% LDT2 0.0326 0.1303 0.1328 0.7043
HHGV 0.2887 0.5622 0.1491 0.0
HDDT 0.0 0.0 0.0 1.0000
• MC i . 0000 0.0 0.0 0.0
EXAMPLE 2 (cont)
CUTOFF =10.00 MICRONS
VMT MIXTURE = NORM
SPECIAL CONDITIONS = NORM
CLASS
VMT FRAC. VMT WT. TOTAL
0.7845 0.0356
0.0827 0.0075
0.0470 0.0045
0.0419 0.0051
0.0354 0.0588
'0.0085 0.0003
-------
ENTER YEAR FOR THIS RUN -- (FORMAT = 14)
1990
ENTER VEHICLE SPEED -- (FORMAT = F4.1)
25,000
ENTER DRIVING CONDITIONS (1= TRANSIENT, 2 = CRUISE)
1
ENTER PARTICLE SIZE CUTOFF (0.1-10.) -FORMAT
VFH. CLASS
i.nv
i. n 1 1
I.DT2
HDGV
HDDT
MC
PARTICIPATE FRACTION BY VEHICLE CLASS
I.F.AVJ
0.0470
0.0547
0.1349
0.3444
0.0
1 . 0000
ORGANIC
0.3432
0 . 2029
0.1313
0.5729
0.0
0 . 0
SULFATE
0.2209
0.1490
0.1471
0.0327
0.0
0.0
DIESEL
0 . 3339
0 . 5934
0 . 5367
0.0
1 .0000
0.0
0.7345
O.V327
0.0470
0.0419
0.0354
0.0035
0.0336
0.0057
0.0036
0.0057
0.0530
ENTER YEAR FOR THIS RUN -- (FORMAT = 14)
!
+Attent ion interrupt at 35335E
if
c
c
c
c
c
c
-------
test example —alternate lead content for 1985 to 1990
ENTER JAMFLC, VMFLAG, -~- ~ "
'TTTTT f ~~
EXAMPLE 3
ALTERNATE LEAD CONTENT OF GASOLINE BY YEAR
FNTF.R FIRST AND LAST YEAR TO CHANGE - (14, IX,:
1985 1990
ENTER LEAD CONTENT OF GASOLINE FOR LEADED
FOR YEAR 1985
0.010
ij'YEAR 1986
1.05O. 0.010
FOR YEAR 1987
1.050 0.005
FOR YEAR 1988
1.050 0.005
FOR YEAR 1989
1.050 0.005
FOR YEAR 1990
1.050 0.005
AND UNLEADED -- FORMAT=F5.3,IX,F5.3
ENTER YEAR FOR THIS RUN -- (FORMAT = 14)
1990
ENTER VEHICLE SPEED -- (FORMAT = F4.1)
25.00
ENTER DRIVING CONDITIONS (1= TRANSIENT, 2= CRUISE)
1
ENTER PARTICLE SIZE CUTOFF (0.1-10.) -FORMAT(F5.2>
10.00
i
-vl
TEST EXAMPLE -- ALTERNATE LEAD CONTENT FOR 1985 TO 1990
CALCULATIONS FOR YEAR 1990 SPEED £5.0 MPH PARTICLE CUTOFF =10.00 MICRONS
FLAG SF.TTINCS:
TAMPERING RATES = NORM
MILEAGE H, REGIST. = NOKM
DRIVING CONDITION?: = TRAN
MISFUELING RATES = NORM
I /M PROGRAM = NO
VMT MIXTURE = NORM
SPECIAL CONDITIONS = USER
PARTICIPATE EMISSIONS (GRAMS/MILE) BY VEHICLE CLASS
VFH.CLASS LEAP
ORGANIC SUl.FATE DIESEL TOTAL
VMT FRAC. VMT WT. TOTAL
LDV
LDT1
LDT2
HHGV
HDDT
MC
0 . 005?
0.0139
0 . 0237
0 . 0559
0.0
0 . 0294
0.0195
0.0234
0 . 0228
0 . 0988
0.0
0.0
0 . 0085
0.0071
0.0078
0.0)00
0.0
0 . 0
0.0174
0.0434
0.0434
0 . 0
1.7319
0 . 0
0.0506
0.0879
0.0977
0.1647
1.7319
0.0294
0.7845
0.0827 .
0.0470
0.0419
0.0354
0.0085
0.0397
0.0073
0.0046
0.0069
0.0613
0 . 000.3
TOTAL BRftKF PARTICIPATES = 0.0125
TOTAL TIRE PARTICIPATES = 0.0020
TOTAL VMT WEIGHTED PARTICIPATE EMISSIONS = 0.1346
VEH.CI.ftS8 I.FAti
PAfaiCULATE FRACTION BY VEHICLE CLASS
ORGANIC SULFATE DIESEL
LDV
LDT1
LDT2
HDCV
HDDT
'-. MC
0.1027
0.1584
0.2426
0.3397
0.0
1.0000
0 . 3856
0 . 2666
0 . 2333
0.5999
0.0
0.0
0.1688
0.0809
0.0795
0.0604
0.0
0.0
0.3429
0.4941
0.4445
0.0
1 . 0000
0.0
-------
p *f *1
1
£
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
t e d operand at end of command
TEST EXAMPLE -- VMT MIK FOR EACH SCENARIO
) ?. ) 1 I 1
1935
30 . 0
1
10. 00
0.750 0.100 0.100 0.050 0.000 0.000
1 985
30.0
1
10.00
0.600 0.100 0.100 0.050 0.150 0.000
1 985
30.0
1
5.00
0.600 0.100 0.100 0.050 0.150 0.000
EXAMPLE 4
:*run -load 4=parti.c.ef.da 5=examplel
TEST EXAMPLE -- VMT MIX FOR EACH SCENARIO
CALCULATIONS FOR YEAR i 985 SPEED 30.0 MPH
PARTICLE CUTOFF =10.00 MICRONS
FLAG SETTINGS:
TAMPERING RATES = NORM
MILEAGE S REGIST. ~ NORM
DRIVING CONDITIONS -- TRAN
MISFUELING RATES = NORM
I/M PROGRAM = NO
VMT MIXTURE = SCEN
SPECIAL CONDITIONS = NORM
PARTICULATE EMISSIONS (GRAMS/MILE) BY
00 VEH. CLASS LEAH ORGANIC SHI. FATE DIESEL TOTAL
I.DV 0.0120 0.0234 0.0096 0.0)28 0.0578
LPT) 0.0235 0.0303 0.0069 0.0)95 0.0803
LOT? 0.04)8 0.0286 0.0062 0.0)95 0.0962
HDCV 0.0987 0.1422 0.0023 0.0 0.2432
Hilbt" " " 0.0 0.0 0.0 1.9)9F! 0.'0
MC 0.03)1 0.0 0.0 0.0 0.0
TOTAL BRAKt: PARTICUI.ATFS = 0.0)25
TOTAL TIRE PARTICULATES = 0.0020
TOTAL VMT WEIGHTED PARTICULATE EMISSIONS = 0.0877
PARUCULrtTE FRACTION BY VEHICLE CLASS
VEH. CLASS LEAD ORGANIC SULFATE DIESEL
LDV 0.2077 0.4051 0.1655 0.2217 W
LDT1 0.2932 0.3776 0.0859 0.2432 A
LDT2 0.4350 0.2978 0.0643 0.2030
HHGV 0.4059 0.5846 0.0095 0.0
HDD7 0.0 0.0 0.0 0.0
MC 0.0 0.0 0.0 0.0
VMT
0
0
0
0
0
0
FRAC .
. 7500
. 1 000
. 1000
.0500
."6" " " "
.0
VMT
0.
0.
0.
0.
0.
0.
WT. TOTAL
0434
0080
0096
0122
0
0
-------
'u>
TEST EXAMPLE -- VMT MIX FOR EACH SCENARIO
CALCULATIONS FOR YEAR i905 SPEED 30.0 MPH PARTICLE CUTOFF =10.00 MICRONS
FLAG SETTINGS:
TAMPERING RATES = NORM MISFUELING RATES
MItflAGF. fc REGIST. = NORM I/M PROGRAM
DRIVING COM IT IONS - 7 RAN
PARTICULATE EMISSIONS (GRAMS/MILE) BY
VEH. CLASS LEAD
LPV
1. DTI
LPT2
HDGV
HDDT
MC
0 . 0 1 20
0 . 0235
0 .04)8
0.0937
0 . 0
0.0311
ORGANIC
0 . 0234
0.0303
0 . 0286
0.1422
0 . 0
0.0
SULFATE
0 . 0096
0.0069
0 . 0062
0.0023
0.0
0.0
DIESEL
0.01 28
0.0195
0 . 0 1 95
0.0
1 .9)98
0.0
= NORM VMT MIXTURE = SCEN
= NO SPECIAL CONDITIONS = NORM
VEHICLE CLASS
TOTAL
0.0573
0.0803
0.0962
0.2432
1 .9198
0.0
VMT FRAC.
0.6000
0.1000
0.1000
0.0500
0 . 1 500
0.0
VMT WT. TOTAL
0.0347
0 . 0080
0.0096
0.0122
0 . 28SO
0.0
TOTAI PRAKF. PART ICLH.A7ES - 0.0)25
TOTAL TIRE PARTICULATES = 0.0020
TOTAI. VMT WEIGHTED PARTICULATE EMISSIONS = 0.3670
VFH. CLASS I.FAP
PARTICULATE FRACTION BY VEHICLE CLASS
ORGANIC SULFATE DIESEL
LDV
I. DTI
LDT2
HDGV
HDDT
MC
0.2077
0 . 2932
0.4350
0 . 4059
0.0
0.0
0.4051
0 . 3776
0.2978
0.5846
0.0
0 . 0
0.1655
0.0859
0.0643
0.0095
0.0
0.0
0.2217
0.2432
0 . 2030
0.0
1 .0000
0.0
TEST EXAMPLE -- VMT MIX FOR EACH SCENARIO
CALCULATIONS FOR YEAR 1985 SPEED 30.0 MPH PARTICLE CUTOFF = 5.00 MICRONS
FLAG SETTINGS:
TAMPERING RATES = NORM MISFUELING RATES
Mil. FACE & REGIST. - NORM I/M PROGRAM
DRIVING CONDITIONS = TRAN
PARTICULATE EMISSIONS (GRAMS/MILE) BY
VEH. CLASS
LDV
LP71
LDT2
HDGV
HDDT
MC
t EAM
0 . 0096
0.0188
0.0333
0 . 0785
0.0
0.0247
ORGANIC
0 . 0204
0 . 0253
0.0238
0..1130
0.0
0.0
SULFATE
0 . 0090
0 . 0065
0 . 0050
0.0018
0.0
0.0
DIESEL
0.0)2)
0.0185
0.0135
0.0
1.8174
0.0
= NORM VMT MIXTURE = SCEi
= NO SPECIAL CONDITIONS = NORI
VEHICLE CLASS
TOTAL
0.0512
0.0691
0.0814
0.1933
1.8174
0.0
VMT FRAC.
0 . 6000
0.1000
0.1000
0 . 0500
0 . 1 500
0 . 0
VMT WT. TOTAL
0.0307
0.0069
0.0031
0.0097
0.2726
0.0
TOTAI. BRAKE PART ICUI. At E S == 0.0106
TOTAL TIRE PARTICULATES = 0.0010
TOTAL VMT WEIGHTED PARTICULATE EMISSIONS = 0.3397
VEH. CLASS LEAP
PARTICULATE FRACTION BY VEHICLE CLASS
ORGANIC SULFATE DIESEL
LDV
LDTJ
LDT2
HDGV
HDDT
MC
0.1S31
0.27)9
0.4096
0.4059
0.0
0.0
0.3933
0 . 3666
0 . 2920
0.5846
0.0
0 . 0
0.1766
0.0940
0.0713
0.0095
0.0
0.0
0.2369
0.2675
0.2271
0.0
1 .0000
0.0
-------
EXAMPLE 5
c
c
Sedit examp
p *f *l
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
Ie2
3
TEST EXAMPLE - ALT. TAMPERING*
214211
.100 .100 .150
1 1
.100 .100 .100 .100 .050 .050
.030 .030 .030 .030 .030 .030
0 0
1 1
15000 14000 13000 12000 11000
9000 8000 7000 6000 5000
0 0
0.050000 0.000100
0.060000 0.000200
0.070000 0.000200
0.070000 0.000200
1987
25.00
2
5.00
LDV REG.&MILEAGE» MISFUELING
.050 .050 .050 .050
.030 .030 .030 .030
10000 10000 10000 10000 10000
4000 4000 4000 4000 4000
-------
EXAMPLE 5 (cont)
TEST EXAMPLE - ALT. TAMPERING, LDV REG.S.MILEAGE, MISFUELING
CALCULATIONS FOR YEAR J 987 SPEED 25.0 MPH PARTICLE CUTOFF = 5.00 MICRONS
FLAG SETTINGS:
TAMPERING RATES = USER
MILEAGE & REG I ST. =• M-»R
DRIVING CONDITIONS = CRSE
MISFUELING RATES = USER
I/M PROGRAM = NO
VMT MIXTURE = NORM
SPECIAL CONDITIONS = NORM
PARTICIPATE EMISSIONS
-------
APPENDIX A
SOURCE CODE LISTING
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CC1d=SHJD Page
1 C THIS IS PROGRAM CALCULATES LEAD, ORGANIC. SULFATE AND
2 C PARTICULATE EMISSIONS. IT IS DESIGNED TO BE SIMILIAR IN
3 C FORMAT TO MOBILES.
4 C
5 C PREPARED FOR E.P.A. BY: ENERGY AND ENVIRONMENTAL ANALYSIS
6 C 1655 NORTH FORT MEYER DRIVE
7 C ARLINGTON VA 22209
8 C
9 C
10 C MAIN PROGRAM
11 COMMON/IOUCOM/IOUREP,IOUERR,IOUASK,IOUDAT,IOUGEN
12 C
13 C DEFAULT UNIT ASSIGNMENTS
14 C DATA FILE (DEFAULT EMISSION FACTORS)
15 IOUDAT=4
16 C USER SUPPLIED DATA
17 IOUGEN=5
18 C OUTPUT REPORT
19 IOUREP=6
20 C C ERROR REPORT
21 IOUERR=6
22 C PROMPTS
23 IOUASK=6
24 C
25 C READ CONTROL SECTION
26 CALL CONSECU99)
27 C READ DATA TABLES
28 CALL EMFIN
29 INERR=0
30 C READ ONETIME DATA
31 CALL ONESECUNERR.&99)
32 IF(INERR.GT.50)GOTO 20
33 C READ PARAMETER SECTION (SCENARIOS)
34 10 CALL PARSEC(INERR,&99)
35 IF(INERR.GT.50)GOTO 20
36 C CALCULATIONS
37 CALL CALC
38 C OUTPUT
39 CALL OUT
40 GOTO 10
41 99 STOP
42 20 WRITEUOUERR, 1000)
43 1000 FORMAT('0','EXCESS DATA ERRORS PREVENT FURTHER ANALYSIS')
44 END
45 C
46 c ***************************************
47 C * INPUT ROUTINES *
48 c **«******»*******»********»*«»*»«*»»»*»
49 C
50 C
51 SUBROUTINE CONSEC(')
52 C
53 C CONSEC reads In and validates the Control Section flags. It also sets the
54 C input mode for the run, resets the output units (if requested) and reads in
55 C and stores the run title.
56 C
57 C Cat led by MAIN.
58 C
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
8 1
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
1 10
1 1 1
12
13
14
15
16
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
Calls MOD (FORTRAN library function), OUTHDR and QUITER.
Input on ca 1 1 :
common blocks:
Output on return:
common blocks:
/IOUCOM/ IOUREP, IOUERR, IOUASK
/PROJCM/ PROJID
Local array subscripts:
IOUCHGO) - IOUCHG ( JU )
IOUOLD(3) - IOUOLD ( JOU )
IOUREF(9) - IOUREF ( JU )
NAMFLG(13) - NAMFLG ( I FG )
Local variable / array dictionary:
Name Type Description
FLGERR I number of flag values found to be out of range
INFORM I dummy variable: QUITER call purpose is to write out run
' inform' at i on , not an error or a warning message
IOUCHG I user entered changes to output unit assignments
IOUOLD I output unit numbers to be used this run
IOUREF I character representations of i/o unit numbers 1 - 9
IOUSUM I sum of acceptable changed i/o unit numbers
IPRSAV I backup of entered IPROMT - used in QUITER if out of range
JUNIT I potential numeric match to entered character for unit change
MDEFLG I decides prompt & read format for Control Section flags:
1 = vertical format 2 = horizontal format
NAMFLG R*8 flag names for prompts
INTEGER TAMFLG, SPCFLG, VMFLAG '-
COMMON/ FLAGS 1 /I PROMT, TAMFLG, SPCFLG, VMFLAG, MYMRFG, IMFLAG , MI SFLG
COMMON/ IOUCOM/ IOUREP, IOUERR. IOUASK , IOUDAT , IOUGEN
COMMON /PROJCM/ PROJID(20)
INTEGER FLGERR
DIMENSION IOUREF (9) ,IOUOLD(3) , IOUCHGO)
REAL*8 NAMFLG(13) ,NAMFNC(2)
EQUIVALENCE
* (lOUOLD(l) , IOUREP) , (IOUOLD(2) , IOUERR) , ( IOUOLD( 3) , IOUASK)
DATA FLGERR/0/
DATA IOUREF/
* '1 ','2 ','3 '.'4 ','5 ','6 ','7 ','8 ','9 '/
DATA NAMFLG/' TAMFLG:',' VMFLAG:',' MYMRFG:',
& ' MISFLG:',' SPCFLG:',' IMFLAG:'/
&.NAMFNC/' FREMCT ',' FRFLT '/
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD Page
117 C Read In prompt/mode flag and output unit changes (if any).
1 18 C
119 READdOUGEN. 100,END=99) I PROMT , IOUCHG
120 100 FORMAT(I1.3A1)
121 C
122 C Output unit numbers parsed here.
123 C
124 C Change alphanumeric to numeric values.
125 C
126 DO 20 JOU=1 ,3
127 C
128 C Test for allowed values (no input units).
129 C
130 DO 10 JU=1 ,9
131 IF(JU.EQ.IOUGEN)GOTO 10
132 IF(JU.EQ.IOUDAT)GOTO 10
133 JUNIT=JU
134 IF(IOUCHG(JOU).EQ.IOUREF(JU)) GOTO 15
135 10 CONTINUE
136 C
137 C Set default for unrecognized value.
138 C
139 JUNIT=0
140 15 IOUCHG(JOU)=JUNIT
141 20 CONTINUE
142 C
143 C Assign output units.
144 C
145 IOUSUM=0
146 DO 25 JOU=1,3
147 IOUSUM=IOUCHG(JOU)+IOUSUM
148 IF(IOUCHG(JOU]I.GT.O) IOUOLD(JOU) = IOUCHG(JOU)
149 25 CONTINUE
150 C
151 C Parse prompt/mode flag here.
152 C
153 C Save input flag for future possible abort.
154 C
155 IPRSAV=IPROMT
156 IF(IPROMT.GE.1.AND.IPROMT.LE.4) GOTO 30
157 C
158 C Abort will happen due to prompt/mode flag error.
159 C Collect title for id of possible multiple units.
160 C
161 IPROMT=1
162 WRITE(IOUASK,210)
163 READ(IOUGEN,110,END=99) PROJID
164 GOTO 60
165 C
166 C Parse prompt/mode flag.
167 C
168 30 MDEFLG=1
169 IF(IPROMT.GE.3)MDEFLG=2
170 IPROMT=MOD(IPROMT-1,2)+1
171 C
172 C Read title.
173 C
174 55 IF(IPROMT.EQ.2) WRITE(IOUASK,210)
-------
Listing of PART.LEAD at 13:40:01 on DEC 2. 1985 for CCid=SHJD
Page
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
21 1
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
210 FORMAT('0'.'Enter project id:')
READUOUGEN, 1 10, END=99) PROJID
110 FORMAT(20A4)
If error unit no longer same as report unit, write run title to it.
60 IF(IOUREP.NE.IOUERR) CALL OUTHDR(IOUERR)
Display unit number values if changed.
IF(IOUSUM.GT.O) CALL QUITERCO.,0,71.INFORM)
IF(IPROMT.EQ.2.AND.IOUASK.NE.IOUREP.AND.IOUASK.NE.IOUERR)
* CALL OUTHDR(IOUASK)
Reset INITPR so that subsequent branch check in OUTPUT will succeed.
INITPR=1
Test IPRSAV (= entered IPROMT value). Out-of-range => run stops.
IF(IPRSAV.LT.1.OR.IPRSAV.GT.4) CALL QUITER(0.,IPRSAV,2,FLGERR)
IF(MDEFLG.EQ.2) GOTO 65
Read flags via vertical input mode.
IF(IPROMT.EQ.2) WRITE(IOUASK,220)NAMFLG(1)
READ(IOUGEN,120,END=99) TAMFLG
IF(TAMFLG.LT.1.OR.TAMFLG.GT.2) CALL QUITER(0.,TAMFLG,3,FLGERR)
IF(IPROMT.EQ.2)WRITE(IOUASK,220)NAMFLG(2)
READUOUGEN, 1 20 , END=99 ) VMFLAG
IF(VMFLAG.LT.1.OR.VMFLAG.GT.3) CALL QUITERCO..VMFLAG,5,FLGERR)
IF(IPROMT.EQ.2) WRITE(IOUASK,220)NAMFLGC3)
READ(IOUGEN,120,END=99) MYMRFG
IF(MYMRFG.LT.1.OR.MYMRFG.GT.4) CALL QUITER(0..MYMRFG,6,FLGERR)
IF(IPROMT.EQ.2) WRITE(IOUASK,220)NAMFLG(4)
READdOUGEN, 1 20 , END=99)MI SFLG
IF(MISFLG.LT.1.OR.MISFLG.GT.3) CALL QUITER(0.,MISFLG,4.FLGERR)
IF(IPROMT.EQ.2) WRITE(IOUASK,220)NAMFLG(5)
READdOUGEN, 120,END=99) SPCFLG
IF(SPCFLG.LT.1.OR.SPCFLG.GT.3JCALL QUITER(0.,SPCFLG.7,FLGERR)
IF(IPROMT.EQ.2) WRITE(IOUASK,220)NAMFLG(6)
READdOUGEN, 120,END=99) IMFLAG
IF(IMFLAG.LT.1.OR.IMFLAG.GT.2) CALL QUITER(0.,IMFLAG,8,FLGERR)
120 FORMAT(II)
220 FORMATC ',' Enter ' ,A8)
GOTO 70
Read flags in via horizontal mode.
65 IF(IPROMT.EQ.2) WRITE(IOUASK,230)
230 FORMATC '.'ENTER TAMFLG, VMFLAG, MYMRFG, MISFLG, SPCFLG, IMFLAG')
READdOUGEN,130,END=99) TAMFLG.VMFLAG,MYMRFG,MISFLG.SPCFLG,IMFLAG
130 FORMAT(I 1,12(IX,I 1 ))
IF(TAMFLG.LT.1.OR.TAMFLG.GT.2) CALL QUITER(0..TAMFLG,3.FLGERR)
IF(MISFLG.LT.1.OR.MISFLG.GT.3) CALL QUITER(0.,MISFLG.4,FLGERR)
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page
233
234
235
236
237
23S
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
O"7K
Z /D
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
IF(VMFLAG.LT. 1 . OR . VMFLAG .GT . 3) CALL QUITER(0 . , VMFLAG. 5 , FLGERR )
IF(MYMRFG.LT. 1 . OR . MYMRFG . GT . 4 ) CALL QUITER (0 ., MYMRFG , 6 , FLGERR )
IF(SPCFLG.LT. 1 .OR.SPCFLG.GT.2) CALL QUITER(0 ., SPCFLG, 7 . FLGERR)
IF(IMFLAG.LT. 1 . OR . IMFLAG . GT . 2) CALL QUITER(0 ., IMFLAG, 8 , FLGERR )
1 or more flag values out of range => warn user & then stop run.
70 IF(FLGERR.GT.O) CALL QUITER(0 ., 0 , 53 , FLGERR)
RETURN
EOF on any attempted read => take alternate return 1 => run aborts.
99 RETURN!
END
SUBROUTINE QUITER ( RVALUE , I VALUE , MESSAG , I NERR )
QUITER prints error and warning messages on unit IOUERR.
CALLED BY CONSEC . PARSEC, ONSEC, GETVMT. GETLCG, FREMCT, FRFLT
Input on cal 1 :
parameter list: RVALUE , IVALUE .MESSAG, INERR
common blocks:
/IOUCOM/ IOUREP, IOUERR, IOUASK
Output on return:
parameter list: INERR
Warning: there is a conditional branch to STOP at statement number
Local array subscripts:
IOBMSG(5) - IOBMSG ( ICH )
NAMFLG(15) - NAMFLG ( ICH )
NAMFNCU6) - NAMFNC ( ICH )
NAMVEH(8) - NAMVEH ( IV )
Local variable / array dictionary:
Name Type Description
IOBMSG I character string vector: 'out of bounds for'
MAXMES I maximum message code value ( = 83 in MOBILE3 release
MCODE I = MESSAG = message code, unless < 0 or > MAXMES => =
NAMFLG R*8 control flag names + 'Unknown'
98.
vers i on)
1
NAMFNC R*8 function names (to identify subprogram source of error)
NAMVEH R*8 vehicle type names
COMMON/ IOUCOM/ IOUREP, I OUERR , I OUASK , IOUDAT , I OUGEN
DIMENSION IOBMSG(5)
REAL*8 NAMFLG(S) ,NAMVEH(6) ,NAMFNC(2)
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page
291
292
293
294
295
296
297
29B
299
300
301
302
303
304
305
306
307
308
309
310
31 1
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
DATA NAMFLG/
' Unknown ' ,
' I PROMT ' ,
' TAMFLG ' ,
' MISFLG ',
' VMFLAG ' ,
' MVMRFG ' .
' SPCFLG ' ,
' IMFLAG '/
C
DATA NAMVEH/
'LOV
'LOT 1
' LOT 2
'HDGV
'HDDT
•MC '/
C
DATA NAMFNC/ 'FREMCT ','FRFLT '/
C
DATA IOBMSG/
* ' ou ' , ' t of',' bou'.'nds '.'for '/
C
DATA MAXMES/40/
C
MCODE=MESSAG
C
C Trap unknown errors.
C
IF(MCODE.LT. 1 . OR .MCODE . GT .MAXMES) MCODE=1
C
C Test for warnings.
C
I F (MCODE . NE . 9 . AND . MCODE . NE . 1 4 . AND . MCODE . NE . 1 5 . AND . MCODE
* . AND. MCODE. NE. 20. AND. MCODE. NE. 21 .AND. MCODE
* . AND . MCODE . NE . 26 . AND . MCODE . NE . 27 . AND . MCODE
* . AND . MCODE . NE . 29 . AND . MCODE . NE . 33 . AND . MCODE
* )WRITE(IOUERR.300)MCODE
300 FORMAT( '-' , '**ERROR ',12,':')
C
C Branch to the appropriate error / warning message.
C
GOTO(
10,11,11,11.11,11,11,11,16,17,
18, 19,20,21 ,22,29,30,31 ,39,49,
52,53,60,61 ,71 ,72,73,74,75,76,
77,78,79.80,81 ,82,83,84.85.86
) .MCODE
C
10 WRITEUOUERR.310) MESSAG .NAMFLG(MCODE)
310 FORMAT( '+', 10X. 'Message code '.14. AS)
GOTO 98
11 WRITE(IOUERR,31 1 ) I VALUE , IOBMSG . NAMFLG(MCODE)
31 1 FORMAT ( '+' , 10X.I4.5A4, 'flag' ,A8)
IF(MCODE.EQ.2) WRITE( IOUERR , 3 1 2)
312 FORMAT(' ',' Prompt /mode type (1 to 4) must be corrected
* ' your next run . ' )
GOTO 98
i
NE. 18
NE.25
NE. 28
NE.34
before
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD Page
349 16 WRITE(IOUERR,316) RVALUE
350 316 FORMATC ','Warning: ',610.3,' speed reduced to 55 m.p.h.',
351 * ' max imum')
352 GOTO 99
353 17 WRITE(IOUERR,317) RVALUE,IOBMSG
354 317 FORMATC'+',10X.G10.3,5A4,'VMT MIX (0. to 1.)')
355 GOTO 98
356 18 WRITECIOUERR.318) RVALUE
357 318 FORMATC'+',10X.G10.3,' sum of VMT MIX is not equal to 1.')
358 GOTO 98
359 19 WRITECIOUERR.319} IVALUE
360 319 FORMATC'+',10X.G10.3,' negative model year mileage')
361 GOTO 98
362 20 WRITECIOUERR.320) RVALUE
363 320 FORMAT(' + ',10X.G10.3,' negative model year registration')
364 GOTO 98
365 21 WRITECIOUERR.321) RVALUE
366 321 FORMATC' '.'Warning: ',610.3,' registration with zero mileage')
367 GOTO 99
368 22 WRITECIOUERR.322) RVALUE
369 322 FORMATC' '.'Warning: ',610.3,' mileage with zero registration')
370 GOTO 99
371 29 WRITECIOUERR.329)
372 329 FORMATC' + ',10X,' RANGE OF YEARS FOR LEAD CONTENT CANNOT BE <= 0')
373 GOTO 98
374 30 WRITECIOUERR.330) RVALUE
375 330 FORMATC' + ',10X.G10.3, ' intercept must be positive')
376 GOTO 98
377 31 WRITECIOUERR.331) RVALUE
378 331 FORMATC' '.'Warning: '.G10.3,' negative slope for ageing vehicle')
379 GOTO 99
380 39 WRITECIOUERR.339) RVALUE
381 339 FORMATC'+'.10X.G10.3,' speed must be positive')
382 GOTO 98
383 49 WRITECIOUERR.349) RVALUE
384 349 FORMATC' '.'Warning: '.G10.3.' MYR sum not = 1. (will normalize)')
385 GOTO 99
386 52 WRITECIOUERR.352) RVALUE
387 352 FORMATC' '.'Warning: '.G10.3,' speed increased',
388 * ' to 5 m.p.h. minimum')
389 GOTO 99
390 53 WRITECIOUERR.353)
391 353 FORMATC'0','You must correct all flags before you run again.')
392 GOTO 98
393 60 WRITECIOUERR.360)
394 360 FORMATC'+',10X,' END OF FILE IN IOUGEN DATA FILE')
395 GOTO 98
396 61 WRITECIOUERR.360) INERR.NAMFNCCI VALUE)
397 361 FORMATC'+',10X,' DEFAULT USED FOR',15.' IN INDEX FUNCTION' ,A8)
398 GOTO 98
399 71 WRITECIOUERR.371) IOUREP,IOUERR,IOUASK
400 371 FORMATC' '.'COMMENT: CURRENT OUTPUT UNIT NUMBERS ARE ',
401 *' IOUREP=' , I 1 , ' IOUERR=' , I 1 , ' IOUASK=',I1)
402 GOTO 99
403 72 WRITECIOUERR.372) RVALUE
404 372 FORMATC'+',10X.610.3,' PARTICLE SIZE CUTOFF INCREASED TO 0.1 MICRONS')
405 GOTO 98
406 73 WRITE(IOUERR,373) RVALUE
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page
407
408
409
410
41 1
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
373 FORMATC/+',10X.G10.3,' PARTICLE SIZE CUTOFF IS DECREASED TO 10 MICRONS')
GOTO 98
74 WRITE(IOUERR.374)RVALUE,IVALUE
374 FORMATC' + ',10X,' LEAD CONTENT OF LEADED GAS CANNOT BE '.G10.3,
&' FOR YEAR',15)
GOTO 98
75 WRITE(IOUERR,375)RVALUE,IVALUE
375 FORMATC'->-', 10X, ' LEAD CONTENT OF UNLEADED GAS CANNOT BE '.G10.3,
&' FOR YEAR' ,15)
GOT098
76 WRITE(IOUERR,376) RVALUE
376 FORMATC'+',10X.G10.3,' TAMPERING RATE MUST BE POSITIVE OR ZERO —'.G10.3)
GOTO 98
77 WRITEUOUERR.377)
377 FORMATC'+',10X,' END-OF-FILE IN IOUDAT DATA FILE')
GOTO 99
78 WRITE(IOUERR.378)NAMFNC(IVALUE)
378 FORMATC'+','ARGUMENT OUT OF BOUNDS IN FUNCTION ',A8)
GOTO 98
79 WRITE(IOUERR,379) RVALUE
379 FORMATC '.'WARNING: VEHICLE SPEED =',G10.3,' LESS THAN 19.6 MPH FOR SULFATE CAL.')
GOTO 98
80 WRITE(IOUERR,380) RVALUE
380 FORMATC ','WARNING: VEHICLE SPEED =',G10.3,' MORE THAN 34.8 MPH FOR SULFATE CAL.')
GOTO 98
81 WRITEUOUERR, 381 ) IVALUE
381 FORMATC ' +' , 10X, 'DRIVING CONDITION FLAG MUST BE 1 OR 2 CANNOT BE ',12)
GOTO 98
82 WRITE(IOUERR,382) IVALUE
382 FORMAT(' + ', 10X, 'YEAR CANNOT BE ',14,
GOTO 98
83 WRITE(IOUERR,383) IVALUE
383 FORMATC'+',10X.'NO SUCH VEH. CLASS '
GOTO 98
84 WRITE(IOUERR,384) IVALUE
384 FORMATC'+',' NO SUCH FUEL TYPE ',14,
GOTO 98
85 WRITECIOUERR.385) IVALUE
385 FORMATC' + ',10X,'NO VEH. CLASS ',14,'
GOTO 98
86 WRITECIOUERR.386) IVALUE
386 FORMATC' + ',10X,'NO FUEL TYPE ',14,' FOR MILEAGE ACCUMULATION DATA')
GOTO 98
C
98 IF(MCODE.EQ.1.OR.MCODE.EQ.2.OR.MCODE.EQ.10.OR.MCODE.EQ.12
& .OR.MCODE.EQ.16.OR.MCODE.EQ.17.OR.MCODE.EQ.19
& .OR.MCODE.EQ.22.OR.MCODE.EQ.23
& .OR.MCODE.EQ.30.OR.MCODE.EQ.31.OR.MCODE.EQ.35
8. . OR . MCODE . EQ. 36 . OR . MCODE . EQ. 37 . OR . MCODE . EQ . 38
8, .OR.MCODE.EQ.39.OR.MCODE.EQ.40) STOP
INERR=INERR+1
C
99 RETURN
END
C
SUBROUTINE EMFIN
C
C READ DATA TABLES FROM UNIT IOUDAT
' — MUST BE LATER THAN 1979')
,12,' FOR REGISTRATION DATA')
' FOR REGISTRATION DATA')
FOR MILEAGE ACCUMULATION DATA')
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD Page
465 C
466 COMMON/IOUCOM/IOUREP.IOUERR,IOUASK,IOUDAT,IOUGEN
467 COMMON/FLAGS 1/IPROMT,TAMFLG,SPCFLG,VMFLAG.MYMRFG,IMFLAG.MISFLG
468 COMMON/EMORG/IORGYR(20),IORGFT(20),IORGCT(20),IORGVH(20),
469 &EFO(20).NORTOT
470 COMMON/EMSUL/ISULYRC20),ISULFT(20),ISULCT(20),ISULVH(20).
471 &EFS(20,2).NSUTOT
472 COMMON/EMDIE/IDIEYR(20),IDIEVH(20),EFDP(20).NDITOT
473 COMMON/EMPB/EFLEAO(20),IPBVR(20),IPBVH(20),IPBFT(20).NPBTOT
474 COMMON/GAS/LCGSIZ,LCGYR(20),PBCG(20.2),NCFEYR(27),CFE(27,6),
475 8.ICFESZ,IFECSZ,FECVS(13),FEC(13,2)
476 COMMON/MSFL/RMISFL(20,4.2),TMISFL(6,2)
477 COMMON/PSTEMP/PSDIST(5,2,6)
478 COMMON/SCENR/IVEAR.VEHSPD,IDRIV
479 COMMON/CONST/FPB.FPBE.FPBL
480 COMMON/TWF/TWFREG(20,6.3),ITWFMA(20,6,3)
481 COMMON/HDDT/IHTRG,IHDDMA(20,4),IHDDYR(21),HDDTRG(21.4)
482 COMMON/HDDV/IHDDMYC13).HDDVCFC13)
483 COMMON/EMCT/FPBEMCC15,6),FNLEMCC15,6,8),FDIEMC(15,6),IEMCYR(15),
484 &IEMCT,TAMPER(3,2)
485 COMMON/FLEET/FSF(26.6,3),IFSFYR(26.6),IFSF(6),VMT(6)
486 C
487 C LDV LDT1 LDT2 ORGANIC - 6 RECORDS
488 1001 FORMAT(I2)
489 READ(IOUDAT,1001)NORG1
490 NORTOT=NORG1
491 1000 FORMAT(I4,3X,I1.IX,I 1, IX.F5.3)
492 DO 10 1=1.NORG1
493 READ(IOUDAT.1000.END=999)IORGYR(I).IORGCT(I).IORGFT(I),EFO(I)
494 10 IORGVHd) = 1
495 C
496 C LDV LDT1 LDT2 SULFATES 13 RECORDS
497 C
498 1010 FORMAT(I4,3X,I1 .1X,I 1,1X.F4. 1,1X.F5.3)
499 READdOUDAT. 100DNSUL1
500 NSUTOT=NSUL1
501 DO 20 1=1.NSUL1
502 READdOUDAT.1010,END=999)ISULYR(I),ISULCT(I).ISULFT(I),
503 &EFS(I,2),EFS(I,1 )
504 20 ISULVH(I)=1
505 C
506 C LDDV DIESEL 3 RECORDS
507 C
508 READdOUDAT, 1001 )NDIES1
509 NDITOT=NDIES1
510 1020 FORMAT(I4,1X.F5.3)
511 DO 30 1=1.NDIES1
512 READdOUDAT, 1 020 , END=999) IDIEYRd ) . EFDP(I)
513 IDIEVH(I)=1
514 30 CONTINUE
515 C
516 C LOOT DIESEL 3 RECORDS
517 C
518 READdOUDAT, 1001 )NDIES2
519 NDITOT=NDITOT+NDIES2
520 NDINXT=NDIES1-H
521 DO 40 I=NDINXT,NDITOT
522 READdOUDAT, 1020 , END=999 ) IDIEYR (I) ,EFDP(I)
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD Page 10
523 IDIEVH(I)=2
524 40 CONTINUE
525 C
526 C HDGV ORGANIC 4 RECORDS
527 C
528 READ(IOUDAT,1001)NORG2
529 NORNXT=NORTOT+1
530 NORTOT=NORTOT+NORG2
531 DO 50 I=NORNXT.NORTOT
532 READdOUDAT, 1 000 , END=999) IORGVR (I ) , IORGCTU ) , IORGFT(I ) ,EFO(I)
533 50 IORGVH(I)=4
534 C
535 C HDGV SULFATE 3 RECORDS
536 C
537 READdOUDAT, 1001 )NSUL2
538 NSUNXT=NSUTOT+1
539 NSUTOT=NSUTOT+NSUL2
540 DO 60 I=NSUNXT,NSUTOT
541 READdOUDAT,1025,END=999)ISULVR(I),ISULCT(I),ISULFT(I),EFS(I , 1 )
542 1025 FORMAT(I4,3X,I 1,IX,I 1 ,1X.F5.3)
543 ISULVH(I)=4
544 60 EFS(I.2)=0.0
545 C
546 C HDDT DIESEL 1 RECORD
547 C
548 READdOUDAT, 1001 )NDIES3
549 NDINXT=NDITOT+1
550 NDITOT=NDITOT+NDIES3
551 DO 70 I=NDINXT,NDITOT
552 READdOUDAT, 1020, END=999) IDIEVR (I) ,EFDP(I )
553 70 IDIEVH(I)=4
554 C
555 C LDV LDT1 LEAD 6 RECORDS
556 C
557 READdOUDAT, 1001 ,END=999)IPB
558 1030 FORMAT(I4, IX,I 1,1X.F5.3)
559 NPBTOT=IPB
560 DO 80 1=1.NPBTOT
561 IPBVH(I)=1
562 80 READdOUDAT,1030,END=999)IPBYR(I),IPBFT(I),EFLEAD(I)
563 C
564 C LDT2 LEAD 6 RECORDS
565 C
566 READdOUDAT. 1001 ,END=999)IPB
567 NPBNXT=NPBTOT+1
568 NPBTOT=NPBTOT+IPB
569 DO 90 I=NPBNXT,NPBTOT
570 READdOUDAT,1030,END=999)IPBYR(I),IPBFT(I).EFLEAD(I)
571 90 IPBVH(I)=3
572 C
573 C MC LEAD 2 RECORDS
574 C
575 READdOUDAT, 1001 ,END=999)IPB
576 NPBNXT=NPBTOT+1
577 NPBTOT=NPBTOT+IPB
578 DO 95 I=NPBNXT,NPBTOT
579 READdOUDAT,1030,END=999)IPBVR(I),IPBFT(I),EFLEAD(I )
580 95 IPBVH(I)=6
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD Page 11
581 C
582 C READ TABLE 2-2 LEAD CONTENT OF GASOLINE BY YEAR
583 C
584 READ(IOUDAT,1001,END=999)LCGSIZ
585 1040 FORMAT(I4,1X.F4.2.1X.F5.3)
586 DO 300 1=1.LCGSIZ
587 300 READdOUDAT.1040,END=999)LCGYR(I),(PBCG(I,IFT),IFT=1,2)
588 C
589 C READ TABLE 2-3 LDV EMCT
590 C
591 CALL GTEMCT(I)
592 C
593 C READ TABLE 2-4 LDV FLEET SALES FRACTIONS
594 C
595 CALL FLTSLE(1)
596 C
597 C READ TABLE 2-5 LDV TRAVEL WEIGHTING
598 C
599 CALL TRAVWT(1,1)
600 C
601 C READ TABLE 2-6 CITY/HIGHWAY FUEL ECONOMY
602 C
603 READdOUDAT, 1001 , END=999) I CFESZ
604 1050 FORMAT(I4,1X,4(F4.1,IX),2(F3.1,IX))
605 DO 310 1=1.ICFESZ
606 READdOUDAT, 1050)NCFEYR (I ) . (CFE(I . JVC) , IVC=1 ,6)
607 IF(NCFEVR(I) .LE.1986)CFE(I,4)=CFE(I,6)
608 310 CONTINUE
609 C
610 C READ TABLE 2-7 FUEL ECONOMY SPEED CORRECTIONS
61 1 C
612 READdOUDAT, 1001 , END=999) I FECSZ
613 1060 FORMAT(F4.1.1X.F5.3,1X.F5.3)
614 DO 320 1=1.IFECSZ
615 320 READdOUDAT, 1060, END=999) FECVSd ) , (FECCI .ID) ,10=1 ,2)
616 C
617 C READ TABLES FOR LDT1 AND LDT2
618 C
619 DO 330 IVHCLS=2,3
620 CALL GTEMCT(IVHCLS)
621 CALL FLTSLE(IVHCLS)
622 CALL TRAVWTdVHCLS, 1 )
623 330 CONTINUE
624 C
625 C READ TABLE 2-14 LOOT TRAVEL WEIGHTING
626 C
627 CALL TRAVWT(2,3)
628 C
629 C READ TABLE 2-15 HDGV FLEET SALES
630 C
631 CALL FLTSLE(4)
632 C
633 C READ TABLE 2-16 HDGV TRAVEL WEIGHTING
634 C
635 CALL TRAVWT(4,1)
636 C
637 C READ TABLE 2-17 HDDT REGISTRATION VS AGE
638 C
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page 12
639 READUOUDAT, 1001 , END=999) IHDDT
640 1065 FORMAT(3X.F5.3)
641 DO 335 IAGE=1,IHDDT
642 335 READ(IOUDAT,1065)TWFREG(IAGE,4,3)
643 C
644 C READ TABLE 2-18 MOTORCYCLE TRAVEL WEIGHTING
645 C
646 CALL TRAVWT(6,1)
647 C
648 C READ TABLE 2-19 MISFUELING
649 C
650 1070 FORMAT(5X,F4.2.1X.F4.2)
651 DO 350 IVHCLS=1,6
652 IF (IVHCLS.EQ.5) GOTO 340
653 READdOUDAT, 1070 , END=999)WRGIM , WRGNIM
654 340 TMISFLdVHCLS, 1 )=WRGNIM
655 TMISFLCIVHCLS.2)=WRGIM
656 350 CONTINUE
657 C
658 C READ TABLE 2-19A MISFUELING
659 C
660 READdOUDAT, 1001 , END=999) IRMI S
661 1080 FORMAT(3X,8F4.2)
662 DO 360 IAGE=1,IRMIS
663 READdOUDAT, 1 080 , END=999) ( (RMI SFL(I AGE , J , K) , K= 1 , 2) , J= 1 , 4)
664 360 CONTINUE
665 C
666 C READ TABLE 2-20 PARTICLE SIZE ML. MNL.C,MNL.NC,MD,MB
667 C
668 DO 365 1=1,5
669 DO 365 J=1,2
670 DO 365 K=1,6
671 365 PSDIST(I,J,K)=0.0
672 2000 FORMAT(5(F5.2,IX))
673 DO 370 J=1 ,3
674 READdOUDAT, 2000, END=999)(PSDISTd ,1.J),I = 1,3)
675 READdOUDAT,2000,END=999)(PSDIST(I,2.J),1=1.3)
676 370 CONTINUE
67/ DO 380 J=4,5
678 READdOUDAT,2000,END=999)(PSDIST(I,1,J),I=1,5)
679 READdOUDAT, 2000, END=999)(PSDISTd .2.J),1 = 1.5)
680 380 CONTINUE
681 READdOUDAT. 2000, END=999)(PSDISTd ,1,6),1 = 1,2)
682 READdOUDAT, 2000, END=999)(PSDISTd .2.6). 1 = 1.2)
683 C
684 C READ TABLE 2-21 HDDV CONVERSION FACTORS
685 C
686 READdOUDAT. 1001 . END=999) IHDDT
687 2010 FORMAT(I4.1X.F6.4)
688 READdOUDAT,2010,END=999)(IHDDMY(I).HDDVCF(I),1=1.IHDDT)
689 C
690 C READ TABLE 2-22 FRACTION OF LEAD BURNED
691 C 1. NON-MISFUELED 2. 1975-1980 MISFUELED 3. 1981+ MISFUELED
692 C
693 2020 FORMAT(3(F4.2,IX))
694 800 READdOUDAT, 2020, END=999) FPB , FPBE , FPBL
695 C
696 C READ TABLE 2-23 CATALYST REMOVAL RATE (TAMPERING)
-------
Listing of PART.LEAD at 13:40:01 on DEC 2. 1985 for CCid=SHJD Page 13
697 C
698 2030 FORMAT(3(F5.3.1X))
699 READUOUDAT,2030.ENO=999)((TAMPER(I.J),I=1,3),J=1,2)
700 C
701 C READ TABLE A-1 AVE ANNUAL MILEAGE FOR HDDT
702 C
703 READ(IOUDAT.1001.END=999)IHDTMA
704 2040 FORMAT(3X,3(I5.IX),15)
705 DO 382 IAGE=1,IHDTMA
706 382 READ(IOUDAT.2040,END=999)(IHDDMA(IAGE.IWT),IWT=1,4)
707 C
708 C READ TABLE A-2 PROJECTIONS OF HEAVY DUTY VEHICLES IN OPERATION
709 C
710 READdOUDAT, 1001 , END=999) IHTRG
711 2050 FORMATU4. 1X.4(F5.3, IX))
712 DO 385 IYR=1,IHTRG
713 READ(IOUDAT,2p50,END=999)IHDDYR(IVR).(HDDTRG(IYR,IWT),IWT=1,4)
714 385 CONTINUE
715 C
716 C READ DEFAULT VMT MIX
717 C
718 READ(IOUDAT,2060)(VMT(IVHCLS),IVHCLS=1,6)
719 2060 FORMAT(6(F6.4,IX))
720 C
721 C SET OTHER DEFAULT FLEET SALES FRACTIONS
722 C
723 IFSF(6)=20
724 DO 405 1 = 1 ,20
725 C ALL HDDT USE DIESEL FUEL
726 FSF(I,4,3) = 1 .
727 C ALL MOTORCYCLES USE LEADED GAS
728 FSF(I,6,1) = 1 .
729 FSF(I,6.2)=0.0
730 FSF(I,6,3)=0.0
731 IFSFYR(I,6)=1900
732 405 CONTINUE
733 RETURN
734 999 CALL QUITER(0.,0,31,INERR)
735 RETURN
736 END
737 C
738 C
739 SUBROUTINE TWFCAL(IYEAR)
740 COMMON/TWF/TWFREG(20,6,3),ITWFMA(20.6,3)
741 COMMON/TWFRAC/TWFAC(20,6,3)
742 COMMON/HDDT/IHTRG,IHDDMA(20,4),IHDDYR(21),HDDTRG(21,4)
743 COMMON/FLEET/FSF(26.6,3),IFSFVRC26.6),IFSF(6),VMT(6)
744 DIMENSION TWFSUM(6,3)
745 C CALCULATE HDDT MILEAGE ACCUMULATION
746 C
747 IYRPTR=1
748 DO 390 1=1.IHTRG
749 390 IF(IHDDYR(I).GE.I YEAR.AND.IYRPTR.EQ.1)IYRPTR=I
750 ITWFMA(1,4,3)=0
751 DO 400 IAGE=2,20
752 HDDTWT=0.0
753 HDDSUM=0.0
754 DO 395 IWT=1,4
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD Page 14
755 HDDTWT=HDDTWT+HDDTRG(IYRPTR,IWT)
756 395 HDDSUM=HDDSUM-HHDDMA(IAGE-1,IWT)*HDDTRG(IYRPTR,IWT)
757 400 ITWFMA(IAGE,4,3)=HDDSUM/HDDTWT
758 C
759 C CALCULATE TRAVEL FRACTION BY AGE
760 C
761 DO 410 IVHCLS=1,6
762 DO 410 IFT=1,3
763 TWFSUM(IVHCLS,IFT)=0.0
764 DO 410 IAGE=1,20
765 410 TWFSUM(IVHCLS,IFT)=TWFREG(IAGE,IVHCLS,IFT)*
766 CITWFMA(IAGE,IVHCLS,IFT)+TWFSUM(IVHCLS,IFT)
767 DO 420 IVHCLS=1,6
768 DO 420 IFT=1,3
769 DO 420 IAGE=1,20
770 IF(TWFSUM(IVHCLS,IFT).NE.0.0)
771 8.TWFAC(IAGE,IVHCLS,IFT)=TWFREG(IAGE.IVHCLS.IFT)*
772 &ITWFMA(IAGE,IVHCLS,I FT)/TWFSUM(IVHCLS,I FT)
773 IF(TWFSUM(IVHCLS,IFT).EQ.0.0)TWFAC(I AGE,IVHCLS,I FT)=0.0
774 420 CONTINUE
775 RETURN
776 END
777 C
778 C
779 SUBROUTINE TRAVWT(IVHCLS,I FT)
780 C READ TRAVEL WEIGHTING FACTOR BY VEH.CLASS AND FUEL TYPE
781 C
782 COMMON/TWF/TWFREG(20,6,3),ITWFMA(20,6,3)
783 COMMON/IOUCOM/IOUREP,IOUERR,IOUASK,IOUDAT,IOUGEN
784 1300 FORMATU2. 1X.F5.3, 1X.I5)
785 DO 10 1=1,20
786 10 READ(IOUDAT,1300)IAGE,TWFREG(I.IVHCLS,I FT),ITWFMA(I.IVHCLS,I FT)
787 IF(IFT.EQ. DGOTO 40
788 IF(IVHCLS.EQ.2.AND.IFT.EQ.3)GOTO 20
789 GOTO 60
790 C SOME CLASSES HAVE EQUIV. WEIGHTING
791 C LDDT I AND LDDT II
792 20 DO 30 1=1,20
793 TWFREGU,3,IFT)=TWFREG(I.IVHCLS.I FT)
794 30 ITWFMAd,3,IFT)=ITWFMA(I,IVHCLS,IFT)
795 GOTO 60
796 40 DO 50 1=1,20
797 C LDDV AND LDV
798 IF(IVHCLS.EQ.1)TWFREGU.IVHCLS,3)=TWFREG(I.IVHCLS,1)
799 IF(IVHCLS.EQ.1)ITWFMA(I,IVHCLS,3)=ITWFMA(I,IVHCLS,1)
800 C HDGV1 AND HDGV2
801 IF(IVHCLS.EQ.4)TWFREG(I.5,1)=TWFREG(I.IVHCLS,1)
802 I F(IVHCLS.EQ.4) ITWFMAd,5,1)=ITWFMA(I,IVHCLS,1)
803 C LEAD AND UNLEAD
804 TWFREG(I.IVHCLS.2)=TWFREG(I.IVHCLS.1)
805 ITWFMAd,IVHCLS,2)=ITWFMA(I.IVHCLS,1)
806 50 CONTINUE
807 60 RETURN
808 END
809 C
810 C
811 SUBROUTINE GTEMCT(IVHCLS)
812 C READ FRACTION OF VEHICLES WITH VARIOUS EMISSION CONTROL DEVICES
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD Page 15
813 C
814 COMMON/EMCT/FPBEMC(15,6),FNLEMC(15.6,8},FDIEMC(15,6),IEMCVR(15),
815 &IEMCT,TAMPER(3,2)
816 COMMON/IOUCOM/IOUREP,IOUERR,IOUASK.IOUDAT,IOUGEN
817 1101 FORMAT(I 2)
818 READdOUDAT, 1101 HEMCT
819 C READ LEADED FRACTION
820 1100 FORMAT(15(F5.3,IX))
821 READdOUDAT, 1 100) (FPBEMC( I , IVHCLS) ,1 = 1 .IEMCT)
822 C UNLEADED FRACTIONS
823 DO 20 1=1,5
824 20 READdOUDAT, 1 100) (FNLEMC( J , I VHCLS, I) ,J=1 , IEMCT)
825 C DIESEL FRACTION
826 READdOUDAT, 1 100) (FDIEMC (J . I VHCLS) , J=1 , IEMCT)
827 C MORE UNLEADED FRACTIONS
828 DO 30 1=6,8
829 30 READdOUDAT, 1 1 00) (FNLEMC( J , I VHCLS, I) ,J=1 , IEMCT)
830 IEMCYRC1)=1900
831 DO 40 1=2,15
832 40 IEMCVRd) = 1973+1
833 RETURN
834 END
835 C
836 C
837 SUBROUTINE FLTSLE(IVHCLS)
838 C READ FLEET SALES
839 COMMON/IOUCOM/IOUREP.IOUERR,IOUASK,IOUDAT, IOUGEN
840 COMMON/FLEET/FSF(26.6,3) , IFSFYR(26,6),IFSF(6),VMT(6)
841 READdOUDAT, 1200)NFSF
842 IFSF(IVHCLS)=NFSF
843 1200 FORMAT(I2)
844 1210 FORMAT(I4,1X,3(F5.3,IX))
845 DO 10 I=1,NFSF
846 IF(IVHCLS.LE.3)READ(IOUDAT,1210)IFSFYR(I,IVHCLS),
847 &FSF(I,IVHCLS,2),FSF(I,IVHCLS,1),FSF(I,IVHCLS,3)
848 IF (I VHCLS. EQ. 4) READdOUDAT, 1 2 10) IFSFYR( I , I VHCLS) ,
849 &FSF(I,IVHCLS,2), FSF(I ,IVHCLS,1)
850 10 CONTINUE
851 40 RETURN
852 END
853 C
854 C
855 SUBROUTINE ONESEC(INERR,*)
856 C ONE TIME DATA SECTION
857 INTEGER TAMFLG,SPCFLG,VMFLAG
858 COMMON/FLAGS 1/1 PROMT,TAMFLG.SPCFLG,VMFLAG,MYMRFG,IMFLAG,MISFLG
859 COMMON/EMCT/FPBEMCC15.6),FNLEMC(15,6,8),FDIEMCC15,6),IEMCYR(15).
860 &IEMCT,TAMPER(3,2)
861 COMMON/MSFL/RMISFL(20.4.2),TMISFL(6.2)
862 COMMON/IOUCOM/IOUREP,IOUERR,IOUASK,IOUDAT,IOUGEN
863 COMMON/PSIZE/PSZCUT.PSL.PSNL.PSNLCT.PSDIE,PSBRK,PSTIRE
864 COMMON/TWF/TWFREG(20,6,3),ITWFMA(20.6,3)
865 C
866 DIMENSION CLSNM(6)
867 C
868 DATA CLSNM/' LDV','LDT1'.'LDT2','HDT1'.'HDT2',' MC '/
869 C
870 C ALTERNATE TAMPERING RATES
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD Page 16
871 IF(TAMFLG.EQ.DGOTO 60
872 IF(IPROMT.EQ.2)WRITE(IOUASK,200)
873 200 FORMAT(-0','ENTER FRACTION OF CATALYST EQUIPPED VEHICLES WITH'/
874 &'CATALYST REMOVED FOR LDV, LDT1 AND LDT2 FORMAT=3(F4.3,1X)')
875 IF(IMFLAG.EQ. 1.AND.I PROMT.EQ.2. )WRITE(IOUASK,210)
876 IF(IMFLAG.EQ.2.AND.I PROMT.EQ.2.)WRITE(IOUASK,220)
877 210 FORMAT('NO I/M PROGRAM IN IN EFFECT')
878 220 FORMATCI/M PROGRAM IS IN EFFECT')
879 READ(IOUGEN,225)(TAMPER(I,IMFLAG),1=1.3)
880 225 FORMAT(3(F4.3,IX))
881 DO 5 1=1,3
882 IF(TAMPER(I.IMFLAG).LT.0.0)CALL QUITER(TAMPER(I,IMFLAG),0,30,INERR)
883 5 IF(TAMPER(I,IMFLAG).LT.0.0)TAMPER(I,IMFLAG)=0.0
884 C
885 C ALTERNATE VEHICLE MILEAGE DISTRIBUTION
886 C
887 60 IF(VMFLAG.NE.3)GOTO 10
888 IF(I PROMT.EQ.2)WRITE(IOUASK,272)
889 272 FORMATC ENTER VMT FRACTIONS FOR LDV, LDT1, LDT2, HDGV, HDDT, MC
890 C -6(F5.3,1X)')
891 CALL GETVMT(INERR,&99)
892 C ALTERNATE REGISTRATION AND MILEAGE ACCUMULATION
893 C
894 10 IF(MYMRFG.EQ.1)GOTO 65
895 IF(MYMRFG.LT.3)GOTO 30
896 IF(IPROMT.EQ.2)WRITE(IOUASK,230)
897 230 FORMAT('0','ALTERNATE REGISTRATION '/
898 &' ENTER VEHICLE CLASS(1-6) & FUEL TYPE(1-3) FORMAT — I1,1X,I1'/
899 &• ENTER 0,0 AS LAST VEHICLE CLASS AND FUEL TYPE')
900 READ(IOUGEN,235)IVHCLS,IFT
901 235 FORMATUl.1X.il)
902 IF(IVHCLS.LE.O)GOTO 30
903 IF(IVHCLS.LT.1.OR.IVHCLS.GE.6)CALL QUITER(0.,IVHCLS,37.INERR)
904 IF(IFT.LT. 1 .OR.IFT.GT.3)CALL QUITER(0.,IFT,38,INERR)
905 IF(IVHCLS.EQ.6)IFT=1
906 IF(IVHCLS.EQ.4)IFT=1
907 IF(IVHCLS.EQ.5)IFT=3
908 IF(IVHCLS.EQ.5)IVHCLS=4
909 IF (I PROMT.EQ.2)WRITE(IOUASK,240)CLSNM(IVHCLS)
910 240 FORMATC ENTER ',A4,' FRACTION OF REGISTRATION FOR AGES 1-10',
911 8.' FORMAT — 10F5.3')
912 READ(IOUGEN,260,END=99)(TWFREGdAGE,IVHCLS,I FT),IAGE=1,10)
913 IF(IPROMT.EQ.2)WRITE(IOUASK,250)CLSNM(IVHCLS)
914 250 FORMATC ENTER ',A4,' FRACTION OF REGISTRATION FOR AGES 11-20+',
915 &' FORMAT — 10F5.3')
916 READ(IOUGEN,260,END=99)(TWFREGCIAGE,IVHCLS,IFT),IAGE=11,20) '
917 260 FORMAT(10F5.3)
918 CHKREG=0.0
919 DO 25 IAGE=1 ,20
920 IF (TWFREGd , IVHCLS, IFT) . LT . 0 . )
921 * CALL QUITER(TWFREG(I.IVHCLS,IFT),0,13,INERR)
922 25 CHKREG=CHKREG+TWFREG(IAGE,IVHCLS, IFT)
923 IF (ABS(CHKREG-1.).GE.0.001)CALL QUITER(CHKREG,0,20,INERR)
924 . IF(ABS(CHKREG-1 . ) .LT.OODGOTO 10
925 DO 27 IAGE=1,20
926 TWFREGdAGE,IVHCLS.IFT)=TWFREG(IAGE,IVHCLS,IFT)/CHKREG
927 IFCIFT.EQ.1)TWFREG(IAGE,IVHCLS,2)=TWFREG(IAGE,IVHCLS.IFT)
928 27 IF(IFT.EQ.2)TWFREG(IAGE,IVHCLS,1)=TWFREG(IAGE,IVHCLS,IFT)
-------
Listing of PART.LEAD at 13:40:01 on DEC 2. 1985 for CCid=SHJD . Page 17
929 GOTO 10
930 C
931 C ALTERNATE MILEAGE ACCUMULATION RATES
932 C
933 30 IF(MVMRFG.NE.2.AND.MYMRFG.NE.4)GOTO 65
934 IF(IPROMT.EQ.2)WRITE(IOUASK,265)
935 265 FORMAT('0','ALTERNATE MILEAGE ACCUMULATION RATES'/
936 &' ENTER VEHICLE CLASS(1-6) & FUEL TYPE(1-3) FORMAT ~ I1.1X.I1'/
937 &' ENTER 0,0 AS LAST VEHICLE CLASS AND FUEL TYPE')
938 READ(IOUGEN,235)IVHCLS,IFT
939 IF(IVHCLS.LE.O)GOTO 65
940 IF(IVHCLS.GT.6)CALL QUITER(0..IVHCLS,39,INERR)
941 IF(IFT.LT.1.OR.IFT.GT.3)CALL QUITER(0.,IFT,40.INERR)
942 IF(IVHCLS.EQ.6)IFT=1
943 IF(IVHCLS.EQ.4)IFT=1
944 IF(IVHCLS.EQ.5)IFT=3
945 IF(IVHCLS.EQ.5)IVHCLS=4
946 IF(IPROMT.EQ.2)WRITE(IOUASK,270)CLSNM(IVHCLS)
947 267 FORMAT(9(I5,IX).15)
948 270 FORMATC'0','ENTER ANNUAL MILEAGE ACCUMULATION FOR ',A4.
949 &' AGES 1-10 FORMAT -- 10(15.IX)')
950 271 FORMAT('0'.'ENTER ANNUAL MILEAGE ACCUMULATION FOR ',A4.
951 8.' AGES 11-20+ FORMAT 10(15, IX)')
952 READ(IOUGEN,267.END=99)(ITWFMA(IAGE,IVHCLS,IFT),IAGE=1.10)
953 IF(IPROMT.EQ.2)WRITE(IOUASK,271)CLSNM(IVHCLS)
954 READ(IOUGEN,267,END=99)(ITWFMA(I AGE,IVHCLS,IFT),IAGE=11,20)
955 DO 50 IAGE=1 ,20
956 IF(ITWFMA(IAGE,IVHCLS,IFT).EQ.0.AND.TWFREG(IAGE,IVHCLS,I FT).NE.0.0)
957 & CALL QUITER(TWFREG(IAGE,IVHCLS,IFT),0, 14,INERR)
958 IF(ITWFMA(IAGE,IVHCLS,I FT).NE.0.AND.TWFREGCIAGE,IVHCLS.I FT).EQ.0.0)
959 & CALL QUITER(TWFREG(IAGE,IVHCLS.IFT),0,15,INERR)
960 50 IF(ITWFMA(IAGE,IVHCLS.IFT).LT.O) CALL
961 * QUITERCO.,ITWFMA(IAGE.IVHCLS,IFT),12,INERR)
962 GOTO 30
963 C
964 C
965 C MISFUELING BY VEHICLE CLASS MILEAGE ACCUMULATION
966 C
967 65 IF(MISFLG.EQ.1.OR.MISFLG.EQ.3) GOTO 80
968 DO 70 IVHCLS=1,4
969 IF(I PROMT.EQ.2)WRITE(IQUASK,275)CLSNM(IVHCLS)
970 275 FORMAT('ENTER FOR ',A4,' ZERO MILE MISFUELING RATE AND SLOPE
971 & —FORMAT=F8.6,1X.F8.6')
972 READ(IOUGEN,280)RMISZM.RMISSL
973 280 FORMAT(F8.6,1X.F8.6)
974 IF(RMISZM.LT.O)CALL QUITER(RMISZM,0,17,INERR)
975 IF(RMISSL.LT.O)CALL QUITER(RMISSL,0,18,INERR)
976 ITMSUM=0.0
977 DO 70 IAGE=1 ,20
978 ITMSUM=ITMSUM+ITWFMA(IAGE,IVHCLS,2)
979 RMISFL(IAGE,IVHCLS,IMFLAG)=RMISZM+RMISSL*ITMSUM
980 70 IF(RMISFL(IAGE,IVHCLS,IMFLAG).GT.1.0)RMISFL(IAGE,IVHCLS,IMFLAG)=1.0
981 C
982 C
983 C ALTERNATE LEAD CONTENT OF GASOLINE
984 C
985 80 IF(SPCFLG.EQ.2)CALL GETLCG(INERR,&99)
986 90 RETURN
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page
18
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
1001
1002
1003
1004
1005
1006
1007
1008
1009
1010
101 1
1012
1013
1014
1015
1016
1017
1018
1019
1020
1021
1022
1023
1024
1025
1026
1027
1028
1029
1030
1031
1032
1033
1034
1035
1036
1037
1038
1039
1040
1041
1042
1043
1044
99 RETURN1
END
SUBROUTINE PARSEC(INERR,*)
C PARAMETER SECTION AND NORMAL USER INPUTS
C
INTEGER TAMFLG.SPCFLG.VMFLAG
COMMON/IOUCOM/IOUREP,IOUERR,IOUASK,IOUDAT,IOUGEN
COMMON/SCENR/IYEAR.VEHSPD,IDRIV
COMMON/PSIZE/PSZCUT,PSL,PSNL,PSNLCT,PSDIE,PSBRK,PSTIRE
COMMON/FLAGS 1/I PROMT,TAMFLG,SPCFLG.VMFLAG.MYMRFG,IMFLAG.MISFLG
C READ YEAR FOR WHICH CALULATION IS TO BE DONE
80 IF(I PROMT.EQ.2)WRITE(IOUASK,290)
290 FORMAT('-','ENTER YEAR FOR THIS RUN — (FORMAT = 14)')
READ(IOUGEN.300.END=99)IYEAR
300 FORMAT(I4)
IF(IYEAR.LT.1980)CALL QUITER(0. , IYEAR,36,INERR)
CALL TWFCAL(IVEAR)
C READ VEHICLE SPEED
IF(I PROMT.EQ.2)WRITE(IOUASK,310)
310 FORMAT('ENTER VEHICLE SPEED — (FORMAT = F4.1)')
READ(IOUGEN,320,END=99)VEHSPD
320 FORMAT(F4.1)
IF(VEHSPD.LT.0.0)CALL QUITER(VEHSPD,0,19,INERR)
IF(VEHSPD.LT.0.0)VEHSPD=-1.*VEHSPD
IF(VEHSPD.LT.5.0)CALL QUITER(VEHSPD,0,21,INERR)
IF(VEHSPD.LT.S.O) VEHSPD=5.0
IF(VEHSPD.GT.55.0)CALL QUITER(VEHSPD,0.9,INERR)
IF(VEHSPD.GT.55.0) VEHSPD=55.0
IF(VEHSPD.GT.34.8)CALL QUITER(VEHSPD,0,34,INERR)
IF(VEHSPD.LT.19.6)CALL QUITER(VEHSPD,0,33,INERR)
C READ DRIVING CONDITIONS
90 IF(I PROMT.EQ.2)WRITE(IOUASK,330)
330 FORMATC ENTER DRIVING CONDITIONS (1= TRANSIENT, 2= CRUISE)')
READ(IOUGEN,340,END=99)IDRIV
340 FORMAT(I1 )
IF(IDRIV.LE.O.OR.IDRIV.GT.2) CALL QUITER(0.,IDRIV,35,INERR)
C READ PARTICLE SIZE CUTOFF
100 IF(IPROMT.EQ.2)WRITE(IOUASK,350)
350 FORMATC ENTER PARTICLE SIZE CUTOFF (0.1-10.) -FORMAT(F5.2) ' )
READ(IOUGEN,360,END=99)PSZCUT
360 FORMAT(F5.2)
IF(PSZCUT.LT.O.1)CALL QUITER(PSZCUT,0,26,INERR)
IF(PSZCUT.LT.O.1)PSZCUT=0.1
IF(PSZCUT.GT.10.)CALL QUITER(PSZCUT,0.27.INERR)
IF(PSZCUT.GT.10.)PSZCUT=10.
CALL PARSIZ
IF(VMFLAG.NE.2) GOTO 95
IF(I PROMT.EQ.2)WRITE(IOUASK,272)
272 FORMATC ENTER VMT FRACTIONS FOR LDV, LDT1, LDT2, HDGV. HDDT, MC
C -6(F5.3,IX)')
CALL GETVMT(INERR,&99)
95 RETURN
99 CALL QUITER(0.,0,23,INERR)
RETURN
END
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page
19
1045
1046
1047
1048
1049
1050
1051
1052
1053
1054
1055
1056
1057
1058
1059
1060
1061
1062
1063
1064
1065
1066
1067
1068
1069
1070
1071
1072
1073
1074
1075
1076
1077
1078
1079
1080
1081
1082
1083
1084
1085
1086
1087
1088
1089
1090
1091
1092
1093
1094
1095
1096
1097
1098
1099
1 100
1 101
1 102
C
C
SUBROUTINE PARSIZ'
C INTERPOLATES PARTICLE SIZE DISTRIBUTION
C
COMMON/ PSIZE/PSZCUT.PSL.PSNL.PSNLCT.PSDIE.PSBRK.PSTI RE
COMMON/PSTEMP/PSDIST(5,2,6)
C
DO 60 ISRC=1 .6
PS=0.0
MXSIZE=5
IF (PSZCUT. LE.PSDISTf 1 . 1 .ISRC)) GOTO 40
DO 10 1=1.5
10 IF(PSDIST(I , 1 .ISRC) .EQ.O.O)MXSIZE=I
IF(MXSIZE.EQ. 1 )PS=PSDIST( 1 ,2, 1 SRC)
IF(PSZCUT.GE.PSDIST(MXSIZE, 1 , ISRC) )PS=PSDIST(MXSIZE , 2 , ISRC)
IF(PS.NE.O.O)GOTO 50
SLOPE=0.0
NPTR=1
DO 20 I=2,MXSIZE
20 IF(PSDIST(I , 1 .ISRC) .GE . PSZCUT . AND. NPTR . EQ. 1)NPTR=I
SLOPE=(PSDIST(NPTR,2, I SRC )-PSDIST( NPTR- 1 ,2, ISRC))/
&(PSDIST(NPTR, 1 ,ISRC)-PSDIST(NPTR-1 , 1 .ISRC))
PS=SLOPE*(PSZCUT-PSDIST(NPTR-1 . 1 , I SRC) ) + PSDI ST(NPTR- 1 ,2,
GOTO 50
40 SLOPE=(PSDIST(2,2,ISRC)-PSDIST(1 ,2, ISRC))/
8,(PSDIST(2. 1 ,ISRC)-PSDIST(1 . 1 ,ISRC))
PS=PSDIST( 1,2. ISRC)-SLOPE*(PSDIST(1 , 1 . I SRC) -PSZCUT)
50 IFUSRC.EQ. 1)PSL=PS
IF(ISRC.EQ.2)PSNLCT=PS
IF(ISRC.EQ.3)PSNL=PS
IF(ISRC.EQ.4)PSDIE=PS
IF(ISRC.EQ.5)PSBRK=PS
IF(ISRC.EQ.6)PSTIRE=PS
60 CONTINUE
RETURN
END
C
SUBROUTINE GETLCGdNERR,*)
C CHANGE LEAD CONTENT OF GASOLINE
COMMON/GAS/LCGSIZ.LCGYR(20) ,PBCG(20,2) ,NCFEYR(27) ,CFE(27
&ICFESZ,IFECSZ,FECVS( 1 3) , FEC( 1 3 , 2)
COMMON/ IOUCOM/ IOUREP , IOUERR . IOUASK , IOUDAT , IOUGEN
COMMON/ FLAGS 1 / 1 PROMT , TAMFLG . SPCFLG , VMFLAG , MYMRFG , IMFLAG .
IF(IPROMT.EQ.2)WRITE(IOUASK, 190)
190 FORMAT( '0' , 'ALTERNATE LEAD CONTENT OF GASOLINE BY YEAR')
IF (I PROMT. EQ. 2) WRITE (IOUASK. 195)
ISRC)
.6).
MISFLG
195 FORMATC ENTER FIRST AND LAST YEAR TO CHANGE - (I4.1X.I4)')
5 READ (IOUGEN, 196 . END=99 ) ISTART , I FIN
196 FORMAT (14. IX, 14)
IF( ISTART. LT. 1974 .OR . I START .GT . I FIN) CALL QUITER(0. ,0,16,
IPTR=1
DO 6 1=1 .LCGSIZ
6 IF(LCGYR(I) .LT. ISTART) I PTR=I
LCGSIZ=IPTR
IF (I PROMT. EQ. 2) WRITE (IOUASK. 200)
200 FORMAT( '0' . 'ENTER LEAD CONTENT OF GASOLINE FOR LEADED
8. AND UNLEADED — FORMAT=F5 . 3 , IX , F5 . 3 ' )
INERR)
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page 20
1 103
1 104
1 105
1 106
1 107
108
109
10
1 1
12
13
14
15
16
17
18
19
120
121
122
123
124
125
126
1 127
1 128
1 129
1 130
1 131
132
133
134
135
136
137
138
139
1 yl n
1 4U
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
DO 20 IVR=ISTART, IFIN
10 I F ( I PROMT . EQ . 2 ) WR I TE ( I OUASK , 2 1 0 ) I YR
210 FORMATC FOR YEAR '.15)
LCGSIZ=LCGSIZ+1
LCGYR(LCGSIZ)=IYR
READ(IOUGEN,220.END=99)(PBCG(LCGSIZ,IFT) ,IFT=1 ,2)
220 FORMAT(F5.3, 1X.F5.3)
IF(PBCG(LCGSIZ, 1) . LT. 0.0) CALL QUITER(PBCG(LCGSI Z , 1) ,IYR.
IF(PBCG(LCGSIZ, 1) . LT . 0 . 0)PBCG( LCGSI Z , 1 )=-1 . *PBCG( LCGSI Z ,
IF(PBCG(LCGSIZ.2) .LT. 0.0) CALL QUITER (PBCG( LCGSI Z , 2) , IYR,
IF(PBCG(LCGSIZ,2) . LT . 0 . 0) PBCG( LCGSIZ , 2)=- 1 . *PBCG( LCGSI Z ,
20 CONTINUE
RETURN
99 RETURN 1
END
SUBROUTINE GETVMT( INERR ,*)
GETVMT gets the fleet vehicle miles travelled (vmt) mix by
Called by ONESEC and PARSEC.
Cal Is QUITER.
Input on cal 1 :
parameter list: INERR
common blocks:
/FLAGS1/ IPROMT
/IOUCOM/ IOUGEN.IOUASK
Output on return:
Local variable / array dictionary:
Name Type Description
VMTCHK R vmt mix check: sum of the vmt fractions = 1.0
(i.e., a normalized vmt distribution has been
Notes:
GETVMT is a new subprogram (added for MOBILES).
INTEGER TAMFLG.SPCFLG.VMFLAG
COMMON/ FLAGS 1 / I PROMT , TAMFLG , SPCFLG , VMFLAG , MYMRFG , IMFLAG ,
COMMON/ IOUCOM/ I OUREP.IOUERR, I OUASK, I OUDAT.IOUGEN
COMMON/FLEET/FSF(26,6,3) , IFSFYR ( 26, 6) ,IFSF(6) ,VMT(6)
I F ( I PROMT . EQ . 2 ) WRI TE ( IOUASK . 200 )
200 FORMAT( '0' . 'Enter VMT split:')
READ(IOUGEN, 100,END=99) VMT
100 FORMAT(6(F5.3, IX))
VMTCHK=0.0
28, INERR)
1)
29, INERR)
2)
vehicle type
+/- 1.E-6
entered)
MISFLG
-------
Listing of PART.LEAD at 13:40:01 on DEC 2. 1985 for CCid=SHJD
Page 21
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
1200
1201
1202
1203
1204
1205
1206
1207
1208
1209
1210
121 1
1212
1213
1214
1215
1216
1217
1218
DO 10 IV=1 ,6
VMTCHK= VMT (IV) + VMTCHK
IF(VMT(IV) .LT.O.O.OR.VMT(IV) .GT. 1 . )
* CALL QUITER(VMTUV) ,0, 10.INERR)
10 CONTINUE
IF(VMTCHK+1 .E-6.LT. 1 . .OR. VMTCHK- 1 .E-6.GT. 1 . )
* CALL QUITER(VMTCHK,0, 1 1 .INERR)
C
RETURN
C
C EOF on any attempted read => take alternate return 1 => run
C
99 RETURN 1
END
C
aborts .
C ****************************************************************
C *
C * LEAD EMISSION FACTOR SUBROUTINES
C *
C *************************************************************
SUBROUTINE LEAD( I YEAR , I VRPTR , IVHCLS , I FT , VEHSPD , IDRIV.EF)
C CALCULATES LEAD EMISSION FACTOR FOR GIVEN VEHICLE. CLASS , YEAR
C FUEL TYPE AND VEHICLE SPEED
C
COMMON/ PSIZE/PSZCUT , PSL , PSNL , PSNLCT , PSDI E , PSBRK , PSTIRE
COMMON/CONST/FPB , FPBE , FPBL
C
EF=0.0
PBRAT=1 .557
IF(IFT.GT.2)GOTO 99
AS=FPB
AMIS=FPBE
IF(IYRPTR.GE. 198 1 ) AMIS=FPBL
CALL GTPBCGUYEAR, 1 ,PBL)
CALL GTPBCG(IYEAR, 2, PBNL)
CALL GETCS(IVHCLS, VEHSPD, IDRIV.CS)
CALL GETCFE(IVRPTR.IVHCLS.FLECON)
IF ( IVHCLS. GT. 3 )GOTO 20
C LEADED LDV, LDT1, LDT2
IF(IFT.EQ.2)GOTO 10
CALL GETPB(IYRPTR, IVHCLS, 1 ,EFPB1 )
CALL GETPB(IYRPTR,IVHCLS,2,EFPB2)
EF=(PBL*EFPB1*PSL+PBNL*EFPB2*PSNL)*AS*PBRAT/(FLECON*CS)
GOTO 99
C UNLEADED LDV, LDT1, LDT2
10 IF(IYRPTR.LT. 1975)GOTO 99
CALL NLCONS(IYEAR,IVRPTR, I VHCLS .TAMFRC , RMIS)
FNOCAT=FREMCT( I VRPTR, IVHCLS, 2, 1)
FCAT=FREMCT(IYRPTR, IVHCLS, 2, 6)
*
*
*
***
,
EF=PBNL*( 1 .-RMIS)*PSNLCT*AS+(PBL*RMIS*PSL*AS)*(FNOCAT+TAMFRC*FCAT)
&+PBL*RMIS*PSL*( 1 .-TAMFRC)*FCAT*AMIS
EF=EF*PBRAT/(FLECON*CS)
GOTO 99
C HDGV
20 IF(IVHCLS.GT.5)GOTO 40
IF(IFT.EQ.2)GOTO 30
C LEADED HDGV
EF=AS*PBL*PBRAT*PSL/FLECON
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page 22
1219
1220
1221
1222
1223
1224
1225
1226
1227
1228
1229
1230
1231
1232
1233
1234
1235
1236
1237
1238
1239
1240
1241
1242
1243
1244
1245
1246
1247
1248
1249
1250
1251
1252
1253
1254
1255
1256
1257
1258
1259
1260
1261
1262
1263
1264
1265
1266
1267
1268
1269
1270
1271
1272
1273
1274
1275
1276
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
GOTO 99
UNLEADED HDGV 1987+
30 IF(IYRPTR.LT.1987.OR.IVHCLS.EQ.5)GOTO 99
CALL NLCONSUYEAR, I VRPTR . VEHSPD . TAMFRC , RMI S)
EF=((1.-RMIS)*AS*PBNL*PBRAT*PSNLCT/FLECON)+(RMIS*AMIS*PBL*
&PBRAT»PSL/FLECON)
GOTO 99
MOTORCYCLES
40 CALL GETPBUYRPTR, IVHCLS, 1 ,EFL)
EF=EFL*PSL
99 RETURN
END
SUBROUTINE GTPBCGCIYRPTR.I FT,PBCON)
GET LEAD CONTENT OF LEADED AND UNLEADED GASOLINE
COMMON/GAS/LCGSIZ.LCGYR(20),PBCG(20,2),NCFEYR(27),CFE(27,6),
&ICFESZ,IFECSZ,FECVS(13),FEC(13,2)
PBCON=0.0
NPTR=1
DO 10 1=1.LCGSIZ
10 IF(IYRPTR.GE.LCGYR(I))NPTR=I
PBCON=PBCG(NPTR,IFT)
RETURN
END
SUBROUTINE GETPB(IYRPTR,IVHCLS,IFT.EFPB)
FIND LEADED EMISSION FACTORS
COMMON/EMPB/EFLEAD(20),IPBYR(20),IPBVH(20),IPBFT(20).NPBTOT
EFPB=0.0
IVH=IVHCLS
IF(IVHCLS.EQ.2)IVH=1
DO 10 1=1.NPBTOT
IF(IVH.NE.IPBVH(I))GOTO 10
IF(IFT.NE.IPBFT(I))GOTO 10
IF(IYRPTR.GE.IPBYR(I))EFPB=EFLEAD(I)
10 CONTINUE
30 RETURN
END
SUBROUTINE GETCS(IVHCLS.VEHSPD.IDRIV.CS)
GET SPEED DEPENDENT CORRECTION FACTOR
COMMON/GAS/LCGSIZ,LCGYR(20),PBCG(20,2),NCFEYR(27),CFE(27,6) ,
8.ICFESZ.IFECSZ,FECVS(13) , FEC( 13.2)
CS=0.0
IF(IVHCLS.GT.3)GOTO 30
IF(VEHSPD.LE.FECVS(1))CS=FEC( 1 .IDRIV)
IF(VEHSPD.GE.FECVS(IFECSZ))CS=FEC(IFECSZ,IDRIV)
IF(CS.NE.O)GOTO 30
NPTR=1
DO 10 I=2,IFECSZ
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page 23
1277
1278
1279
1280
1281
1282
1283
1284
1285
1286
1287
1288
1289
1290
1291
1292
1293
1294
1295
1296
1297
1298
1299
1300
1301
1302
1303
1304
1305
1306
1307
1308
1309
1310
131 1
1312
1313
1314
1315
1316
1317
1318
1319
1320
1321
1322
1323
1324
1325
1326
1327
1328
1329
1330
1331
1332
1333
1334
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
10 IF(FECVS(I).GE.VEHSPD.AND.NPTR.EQ.1)NPTR=I
SLOPE=(FEC(NPTR.IDRIV)-FEC(NPTR-1,IDRIV))/
&(FECVS(NPTR)-FECVS(NPTR-1))
CS=FEC(NPTR-1.IDRIV)+SLOPE*(VEHSPD-FECVS(NPTR-1
30 RETURN
END
SUBROUTINE GETCFE(IVRPTR.IVHCLS.FLECON)
GET COMBINED CITY/HIGHWAY FUEL ECONOMY
COMMON/GAS/LCGSIZ.LCGYRC20),PBCG(20.2),NCFEYR(27),CFE(27,6),
&ICFESZ,IFECSZ,FECVS(13),FEC(13,2)
FLECON=0.0
NPTR=1
IVH=IVHCLS
IF(IVHCLS.EQ.5.AND.IYRPTR.LT.1987)IVH=4
IF(IVHCLS.EQ.6)GOTO 40
DO 10 1=1.ICFESZ
10 IF(IVRPTR.GE.NCFEYR(I))NPTR=I
FLECON=CFE(NPTR.IVH)
40 RETURN
END
SUBROUTINE NLCONSCI YEAR,IYRPTR,IVHCLS.TAMFRC,RMIS)
UNLEADED EMISSIONS FACTORS CONSTANTS
INTEGER TAMFLG.SPCFLG.VMFLAG
COMMON/FLAGS 1/IPROMT,TAMFLG,SPCFLG,VMFLAG.MYMRFG,IMFLAG.MISFLG
COMMON/MSFL/RMISFL(20,4,2),TMISFL(6,2)
COMMON/EMCT/FPBEMC(15,6).FNLEMC(15,6,8),FDIEMC(15.6),IEMCYR(15),
&IEMCT,TAMPER(3.2)
TAMFRC=0.0
RMIS=0.0
IAGE=IYEAR-IYRPTR
IFUYRPTR.LT. 1975)GOTO 20
UNLEADED LDV, LDT1.LDT2, HDGV
IF(IVHCLS.LE.3)TAMFRC=TAMPER(IVHCLS.IMFLAG)
10 IF(IVHCLS.GT.4)RMIS=0.0
IF(IVHCLS.LE.4)RMIS=RMISFL(IAGE,IVHCLS,IMFLAG)
IF(MISFLG.EQ.3.AND.IVHCLS.LE.3)RMIS=TMISFL(IVHCLS,IMFLAG)
IF(IVHCLS.EQ.4.AND.IYRPTR.LT.19B7)RMIS=0.0
IF(RMIS.LT.0.0)RMIS=0.0
20 RETURN
END
********«******<
* ORGANIC EMISSION FACTORS SUBROUTINES
*
SUBROUTINE ORGANC(I YEAR,IYRPTR,IVHCLS,I FT,EFORG)
FIND ORGANIC EMISSIONS FACTORS
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page 24
1335
1336
1337
1338
1339
1340
1341
1342
1343
1344
1345
1346
1347
1348
1349
1350
1351
1352
1353
1354
1355
1356
1357
1358
1359
1360
1361
1362
1363
1364
1365
1366
1367
1368
1369
1370
1371
1372
1373
1374
1375
1376
1377
1378
1379
1380
1381
1382
1383
1384
1385
1386
1387
1388
1389
1390
1391
1392
COMMON/PSIZE/PSZCUTtPSL,PSNL,PSNLCT,PSDIE.PSBRK,PSTIRE
EFORG=0.0
IF(IFT.GT.2)GOTO 99
IF(IVHCLS.GT.3)GOTO 20
IF(IFT.EQ.2)GOTO 10
LEADED LDV, LDT1, LDT2
CALL GETEFOdYRPTR,IVHCLS,1,1,EF)
EFORG=EF*PSL
GOTO 99
UNLEADED LDV, LOT 1 , LDT2
10 IF(IYRPTR.LT.1975)GOTO 99
CALL NLCONS(IYEAR,IVRPTR.IVHCLS,TAMFRC,RMIS)
FCAT=FREMCT(IYRPTR.IVHCLS.2,6)
FNOCAT=FREMCT(IVRPTR,IVHCLS,2.1)
CALL GETEFOdYRPTR,IVHCLS,6,2.EFNLCT)
CALL GETEFO dVRPTR. IVHCLS, 6,1.EFLCAT)
CALL GETEFO(IYRPTR,IVHCLS,1,2,EFNL)
EFORG=(1.-RMIS)*FCAT*EFNLCT*PSNLCT+RMIS*FCAT*EFLCAT*PSL+
&FNOCAT*EFNL*PSNL
GOTO 99
HDGV
20 IF(IVHCLS.GT.5)GOTO 99
CALL GETCFE(IYRPTR,IVHCLS,FLECON)
IF(IVHCLS.EQ.4.AND.IYRPTR.GE.1987.AND.I FT.EQ.2)GOTO 30
IF(IVHCLS.EQ.4.AND.IYRPTR.GE.1987.AND.IFT.NE.2)GOTO 99
IF(IVHCLS.EQ.5.AND.IYRPTR.LT.1987)GOTO 99
IF(IVHCLS.EQ.4.AND.IYRPTR.GE.1987.AND.I FT.NE.2)GOTO 99
ALL LEADED HDGV
CALL GETEFO(IYRPTR,IVHCLS,1,1.EFL)
EFORG=EFL*PSL*5.0/FLECON
GOTO 99
UNLEADED HDGV 1987+
30 IFdYRPTR.LT. 1987. OR. IVHCLS. EQ.5)GOTO 99
CALL NLCONSdVEAR, I YRPTR , IVHCLS, TAMFRC , RMI S)
CALL GETEFO(IYRPTR,IVHCLS,8,2,EFNLCT)
CALL GETEFO(IYRPTR,IVHCLS,8,1.EFLCAT)
EFORG=((1.-RMIS)*EFNLCT*PSNLCT+RMIS*EFLCAT*PSL)*5.0/FLECON
99 RETURN
END
SUBROUTINE GETEFO(IYRPTR.IVHCLS,I CAT,I FT,EF)
GET ORGANIC EMISSION FACTOR
COMMON/EMORG/IORGVR(20),IORGFT(20),IORGCT(20),IORGVH(20),
&EFO(20).NORTOT
EF=0.0
IVH=IVHCLS
IF(IFT.GT.2)GOTO 10
IF(IVHCLS.EQ.5)IVH=4
IF(IVHCLS.LE.3)IVH=1
DO 10 1=1.NORTOT
IF(IVH.NE.lORGVH(I))GOTO 10
IF(ICAT.NE.IORGCT(I))GOTO 10
IF(IFT.NE.lORGFT(I))GOTO 10
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page
25
1393
1394
1395
1396
1397
1398
1399
1400
1401
1402
1403
1404
1405
1406
1407
1408
1409
1410
141 1
1412
1413
1414
1415
1416
1417
1418
1419
1420
1421
1422
1423
1424
1425
1426
1427
1428
1429
1430
1431
1432
1433
1434
1435
1436
1437
1438
1439
1440
1441
1442
1443
1444
1445
1446
1447
1448
1449
1450
C
C
C '
C*
C '
C '
C
C
C 1
C
C 1
C 1
C
C
C
C 1
IF(IYRPTR.GE.IORGVR(I))EF=EFO(I)
10 CONTINUE
RETURN
END
**********************************
*
* SULFATE EMISSION FACTORS SUBROUTINES *
* *
*********************************************************
SUBROUTINE SULFAT(I YEAR,IYRPTR,IVHCLS,VEHSPD,I FT,EF)
FIND SULFATE EMISSIONS BY YEAR, VEHICLE.CLASS, AND FUEL TYPE
COMMON/EMCT/FPBEMC(15.6),FNLEMC(15,6,8),FDIEMC(15,6),IEMCYR(15),
8.IEMCT,TAMPER(3,2)
COMMON/PSIZE/PSZCUT,PSL.PSNL,PSNLCT,PSDIE,PSBRK,PSTIRE
EF=0.0
IF(IFT.GT.2)GOTO 99
IF(IVHCLS.GT.3)GOTO 20
IF(IFT.EQ.2)GOTO 10
LEADED NOCAT LDV, LDT1, LDT2
CALL GETEFSUYRPTR, I VHCLS, 1 , 1 , VEHSPD, EF)
EF=EF*PSL
GOTO 99
UNLEADED 1975+ LDV, LDT1, LDT2
10 IFUYRPTR.LT. 1975)GOTO 99
CALL GETEFSUYRPTR, I VHCLS, 7, 2, 19.6.EFNOAR)
CALL GETEFSUYRPTR, I VHCLS, 8, 2, 19.6.EFAIR)
CALL GETEFSUYRPTR, I VHCLS, 1,2, 19.6.EFNOCT)
CALL GETEFSUYRPTR. I VHCLS. 6, 1 , 19.6.EFMIS)
CALL NLCONSUVEAR, I YRPTR, IVHCLS , TAMFRC , RMIS)
FNOAR=FREMCT(IYRPTR,IVHCLS.IFT,7)
FAIR = FREMCT(IYRPTR,IVHCLS.I FT,8)
FNOCAT = FREMCT(IYRPTR,IVHCLS,I FT, 1)
EFLO=(1.-RMIS)*(FNOAR*EFNOAR*PSNLCT+FAIR*EFAIR*PSNLCT+
&FNOCAT*EFNOCT*PSNL)+RMIS*EFMIS*PSL
CALL GETEFSUYRPTR, IVHCLS. 1 , 2 . 34 . 8 . EFNOCT)
CALL GETEFSUYRPTR , IVHCLS . 2 , 2 , 34 .8 , EFOX)
CALL GETEFSCIVRPTR,IVHCLS,4,2,34.8,EF3W)
CALL GETEFSUYRPTR , IVHCLS , 3 , 2 , 34 . 8 , EFOXAR)
CALL GETEFSUVRPTR, I VHCLS , 5 , 2 , 34 . 8 , EF3WOX)
CALL GETEFSUYRPTR, I VHCLS, 6, 1 .34.8.EFMIS)
FOX=FREMCT(IYRPTR,IVHCLS.IFT,2)
FOXAR=FREMCT(IYRPTR,IVHCLS,IFT,3)
F3W=FREMCT(IYRPTR,IVHCLS,IFT,4)
F3WOX=FREMCT(IYRPTR,IVHCLS,I FT,5)
EFHI=(1.-RMIS)*(FNOCAT»EFNOCT*PSNL+FOX*EFOX*PSNLCT+
&F3W*EF3W*PSNLCT+FOXAR*EFOXAR*PSNLCT+F3WOX*EF3WOX*PSNLCT)+
&RMIS*EFMIS*PSL
SLOPE=(EFHI-EFLO)/(34.8-19.6)
EF=SLOPE*(VEHSPD-19.6J+EFLO
GOTO 99
HDGV
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page 26
1451
1452
1453
1454
1455
1456
1457
1458
1459
1460
1461
1462
1463
1464
1465
1466
1467
1468
1469
1470
1471
1472
1473
1474
1475
1476
1477
1478
1479
1480
1481
1482
1483
1484
1485
1486
1487
1488
1489
1490
1491
1492
1493
1494
1495
1496
1497
1498
1499
1500
1501
1502
1503
1504
1505
1506
1507
1508
20 IF(IVHCLS.GT.5)GOTO 40
CALL GETCFEdVRPTR,IVHCLS,FLECON)
CALL GETEFS(IYRPTR. IVHCLS, 1 , 1 , VEHSPD , EFNOCT)
IF(IFT.EQ.2)GOTO 30
EF=EFNOCT*PSL*5.0/FLECON
GOTO 99
UNLEADED HDGV 1987+
30 IF(IYRPTR.LT.1987.OR.IVHCLS.EQ.5)GOTO 99
CALL NLCONSd YEAR,IYRPTR.IVHCLS.TAMFRC,RMIS)
CALL GETEFSCIYRPTR,IVHCLS,8,2,VEHSPD,EFCAT)
EF=((1.-RMIS)*EFCAT*PSNLCT+RMIS*EFNOCT*PSL)*5.0/FLECON
GOTO 99
MOTORCYLES
40 EF=0.0
99 RETURN
END
SUBROUTINE GETEFSCIYRPTR,IVHCLS,ICAT.IFT,VEHSPD,EF)
GETS SULFATE EMISSION FACTOR
COMMON/EMSUL/ISULYR(20),ISULFT(20),ISULCT(20),ISULVH(20),
8.EFSC20.2) .NSUTOT
EF=0.0
EF1=0.0
EF2=0.0
IVH=IVHCLS
IF(IFT.GT.2)GOTO 99
IF(IVHCLS.GE.6)GOTO 99
IF(IVHCLS.LE.3)IVH=1
IF (IVHCLS.EQ.5)IVH=4
DO 10 1 = 1, NSUTOT
IF(IVH.NE.ISULVH(I))GOTO 10
IF(IFT.NE.ISULFT(I))GOT010
IFdCAT.NE. ISULCTd))GOTO 10
IF(IVH.EQ.4.AND.IYRPTR.GE.ISULYR(I))EF=EFS(I,1)
IF(EF.NE.O.O)GOTO 99
IF(IYRPTR.GE.ISULYR(I).AND.EFS(I.2).LT.20.)EF1=EFS(I.1)
IF(I YRPTR. GE.ISULYRd).AND.EPS(I,2).GT.20.)EF2=EFS(I,1)
10 CONTINUE
20 IF(EF1.EQ.O.O.AND.EF2.EQ.O.O)GOTO 99
SLOPE=(EF2-EF1)/(34.8-19.6)
EF=EF1+SLOPE*(VEHSPD-19.6)
99 RETURN
END
SUBROUTINE DIESEL(IYRPTR,IVHCLS,EF)
COMMON/PSIZE/PSZCUT,PSL,PSNL,PSNLCT,PSDIE,PSBRK,PSTIRE
COMMON/EMDIE/IDIEYR(20),IDIEVH(20),EFDP(20).NDITOT
COMMON/HDDV/IHDDMY(13),HDDVCF(13)
EF=0.0
IVEH=IVHCLS
IF(IVHCLS.EQ.3)IVEH=2
IF(IVHCLS.GE.5)GOTO 99
DO 10 1=1,NDITOT
IF(IVEH.NE.IDIEVH(I))GOTO 10
IF(IYRPTR.GE.IDIEYR(I))EF=EFDP(I)
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page 27
1509
1510
1511
1512
1513
1514
1515
1516
1517
1518
1519
1520
1521
1522
1523
1524
1525
1526
1527
1528
1529
1530
1531
1532
1533
1534
1535
1536
1537
1538
1539
1540
1541
1542
1543
1544
1545
1546
1547
1548
1549
1550
1551
1552
1553
1554
1555
1556
1557
1558
1559
1560
1561
1562
1563
1564
1565
1566
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
10 CONTINUE
IF(IVHCLS.EQ.4)GOTO 20
EF=EF*PSDIE
GOTO 99
HDDT
20 DO 30 1=1.13
30 IF(IYRPTR.GE.IHDDMV(I))CF=HDDVCF(I)
EF=EF*PSDIE*CF
99 RETURN
END
****
* FUNCTIONS
*********
*
*********
' FUNCTION FREMCT(IVRPTR,IVHCLS,IFT,ICAT)
COMMON/EMCT/FPBEMC(15,6),FNLEMC(15,6.8),FDIEMC(15,6),IEMCVR(15)
8.IEMCT,TAMPER(3.2)
IFdVHCLS.LT.1.OR.IVHCLS.GT.6)CALL QUITER(0..1.32.INERR)
IF(IFT.LT.1.OR.I FT.GT.3)CALL QUITER(0. , 1 ,32,INERR)
IF(ICAT.LT.1.OR.ICAT.GT.8)CALL QUITER(0.,1,32,INERR)
FREMCT=0.0
NPTR=1
DO 10 1=1,IEMCT
10 IF(IVRPTR.GE.IEMCVR(I))NPTR=I
IF(IFT.EQ.1)FREMCT=FPBEMC(NPTR,IVHCLS)
IF(IFT.EQ.2)FREMCT=FNLEMC(NPTR,IVHCLS,ICAT)
IF(IFT.EQ.3)FREMCT=FDIEMC(NPTR,IVHCLS)
RETURN
END
FUNCTION FRFLT(IVRPTR,IVHCLS,I FT)
COMMON/FLEET/FSF(26,6,3),IFSFYR(26,6),IFSF(6),VMT(6)
IF(IVHCLS.LT.1.OR.IVHCLS.GT.6)CALL QUITER(0.,2,32,INERR)
IF(IFT.LT. 1.OR.I FT.GT.3)CALL QUITER(0.,2,32,INERR)
FRFLT=0.0
NPTR=1
NFSF=IFSF(IVHCLS)
DO 10 1=1.NFSF
10 IF(IYRPTR.GE.IFSFYR(I,IVHCLS))NPTR=I
FRFLT=FSF(NPTR,IVHCLS,IFT)
RETURN
END
* CALCULATIONS
***<
***********
*
**********************
SUBROUTINE CALC
COMMON/TWFRAC/TWFAC(20,6 , 3)
COMMON/PSIZE/PSZCUT,PSL,PSNL.PSNLCT,PSDIE,PSBRK,PSTIRE
COMMON/FLEET/FSF(26,6,3),IFSFYR(26,6),IFSF(6),VMT(6)
COMMON/ANSWER/PBSUM(6),ORGSUM(6),SULSUM(6),DIESUM(6),CLSSUM(6),
8.BRKSUM.TIRSUM.SUMALL
COMMON/SCENR/IVEAR.VEHSPD.IDRIV
-------
Listing of PART.LEAD at 13:40:01 on DEC 2. 1985 for CCid=SHJD
Page
28
1567 C
1568 C START VEHICLE CLASS LOOP
1569 DO 10 IVHCLS=1,6
1570 10 CLSSUM(IVHCLS)=0.0
1571 DO 30 IVHCLS=1,6
1572 PBSUM(IVHCLS)=0.0
1573 ORGSUM(IVHCLS)=0.0
1574 SULSUM(IVHCLS)=0.0
1575 DIESUM(IVHCLS)=0.0
1576 C START 20-YEAR LOOP
1577 IFIRVR=IYEAR-19
1578 DO 20 IYRPTR=IFIRYR,IYEAR
1579 IAGE=IVEAR-IYRPTR+1
1580 IF(IAGE.GT.20)IAGE=20
1581 C FUEL TYPE LOOP (LEADED, UNLEADED)
1582 DO 15 IFT=1,2
1583 FR=FRFLT(IVRPTR,IVHCLS,IFT)
1584 TW=TWFAC(IAGE,IVHCLS,I FT)
1585 IFCIVRPTR.LT.1987.AND.IVHCLS.EQ.5)GOTO 15
1586 IF(IVRPTR.GE.1987.AND.IVHCLS.EQ.5.AND.IFT.EQ.2)GOTO 15
1587 IF(IYRPTR.GE.1987.AND.IVHCLS.EQ.4.AND.IFT.EQ.1)GOTO 15
1588 IF(IYRPTR.GE.1987.AND.IVHCLS.EQ.5.AND.IFT.EQ.1)FR = FRFLT(IYRPTR,4 , 1 )
1589 CALL LEAD(IYEAR,IYRPTR,IVHCLS,IFT,VEHSPD,IDRIV.EFPB)
1590 PBSUMdVHCLS)=PBSUM(IVHCLS)+EFPB*FR*TW
1591 CALL ORGANCCIYEAR,IYRPTR,IVHCLS,IFT,EFORG)
1592 ORGSUM(IVHCLS)=ORGSUM(IVHCLS)+EFORG*FR*TW
1593 CALL SULFATUYEAR,IVRPTR,IVHCLS,VEHSPD,IFT.EFSUL)
1594 SULSUM(IVHCLS)=SULSUM(IVHCLS)+EFSUL*FR*TW
1595 FRTW=FR*TW
1596 15 CONTINUE
1597 IFT=3
1598 CALL DIESELUYRPTR, IVHCLS, EFDIE)
1599 FR=FRFLT(IVRPTR,IVHCLS,IFT)
1600 DIESUM(IVHCLS)=DIESUM(IVHCLS)+EFDIE*FRFLT(IYRPTR,IVHCLS,IFT)*
1601 CTWFACdAGE, IVHCLS. IFT)
1602 20 CONTINUE
1603 C WEIGHT BY AREA TRAVEL FRACTION
1604 IF(IVHCLS.LE.3)CLSSUM(IVHCLS)=(PBSUM(IVHCLS)+ORGSUM(IVHCLS)*
1605 &SULSUM(IVHCLS)+OIESUM(IVHCLS))*VMT(IVHCLS)
1606 IF (IVHCLS. EQ.4.0R. I VHCLS . EQ. 5) CLSSUM(4) =CLSSUM(4)-i-(PBSUM( IVHCLS)
1607 &+ORGSUM(IVHCLS)+SULSUM(IVHCLS))*VMT(4)
1608 IF(IVHCLS.EQ.4)CLSSUM(5)=DIESUM(4)*VMT(5)
1609 IF(IVHCLS.EQ.6)CLSSUM(6)=PBSUM(6)*VMT(6)
1610 30 CONTINUE
1611 C SHIFT HDGV2 TO HDGV
1612 PBSUM(4)=PBSUM(4)+PBSUM(5)
1613 ORGSUM(4)=ORGSUM(4)+ORGSUM(5)
1614 SULSUM(4)=SULSUM(4)+SULSUM(5)
1615 PBSUM(5)=0.0
1616 ORGSUM(5)=0.0
1617 SULSUM(5)=0.0
1618 C SHIFT HDDT TO CLASS FIVE
1619 DIESUM(5)=DIESUM(4)
1620 DIESUM(4)=0.0
1621 C
1622 C FIND TOTAL EMISSIONS
1623 C
1624 SUMALL=0.0
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page 29
1625
1626
1627
1628
1629
1630
1631
1632
1633
1634
1635
1636
1637
1638
1639
1640
1641
1642
1643
1644
1645
1646
1647
1648
1649
1650
1651
1652
1653
1654
1655
1656
1657
1658
1659
1660
1661
1662
1663
1664
1665
1666
1667
1668
1669
1670
1671
1672
1673
1674
1675
1676
1677
1678
1679
1680
1681
1682
DO 40 1=1.6
SUMALL=SUMALL+CLSSUMCI)
40 CONTINUE
BRKSUM=PSBRK*0.0128
TIRSUM=PSTIRE*0.002
SUMALL=SUMALL+BRKSUM+TIRSUM
RETURN
END
C
C
C *
C
C
OUTPUT ROUTINES
>»*
*
*************
SUBROUTINE OUT
INTEGER TAMFLG,SPCFLG,VMFLAG
COMMON/PROJCM/PROJIDC20)
COMMON/IOUCOM/IOUREP,IOUERR,IOUASK,IOUDAT.IOUGEN
COMMON/FLAGS 1/I PROMT,TAMFLG.SPCFLG,VMFLAG,MYMRFG,IMFLAG,MISFLG
COMMON/PSIZE/PSZCUT,PSL,PSNL.PSNLCT,PSDIE,PSBRK,PSTIRE
COMMON/ANSWER/PBSUMC6),ORGSUM(6),SULSUM(6),DIESUM(6),CLSSUM(6),
&BRKSUM.TIRSUM,SUMALL
COMMON/FLEET/FSF(26,6,3),IFSFVR(26,6),IFSF(6),VMT(6)
COMMON/SCENR/IVEAR.VEHSPD.IDRIV
DATA NAMTAM/'NORM'
NAMIM/' NO '.
NMYMR/'NORM',
NAMVM/'NORM',
NAMSPC/'NORM'
NAMORW 'TRAN'
NAMMIS/'NORM'
CLSNM/'LDV ',
&
&
&
&
&
&
&
DIMENSION NAMTAM(2),NAMIM(2),NMYMR(4),NAMVM(3),NAMSPC(3) ,
&NAMDRV(2),CLSNM(6).NAMMISC3)
,'USER'/ ,
' YES'/,
'MILE',' REG',' M+R'/,
'SCEN','USER'/,
,'USER'/,
,'CRSE'/,
,'USER','SMPL'/,
'LDT1','LDT2','HDGV,'HDDT',' MC '/
CALL OUTHDR(IOUREP)
WRITEUOUREP, 1030)IYEAR,VEHSPD,PSZCUT
1030 FORMAT('OCALCULATIONS FOR YEAR ',14,' SPEED '.F5.1.' MPH',
8.' PARTICLE CUTOFF = ' , F5 . 2 , ' MICRONS')
WRITE(IOUREP,lOOO)NAMTAM(TAMFLG).NAMMIS(MISFLG),NAMVM(VMFLAG)
1000 FORMAT('OFLAG SETTINGS: '/' TAMPERING RATES = '.A4,
&' MISFUELING RATES = ',A4,' VMT MIXTURE = ',A4)
WRITE(IOUREP.lOIO)NMVMR(MYMRFG),NAMIM(IMFLAG),NAMSPC(SPCFLG)
1010 FORMATC MILEAGE & REGIST. = ',A4,' I/M PROGRAM = ',
&A4,' SPECIAL CONDITIONS = ',A4)
WRITE(IOUREP,1020)NAMDRVCIDRIV)
1020 FORMATC DRIVING CONDITIONS = ',A4)
C PRINT BREAKDOWN
WRITE(IOUREP,1025)
1025 FORMATC'0',9X,'PARTICULATE EMISSIONS (GRAMS/MILE) BY VEHICLE CLASS')
WRITE(IOUREP.1040)
1040 FORMATC'0','VEH.CLASS',2X,' LEAD ',2X,'ORGANIC',2X,'SULFATE',2X,
&' DIESEL ' ,2X, 'TOTAL' ,4X,'VMT FRAC. ' ,2X.'VMT WT. TOTAL')
DO 10 IVHCLS=1,6
CLS=0.0
IFCVMTCIVHCLS).NE.O.O)CLS=CLSSUM(IVHCLS)/VMT(IVHCLS)
10 WRITECIOUREP.1050)CLSNM(IVHCLS),PBSUMCIVHCLS),ORGSUM(IVHCLS).
&SULSUM(IVHCLS),DIESUM(IVHCLS),CLS.VMTCIVHCLS),CLSSUMCIVHCLS)
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page 30
1683
1684
1685
1686
1687
1688
1689
1690
1691
1692
1693
1694
1695
1696
1697
1698
1699
1700
1701
1702
1703
1704
1705
1706
1707
1708
1709
1710
171 1
1712
1713
1714
1715
1716
1717
1718
1719
1720
1721
1722
1723
1724
1725
1726
1727
1728
1729
1730
1731
1732
1733
1734
1735
1736
1737
1738
1739
1740
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
1 050 FORMATC3X , A4 , 3X . 2X , 4C F7 . 4 , 2X) , F7 . 4 , 4X , F7 . 4 , 4X, F7 . 4)
WRITECIOUREP.1052)BRKSUM,TIRSUM
1052 FORMATC TOTAL BRAKE PARTICIPATES = '.F6.4/
*' TOTAL TIRE PARTICULATES = '.F6.4)
WRITECIOUREP,1055)SUMALL
1055 FORMATC TOTAL VMT WEIGHTED PARTICIPATE EMISSIONS = '.F7.4)
OUTPUT FRACTIONS
WRITEUOUREP, 1060)
1060 FORMATC'0',18X.'PARTICULATE FRACTION BY VEHICLE CLASS')
WRITEUOUREP. 1070)
1070 FORMATC1X.'VEH.CLASS',2X,' LEAD ',2X,'ORGANIC',2X,'SULFATE'
&2X,'DIESEL ')
DO 20 IVHCLS=1,6
PB=0.0
ORG=0.0
SUL=0.0
DIE=0.0
IF(CLSSUMCIVHCLS).EQ.O.OJGOTO 20
PB=PBSUM(IVHCLS)*VMT(IVHCLS)/CLSSUM(IVHCLS)
ORG=ORGSUM(IVHCLS)*VMT(IVHCLS)/CLSSUM(IVHCLS)
SUL=SULSUM(IVHCLS)*VMT(IVHCLS)/CLSSUM(IVHCLS)
DIE=DIESUM(IVHCLS)*VMT(IVHCLS)/CLSSUMCIVHCLS)
20 WRITECIOUREP,1050)CLSNM(IVHCLS),PB,ORG.SUL.DIE
RETURN
END
SUBROUTINE OUTHDR(IOUOUT)
OUTHDR writes the run title on unit IOUOUT.
CALLED BV CONSEC. OUT
Input on cal1 :
parameter list: IOUOUT
common block: /PROJCM/ PROJID
Output on return:
common block: /SVSCOM/ INITPR
COMMON /PROJCM/ PROJID(20)
WRITE(IOUOUT,200) PROJID
200 FORMATC'1',20A4)
RETURN
END
BLOCK DATA
COMMON/EMORG/IORGYRC20),IORGFT(20),IORGCTC20),IORGVHC20),
&EFOC20),NORTOT
COMMON/EMSUL/ISULYRC20),ISULFT(20).ISULCTC20), ISULVHC20),
&EFSC20.2).NSUTOT
COMMON/EMDIE/IDIEYRC20).IDIEVHC20).EFDPC20).NDITOT
COMMON/EMPB/EFLEADC20).IPBVRC20),IPBVH(20).IPBFTC20).NPBTOT
COMMON/GAS/LCGSIZ.LCGYRC20),PBCGC20.2).NCFEYRC27),CFE(27,6),
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page 31
1741
1742
1743
1744
1745
1746
1747
1748
1749
1750
1751
1752
1753
1754
1755
1756
1757
1758
1759
1760
1761
1762
1763
1764
1765
1766
1767
1768
1769
1770
1771
1772
1773
1774
1775
1776
1777
1778
1779
1780
1781
1782
1783
1784
1785
1786
1787
1788
1789
1790
1791
1792
1793
1794
1795
1796
1797
1798
C
C
C
C
C
C
C
C
C
&ICFESZ.IFECSZ,FECVS(13),FEC(13,2)
COMMON/MSFL/RMISFL(20,4.2),TMISFL(6.2)
COMMON/PSTEMP/PSDIST(5,2,6)
COMMON/TWF/TWFREG(20,6.3).!TWFMA(20.6,3)
COMMON/HDDT/IHTRG,IHDDMA(20,4),IHDDYR(21),HDDTRG(21,4)
COMMON/HDDV/IHDDMV(13),HDDVCF(13)
COMMON/EMCT/FPBEMC(15.6),FNLEMC(15,6,8).FDIEMC(15.6),IEMCYR(15),
&IEMCT,TAMPER(3.2)
COMMON/FLEET/FSF(26,6.3).IFSFYR(26.6),IFSF(6),VMT(6)
DATA IORGVR/20*0/,
& IORGFT/20*0/,
& IORGCT/20*0/,
8. IORGVH/20*0/,
& EFO/20*0.0/.
& NORTOT/0/
DATA
&
&
&
&
&
&
DATA
&
&
&
DATA
&
&
&
&
DATA
&
&
&
&
&
&
&
ISULVR/20*0/.
ISULFT/20*0/,
ISULCT/20*0/,
ISULFT/20*0/,
ISULVH/20*0/.
EFS/40*0.0/,
NSUTOT/0/
IDIEYR/20*0/,
IDIEVH/20*0/,
EFDP/20*0.0/,
NDITOT/0/
EFLEAD/20*0.0/.
IPBVR/20*0/ ,
IPBVH/20»0/,
IPBFT/20*0/,
NPBTOT/0/
LCGSIZ/0/,
LCGYR/20*0/,
PBCG/40*0.0/,
NCFEYR/27*0/,
CFE/162»0.0/,
ICFESZ/0/.
FECVS/13*0.0/.
FEC/13*0.0/
DATA RMISFL/160*0.0/,
TMISFL/12*0.0/
PSDIST/60*0.0/
DATA
DATA
TWFREG/360*0.0/,
ITWFMA/360*0/
DATA IHTRG/0/.
& IHDDMA/80*0/,
& IHDDYR/21*0/,
8. HDDTRG/84*0.0/
DATA IHDDMY/13*0/,
-------
Listing of PART.LEAD at 13:40:01 on DEC 2, 1985 for CCid=SHJD
Page 32
1799
1800
1801
1802
1803
1804
1805
1806
1807
1808
1809
1810
181 1
1812
1813
& HDDVCF/13*0.0/
DATA FPBEMC/90*0.0/,
& FNLEMC/720*0.0/,
& FDIEMC/90*0.0/,
& IEMCYR/15/,
& IEMCT/0/.
& TAMPER/6*0.0/
DATA FSF/468*0.0/,
& IFSFVR/156*0/f
& IFSF/6*0/,
8. VMT/6*0.0/
END
-------
APPENDIX B
PROGRAM DATA FILE
-------
Listing of PART.DATA at 13:39:45 on DEC 2, 1985 for CCid=SHJD Page
1
2
3
4
5
6
7
B
9
10
1 1
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
6
1900
1970
1975
1975
1975
1975
13
1900
1975
1974
1975
1975
1975
1975
1975
1975
1975
1975
1975
1975
3
1900
1981
1987
3
1900
1981
1987
4
1900
1987
1987
1987
3
1900
1987
1987
1
1900
6
1900
1900
1971
1971
1975
1975
6
1900
1900
1971
1971
1979
1979
2
1900
1978
17
0
0
0
0
0
0
S
S
S
S
S
S
S
S
S
S
S
S
S
0.
0.
0.
0.
0.
0.
0
0
0
0
S
S
S
0.
1
2
1
2
1
2
1
2
1
2
1
2
1
1
1 1 0.
1 1 0.
1 1 0.
620.
610.
1 2 0.
1119
6 1 19
1 1 34
7 2 19
8 2 19
1 2 19
1 1 34
1 2 34
2 2 34
4 2 34
3 2 34
5 2 34
6 1 34
700
300
200
BOO
300
260
1 1 0.
820.
8 1 0.
1 1 0.
1 1 0.
820.
1 1 0.
700
0.887
0.113
0.916
0.084
0.724
0.276
0.887
0.113
0.916
0.084
0.724
0.276
0.198
0.046
193
068
030
017
068
030
.6
.6
.8
.6
.6
.6
.8
.8
.8
.8
.8
.8
.8
370
054
163
370
006
048
006
0
0
0
0
0
0
0
0
0
0
0
0
0
.002
.002
.001
.005
.016
.002
.001
.001
.005
.001
.020
.025
.001
-------
Listing of PART.DATA at 13:39:45 on DEC 2, 1985 for CCid=SHJD
Page
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
1 10
1 1 1
1 12
1 13
1 14
1 15
1 16
1900 1 .79 0.014
1975 1 .82 0.014
1976 2.02 0.014
1977 2.03 0.014
1978 1 .94 0.014
1979 1 .85 0.014
1980 1 .38 0.014
1981 1.15 0.014
1982 .24 0.014
1983 .14 0.014
1984 .10 0.014
1985 . 10 0.014
1986 .10 0.014
1987 . 10 0.014
1988 1.10 0.014
1989 1.10 0.014
1990 1 . 10 0.014
15
1 .000 0. 128 0. 134 0. 158
0.000 0.081 0.020 0.000
0.000 0.597 0.637 0.650
0.000 0.322 0.343 0.350
0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000
0.000 0.003 0.003 0.004
0.000 0.919 0.980 1.000
0.000 0.597 0.637 0.650
0.000 0.322 0.343 0.350
22
1900 0.000 1.000 0.000
1975 0.869 0.128 0.003
1976 0.863 0. 134 0.003
1977 0.838 0.158 0.004
1978 0.865 0. 126 0.009
1979 0.875 0.097 0.028
1980 0.966 0.000 0.034
1981 0.939 0.000 0.061
1982 0.954 0.000 0.046
1983 0.947 0.000 0.053
1984 0.940 0.000 0.060
1985 0.934 0.000 0.066
1986 0.927 0.000 0.073
1987 0.920 0.000 0.080
1988 0.910 0.000 0.090
1989 0.900 0.000 0.100
1990 0.887 0.000 0. 1 13
1991 0.887 0.000 0.113
1992 0.886 0.000 0. 1 14
1993 0.886 0.000 0. 1 14
1994 0.885 0.000 0.115
1995 0.885 0.000 0. 1 15
01 0.028 12818 0358.9 0
02 0. 107 12639 1352.4 0
03 0. 100 1 1933 1 193.3 0
04 0.094 11268 1059.2 0
05 0.088 10639 0936.2 0
06 0.080 10045 0803.6 0
07 0.075 09485 071 1 .4 0
0. 126 0.097
0.000 0.000
0.650 0.640
0.350 0.345
0.000 0.015
0.000 0.000
0.009 0.028
1.000 1.000
0.650 0.655
0.350 0.345
.038
. 142
. 125
.111
.098
.084
.075
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.421 0.033 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.514 0.099 0.142 0.109 0.000 0 .-000 0.000 0.000 0.000
0.053 0.263 0.313 0.244 0.396 0.532 0.587 0.641 0.704
0.012 0.605 0.545 0.647 0.604 0.468 0.414 0.359 0.296
0.034 0.061 0.046 0.053 0.060 0.066 0.073 0.080 0.090
1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
0.474 0.296 0.149 0.194 0.237 0.381 0.450 0.518 0.554
0.526 0.704 0.851 0.806 0.763 0.691 0.587 0.482 0.446
-------
Listing of PART.DATA at 13:39:45 on DEC 2, 1985 for CCid=SHJD
Page
1 17
1 18
1 19
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
08 0.
09 0.
10 0.
11 0.
12 0.
13 0.
14 0.
15 0.
16 0.
17 0.
18 0.
19 0.
20 0.
27
1900
1970
1971
1972
1973
1974
1975
1976
1977
1978
1979
1980
1981
1982
1983
1984
1985
1986
1987
1988
1989
1990
1991
1992
1993
1994
1995
13
05.0
10.0
15.0
20.0
25.0
30.0
32.7
35.0
40.0
45.0
50.0
55.0
60.0
15
1 .000
0.000
069 08955 0617
062 08455 0524
056 07983 0447
050 07538 0376
043 071 17 0306
037 06720 0248
031 06345 0196
024 05991 0143
018 05657 0101
012 05341 0064
008 04043 0032
006 04762 0028
008 04496 0036
13
13
13
13
12
12
13
14
15
16
17
20
21
22
22
22
23
23
24
24
25
25
26
26
27
27
29
0.
0.
0.
0.
0.
0.
t
t
0
0
.9
.9
.2
. 1
.9
.6
.5
.8
.5
.8
.2
.0
.4
.2
.2
.8
.2
.8
.3
.8
.2
.7
.2
.6
.2
.6
.0
323
553
692
790
885
963
000
022
053
073
078
063
023
10
10
10
10
09
09
1 1
12
13
13
14
16
17
18
19
19
20
21
23
23
23
24
24
24
25
25
26
0
0
0
1
1
1
1
1
1
1
1
1
1
. 188
. 123
.6 07
.6 07
.4 07
.2 07
.9 07
.6 06
.6 08
.3 09
.0 09
.4 09
.2 09
. 1 11
.7 13
.6 13
.2 13
.9 13
.7 14
.4 14
.0 14
.3 14
. 1 14
.0 15
.5 15
.4 15
.3 15
.8 16
.2 16
.467
.709
.997
. 153
.248
.294
.303
.303
.288
.256
.210
. 159
. 104
0.088
0.225
.9
.9
.7
.4
.0
.9
.8
.7
.4
.6
.8
.5
.3
.6
.7
.9
.0
.3
.5
.7
.9
.2
.4
.7
.9
.2
.4
0
0
.9 0
.2 0
.0 0
.9 0
.0 0
.6 0
.7 0
.8 0
.8 0
. 1 0
.3 0
.6 0
.0 0
00.
00.
00.
00.
00.
00.
00.
00.
00.
00.
00.
00.
00.
00.
00.
00.
00.
00.
09.
09.
09.
09.
09.
09.
09.
09.
10.
.038
.083
.065
.055
.047
.040
.032
.026
.021
.015
.01 1
.007
.003
.003
.004
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
5
5
6
7
7
8
8
9
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
5
5
5
5
5
5
5
5
5
0.
0.
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.0
.6
.6
.6
.6
.7
.7
.7
.7
.8
027
069
6
6
6
6
6
6
6
7
7
8
8
8
8
8
8
8
9
9
9
9
9
9
9
9
9
9
9
.5
.4
.4
.4
.5
.7
.8
.3
.7
.0
.2
.4
.6
.8
.9
.9
.0
.0
.0
. 1
.2
.2
.3
.4
.4
.5
.6
0.030 0.021 0.026 0.021 0.022 0.000 0.000 0.000 0.000 0.000
0.034 0.027 0.011 0.000 0.000 0.000 0.000 0.000 0.000 0.000
-------
Listing of PART.DATA at 13:39:45 on DEC 2, 1985 for CCid=SHJD Page
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
21 1
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
0.000 0.570 0.504 0.596
0.000 0.307 0.271 0.321
0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000
0.000 0.002 0.003 0.005
0.000 0.877 0.775 0.917
0.000 0.570 0.504 0.596
0.000 0.307 0.271 0.321
22
1900 0.000 1.000 0.000
1975 0.810 0. 188 0.002
1976 0.909 0.088 0.003
1977 0.957 0.038 0.005
1978 0.964 0.027 0.009
1979 0.942 0.030 0.028
1980 0.945 0.021 0.034
1981 0.914 0.026 0.060
1982 0.899 0.021 0.080
1983 0.878 0.022 0. 100
1984 0.870 0.000 0. 130
1985 0.840 0.000 0. 160
1986 0.820 0.000 0. 180
1987 0.790 0.000 0.210
1988 0.760 0.000 0.240
1989 0.730 0.000 0.270
1990 0.706 0.000 0.294
1991 0.697 0.000 0.303
1992 0.688 0.000 0.312
1993 0.679 0.000 0.321
1994 0.670 0.000 0.330
1995 0.661 0.000 0.339
01 0.023 17394 0400.1 0
02 0.089 17079 1520.0 0
03 0.085 15839 1346.3 0
04 0.081 14690 1 189.9 0
05 0.076 13624 1035.4 0
06 0.072 12636 0909.8 0
07 0.068 1 1719 0796.9 0
08 0.064 10868 0695.6 0
09 0.060 10080 0604.8 0
10 0.055 09348 0514. 1 0
11 0.050 08670 0433.5 0
12 0.046 08041 0369.9 0
13 0.042 07457 0313.2 0
14 0.038 06916 0262.8 0
15 0.034 06415 0218. 1 0
16 0.029 05949 0172.5 0
17 0.025 05517 0137.9 0
18 0.021 051 17 0107.5 0
19 0.017 04746 0080.7 0
20 0.025 04402 01 10. 1 0
15
1 .000 1 .000 1 .000 1 .000
0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000
0.605
0.325
0.000
0.000
0.009
0.930
0.605
0.325
.036
. 135
. 120
. 106
.092
.081
.071
.062
.054
.046
.039
.033
.028
.023
.019
.015
.012
.009
.007
.010
0.988
0.000
0.000
0.000
0.000
0.000
0.561 0.564 0.574 0.066 0.092 0.003 0.000 0.000 0.000 0.000
0.405 0.409 0.415 0.887 0.687 0.595 0.550 0.550 0.150 0.150
0.000 0.000 0.000 0.027 0.054 0.126 0.150 0.150 0.350 0.350
0.000 0.000 0.000 0.021 0.167 0.276 0.300 0.300 0.500 0.500
0.028 0.034 0.060 0.080 0.100 0.130 0.160 0.180 0.210 0.240
0.966 0.973 0.989 1.000 1.000 1.000 1.000 1.000 1.000 1.000
0.561 0.564 0.574 0.069 0.073 0.056 0.080 0.080 0.200 0.200
0.405 0.409 0.415 0.931 0.927 0.944 0.920 0.920 0.800 0.800
,
0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.008 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
0.496 0.500 0.500 0.222 0.284 0.000 0.000 0.000 0.000 0.000
0.496 0.500 0.500 0.704 0.577 0.823 0.800 0.800 0.000 0.000
0.000 0.000 0.000 0.000 0.012 0.003 0.000 0.000 0.020 0.020
0.000 0.000 0.000 0.074 0.127 0.174 0.200 0.200 0.980 0.980
-------
Listing of PART.DATA at 13:39:45 on DEC 2, 1985 for CC1d=SHJD Page
0.025 0.050 0.050 0.080 0.113 0.147 0.180 0.194 0.208 0.222
0.992 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
0.496 0.500 0.500 0.222 0.293 0.000 0.000 0.000 0.020 0.020
0.496 0.500 0.500 0.778 0.707 1.000 1.000 1.000 0.980 0.980
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
0.000 0.000 0.000 0.000 0.012
0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000
0.000 0.000 0.000 0.000 0.000
22
1900 0.000 1 .000 0.000
1975 0.000 0.998 0.002
1976 0.000 0.997 0.003
1977 0.000 0.995 0.005
1978 0.000 0.991 0.009
1979 0.972 0.000 0.028
1980 0.966 0.000 0.034
1981 0.940 0.000 0.060
1982 0.920 0.000 0.080
1983 0.900 0.000 0. 100
1984 0.870 0.000 0. 130
1985 0.840 0.000 0 . 1 60
1986 0.820 0.000 0. 180
1987 0.790 0.000 0.210
1988 0.760 0.000 0.240
1989 0.730 0.000 0.270
1990 0.706 0.000 0.294
1991 0.697 0.000 0.303
1992 0.688 0.000 0.312
1993 0.679 0.000 0.321
1994 0.670 0.000 0.330
1995 0.661 0.000 0.339
01 0.023 18352 0422.1 0.036
02 0.089 18001 1602.1 0.138
03 0.085 16622 1412.9 0.122
04 0.081 15348 1243.2 0.107
05 0.076 14172 1077.1 0.093
06 0.072 13087 0942.3 0.081
07 0.068 12084 0821.7 0.071
08 0.064 11158 0714.1 0.062
09 0.060 10303 0618.2 0.053
10 0.055 09514 0523.3 0.045
11 0.050 08785 0439.3 0.038
12 0.046 08112 0373.2 0.032
13 0.042 07491 0314.6 0.027
14 0.038 06917 0262.8 0.023
15 0.034 06386 0217.1 0.019
16 0.029 05897 0171.0 0.015
17 0.025 05446 0136.2,0.012
18 0.021 05028 0105.6 0.009
19 0.017 04643 0078.9 0.007
20 0.025 04287 0107.2 0.009
01 0.023 17552 0403.7 0.035
02 0.087 17230 1499.0 0.129
03 0.083 15964 1325.0 0.114
04 0.079 14791 1168.5 0.101
05 0.075 13705 1027.9 0.088
06 0.071 12699 0901.6 0.078
07 0.066 11766 0776.6 0.067
08 0.062 10901 0675.9 0.058
09 0.058 10101 0585.9 0.050
10 0.054 09359 0505.4 0.043
11 0.049 08671 0424.9 0.037
-------
Listing of PART.DATA at 13:39:45 on DEC 2, 1985 for CCid=SHJD Page
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
31 1
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329 09 0.049 09952 0487.6 0.037
330 10 0.041 09010 0369.4 0.028
331 11 0.036 08156 0293.6 0.023
332 12 0.030 07384 0221.5 0.017
333 13 0.026 06685 0173.8 0.013
334 14 0.022 06052 0133.1 0.010
335 15 0.020 05479 0121.0 0.009.
336 16 0.016 04960 0079.4 0.006
337 17 0.014 04490 0062.9 0.005
338 18 0.012 04065 0048.8 0.004
339 19 0.010 03680 0036.8 0.003
340 20 0.049 03332 0163.3 0.013
341 20
342 01 0.000 00000 00000.0 0.000
343 02 0.166 67910 11273.1 0.241
344 03 0.138 61749 08521.4 0.182
345 04 0.115 56155 06457.8 0.138
346 05 0.097 51073 04954.1 0.106
347 06 0.080 46457 03716.6 0.079
348 07 0.067 42260 02831.4 0.060
12 0.
13 0.
14 0.
15 0.
16 0.
17 0.
18 0.
19 0.
20 0.
20
1900
1977
1978
1979
1980
1981
1982
1983
1984
1985
1986
1987
1988
1989
1990
1991
1992
1993
1994
1995
01 0.
02 0.
03 0.
04 0.
05 0.
06 0.
07 0.
08 0.
045 08035 0361
041 07444 0305
037 06897 0255
033 06391 0210
029 05921 0171
025 05487 0137
020 05084 0101
016 04710 0075
025 04364 0109
0.000 1.000
0.000 1.000
0.000 1.000
0.000 1.000
0.000 1.000
0.000 1.000
0.000 1 .000
0.000 1.000
0.000 1.000
0.000 1.000
0.000 1.000
0.823 0. 177
0.824 0. 176
0.825 0. 175
0.826 0. 174
0.828 0. 172
0.829 0.171
0.833 0. 167
0.837 0. 163
0.840 0. 160
000 00000 0000
148 19967 2955
126 18077 2277
107 16365 1751
092 14815 1363
078 13413 1046
067 12143 0813
058 10993 0637
.6
.2
.2
.9
.7
.2
.7
.4
. 1
.0
. 1
. 7
. 1
.0
.2
.6
.6
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
.031
.026
.022
.018
.015
.012
.009
.006
.009
.000
.227
. 175
. 134
. 105
.080
.062
.049
-------
Listing of PART.DATA at 13:39:45 on DEC 2, 1985 for CCid=SHJD Page
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
363
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
08
09
10
1 1
12
13
14
15
16
17
18
19
20
01
02
03
04
05
06
07
08
09
10
1 1
12
13
14
15
16
17
18
19
20
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
LDV
LDT1
LDT2
HDGV
MC
20
01
02
03
04
05
06
07
08
09
10
1 1
12
13
14
15
16
17
18
19
.056
.047
.040
.033
.027
.023
.019
.015
.013
.01 1
.009
.008
.034
.000
. 167
. 159
. 134
. 142
.131
.080
.051
.028
.010
.098
.000
.000
.000
.000
.000
.000
.000
.000
.000
38447
34982
31832
28968
26363
23995
21843
19883
18101
16481
15007
13665
12444
0000
4100
2800
2100
1600
1200
0800
0600
0400
0200
0200
0200
0000
0000
0000
0000
0000
0000
0000
0000
02153.0
01644.2
01273.3
00955.9
0071 1 .8
00551 .9
00415.0
00298.2
00235.3
00181 .3
00135. 1
00109.3
00423. 1
000.
685.
445.
281 .
227.
157.
063.
030.
01 1 .
002.
019.
000.
000.
000.
000.
000.
000.
000.
000.
000.
0
7
7
0
0
B
7
4
1
1
6
0
0
0
0
0
0
0
0
0
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.046
0.035
0.027
0.020
0.015
0.012
0.009
0.006
0.005
0.004
0.003
0.002
0.009
000
356
232
146
1 18
082
033
016
001
010
000
000
000
000
000
000
000
000
000
000
0.09 0.20
0.20 0.46
0.21 0.47
0. 19 0.40
0.00 0.00
.
.
.
.
.
.
.
t
04 .
07 .
10 .
13 .
16 .
18 .
21 .
23 .
25 .
27 .
29 .
31 .
33 .
34 .
36 .
37 .
39 .
40 .
41 .
04 .22
05 .27
06 .31
07 .35
08 .38
09 .42
09 .45
10 .47
11 .50
1 1 .52
12 .55
12 .57
13 .59
13 .60
14 .62
14 .64
15 .65
. 13
. 14
. 16
. 17
. 18
.19
.20
.21
.21
.22
.23
.24
.24
.25
.25
.26
.26
15 .66 .26
15 .68
.27
.23
.27
.32
.36
.39
.43
.46
.49
.51
.54
.56
.58
.60
.62
.63
.65
.66
.68
.69
. 13
. 15
. 16
. 17
. 18
. 19
.20
.21
.22
.23
.23
.24
.25
.25
.26
.26
.26
.27
.27
. 18
.23
.28
.32
.36
.39
.42
.45
.48
.50
.52
.54
.56
.57
.59
.60
.61
.62
.63
. 12
. 13
. 15
. 16
. 17
.18
. 19
.20
.21
.22
.22
.23
.23
.24
.24
.25
.25
.25
.25
-------
20 .42
0.20
0.23
0.20
0.87
0.20
0.42
0.20
0.73
0.43
0.09
0. 10
0.01
13
1951 2
1963 2
1966 2
1969 3
1972 3
1975 3
1980 2
. 16 .
2.00
0.43
2.00
0.89
2.00
0.66
1 .00
0.86
1 . 10
0. 16
10.00
1 .00
.7420
.7307
.8267
.0080
.1917
. 1420
.7780
69
10
0
10
0
10
0
2
0
4
0
.27 .70 .
.00
.64
.00
.97
.00
.90
.00 2.50
.90 0.92
.70 7.00
.82 0.90
,28 .6
10.00
1 .00
10.00
0.98
Listing of PART.DATA at 13:39:45 on DEC 2, 1985 for CCid=SHJD Page
407 20 .42 .16 .69 .27 .70 .28 .64 .26
408
409
410
41 1
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428 1982 2.5580
429 1985 2.4700
430 1986 2.4260
431 1987 2.3600
432 1993 2.3175
433 1997 2.2970
434 0.75 0.40 0.44
435 0.045 0.195 0.195
436 0.017 0.050 0.050
437 20
438 01 18352 45544 53370 82288
439 02 16946 39671 46901 74984
440 03 15648 34558 41190 68328
441 04 14449 30092 36206 62263
442 05 13342 26213 31812 56737
443 06 12320 22834 27948 51700
444 07 11376 19898 24556 47111
445 08 10504 17332 21575 42930
446 09 9700 15098 18956 39119
447 10 8956 13152 16655 35647
448 11 8270 11456 14632 32483
449 12 7637 9979 12856 29599
450 13 7052 8693 11296 26972
451 14 6511 7572 9925 24578
452 15 6012 6596 8719 22396
453 16 5552 5746 7661 20408
454 17 5126 5005 6728 18597
455 18 4734 4360 5913 16946
456 19 4371 3798 5196 15442
457 20 4036 3308 4565 14071
458 21
459 1980 0.000 0.006 0.112 1.521
460 1981 0.000 0.006 0.124 1.581
461 1982 0.049 0.009 0.135 1.599
462 1983 0.104 0.013 0.141 1.592
463 1984 0.185 0.022 0.153 1.641
464 1985 0.274 0.037 0.166 1.719
-------
Listing of PART.DATA at 13:39:45 on DEC 2. 1985 for CC1d=SHJD
Page
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
1986 0.370
1987 0.475
1988 0.588
1989 0.707
1990 0.831
1991 0.960
1992
1993
1994
1995
1996
1997
1998
1999
2000
.092
.225
.354
.480
.600
.712
.816
.912
.999
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.
0.7845 0.0827
053
071
089
106
122
137
151
165
178
190
202
212
222
230
238
0.
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
. 177
. 185
. 193
.201
.208
.215
.222
.229
.237
.245
.253
.261
.269
.276
.283
0470 0.
1 .
1 .
2.
2.
2.
2.
2.
2.
2.
2.
2.
3.
3.
3.
3.
816
927
041
151
258
362
471
581
693
807
914
015
108
194
273
0419 0.0354 0.0085
------- |