&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

-------