Volume 2: Users' Manual for The Revised GAP Regional Econometric Model
THE OAP REGIONAL ECONOMETRIC MODEL:
A REVISED VERSION
Prepared for:
Environmental Protection Agency
Office of Air Programs
Research Triangle Park
North Carolina
September 25,1972
-------
FINAL REPORT
CONTRACT NUMBER 68-04-0039
VOLUME 2: USERS' MANUAL
THE OAP REGIONAL ECONOMETRIC
MODEL: A REVISED VERSION
Prepared for:
Environmental Protection Agency
Office of Air Programs
Research Triangle Park, North Carolina
Prepared by:
CONSAD Research Corporation
121 North Highland Avenue
Pittsburgh, Pennsylvania 15206
Fu-chen Lo and Kathryn Mason
September 25, 1972
-------
I
1----
I
T ABLE OF CONTENTS
Page
1. O. USERS' GUIDE: A GENERAL
DESCRIPTION. . . . . . . . . . .
. . . .
. . . .
. . . . .
1
1. 1. Introduction. . . . . . . . . . . . . . . . . . . . . .
1.2. A General Description. . . . . . . . . . . . . . . . .
1
3
2. O. PROGRAM RMS (Regional Model
Sim'l11ation) . . . . . . . . . . . . . . . . . . . .
.....
9
2. 1. Introduction. . . . . . . . . . . . . . . . . . . . . .
2.2. Job Deck Set-up... . . . . . . . . . . . . . . . . . . .
2.3. Data Input Deck. . . . . . . . . . . . . . . . . . . .
2.4. Description of the Program Deck. . . . . . . . . . .
2. 5. Program RMS (Regional Model Simulation) . . . . . .
9
10
11
32
38
3. O. PROGRAM lOA (INPUT-OUTPUT ANALYSIS). . . . . . . . 112
3. 1. Int r od ucti on . . . . . . . . . . . . . . . . . . . . . . . 112
3.2. Job Deck Setup. . . . . . . . . . . . . . . . . . . . . 113
3. 3. Data Input Deck. . . . . . . . . . . . . . . . . . . . 114
3.4. Description of the Program Deck. . . . . . . . . . . 116
3. 5. Program lOA (Input-Output Analysis) . . . . . . . . . 120
4.0
PROGRAM FBE (Feedback Effects) . . .
.......
. . . 133
4. 1. Introduction. . . . . . . . . . . . . . . . . . . . . . 133
4.2. Job Deck Set-Up. . . . . . . . . . . . . . . . . . . . 134
4. 3. Data Input Deck. . . . . . . . . . . . . . . . . . . . 135
4.4. Description of the Program Deck. . . . . . . . . . . 136
4. 5. Program FBE (Feedback Effect) .. . . . . . . . . . 140
APPENDIX: EPA Control Cost Estimates. . .
........
. . 144
o
ii
-------
I
I . -
LIST OF FIGURES AND TABLES
Figure 1. 1 Block Flow Chart of the Program RMS
Figure 3.1 Flow Chart of Program lOA
. . . .
Figure 4. 1 Flow Chart of Program FBE
. . . . .
Table 1. 1
Relations Between Simulation Option and
Equation Blocks in Program RMS. . .
iii
. . . . .
. . . . .
~
Page
6
117
138 '
7
-------
1.0
USERS' GUIDE: A GENERAL DESCRIPTION
1. 1 ,
Introduction
This is a revised users I guide of the OAP Regional Econometric
Model;:' It is intended to facilitate an efficient simulation of var ious
control strategies.
As described and demonstrated in Volume 1, air
pollution control strategies are first transformed into a set of appro-
priate model inputs. >:":< The computer program accepts thes e inputs
as "exogenous shocks" to the model system and generates as output,
from the model system, measurements of the "changes" in the key
economic variables in the Air Quality Control Regions (AQCRs) under
stud y .
The computer simulation program of the OAP Regional Econo-
metric Model includes three programs:
Program RMS (Regional Model Simulation),
Program lOA (Input-Output Analysis), and
Program F BE (Feedback Effects).
As indicated in Volume 1, the input-output model and the interregional
feedback effects were introduced in order to measure the external
. >:'See CONSAD Research Corporation, OAP Regional Econometric
Model: A Revised Version, Vol. 1, September 25, 1972.
>:<>:
-------
market of an AQCR.
The main part of the simulation program, con-
sequently, is the Prog'ram RMS - - the simulation tool for the Regional
Econometric Model.
Program RMS has been revised so that the growth rates by
AQCRs could be included to provide economic projection of future
year s.
In strategy simulation, cumulative effects over time of any
given strategy are also included.
At the san~e time, a new set of
control costs which are the preliminary estimates corresponding to
the control implied in the Federal Register of August 14, 1971, as
promulgated by EP A.
These estimates are termed as EPA cost esti-
mates in the rest of the report.
2
-------
1.2 A Get;leral Description
PROGRAM RMS
The Regional Econometric Model developed in Volume 1 consists
of four major modules or blocks of equations, namely:
manufacturing industry equations
income and regional economy equations
other industry equations
labor market equations
In the development of the Program RMS, linkages between the
modules have been rearranged in order to provide an efficient and
broad range of strategy simulations.
The arrangement of the modules
(or block of equations) in the program provides a number of l'simula-
tion options II which are related to a number of control inputs corres-
ponding to the control strategy under consideration.
At present, eight simulation options have been included, al-
though additional simulation options can be added, depending upon the
need of the user and the full exploration of the model capability.
The
currently available simulation options are as follows.
Option l: Industry control cost.
This option takes direct control costs of the manufacturing industry
in terms of reduction of profit and investment which increase the cost
of production. Control cost data, at present, is the EPA cost estimates
by AQCR by year.
3
-------
Option 2:
Benefit by AQCR.
This option takes regional benefits in terms of increase in the regional
disposable income by AQCR, resulting from changes in property values,
savings in health expenditures, etc., from cleaner air.
Option 3:
Local taxation.
This option takes increases in local taxation for possible subsidies for
industrial control costs in the implementation plan.
Option 4:
Electricity price inc'rease to industry.
This option takes the p'rice increas es for the electric power which is
being passed to industries.
Option 5: Electricity price increase to the residents.
This option takes price
are being pas sed to the
come.
increases for the electric power industry which
residents as a reduction in their disposable in-
Option 6:
Control cost by two-digit SIC.
This option takes control costs by two-digit SIC manufacturing indus-
tries rather than the aggregate control cost by AQCR used in option 1.
Option 7: Increase in property value (due to improved air quality).
This option takes a percentage increase in property value as control
input.
Option 8: Feedback effects.
This option takes the output from Program FBE to measure the inter-
regional feedback effects.
Each option consists of modules A, B, C, and D.
However,
depending upon whether aggregate o'r two-digit SIC manufacturing sec-
taring and other simulation specifications are available, ten blocks of
4
-------
equations, each related to one module, were included in this Program
RMS.
A block flow chart for the simulation options available in Pro-
gram RMS is shown in Figure 1. 1, and Table 1. 1 shows the relation
between simulation option module of the Regional Econometric Model
and equation blocks included in the block flow chart.
For example, in simulation option 1, the industry control costs
are read in from EPA cost estimates control input data.
This leads
to a reduction in both profits and investment of the manufactur ing in-
dustry, thus affecting the entire manufacturing industry module, and
in turn regional income, consumption, government sector of module
B and other industry of module C, and finally labor market of module
D.
If other simulation options were tried, the sequence will refer to
corres ponding blocks in the Program and "net" effects will be summa-
rized in the program output.
Usually each strategy simulation will
involve more than one simulation option.
An example of the simu1a-
tion options involved in a given strategy will be included in Section
2.3.2 of this volume (preparation of parameter cards).
PROGRAM lOA
This program has the capability of aggregating a 100-sector
United States input-output table to a desired number of sectors.
Cur-
rently, a 42-sector 1-0 model has been used (see Appendix D, Volume 1
of the May 15 repo'rt).
The program then provides a new 1-0 coefficient
5
-------
j- - - - - .- - - - - - -I
I I
I ,
I Set 2 A I
8+ Manufacturing I
6 I industries
I
I (2 -digit SIC) I
I
I I
I I
I I
I
I I
I Set 1 A '; I
I Manufacturing
I industries
I (2 -digit SIC)
I
I
I
------------
0'
Set3 A
Manufacturing
industries
(2-digit SIC)
(Two-digit SIC detail)
o
= simulation options
FIGURE 1. 1
BLOCK FLOW CHART OF THE
PROGRAM RMS
Set4 A
Manufacturing
industries
Set 8 B
Income deter-
mination
Set 5 B
Income deter-
mination
Set 9 A
Manufacturing
industries
Set 6 C
Other industrie
Set 7 D
Labor market
OJ
-------
Table 1. 1.
~e .a :ions Je :ween Simulation Op :ion anc Equa :ion B .ocks in Program Rl'v:S
-..J
Simulation 1 2 3 4 5 6 7 8
Option of Industry Benefit Local Electricity Electricity Control Change in Feedback
. Prog ram Control by Taxa- price in- price in- Cost by property effects
Equation RMS Cost AQCR tion crease to creas e to 2 -digit value
Blocks in
Program RMS industry res identials SIC
Block No. 1 A
Block No. 2 A
Block No. 3
Bldtk No. 4 A A
Block No. 5 B B B B
Block No. 6 C C C C C C C C
Block No. 7 D D D D D D D D
Block No. 8 B B B B
Block No. 9 A A A A
Mod ule :
A
B
C
D
Manufacturing Industry
Income and Regional Economy
Other Industry
Labor Market
-------
1
r:-
matrix an~ the corresponding Leontief inverse matrix with new sector-
ing.
Finally, the program takes changes in final demand, either change
by sector or change by average propensity to consume by sector and
produces the change in gross domestic product by sector.
Thus, the changes in the final demand, both as a result of cos ts
and benefits of air pollution control, will caus e the changes in the gros s
domestic product (value of shipments by sector).
Output from Program
lOA becomes the input in the Program FBE for distribution.
PROGRAM FBE
This program takes the output from Program lOA and first trans-
fers the value of shipment into value-added by sector and then distri-
butes the feedback effect to each AQCR with the l"egional market share
matr ix.
The output from Program FBE becomes the input deck for
simulation option 8 in the Program RMS.
8
-------
2. 0 PROGRAM RMS (Regional Model Simulation)
2. 1 Introduction
Program RMS is a large scale computer simulation program
with a 2483 statement program which includes 13 subroutines to solve
both linear and non-linear simultaneous equation system.
For the non-
linear system, an iterative algorithu1--an extended Newton1s method>:<
has been introduced.
Since the program has been designed for a broad range of
poli<::y questions, the simulation of any given strategy requiring a
proper preparation not only of control input (see Section 2.3) but also
the choice of the simulation options appropriate to the strategy.
At the present state, EPA cost estimates by AQCR. have been
included as basic control input. >«>:<
Detailed specification of the prograrn prepaJ."ation is given in
the following sections.
>:< See M. K. Evans and L. R. Klein, The Wharton Econometric
Forecasting Model, University of Pennsylvania, 1967, Chapter IV.
>:<>:< See AppendiX..
9
-------
2.2 Job Deck Set- Up
The Program RMS (Regional Model Simulation) is a 2, 100 state-
ment FOR TRAN program which uses approximately 392K core and runs
80 seconds on an IBM 360/75, depending upon the scope of the strategy
to be simulated.
~
n:st
/ p~ogramrJ\
\ Control ca)
FBE Deck-ll
if benf.I$I~
Benefit \
card, .-J
AQCR codes
1\. Years par~aJ
meter card
,
~Local:1
\" nomic Dat I
I
/ );c'
/ /GO
Prog. RMS /
(Reg. Mode
Simulation)
Program
Deck
One Set for
Each Year
Requested
One Set for Each Run
//Job
Job control
lang. (JC L)
10
-------
"
2. 3 Data ;Input Deck
The input section has undergone considerable modification.
What
was then referred to as the STANDARD DA TA INPUT DECK, is now,
after some reordering and a few additions, read from a tape.
The
data stored on this tape is needed for every RMS simulation regardle~s'
of strategies or options to be performe.d.
A detailed description of the
contents of the tape follows.
2.3. 1 Standard Data Input Tape
The tape includes 14 sets of parameter and input data read from
a total of 5,826 cards arranged in the following order.
1.
Regional Model Coefficients: 20 cards, one card for each
two-digit SIC and one for all manufacturing industries
Gols. 5-10
.15, two-digit SIC code
Cols. 11-20
F 1 0.4, labor share from production
function by industry
Gols. 21-30
Cols. 31-40
Cols. 41-50
F 1 0.4 capital share by industry
F 1 0.4, depreciation rate by industry
FI0.4, the coefficient of profit in
investment function by industry
Cols. 51-60
FI0.4, ratio of value-added to capital
by industry
11
-------
AQCR Codes and Names: 92 label cards, one for each of the
91 AQCRs for which data was available and one which is the
label card for the net effect over all AQCRs studied
.2.
Cols. 1-4
Co1s. 5 -80
14, AQCR code number (blank for the 92nd
card)
19A4, name of AQCR associated with the code
in co1s. 1-4
3. Relative Codes Used by the Program for the 100 AQCRs:
5 cards
Card 1 2 014, relative codes for AQCRs 1-20':<
Card 2 2014, relative codes for AQCRs 21-40':'
Card 3 2014, relative codes for AQCRs 41-60':<
Card 4 2014, relative codes for AQCRs 61-80':'
Card 5 2014, relative codes for AQCRs 81-100':<
4.
Fue 1 Submode 1 Coefficients: 58 ca rds
(a) Share coefficients from fuel equation by two -digit SIC
(19 cards)
Cols. 9 - 10
Co1s. 11-20
Cols. 21-30
Co1s. 31-40
Cols. 41-50
Gals. 51-60
12, code of the two -digit SIC for which this
data applies
FI0.5, coefficient of coal
FI0.5, coefficient of coke
F 10. 5, coefficient of oil
F 1 O. 5, coefficient of natural gas
F 1 O. 5, coefficient of electricity
>:, Code is 0 if no data was available for a given AQCR.
12
-------
(b) Miscellaneous fuel coefficie nts by two -digit SIC including
one coefficient for all manufacturing industries (20 cards)
Cols. 9 -10
Cols. 11-20
Cols. 21-30
Cols. 31-40
12, code of the two -digit SIC for which this
data applies
FIO.5, constant determined in the invest-
ment equation
F 1 O. 5, cons tant multipJier from the fuel
equation
F 10. 5, coefficient of the energy demand
equation
(c) Coefficient relating consumption of electrical energy to
value -added by two -digit SIC (19 cards)
Cols. 9-10
Cols. 11-20
5.
12, code of the two-digit SIC for which this
data applies
F10.6, coefficient
PJQCR Economic Data: 455 cards, five cards for each AQCR
for which data was available
First Card:
Cols. 1-10
Gols. 11-20
Cols. 74-76
Second Card
Cols. 1-10
Cols. 11-20
F 10. 1, labor force in each AQCR
F 10. 3, unemployment rate in each AQCR
13, AQCR code
F 1 O. 1, 1967 employment in manufacturing
industries in each AQCR
FIO.l, 1967 wage bill for manufacturing
industries
13
-------
Cols'. 21-30
Cols. 41-50
Cols. 74-76
Third Card:
Cols. 1-10
Cols. 11-20
Co Is . 21-30
Cols. 74-76
Fourth Card:
Cols. 1-10
Cols. 74-76
Fifth Card:
Cols. 1-10
Cols. 11-20
Cols.21-30
Cols. 74-76
~
FIO.l, 1967 value-added in manufacturing
industries
FlO.l, 1967 investment in manufacturing
indus tr ie s
13, AQCR code
FlO.O, 1967 total regional personal income
F 10. 0, 1967 local gove rnment revenue
F 1 O. 0, 1967 local government expenditures
13, AQCR code
FIO.O, 1967 total personal consumption
13, AQCR code
FIO.l, 1963 number of employees in manu-
facturing industries
FIO.l, 1963 wage bill for manufacturing
incomes
FIO. I, 1963 value-added for manufacturing
industries
13, AQCR code
14
-------
6.
Electric Consumption (10 million KWH) by AQCR: 91 cards
Cols. 1- 10
Cols. 11-20 .
Cols. 21-30
Cols. 31-40
Cols. 74-76
F10.0, total electric consumption in each
AQCR
FIO. O~ residential electric power consump-
tions
F 10.0, electric power consumed by manu-
facturing industries
F 1 O. 0, electric po\ver consumed in other
industries
13, code number for the AQCR for which
data applies
7.
Detailed Two-Digit SIC Data by AQCR: 2,089 cards
(a) Data for 1963
Co1s. 1-10
Co1s. 11-20
Cols. 21-30
Cols. 74-76
Cols. 77-78
FIO.1, number of employees
FIO. 1,' wage bill
FIO.l, value-added
13, AQCR code
12, two-digit SIC code
A "9's" card follows the 1963 data.
(b) Data for 1967
Cols. 1-10 F 1 O. 1, number of employees
Co Is. 11-20 FI0.1, wage bill
Cols. 21-30 F10.1, value -added
Cols. 41-50 F 10. 1, investment
15
-------
Cols. 74-76
Cols. 77-78
13, AQCR code
12, two-digit SIC code
A 1191s" card follows the 1967 data.
8.
Detailed Fuel Information by Two-Digit SIC by AQCR
There are two cards for each AQCR/two-digit SIC combina-
tion (2,058 cards)
Card 1
Cols. 1-10
Cols. 21-30
Cols. 31-40
Cols. 41-50
Cols. 51-60
Cols. 74-76
Cols. 77-78
Card 2
Cols. 1-10
Cols. 11-20
. Cols. 21-30
Cols. 31-40
Cols. 41-50
Cols. 51-60
F 1 O. 0, total cost of fuel consumption (all
type s)
F10.0, quantity of coal
F10.0, cost of coal
F 10. 0, quantity of coke
F 10.0, cost of coke
13, code of the AQCR for which this data
applies
12, two-digit SIC code
F 1 O. 0, quantity of oil
F 1 O. Q, cost of oil
F 1 O. 0, quantity of natural gas
F10.0, cost of natural gas
F10.l, quantity of electricity
F 1 O. 0, cost of electricity
16
-------
Cols. 74-76
13, AQCR code
"I
. A "91 s" card follows the fuel data.
9. Feedback/Distribution Data by AQCR/Two-Digit SIC:
364 cards .
(a) Feedback/distribution to manufacturing industries
(3 cards per AQCR = 273 cards)
Card 1
Cols. 1-70
Cols. 74-76
Card 2
Cols. 1-70
Cols. 74-76
Card 3
Cols. 1-70
Cols. 74-76
7F10.4, feedback/distribution to two-digit
SICs 1-7 .
13, Code of the AQCR for which this data
applie s
7FI0.4, feedback/distribution to SICs 8-14
13, code of the AQCR for which this data
applie s
5F10.4, feedback/distribution to SICs 15-19
13, code of the AQCR for which this data
applies
(b) Non-manufacturing industries feedback/distribution by
AQCR (91 cards)
-Cols. 11-20
F 1 O. 3, feedback distribution
17
-------
10.
11.
12.
13.
Growth rate s for 1970: 104 cards, two cards for each of
the 51 AQCR for which detailed two digit SIC growth rates
were available and two codes for national average growth
rate s by two digit SIC
Card 1
Cols. 1-75
Cols. 78-80
Card 2
Cols. 1-30
Cols. 78- 80
F5. 0, growth rates for SIC 20 through
SIC 35 (excluding SIC 21)
13, AQCR code
F5. 0, growth rates for SIC 36 through
SIC 39, plus the average growth rate for
all manufacturing industries
13, AQCR code
Growth rate s for 1975: 104 cards, same format as above.
Growth rate s for 1980:
104 cards, same format as above.
Total personal income: 153 cards, three cards for each
AQCR for which income data was available
Card 1
Cols. 11-20
Cols. 78-80
Card 2
Cols. 11-20
Cols. 78- 80
Card 3
Cols. 11-20
Cols. 78-80
Fla. 0, Total personal income for 1970
13, AQCR code
Fla. 0, total personal income for 1975
13, AQCR code
FlO.O, total personal income for 1980
13. AQCR code
18
-------
14.
Total personal income: 38 cards, one card for each of
the AQCRs with code greater than 54 for which no detailed
GBE growth factors were available
Cols. 11-20
IlO, AQCR code
Cols. 31-40
F 1 O. 2, total per sonal income growth rate,
average for 1959-1968
This conclude s the de scription of the nine data sets stored on tape
required for any run of RMS.
19
-------
';\
. 2.3.2 Regional Economic Data and Coefficients
1., Population Growth Rates for 91 AQCRs (6 cards),
Format: (16F5.2)
2.
Local Economic Data (91 cards)
Cols. 1-10 1967 Population (FI0. 0)
Cols. 11-20 Local Taxation (F 10. 0)
Cols. 21-30 Property Tax (per capita) (FIO. 0)
Cols. 31-40 Property Valuation (FI0. 0)
Cols. 76-78 AQCR Code (B)
3.
New Investment Equation Coefficients- -one set for each
industry (20 cards),
2. 3. 3 Input Parameter Cards
1.
The Years Parameter Card, for the time series run (with
constant modified switch)
Cols. 1-4 Total number of years included
Cols. 6-10 First year (IS)
Cols. 11-15 Second year (IS)
Cols. 16-20 Third year (IS)
and so on until all the years have been stated, up to 11 years.
20
-------
2.
AQCR Codes
Cols. 1-5
15, number of AQCRs for which updated data
is needed
Cols. 6-10
15, AQCR code for the first city
Cols. 11-15
15, AQCR code for the next city
and so on in cols. of 5 until all the city codes required have
been specified.
3.
The Benefit Card
Cols. 1-10
FI0.2, estimated national benefit from
cleaner air in terms of change in total
consumption expenditure; in the fast runs
this number has been set to either $10 or
$15 billion for both the 3 and 5 year runs
(in millions on card)
Gols. 11-20
FI0.2, estimated national benefit less
compact of cost due to government as sis-
tance. If there is not government assis-
tance, this number should equal the one
in the previous ten columns (in millions)
Cols. 21-30
FI0.2, normalizer: if benefit = $10 billion,
it equals 2, 000. If benefit = $15 billion, it
equals 3, 000.
Cols. 31-40
FI0.2, percentage reduction in manufac-
turing costs resulting from government as-
sistance
4.
FBE Deck. The FBE deck (364 cards) is a set of output from
Program FBE where changes in level of benefits were intro-
duced through use of Programs lOA and FBE. This deck is
needed only if benefit is not equal to $10 billion.
21
-------
1- .
I
.- ---- - - -. --. - --~- -----,,-,-- -.. ... ... -- - .~-
5.
Input Parameter Card
A 11 of the following cards are nece s sary for each year
reque sted.
5. a.
Card Column
Variable Name
Description
1-4 (14)
INPT (1)
Information for aggregate
manufacturing industries
des ired.
(0 or 1) - control input
data will be supplied.
(0) - No
(1) - Yes
5-8(14)
NOPE
Number of simulation
options to be purs ued
for aggregate manufac-
turing industries (0, 1,
2.~ 3, 4 , 5 ) .
9-12
+(NOPE-l)
* 4 [NOPE :(-,
14]
LPARA(I),
1= 1, NOPE
Type of simulation options
to be pursued in this run.
Specify "NOPE" options.
o - none
1 - Industry control cost
2 - Benefit by AQCR
3 - Change in local taxa-
tion
4 - Electricity-price in-
~reas e to indus try
5 - Electricity-price in-
crease to residential
us e rs
N through
N + 3
MCNT
If INPT( 1) = 1, MCNT
specifies the number of
cities for which cost
changes a re given or if
INPT(l) = 0, MCNT is
specified as 1.
22
-------
Ca rd Column
(N + 3) through
(N -+ 7)
[14]
(n + 8) through
(N + 11)
[14]
Variable Name
INPT (2)
NOPE1
(N + 12) through LPNT
(N + 1 5 )
[14]
(ri + 16) through LPNT2(I)
(n + 19) +
(NO PE 1 - 1)
*4
[NOPE 1 * 14]
M through
(M + 3)
NUMCIT
2'3
"
Description
Detailed two dig it SIC
infornlation to be included
(0 or 2) and simulated
o - not wanted
2 - des ired
Number of simulation op-,
tions to be purs ued fo r
detailed two-digit SIC
data (0, 1, 3) -
= 1
Types of detailed two-
digit SIC simulation op-
tions to be purs ued,
specify "NOPE I" options.
o - none
1 - control cos t by 2 dig it
SIC. This is s imu-
lation option number
6.
3 - national feedback/
distribution affect.
This indicates simu1
simulation option
number 8 is to be
pursued.
The number of cities for
which this cont ro1 ca rd
specification applies.
-------
r
I
I
Card Column Variable Name Description
The next five variables
control the printout
options.
L through 1PRN1 =0 312{A) summary table
(L+3) requested but delete
[14J printout of options 1,
4 & 5.
=1 Print results of all
strategies
(L+4) through 1PRN2 =0 Do not print 312(A)
(L+7) table (Le., summary
[14J table of 1 st five strate-
gies
(L+8) through IPRN3 =0 overrides all preceed-
(L+ll) ing requests and prints
[14J out only the net sum-
mary table s for all op-
tions requested
(L+ 12) ISB U G =1 Turns on debug print-
(L+ 15) out patche s
[14J =0 Skips option
(L+ 16) ISBUG2 =1 Detailed data tables
(L+ 19) per year are given
14 =0 Skips option
(L+20) 1CTAB =1 Convergence tables for
(L+23) simultaneous systems
[14J design
=0 Skips option
24
-------
Card Column
Variable Name
De scription
(L+24 )
(L+27)
[I4J
IGOV
= 1 Impact of costs with
government assistance
. =0 Impact without govern-
ment assistance
In addition to the above information, the user should realize that
different simulation options may require different control input data.
25
-------
AN EXAMPLE OF.PARAMETER CARD SPECIFICATION
Ca rd Column
Punch
Description
4
1
Aggregate control cost for manufac-
tur ing, etc.
8
2
Two simulation options are to be pur-
. sued in this strategy simulation for
agg regate pa rt of the model system.
12
4
Option 4 is chosen.
16
5
Option 5 is als 0 chos en.
19-20
.1
Control input data for 91, AQCR I s is
included.
N
'"
24
2
The simulation is to pursue at least one
option involving detailed two-digit
SIC information.
28
1 .
1 - by each AQCR
32
1
One option to be pursued for detailed
2 -digit SIC information.
36
..3
Option 8 is to be pursued.
39-40 .
91
Two-digit SIC data available for 91
AQCR's.
';)
-------
I '
I
I
AN EXAMPLE OF,PARAMETER CARD SPECIFICATION
Card Column Punch Des cr iption
41-48 B LA NK The IPRNS 1, IPRNZ'\ printout options
are bypas sed
49-52 0 Final summary table reque sted
53-56 B LA NK No debugging switches needed
57-60 B LA NK Avoide update tables
61-64 B LA NK Delete convergence tables
N 65-68 1 Impact of cost with government
-J assistance requested
-------
5.b.
Title Cards. On the next two cards, specify the heading
for the final summ.ary table. Example:
Card 1
Cols. 1-80
Five Year Extended Implementation
Card 2
Cols. 1-80
Without Government Assistance
5. c.
Percentage Distribution of Benefit Card
Cols. 1-5
F5.0, percentage distribution of total
benefit to first year in the time series
Cols. 6-10
F5.0, percentage distribution for the
second year (example, 15. for 15%)
Cols. 11-15
F5.0, percentage distribution for the
third year
Cols. 35-40
F5. 0, total for all numbers in previous
columns on this card
5. d.
One Blank Card Necessary Here.
The initial parameter card required for each year makes it pos-
sible to request different simulation options.
As a consequence, the
input control data is subject to a corresponding variability.
One of
the most mechanical adjustments occurs if a limited number of AQCRs
are requested in parameter card # 1.
On response, the control cost
input deck should only include data for the specified number of AQCRs.
28
-------
A more complicated problem results from the fact that the \.~ser
must provide control cost data for seven of the eight possible options
for each year of the run.
Furthermore. most changes in option speci-
fication must be accompanied by modifications on the cost input sec-
tions of the program of which these are two: one for the regional
cost options (one through six) and the other for the detailed SIC
options (six through eight).
These changes are left to the discretion
and judgment of the user and are dependent on the type and form of
the cost data that may become available.
The one restriction is that
the regional cost data must preceed the detailed cost data, as can be
seen from the example that follows.
CONTROL COSTS DECK:
The EPA cost estimates by AQCR
has been assembled as input for this program.
Two control decks
exist:
1,
The three year straight implementation deck (1973-75).
2.
The five year extended cost distribution deck (1973-77).
This data was set up to accommodate the following options:
1.
Manufacturing industrie s control costs (in the form of
stationary combustion and solid waste).
29
-------
4.
Electric power price increase passed to manufac-
turing industries.
5.
Residential price increased on electricity.
6.
Industrial cost impact on high emission two-digit SIC.
8.
National feedback effect to the region through detailed
SICs.
Both the three and the five year cost decks have the following format
for each AQCR.
Card 1 (necessary for regional options)
Cols. 1-10
F10.5, investment cost - stationary combus-
tion
Cols. 11-20
Cols. 21-30
Cols. 71-73
Cols. 74-76
F 1 O. 5, annual cost - stationary combustion
F10.5, annual cost - solid waste
13, years (example - 72)
13, AQCR code
Card 2 (impact on detailed SIC-investment cost and necessary
for option number six)
Cols. 1-10 F10.4, investment cost for SIC 20
Cols. 11-20 FIO.4, investment cost for SIC 26
Cols. 21-30 FIO.4, investment cost for SIC 28
Gals. 31-40 FIO.4, investment cost for SIC 29
Cols. 41-50 FIO.4, investment cost for SIC 32
30
-------
Cols. 51-60
F10.4, investment for SIC 33
Cols. 61-70
FIO.4, investment cost for all manufacturing
industries
Cols. 71-73
13, year
Cols. 74-76
13, AQCR code
Card 3 (impact for detailed SICs - annual costs card, necessary
for option numbe r six)
Cols. 1-10
Cols. 11-20
Cols. 21-30
Cols. 31-40
Cols. 41-50
Cols. 51-60
Cols. 61-70
F 1 0.4, annual cost for SIC 20
FIO.4, annual cost for SIC 26
F 1 0.4, annual cost for SIC 28
F 1 0.4, annual cost for SIC 29
FIO.4, annual cost for SIC 32
F 1 0.4, annual cost for SIC 33
F 1 0.4, annual cost for all manufacturing
industries
Cols. 71-73
year
Cols. 74-76
AQCR code
Note that by changing the benefit card, one variable on the parameter
card (IGOV =0 or I). and the second title card.
The program has been
run with and without government assistance and with corresponding
variations in national feedback effects.
31
-------
2.4 Desc::iption of the Program Deck
',\
The program RMS consists of a MAINLINE and 13 subroutines,
PRIN1, VALUE, VALUE1, VALUE2, VALUE3, VALUE4, SIMUL,
GRADN and VERT, OLETS, PHELP, DEBUG, ABCD, a total of 2483
cards.
The program cards are sequenced along the right-hand edge in '
columns 73-80 with the prefix RMS.
The program begins by reading the standard data input deck from
tape and cards.
The program parameters are then read from the input
parameter card(s).
The parameters specify to the program how many
'and which of the various simulation options have been chosen by the
user.
The program now reads from the control input deck any special-
ized data which may be required by the first simulation option which
has been specified on the input parameter card(s).
The program now
develops the specifications for the equation sets which are required for
this simulation option.
Control is next transferred to subroutine
SIMUL which controls the solution of the equation sets.
SIMUL calls
subroutine VALUE which determines which type of set of equations is
to be solved and then calls for the appropriate subroutine (VALUEl,
VALUEZ, VALUE3 or VALUE4), which actually sets up the complete
system of equations in the manner necessary for solution.
Control is
returned to SIMUL which now calls subroutines GRADN and VERT to
32
-------
do the actual solving of the non-linear parts of the equation system.
After the new values for the system have been determined, control is
returned to the MAINLINE.
The solution procedure involving sub-
routines SIMUL, VALUE, GRADN and VERT is repeated until all sets
of eq"uations needed to implement the desired simulation option have
been solved.
The program now calls subroutine PRIN 1 which prints
out the or iginal values, the values afte r iteration and the pe rcent
change.
If all of the simulation options specified by the user on the
input parameter card(s} have been pursued, the program now prints
out a table show ing the total effect of the combined res ults of the s imu-
lation options and the program terminates or begins the next year I s
simulation.
If all specified options have not been pursued, the pro-
gram again reads from the control input deck any specialized data
required for the next simulation option which has been specified and
proceeds as above.
33
-------
Read stan-
dard input
da.ta
Read para-
mete r input
ca rds
Yes
ead air pollu-
tion control cost
data for AQCR's
re ue sted
@
2
Go to (2,3,4)
based on in-
put param.-
etel."S
Cost affects
p!"ofit direct-
ly or through
increase in
cost of elec.
MACRO FLOWCHART
OF RMS
No
8
Yes
8
cp
Cost affects
income
Cost of con-
trol affects
local taxation
34
-------
11
No.
No
Go to. 9, 10,11
based on inpu
parameters
Calculate
effect of
. national feed-
back from
I/O table
Read AQCR
code and no. of
industries (con-
trol in ut)
Read control
cost data for the
no. of industries
specified above
Solve system of
equations using
new costs
-------
Print out
detailed SIC
information
@
36
-------
'olve system of
equations for new
values from
changed costs
Print out
re sults
7
Print out net
results for all
simulation
options
..
37
-------
2.5 Program RMS (Regional Model Silnulation)
~
38
-------
LEV~- 20.1
ISN 0002
ISN 0003
ISN 0001f
ISN 0005
ISN 0006
ISN 0007
ISN 0008
ISN 0009
w. I S N 0 0 1 0
-.!) ISN 0011
ISN 0012
ISN 0013
ISN 0011f
ISN 0015
ISN 0016
ISN 0017
ISN 0018
15N 0019
15N 0020
15N 0021
ISN 0022
ISN 0023
ISN 002"
ISN 0025
ISN 0026
(MAY 7 )
05/360
F'ORTH'.J H
(OMP L~~ O;)'ONS - NA'1E= '1A N,O;)"=OO,LINE(NT=QO,S:z::=OOnOK,
SOURCE,B(D,NOLIST,NODECK,LOAD,MAP,NOEDIT,NOID,NOXREF
c
(
(
(
(
C
(
(
(
(
C
...e...............e.o..~.......O.$.O..OOOO&.~OOO.....eo.....
e
. A SIMULATION OF E(ONOMIC EFFECTS OF AIR POLLUTION
ON AIR QUALITY CONTROL REGIONS
PROGRAM RMS
PREPARED BY CONS AD RESEARCH CORPORATION
REVISED AUGUST 1972
.
.
.
.
.
(ONTROL .
...............................~.............................
REAL KJNDX,LJNDX
DIMENSION A(20),
DIMENSION ALPHA(20),
DIMENSION BEN(8)
DIMENSION C(92), CD(92),
DIMENSION COFI(91),COST(91,S)
DIMENSION 0(20),01(91)
DIMENSION DZ(91),
DIMENSION DNBARN(92),
DIMENSION EOLDH(91,20)
DIMENSION EOLD(91,20,S),
DIMENSION GD(9Z),
DIMENSION IAQCR(92,ZO),
DIMENSION INPT(2),
DIMENSION KEEP(S),
DIMENSION LPARA(10),
DIMENSION PDN(91),
DIMENSION PICON(20),
DIMENSION PPIO(91)
DIMENSION PXD(91),
DIMENSION PZO(91),
DIMENSION QBARN(92),
DIMENSION QCN(92),
DIMENSION QMN(92),
DIMENSION QTN(92),
AB(91),
8(91,20),
ALIT(20)
BETA(ZO)
CN(92)
DELTC(91),
DNR(91),
DNBAR(9Z)
EO(91,5)
G(92),
GN(9Z),
INO(ZO) ,
JKEEP(20),
LABELl 10,ZO),
LPNT2(3),
PI(92,ZO),
PID(9Z) ,
GAMA(20,5)
HOLO(S)
INOCI(92,20)
KCITY(100)
LKEEP(3,S)
MCITY(91)
PICOEF(ZO)
PIN ( 92)
P X I D ( 9 I ) ,
QBAR(92),
QC(92),
IHH92),
rH(92),
RNET(91 ),
PXKD(91)
QBARD(92)
QCO(92)
QMO(92)
QTD(92)
SUMN(92,20)
.
RMS
RMS
RMS
RMS
RMS
RMS
Rt1S
RMS
RHS
RMS
Rt1S
RMS
RMS
RMS
RMS
RMS
Rt~S
RMS
RMS
RMS
RMS
RMS
RMS
RMS
Rt1S
RMS
RMS
Rt1S
RMS
RMS
Rt1S
RMS
RMS
RMS
RMS
RMS
.
.
.
.
.
DATE
72.22
1
2
3
Q
5
6
7
8
9
1 0
1 1
12
1 3
1 If
1 5
1 6
1 7
18
19
20
21
22
23
2..
25
26
27
28
29
30
3 1
32
33
3tf
35
36
-------
C R'S 37
ISN 0027 DIMENSION QGIVENIJ) RMS 38
C RMS 39
ISN 0028 DIMENSION SUMND(92), SUMNN(92), T 1 92) RMS lfO
ISN 0029 DIMENSION TD(92), TN(92), TR(91) RMS .. 1
ISN 0030 DIMENSION U(92), UD(92), UN(92) RMS lf7-
ISN 0031 DINENSION VK(20), WI91,20), WBI91,20) Rr1S lf3
ISN 0032 DI~1ENSION WBD(91), WD(91), WN(92) RMS If..
ISN 0033 DIMENSION WT(91 ,20), XDI92>, XDN(92) RMS ..S
ISN 003" DIMENSION XD1(91), XII92,20), XID(92) RMS "6
ISN 0035 DIMENSION XINI97.), XK(92,20), XKD(92) RMS "7
ISN 0036 DIMENSION XKN(92), XKTI91,20), XL(92) RMS lf8
I ISN 0037 DIMENSION XLD(92), XLNI921, XND(92) RMS lf9
ISN 0038 DIMENSION XNBAR(92), XNN(92), XNT(92) RMS So
ISN 0039 DIMENSION XX(92,20), XX63(91,20), YR(91) RMS S 1
ISN 00"0 DIMENSION YY(92), YYD(92), YYN(92) RMS 52
ISN OOlfl DIMENSION YGROW(38) RMS 53
ISN 00"2 DIMENSION ZERO( 1932) ,ZZ(91 ,20) RMS 5q
C RMS S5
C RMS S6
C RMS 57
~ C RMS 58
o
ISN 00"3 DIMENSION XLH(92) ,UH(92) ,SUMNH(92,20) ,WBHI92,20) ,XXH(92,20) RMS S9
ISN OOlf.. DIMENSION XIH(9Z,20) ,YYHI9Z) ,TH(92),GH(92),eHI92),XNTHI92) RMS 60
ISN 00"5 DIMENSION XNBAH(92),QTH(9Z) ,QMH(92),QeH(92),QAARHI92),XINCI60,3) RMS 6 1
ISN 00lf6 DIMENSION GRFACI100,20,3),IYEAR(11),ICHOI91) ,XXT(92,20) RMS 62
e RMS 63
C RMS 6lf
ISN 00..7 DIMENSION SCel (3) ,5((215) RMS 65
C RMS 66
C RMS 67
C RMS 68
C RMS 69
ISN 00..8 DIMENSION HHIIJ),HHA(7) RMS 70
ISN 00lf9 DIMENSION ISPIN(]) RMS 7 1
ISN 0050 DIMENSION TITLE(20,Z) RMS 72
ISN 0051 DIMENSION DELTR(5) RMS 73
ISN 0052 DIMENSION OFUI (5) ,OFUA(5) ,FU( 100,5) ,FUll 100,5) RMS 7..
-------
:SN 0053
ISN 0051f
lSN 0055
ISN 0056
ISN 0057
ISN 0058
ISN 0059
ISN 0060
ISN 0061
ISN 0062
ISN 0063
ISN 006&f
~ ISN 0065
......
ISN 0066
ISN 0067
ISN 0068
ISN 0069
ISN 0070
I SN 0071
ISN 0072
ISN 0073
D . MENS I ON
DIMENSION
DIMENSION
DI~1ENSION
DIMENSION
DIMENSION
DIMENSION
~:~UN 92 ,QI'V(20,3) ,PIH(92,20 ,OZ:RO(IIOIf)
XKH(92,20)
OTP(92) ,OTA(92) ,OTO(92) ,OPOP(92) ,OAPV(92)
TP(92) ,TO(92) ,TA(92) ,APV(92) ,POP(92) ,TPD(92) ,TOD(92)
TAD(92) .APVH(92) ,POPH(92) ,TPH(92) tTOH(92) ,TAH(92)
TON(92) ,TAN(92) .TPN(921 ,POPGR(91)
APVD(92)
c
~MS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RI1S
RMS
RMS
COMMON IPRITTI TITLE,BENGOV RMS
COMMON/NPRINI LABEL,HCITY,KCITY, LSTRAT,INPT,IKSW.IHELP,IFEEDRMS
COMMON I IPRIN IINTT,LPNTT ,JyEAR,IAQCR,IJCIT.MIN,IA312 RMS
Rf-1S
RMS
RMS
RMS
COMMON I ZEROT I PIN,XIN,XKN,XnN,YYN,GN.TN,CN,DNBARN,XNN,XLN,WN,
1 UN,5UMNN,QBARN,QMN.QCN,QTN
2 ,TPN.TON,TAN
C
C
COMMON IFATI F( 1000) ,LSTORE,KPIST,ARAB( 10) ,LCHECK
c
COMMON ISENSI ICTAB,NBEFOR,NEWCIT
COMMON ICATI Z(20,20)
c
COMMON 151MI Ml ,M2,MCA5ES,Y( 1000) ,XI 1000) ,MARK,LVAR,MKEEP
c
COMMON/ALET/PERXX,PERW,PERWB,PERSN,PERP,PERXI,PERXK,pERYY,PERNB,
1 PERT,PERG,PERC,PERQM,PERQA,PERQC,PERTP,PERTO,PERTA
COMMON I BUG I 158UG,LREAD,LWRIT
COMMON IBUGGI IPRNI,IPRN2,IPRN3
c
COMMON/CPRIN/PXD,PXID,PXKD,PPIO,PDN,EOLD,ED
c
COMMON IAPRIN/PI,PID,XI,XID,XK,XKD,XX,XD,C,CD,YY,YYD,5UMN,SUMND,
1 XNBAR,DNBAR,XNT,XNO,U,UD,XL,XLD,G,GO,T,TD,QT
2 ,QTD.QC,QCD,QM,QMD,QBAR,QBARD
C
c
COMMON/EVERY/OWB(92) ,OW(92) ,OX(92) ,01 (92) ,OY(92) ,OT(92) ,OG(92),
1 OC(92) ,OXBAR(92) ,OQC(92) ,OQM(921 ,OQB(92), GRFAC,XINC,VK,
2 YGROW,IZS
75
76
77
78
79
AO
8 1
82
83
'8 If
85
86
87
88
89
90
9 1
92
93
91f
95
96
97
98
99
100
1 r:J 1
102
103
101f
105
106
107
108
109
1 10
1 1 I
11'1
-------
T
ISN 007&f
ISN 0075
ISN 0076
ISN 0077
ISN 0078
ISN 0079
ISN 0080
ISN 0081
ISN 0082
,.j:>..
N ISN 0083
ISN 008&f
ISN 0085
ISN 0086
ISN 0087
ISN 0088
ISN 0089
ISN 0090
ISN 0091
ISN 0092
ISN 0093
ISN 009&f
ISN 0095
3 ,OTP,OTA,O.O,O;)OP,OA;)V,PO;)GR
COMMONI WNEWWI WAGECOCZO,3)
RMS
RMS
RM5
RMS
RMS
RMS
RM5
RM5
RM5
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RM5
RMS
RMS
RHS
RMS
RMS
RMS
RMS
RM5
RMS
RMS
RMS
RMS
RM5
RMS
RMS
RMS
RM5
c
COMMON /ATAXI TP,TPD,TO,TOD,TA,TAO
C
EQUIVALENCE COZEROCI ),OWS(I»
C
EQUIVALENCE CZEROCl),PINCl»
C
C
EQUIVALENCE CINDCI(I),XX(I»
C
C
C
DATA LKEEPI 3,1,I,Cf,3,2,5,&f,3,O,5,&f,0,0,SI
DATA SCC1/' Y',' G', . C'I
DATA 5CC2/' N-T',' L',' w',' U','SUMN'I
C
C
DATA ISPIN/20,26,28,29,32,33,CfOI
C
DATA JKEEP/20-21
C
ADJUST = 2000.
NOX = ..1
LLLLL = 0
LLKKK = °
LSTRAT = 1
IZS = °
LPNTT = 0
INTT = 0
ICAN = 0
C
C
C
BEWARE '" SIC AT A TIME WILL BE DIFF.
DO 106 I = 1,1932
106 ZEROCI) = 0.0
ISBUG = 1
3
1 1 &f
1 15
1 1 6
1 1 7
1 1 8
1 1 9
120
1 2 1
122
123
12&f
125
126
127
128
129
130
1 3 1
132
133
13&f
135
136
137
138
139
1&+0
1 &f 1
1 &+ 2
1&+3
1 q q
1&+5
1&+6
1 &+ 7
Ilia
lli 9
150
-------
S~ 0096 -~EA) = 5 HS 15
ISN 0097 LWRIT = 6 RMS 152
ISN 0098 A I A = 0./897. RMS 153
ISN 0099 AAA = 960.111Lf5. RMS 15q
JSN 0100 DO £f07 J = 1, 100 RMS 155
JSN 0 1 0 1 DO £fa7 J = 1,5 Rt~S 156
ISN 0102 FUII,J) = o.u RMS 157
ISN 0103 FUIII,J) = 0.0 RMS ISB
JSN 010q Lf07 CONTINUE RMS 159
JSN 0105 RIT=232.96 RMS 160
ISN 0106 RAT = 10S.6 RMS 1 (, 1
ISN 0107 DELTRll)=IRIT..OS+RAT..OS)/683£f62.2 RMS 162
JSN 0108 DELTR(2)=IRIT..10+RATo.1S)/683Q62.2 RI15 1 (,3
ISN 0109 DELTR(3)=(RIT'.35+RAT..SO)/683~62.2 RMS 1 (, q
ISN 01 10 DELTRlq)=IRIT..qO+RAT".90)/683Q62.2 RMS 16S
JSN 01 1 1 OELTRIS)=IRIT'.10+RAT-1.00)/6A3£fb2.2 RMS 166
ISN 01 17. OI(Y = 683£f62.2 RMS 167
JSN 01 13 OFUl(1) = IRITo.OS)/OKY RMS 168
JSN 01 1 q OFUIIZ) = IRITo.IO)/OKY RMS 169
ISN 0 1 IS DFUI(3) = IRITe.35)/Ot
-------
S~ 0130 IAPROP = ( R~ 5 89
JSN 0131 LTAPE = 9 RI1S 190
JSN 0132 JYEAR I: -2 RMS 19 1
ISN 0133 IDYEAR :: 1971 RMS 192
ISN 01314 IFEED :: 0 RMS 193
ISN 0135 'iDS FORMATI20A'i) RMS 1914
ISN 0136 'il0 FORMATI20I'i) RMS 195
ISN 0137 ~'i0 FORMAT(10X,5FIO.'i) RMS 196
ISN 0138 DO £f0'i I = 1,91 RM5 197
ISN 0139 DO £fOq J = 1,20 RMS 198
ISN 01140 XXII,J) :: 0.0 RMS 199
ISN 0 1 14 1 XXHII.J) :: 0.0 RMS 200
ISN 01£f2 00 6 I K=I,5 RMS 201
ISN 01143 61 EOLOII,J,K) = 0.0 RMS 202
JSN 01q£f liQIi CONTINUE 'RI1S 203
C RMS 20'+
ISN OlliS DO li93 I = 1 tll Oli RMS 205
ISN 01£f6 £f93 OZEROII) = 0.0 RMS 206
C RMS 207
C GENERAL INPUT VARIABLES RMS 208
~ C RI1S 209
~ I SN 01li7 REAO ILREADtl78) IPOPGRII), 1=1,91) RMS 210
JSN OlliS 178 FORMAT 116F5.2) RMS 21 1
ISN 01'i9 00 166 I = 1, 9 1 RMS 212
ISN 0150 READ ILREADtl77) POPHII),TL,PTP,APVHII) RMS 213
ISN 0151 177 FORMAT lliFI0.0) RMS 2114
ISN 0152 TPHII) = POPHII).PTP/I000. RMS 215
ISN 0153 TOH(J) :: TL-TPH(I) RMS 216
ISN 01514 POPGRII) = POPGRII) / 100. RMS 217
ISN 0155 166 CONTINUE RMS 218
ISN 0156 00 5300 I = 1,20 RMS 219
JSN 0157 5300 READ 15,5305) IWAGECO( I ,J) ,J=l ,3) RMS 220
ISN 0158 5305 FORMAT (10X,JFI0.3) RMS 221
C RMS 222
ISN 0159 DO li30 I = 1,20 RMS 223
ISN 0160 READ ILTAPE,liliO) ALPHA(I),BETAII) ,O(I),PICOEFII),VKII) RMS 22'i
ISN 01 b 1 IF (ISBUG) Lf30,Lf30,li20 RMS 225
ISN 0162 £fZO WRITEILWRIT.liliO) ALPHA(I) ,BETA I I) .D(J).PICOEFII),VKII) RMS 22'ci
-------
SN 0.63
lSN 016&+
I~N OJ6~
ISN Bt66
ISN 0167
lSN 0168
ISN 0169
ISN 0170
lSN 0171
ISN 0172
ISN 0173
ISN 017&+
ISN 0175
ISN 0176
ISN 0177
lSN 0178
ISN 0179
ISN 0180
~
U1 ISN
. ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
0181
0182
018'1
0185
0186
0187
0188
0189
0190
0191
0192
0193
ISN 019'1
ISN 0195
'f30 CONTINUE
DO '460 I = 1,92
RfAOlbTAP~,ttpl)
'f~8 ~ORMAT (1Hi)
1.\59 FORMAT(IHlt1I1H
lUDY',/)
'f60 CONTINUE
'f61 FORMAT (I&f,19A'f)
WRITE (LWRIT,I.\S9)
KITTY = 0
DO 'f6'f I = 1,91
WRITE (LWRIT,'fb2) ( IAQCR(I,J) ,J=l,20)
KITTY = KITTY + 1
IF (KITTY-31) 'f6'f,'f63,'163
1.\62 FORMAT(IH ,23X,I3,ZX,19A'f)
"t63 KITTY = 0
WRITE (LWRIT,'fS8)
"t6'f CONTINUE
WRITE (LWRIT,"tS81
(IA~C;R( I ,J) ,,);;'1 ,zeJ)
,'AQCR S FOR wHICH DATA WAS AVAILABLE FOR .THIS
C
READ(LTAPE,'IIO) (KCITY(II,I=I,100)
JF(LLKKK.EQ.l) GO TO 720
LLKKK = 1
READ(LTAPE,"t90) «GAMA(IKIT,JKIT), JKIT = 1,5), IKIT = 1,191
"t90 FORMAT(lOX,5FI0.5)
DO '192 KIT = 1,20
READ(LTAPE,&f91) PICON(KJT), A(KITI, ALIT(KIT)
'f91 FORMAT(10X,3FI0.5)
'f92 CONTINUE
DO 799 KIT = 1,19
799 READ(LTAPE,798) AB(KIT)
798 FORMAT(10X,FI0.61
C
C
C
C
1970-1980 INPUT SEGMENT
DO 'ISO I = 1,91
READ(L TAPE,"t&f5) XLH( J) ,UH( J) ,SUMNH( I ,20) ,WaH( 1,20) ,XXH( J ,20),
R~ S
RMS
RM9
RMS
STRMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS.
RMS
RMS
RMS
RMS
RMS
RMS
227
226
2?9
230
231
232
233
23q
235
236
237
238
239
2'10
21.\1
2'12
2'13
2'1"t
2'15
21.\6
2'17
2'18
21.\9
250
251
252
253
25'1
255
256
257
258
259
260
261
262
263
26'1
-------
ISN 0196
ISN 0197
ISN 0198
ISN 0199
ISN 0200
ISN 0201
ISN 0202
ISN 0203
ISN 0201f
ISN 0205
ISN 0206
ISN 0207
ISN 0208
ISN 0209
ISN 0210
ISN 0211
ISN 0212
ISN 0213
oj::. ISN 0211f
0' ISN 0215
ISN 0216
ISN 0217
ISN 0218
ISN 0219
ISN 0220
ISN 0221
ISN 0223
ISN 0221f
ISN 0225
ISN 0226
ISN 0227
ISN 0228
ISN 0229
ISN 0230
Xli ,2J),YYHII),THII),G.I ,CI ~'S 265
XKHII,7.0) = XXH(I,20IlVKI20) RMS. 266
WII,20) = WBHII,20)/SUMNHII,20) RMS 267
PIHII,20)=XXHII,20)-WBHII,20) RMS 268
TAHII) = THII)-TPH(J) "TOHII) RMS 269
IF (ISBUG) 1f1f7,qq7,qqq RMS 270
qqq WRITEILWRIT,qlf6)XLHI I) ,UHI I),SUMNHII,20),WBHII,20),XXHII,20), RMS 271
1 XIHII,20),YYH(J),TH(I),GH(I),CH(!) RMS 272
qq5 FORMAT(FI0.1,FIO.3,/,3FIO.I ,IOX,FIO.l,/,3F10.0,/,FI0.0,1,80X) RMS 273
qq6 FORMAT(IH ,FI0.l,IX,FI0.3,lX,IfIFI0.1,IX),Q(FI0.0,IX)) RMS Z7q
qlf7 XNTHI I) = I 1.0"UHI I) ).XLHI I) RMS 275
XNBAHII) = XNTHII)-SUMNHII,20) RMS 276
QSO CONTINUE RMS 277
DO Q70 I = 1,91 RMS 278
READILTAPE,lf66) QTHII) ,QCHII),QMH(I) ,GBARHII) RMS 279
WRITEILWRIT,1f66) QTHII),QCH(I),QMH(I),QBARHII) RMS 280
FORI-1ATIIH ,20IS) RMS 281
CONTINUE RMS 2A7.
FORMATlIfFI0.0) RMS 283
READILTAPE,599Q) XKIT,KIND RMS 281f
FORMAT (FI0.l,66X,I2) RMS 285
IF(KIND..99) 5993,S996,5996 RMS 286
5996 READILTAPE,5997) XKIT,XKIT1,XKIT2,XKIT3,ICT,KIND RMS 287
5997 FORMATI3rl0.1,lOX,Fln.l,7.3X,13,I2) RMS 288
IF (KIND-99) 5998,5999,5999 RMS 2A9
5998 KIND = KIND-20+120/KIND) RMS 290
ICITY = KCITY(ICT) RMS 291
IF (ICITY.EQ.O) GO TO 5996 RMS 292
WBHIICITY,KIND) = XKITI RMS 293
SUMNH(ICITY,KIND) = XKIT RMS 29Q
XXH(ICITY,KIND) = XKIT2 RMS 295
XIHIICITY,KINO) = XKIT3 RMS 296
W(ICITy,KIND) = WBHIICITY,KIND)/SUMNHIICITy,KIND) RMS 297
XKH(ICITY,KIND) = XXHIICITY,KIND)/VK(KIND) RMS 298
PIHIICITY,KINO)=XXHIICITY,KINO)-WBH(ICITY,KIND) RM5 299
GO TO 5996 RMS 300
C RMS 301
-------
ISN 0231
ISN 0232
ISN 0233
ISN 023~
ISN 0235
ISN 0236
ISN 0237
ISN 0238
ISN 0239
ISN 02~0
ISN 02~1
~ ISN 02~7.
-.)" I SN 02~3
ISN 02~~
ISN 02'f5
ISN 02'f6
ISN 02'f7
ISN 02'f8
ISN 02'f9
ISN 0250
ISN 0251
ISN 0252
ISN 0253
ISN 025'f
ISN 0255
ISN 0256
ISN 07.57
ISN 0258
.....
::X"INCT
.....
C
C
C
C
C
C
C...........
C..... EXTINCT
C
5999
1967 INPUT SECTION
.....
~MS
RMS
RMS
RM5
RMS
RMS
RMS
RMS
RMS
RMS
,"KIND,IXII),YII),I=3,5)RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
Rr1S
RMS
RMS
RMS
RMS
RMS
RMS
CONTINUE
REAOILTAPE,6000) CT,IX(I),Y(I) ,I=I,2),ICT
IF (KIND-99) 6002,6006,6006
6002 KIND = KIND - 20 + (20 / KIND)
ICITY = KCITY(ICT)
DO 6003 I = I,S
EOLDIICITY,KIND,I) = XII)
6000 )FORMATIFI0.0.1nX,'fFI0.0.13X,13,I2,I,QFI0.0.FI0.I,FI0.0
6003 CONTINUE
GO TO 5999
6006 CONTINUE
DO 2350 ICK = 1,91
DO 23'f5 JNK = 1,19
KG = JKEEP(JNf()
DO 23'f0 KNK = 1.5
KGO = Lf(EEP(KG,KNK)
IF (KGO) 23~5,23'tS.2310
2310 IF (EOLD(ICK,JNK,KGO» 23't0.232S,23QO
2325 EOLDIICK.JNK.KGO) = .001
23't0 CONTINUE
23'f5 CONTINUE
2350 CONTINUE
C SAVE ORIGINAL DATA FOR LATER USE
DO 800 I = 1,91
OW(II = WBH(I.201/SUMNH(I,20)
OWB(I) = wBH(I,20)
OX(I) = XXH(I.20)
01(1) = XIH(I,20)
OYII) = YYH(I)
303
30Q
30;
306
307
308
309
310
3 1 1
312
313
31 't
315
316
317
318
319
37.0
321
322
323
32't
325
326
327
328
329
330
331
332
333
33'f
335
336
337
338
339
3't0
-------
'SIJ 0259
ISN 0260
ISN 0261
ISN 0262
ISN 0263
ISN 026tt
ISN 0265
lSN 0266
ISN 0267
ISN 0268
ISN 0269
ISN 0270
lSN 0271
ISN 0272
ISN 0273
ISN 027tt
ISN 0275
ISN 0276
ISN 0277
ISN 0278
H::- ISN 0279
00 ISN 0280
ISN 0281
ISN 0282
ISN 0283
ISN 02att
ISN 0285
ISN 0286
ISN 0287
ISN 0288
ISN 0289
ISN 0290
ISN 0291
ISN 0292
ISN 0293
ISN 029tt
ISN 0295
ISN 0296
OT(I) = TH( )
OG(I) = GH(I)
0(1) = (H(I)
OXBAR(I) = XNBAH(I)
OQ(I) = Q(H(I)
OQM(I) = QMH(I)
OQB(I) = QBARH(I)
OTP(J) = TPH(I)
OTA(I) = TAH(I)
OTO(1) c TOH(I)
OPOP(I) = POPH(I)
OAPV(I) = APVH(I)
OW(92) = OW(92) + OW (I)
OWB(92) = OWB(92) + OWB(I)
OX(92) = OX(92) + OX (I)
01(92) = 01(92) + 01(1)
OY(92) = OV(92) + OV(I)
OT(92) = OT(92) + OT(I)
OTP(92) = OTP(92) + OTP(I)
OTA(92) = OTA(92) + OTA(I)
OTO(92) = OTO(92) + OTO(I)
OPOP(92) = OPOP(92) + OPOP(I)
OAPV(92) = OAPV(92) + OAPv(l)
OG(92) = OG(92) + OG(I)
0(92) = 0(92) + 0(1)
OXBAR(92) = OXBAR(92) + OXBAR(I)
OQ(92) = OQ(92) + OQ(I)
OQM(92) = OQM(92) + OQM( I)
OQB(92) = OQB(92) + OQR(I)
800 (ONTINUE
PERUN (92) = .O'f
00 201 I = 1,7
READ (LTAPE,202) K,lK
K = K - 20&(20/K)
201 JKEEP (K) = IK
DO 7008 I = J ,20
HEAD (LREAD,7007) (QINV(I,J),J=I,3)
7008 CONTINUE
~MS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
Rr-lS
RMS
RMS
RMS
RMS
RHS
RMS
RMS
RMS
RMS
Rr-lS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
Rr-IS
RMS
RMS
RMS
RMS
3'fl
3'f2
34)
3'f'f
3'f5
3'+6
3'+7
348
3'19
350
351
352
353
35tt
355
356
357
358
359
360
361
362
363
36'f
365
366
367
368
369
370
371
372
373
37'f
375
376
377
378
-------
S I 0297
ISN 0298
ISN 0299
ISN 0301
ISN 0302
ISN 0303
ISN 030'f
ISN 0305
ISN 0306
ISN 0307
ISN 0308
ISN 0309
lSN 0310
lSN 0311
ISN 0312
ISN 0313
ISN 031'f
ISN 0315
ISN 0316
ISN 0317
of:. ISN 0318
-.!)
lSN 0319
ISN 0320
ISN 0321
ISN 0323
ISN 032'f
ISN 0325
lSN 0327
ISN 0328
ISN 0329
ISN 0330
lSN 0331
ISN 0332
ISN 0333
ISN 0335
7007 FO~MAT(10X,3F 0.3
202 FORMAT (213)
IF(LLLLL.EQ.l) GO TO 38
LLLLL = 1
DO 5991 ICT = 1,91
1 = ICT
PERUN (J) = .O'f
00 5990 MJ=ltl9
5990 EOLDH(ICT,MJ) = EOLO(JCT,MJ,5)
EOLDH(ICT,20) = QTH(ICT)
REAO(LTAPE,S992) (B(I,J), J = 1,19)
5991 CONTINUE
00 298 I = 1,91
READ (LTAPE,299) RNET(I)
298 CONTINUE
DO 'f01 I = 1,8
READ(LTAPE./i05)(LABEL(I,J), J = 1,20)
liOI CONTINUE
299 FORMAT(ZOX,FI0.3)
38 CONTINUE
LPNTI = 0
R~S
RMS
RMS
RM5
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
Rf~S
RMS
RM5
RMS
Rr1S
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
C
C
C
7009
GROWTH FACTOR CALCULATIONS FOR 1970-1980 ONLY
REAO(LREAD,1070) INY,(IYEAR(JRK).JRK=I,INY)
WRITE(6,1315)
IF( ICAN,EQ,O) GO TO 6698
DO 'f08 I = 1.100
MYCITY = KCITY(I)
IF (MYCITy ,EQ,O) GO TO 'foe
INRITE(6,'f09) (FU( I ,J) ,J=1 ,5),I
WRITE(7,/i09) (FU(I,J),J=I,5),I
'f09 FORMAT(IH ,F9,3,'fFI0,3,26X,I/i)
WRITE(6.'f091 (FUI(I,J),J=I,5ItJ
WRITE(7,Lt09) (FUI(I.JI,J=I,5),I
CALL OLETS(MYCITY)
IF (PERXX,EQ,O,O) GO TO 'foe
OIVINV = FUIII,5)/PERXX
379
3BO
381
382
383
38'f
385
3Bb
387
3118
389
390
391
392
393
39'f
395
396
397
398
399
'fDO
'fOI
'f02
'f03
'f0'f
LtOS
'fOb
Lt07
'f08
'f09
LtJO
LI t 1
Ll12
Lll 3
Lll/i
/its
Lt 1'6'
-------
S~ 0336 D VTOT=CFU ,5)+FUICI,5 )/:I~RXX ~MS '17
ISN 0337 I F CPERYY.EQ.O.O) GO TO ~08 RMS tf I 8
ISN 0339 YYINV=FUICI,5)/PERYY RMS tfl9
ISN 03'fO YYTOT=CFUCI,5)+FUICI,5»/PERYY RMS tf20
ISN 03'f1 WRITECb,20) DIVTNV,DIVTOT,YYINV,YYTOT,T RMS 'f /' 1
ISN 03'f2 WRITEC7,20) DIVINV,OIVTOT,YYINV,YYTQT,I RM5 'i22
ISN 03tf3 20 FORt1ATIIH F9.'i,3F10.tf,36X,I'i) RMS tf23
I5N 03tf'i tfOB CONTINUE RMS 'f2'i
ISN 03'i5 6698 lCAN = 1 RMS tf25
lSN 03tf6 WRITECLWRIT,1070) lNy,C1YEARCJRK),JRK=I,INY) RMS '176
ISN 03'f7 I M I I = PEAR C I NY> - 1968 - 110(5 RMS 'i27
ISN 03'fa 1070 FORMATC 1'f,lX,J2IS) RMS tf?8
lSN 03'f9 IF (INV.EQ.O) STOP RMS 'f29
ISN 0351 DO 1020 JRK = 1,Stf RMS tf30
ISN 0352 READ ILTAPE,7010) CGRFACCJRK,KRK,I) ,KRK=1,20) RMS 'f 3 I
ISN 0353 7010 FORMAT CI'fF5.0,I,6F5.0) RMS tf32
lSN 0351f 7015 FORMATC1'fF5.3,1,6F5.3) RMS tf.33
lSN 0355 1020 CONTINUE RMS tf3'f
ISN 0356 DO 1025 JRK = 1,5'1 RMS 'f35
ISN 0357 DO 1025 KRK = 1,20 RMS 'i36
\J1 lSN 0358 GRFACIJRK,KRK,l) = GRFACIJRK,KRK,I)/100. RMS '137
o ISN 0359 1025 CONTINUE RMS '138
ISN 0360 00 1030 JRK = 1,5tf RMS 'f39
ISN 0361 REAOILTAPE,701S) (GRFACCJRK,KRK,2),KRK=1,20) RMS ''fO
ISN 0362 1030 CONTINUE RMS 'I'll
ISN 0363 DO 10'f0 JRK = 1,5 tf Rt1S 'f'l2
ISN 036tf READCLTAPE,7015) CGRFACIJRK ,KRK,31 ,KRK=l ,20) RM5 'f'l3
ISN 0365 10'10 CONTINUE RMS 'q'l
lSN 0366 00 1060 JRK = 1,5tf Rrv1S 'iLlS
ISN 0367 READILTAPE,1050)CXINCCJRK,KRK) ,KRK=l ,3) RM5 tf'f6
ISN 0368 1050 FORMATC10X,FI0.0,1,10X,FIO.0,1,10X,F10.0) RMS 'fLl7
ISN 0369 1060 CONTINUE Rt1S ''f8
( RMS tf'f9
lSN 0370 7533 FORMATCtfFIO.2) RMS 'f50
C RMS '151
C RMS '152
C RMS '153
C FOR AQCRS NOT INCLUDED IN DBE REPORT CCDDE GREATER THAN SIf ) RMS 'Sq
-------
ISN 0371
ISN 0372
ISN 0373
ISN 037'+
ISN 0375
ISN 0376
ISN 0377
ISN 0378
ISN 0380
ISN 0381
ISN 0382
ISN 0383
ISN 038'+
U'1
..... IS N
ISN
ISN
ISN
0385
0386
0387
0388
ISN 0389
ISN 0391
ISN 0392
ISN 0393
ISN 0395
ISN 039/,
ISN 0397
ISN 0398
C
C
C
C
C
1051
1310
C
C
1052
1053
10SLi
C
1055
c
721 1
C
C
C
7225
C
READ(LTAPE,1051) YGROW
FORMAT(30X,FlO.2,,+OX)
FORMAT(IH ,IJ,12F9.1,15)
DO 105'f JRK = 1,20
DO 105'f KRK = 1 t3
DO 1053 KC = 55,100
LCITY = KCITY(KC)
IF(LCITY.EQ.O) GO TO
GRFAC(LCITy,JRK,KRK)
GO TO 1053
GRFAC(LCITY,JRK,KRK)
CONTINUE
CONTINUE
1052
= GRFAC(32,JRK,KRK)
= 0.0
DO lOSS IJ=It38
WRITE(LWRIT,l070) IJ
YGROWIIJ) = YGROW(IJ)1100.
REA I) ( L REA D , I 080) I N C , ( I C H 0 ( J R K ) ,J R K = 1 , I N C )
IF(IA312.EQ.O) GO TO 7213
READ(LREAD,7211) BENGOV,BENEF,THISl,P
FORMATI,+FlO.2I
...
OVERWRITES FEEDBACK MATRIX IF BENEFIT NOT EQ. $10 BILLION
IF(BENGOV.EQ.I0000.) GO TO 7213
DO 7225 ICT = 1,91
READ(LREAD,5992)(B(ICT,J) ,J=l,19)
WRITE(LWRIT,5992) (B( ICT,J) ,J=1 ,19)
CONTINUE
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RM5
RMS
RMS
RMS
RI1S
RMS
RI15
RMS
RMS
RMS
RM5
RMS
RMS
RMS
RMS
RMS
RMS
RMS
..RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
j 55
'+56
/.f57
'fS8
'fS9
'+60
Li61
Li62
'+63
Li6,+
'+65
'+66
Lib]
'+68
'+69
'+70
Li 7 1
'+72
'+73
'f7,+
'+75
'+76
Li77
'+78
'f79
'+80
LiBl
'f82
'fB.3
'f8'+
Li8S
Li86
LiB7
'+88
Li89
'+90
'+91
'+92
-------
'5' 0399
ISN 01400
ISN 0~01
ISN 0~02
ISN O~03
ISN o~o~
I5N O~05
ISN O~06
ISN O~07
ISN O~08
ISN O~09
ISN O~10
ISN 0~11
ISN 0~17.
U1 ISN 0~13
N
ISN 0~1~
ISN O~15
ISN 0~16
ISN O~17
ISN O~18
ISN 0~19
ISN 0~21
ISN 0~22
ISN O~23
ISN O~2Lf
ISN OLf25
ISN OLf26
ISN OLf27
7213
C
C
C
C
C
C
CONTI\JUE
ELIMINATES REQUESTS FOR INFO. ABOUT CITIES FOR WHICH NO DATA IS
AVAILABLE (EXAMPLE ~ NO DATA FOR CITY NO. 13)
NIL = INC
DO 1063 I=ltlNC
106Lf KIT = ICHD(I)
lAB = KCITY(KIT)
IF( lAB) 1063,\n61.1063
1061 NIL = NIL ~ 1
DO 1062 NIPP = I,NIL
ICHO(NIPP) = ICHO(NIPP + 'I)
1062 CONTINUE
WRITE(LWRIT,IOb5) KIT
GO TO 106Lf
1063 CONTINUE
INC = NIL
FORMA T (/'
1065
C
C
1080
1100
1 1 10
C
C
1355
NO DATA AVAILABLE FOR CITY NO. ',12)
FORMAT(1615)
DO 1100 JG = l,INY
IYEARIJG) = IYEARCJG)~1969
JYKNT = 0
CONTINUE
IFCJYKNT.EQ.O) GO TO 1115
U(92) = PERUN(92)
UN(92) = ~XNN(92)/XL(92)
CALL PRINI (92)
006699 I = 1,91
XXCI,20)=XXCI,20)+XDNII)
PICI,20) = PICI,20) + PH./CI)
XKCI,20) = XKCI,20) + XKNCI)
R'S
R"15
RMS
RMS
RMS
RHS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
Rf~S
RMS
RMS
RMS
RMS
RMS
RMS
RMS
. RMS
RMS
RMS
RMS
~93
Lf9Lf
Lf95
Lf96
Lf97
Lf98
Lf99
500
501
507.
SO,)
50Lf
505
506
507
508
509
510
5 1 1
512
513
51/.f
515
516
517
518
519
520
521
522
523
52Lf
525
526
527
528
529
53'U'
-------
ISN 01i28
ISN 0&f29
ISN 0'+30
lSN 0'+31
ISN 0'+32
ISN 0'+33
ISN Q&f3&f
ISN Q&f35
ISN Oli36
ISN 01i37
ISN 01i38
ISN 0&f39
ISN O&f&fQ
ISN O&f&fl
ISN 0&f&+2
ISN 0&+&+3
lSN O'+"t"t
ISN 0'+"t5
ISN 0&f"t6
ISN 0'+&f7
! U1 ISN 0'+'+8
w ISN 0'+'+9
ISN 0&f50
ISN 0&f51
lSN 0'+52
ISN 0'+53
ISN 0&f5'+
ISN 0'+55
ISN 0'+56
ISN 0'+57
ISN 0"t58
ISN 0&fS9
ISN 0"t60
ISN 0&f61
lSN 0'+62
ISN 0'+63
ISN 0'+6';
ISN 0"t65
6699
1220
XC ,201 = X CI,20 + X NCI)
rYCI) = YYCI) + YYNCII
GCI) = G(I ) + GNII)
TCI) = TI1) + TN II)
TPIII = TPCII + TPNCII
TACI) = TAlI I + TANIII
TOCI) = TOCI) + TONCI)
(CI) = (CI) + (NIl)
XNBARCII = XNBARCI) + DNBARNCII
XNTCII = XNTC II + XNNCII
XU I) = XLIII + XLNCII
WCI,201 = WCl,20) + WNCII
SUMNCI,2(l) = SUMNCI,201. + SUMNNCJ)
WBCI,201 = XXCI,201-PICI ,20)
QBARCII = QBARCII + QBARNCII
QMCII = QMCII + QMNCII
QCCII = Q(II) + QCNII)
QTCII = QTCII + QTNCII
CONTINUE
DO 1225 1=1,92
DO 1220 J=1,20
XXHCI,J) = XXCI,JI
~8HCI,JI = WBCI,JI
SUMNHCI,J) = SUMNCI,J)
X I H C I , J I = X I C I , J )
EOLDHCI,JI = EOLDCI,J,SI
XKH(I,J) = XK(l,JI
PIH II,J) = PICI,J)
CONTINUE
YYHCI) = YYCIJ
XNBAHCl) = XNBARCI)
THCII = TCII
TPHCI) = TPCI)
TOHII) = TOCI)
TAHCJ) = Till I)
POPHCI) = POP(I)
APVHCI) = APV(I)
GH(II = G(I)
R~ 5
RMS
RMS
RMS
RMS
RM5
Rt'1S
RM5
RM5
RMS
RMS
RMS
RM5
RM5
RMS
RMS
RMS
RMS
RMS
RM5
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
Rt>IS
RMS
RMS
RMS
RMS
531
532
533
53'f
S3S
536
537
538
539
SliD
5'+ 1
S&f2
5'f3
5'f'+
S'f5
5'+6
5'+7
5'+8
5'f9
550
551
55l
553
55"t
555
556
557
558
559
5'-0
561
562
563
S6'+
56S
566
567
51,8
-------
1--- --
15' Off66 C. ( ) = C ( ) H5 569
15N 0'+67 QBARH( I) = QBAR(I) RMS 570
ISN 01461\ QCH ( I) = QC ( I ) RMS 571
ISN 0'+69 (~MH(I) = QM (I) RMS 572
ISN 0'+70 CHH(I) = Q T ( I ) RM5 573
I5N 0'+71 1225 CONTINuE RMS 571f
C RMS 575
C RMS 576
C RMS 577
C RMS 578
I5N 0'+72 DO 6106 I = 1,1932 RMS 579
ISN Off73 6106 ZERO(I) = 0.0 RMS 5RO
ISN 0'+7'+ 1 1 15 CONTINUE RMS 581
ISN 0'+75 I F (JYKNT.EQ.INY) GO TO 7009 RMS 5RZ
ISN 0'+77 JYKNT = JYKNT + 1 RMS 5R3
ISN 0'+78 JYEAR = IYEAR(JYKNT) RMS 5R'+
C RMS 5R5
C ... TO GENERATE oe5 SET TEST TO 8 - 1977 - NOW 1978 ... Rr~s 586
c RMS 587
ISN 0'+79 IF(JYEAR.EQ.IMII) GO TO 95 RMS 588
\.11 C RM5 5R9
*" I SN 0'+81 KYY = JYEAR+1969 RHS 590
ISN OffB2 XX(92,20) = 0.0 RMS 591
ISN 0'+83 P I ( 9,2 , 20) = 0.0 RMS 592
ISN O'+B'+ XI(92,20) = 0.0 RHS 593
ISN 0'+85 XK(92,20) = 0.0 RMS 59'+
ISN 0'+86 SUMN(92,20) = 0.0 RMS 595
ISN 0'+87 XN8AR(92) = 0.0 Rt-1S 596
ISN 0'+88 YY(92) = 0.0 Rr~s 597
ISN 0'+B9 XNT(92) = 0.0 RMS 598
ISN 0'+90 C ( 92) = 0.0 RMS 599
ISN 0'+91 XL(92) = 0.0 RMS 600
ISN 0'+92 QT(92) = 0.0 RMS 601
ISN 0'+93 QM(92) = 0.0 RMS 602
ISN O1f9'+ Q8AR(92) = 0.0 RMS 603
ISN 01f95 QC(92) = 0.0 RMS 60'+
ISN 0'+96 G(92) = 0.0 RMS 605
ISN 0'+97 T(92) = 0.0 RM5 600
-------
ISN 0'i9~
ISN 0'i99
ISN 0500
ISN 0501
ISN 0502
ISN 0503
ISN 0505
ISN 0507
ISN 050A
ISN 0509
ISN 0510
ISN 0511
ISN 0517-
ISN 0513
ISN 0515
ISN 0517
ISN 0518
ISN 0519
ISN 0520
ISN 0521
U'1 ISN 0522
U'1 ISN 0523
ISN 0521f
ISN 0525
ISN 0526
ISN 0527
ISN 052~
ISN 0530
ISN 0532
ISN 0531f
ISN 0535
ISN 0536
ISN 0537
ISN 0538
ISN 0539
ISN 05'iO
ISN 05'il
ISN 05&+Z
.P(92) = (.0
TO(92) = 0.0
TA192) = 0.0
POP(92) = 0.0
APV(9Z) = 0.0
IF (JYKNT.EQ.l) IEX=JYEAR+3
IF(JYKNT.NE.l) lEX = 1
IKX = JyEAR+2
Ir~x = lEX
DO 310 MPCITY = l,INC
lEX = IRX
KIT = ICHO(MPCITY)
LCITY = KCITY(KIT)
IF (LCITY.GT.S3) GO TO ll'iZ
IF (XINC(LCITY,I).EQ.O.O) GO TO Il'iif
YRAG1=XINC(LCITY,I)/YYH(LCITY)
YRAG2 = XINCILCITY,2J/XINC(LCITY,I)
YRAG3 = XINC(LCITy,3)/XINCCLCITY,2)
YRAMI = ALOGI0(YRAGIJ/3.
YRALI = 10...YRAMI
YRAM2 = ALOG10CYRAG2J/5.
YRAL2 = 10...YRAM2
YRAM3 = ALOGI0(YRAG3)/5.
YRAL3 = 10...YRAM3
GO TO 11'i3
11'i'i RATE = 0.0
IF (KIT.EQ.23) RATE = .0762
IF CKIT.EQ.32) RATE = .0758
IF (KIT.EQ.33) RATE = .06'i5
YCATCH = (1.0+RATE).-IKX
YRON = 1.0
YRAGI = 1.0
YRAG2 = 1.0
YRAG3 = 1.0
GO TO ) l'i3
11'iZ IKX = JYEAR+Z
YCATCH = (1.0+YGROW(LCITY-53J)..IKX
YRON = 1.0
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
Rr..,s
RMS
RMS
RMS
RW5
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
607
60~
609
610
6 1 1
612
613
6\q
615
616
617
618
619
620
621
62?-
623
62'i
625
626
6?7
628
629
630
631
632
633
63'i
635
636
637
638
639
6lfO
6 &+ 1
6&+Z
6&+3
61f'i
-------
5 I OS'i:
ISN OStftf
ISN OSIfS
ISN 05tf6
ISN OStf7
ISN OStf8
ISN 05S0
ISN 05S2
ISN OSStf
ISN OSS6
ISN OS58
ISN 05S9
ISN OS61
ISN 0562
ISN 0563
ISN 056£+
ISN 056S
ISN 0566
ISN 0567
ISN 0568
U1 ISN 0569
0'
ISN OS70
ISN 0571
ISN 0572
ISN 0573
ISN 057£+
ISN 0575
ISN 0576
ISN 0578
ISN 0579
ISN OS80
ISN OS81
ISN 0582
ISN 0583
ISN 058£+
Y~AG = 1.0
YRAG2 = 1.0
YRAG3 :: 1.0
I1tf3 00 309 JGE = 1,20
GRON = 1.0
IF (XXH(LCITY,JGE).EQ.O.O) GO TO 309
IF (GRFACILCITY,JGE,1 ).EQ,O.O) GO TO 309
IF (JY~NT.NE.J) GO TO 1237
IF (JYKNT.EQ.l,ANO.JYEAR.EQ.l) GO TO 1237
IF(GRFACILCITy,JGE,1 I.GT.O,O) AF1=ALOGI0IGRFAC(LCITY,JGE,I))/3.
RAFI = 11).-.AFI
IFIGRFACILCITY,JGE,I).GT.O,0.AND.GRFACILCITY,JGE,2,I,GT,O,O)
1 AF2=(ALOGI0IGRFAC(LCITY,JGE,2))-ALOGI0(GRFACILCITY,JGE,1)))/5,
RAF2 = 10.-.AF2
ZYR = JYEAR-2
XCATCH = (RAF1.-3,).IRAF2.-ZYR)
XXHILCITY,JGE) = XXHILCITy,JGE).XCATCH
WBHILCITY,JGE) = WBHILCITy,JGE).XCATCH
PIHILCITY,JGE)=PIHILCITY,JGEI-XCATCH
XIHILCITY,JGEI = XIHILCITY,JGE)-XCATCH
EOLDHILCITY,JGE) = EOLDH(LCITY,JGE)*XCATCH
1237 GO TO (1120,1130,1130,1130,1130,1130,11£+0,11£+0,11'10,11,+0,11£+0),
1 JYEAR
1120 GRON = GRFAC(LCITY,JGE,I)
YRON = YRAL1..3.
KJNDX = 0
LJNDX = 0
GO TO 1150
1131 FORMAT 1/,1H ,15,5FI0.3)
1130 IF (GRFACILCITY,JGE,2)-GRFAC(LCITY,JGE,I).LE.0.0)GO TO 1133
GRONI = (ALOG11)(GRFAC(LCITY,JGE,?)I-ALOGI0(GRFAC(LCITY,
1 JGE, 1 ) I ) 15. .
GRON = 10,QoGRONI
YRON = yRAL2
KJNOX = JYEAR-2
LJNDX = 0
GO TO 1150
1133 WRITE(6,1131) LCITY,GRFAC(LCITY,JGE,2),GRFAC(LCITY,JGE,I)
~ I 5
RMS
RMS
RMS
RNS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
Rr~s .
Rr-1S
RMS
RMS
RMS
Rr-1S
RMS
RM5
RMS
RMS
RMS
RMS
RMS
RMS
RMS
6lf5
6lfb
6Lf7
6£18
6£19
650
651
652
653
6Stf
65S
656
657
658
659
660
661
662
663
66Lf
665
666
667
668
669
670
671
672
673
67Lf
675
676
677
678
679
680
681
682
-------
SIJ 0585
ISN 0586
ISN 0587
ISN 0589
ISN 0590
ISN 0591
ISN 0597.
ISN 0593
ISN 0591f
ISN 0595
ISN 0596
ISN 0597
ISN 0599
ISN 0600
ISN 06UI
ISN 0602
ISN 0603
ISN 0605
U1 ISN 0606
I -J. ISN 0607
I ISN 0608
ISN 0609
ISN 0610
ISN 0611
ISN 0612
ISN 0613
ISN 0611f
ISN 0615
ISN 0617
ISN 0619
ISN 0620
ISN 0621
ISN 0622
ISN 0621f
ISN 0625
ISN 0627
YRON = CYRAG2-1.: )/5.0+.
GO TO 1150
I1lfO IF CGRFACCLCITY,JGE,31-GRFACCLCITY,JGE,2).LE.Q.0)GO TO I1lfl
GRONI = CALOGI0CGRFACCLCITY,JGE,3))-ALOGI0CGRFACCLCITY,
1 JGE, 2) ) ) /5.
GRON = lo..-GRONl
YRON = yRAL3
KJNDX = 5
LJNOX = JYEAR-7
GO TO 1150
I1lfl WRITEC6,1131) LCITy,GRFACILCITY,JGE,3),GRFACILCITY,JGE,2)
YRON = IYRAG3-1.0)/5.0+1.
1150 IF CGRON.EQ.o.O) GO TO 309
1160 XXILCITY,JGE) = XXHCLCITY,JGE)..GRON
VFAC = GRON
00 1161 JS = 1,IEX
1161 WCLCITy,JGE) = WILCITY,JGE) . 1.018
IFISUMNHCLCITY,JGE).EQ.O.o.OR.WILCITY,JGE).EQ.O.O)
1 WRITECLWRIT,1770) KIT,JGE
1770 FORMATIIHI ,If IS)
PIILCITY,JGE) = PIH(LCITY,JGE)-VFAC
WB(LCITY,JGE) = XXCLCITY,JGE)-PIILCITY,JGE)
XICLCITY,JGE) = XIHCLCITY,JGE).VFAC
XKILCITY,JGE) = XXCLCITY,JGE)/VKCJGE)
EOLOCLCITy,JGE,S) = EOLDHCLCITy,JGE) . VFAC
XXTCLCITY,JGE) = XXHCLCITy,JGE)
XKTCLCITY,JGE) = XKH(LCITY,JGE)
SUMNCLCITY,JGE) = WBCLCITy,JGE) / WCLCITY,JGE)
309 CONTINUE
IF CJYKNT.EQ.I.AND.JYEAR.EQd) GO TO 395
IF CJYKNT.NE.I) GO TO 39S
PCATCH = Il.+POPGRCLCITY))o3.
. PCATCH = PCATCH+Cl.+POPGRCLCITYI ).KJNOX
PCATCH = PCATCH + 11.+POPGRILCITY)).LJNOX
IF CLCITY.GT.S3.0R.XINCCLCITY,1 ).EQ.O.o) GO TO 1151
YCATCH =CYRALl0.3)ol.1qlf
IFCKJNOX.EQ.O.o) GO TO 1151
YCATCH = YCATCH.YRAL2..(KJNDX)
~ ' 5
RMS
RMS
RI-1S
RMS
RMS
RMS
RW:;
RMS
Rr-1S
RMS
RMS
RMS
RMS
Rt"S
RMS
RMS
RMS
RMS
RMS
Rt"S
RMS
RMS
RMS
RMS
Rr1S
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
683
681f
685
686
687
68A
689
690
691
692
693
69q
695
696
697
698
699
700
701
707-
703
70q
70S
706
707
708
709
710
71 1
712
713
71 q
715
716
717
718
719
720
-------
151/ 0628
ISN 0630
ISN 0631
ISN 0632
ISN 0633
ISN 063~
ISN 0635
ISN 0636
ISN 0637
ISN 0638
ISN 0639
ISN 06~0
ISN 06~1
ISN 06~2
[SN 06~3
[SN 06~~
ISN 06~5
ISN 06~6
UlISN 06~8
(X) ISN 0650
ISN 0651
ISN 0653
ISN 0655
[SN 0657
ISN 0658
ISN 0659
ISN 0660
[SN 0661
ISN 0662
[SN 0663
[SN 066~
ISN 0665
ISN 0666
ISN 0667
ISN 0668
IFCLJ~)X.EQ.O.) GO TO 51
YCATCH = YCATCHeYRAL3..LJNDX
1151 YYH(LCITY) = YYHCLCITY).YCATCH
IRX = lEX
lEX = 1
XNBAH(LCITY) = XNBAH(LCITY)~(YYH(LCITY)-XXHCLCITY,20»/
1 (CYYHCLCITY)/YCATCH)-(XXH(LCITY,20)/XCATCH»
THCLCITY) = TH(LCITY).YCATCH
POPH(LCITy) = POPHCLCITY).PCATCH
TPH(LCITY) = TPH(LCIiY)eYCATCH
APVH(LCITY) = APVHILCITY) ~ YCATCH
TOHILCITY) = TOHCLCITY) . YCATCH
TAHILCITY) = TAH(LCITY) . YCATCH
GHCLCITY) = GHCLCITy).eYCATCH
CHCLCITY) = CHILCITY).YCATCH
QMHILCITY) = QMHILCITY)eXCATCH
QCHILCITY) = QCH(LCITY)~YCATCH
QBARHILCITY) = QBARHILCITY)e(YYHCLCITY)-XXHILCITY,20»/
II IYYHILCITY)/YCATCH)-(XXHILCITY,20)/X~ATCH»
395 IF(LCITY.GT.53) GO TO 120~
IF(XINCILCITY,l).NE.O.O) GO TO 1205
RATE = 0.0
IF(KIT.EQ.23) RATE = .0762
IFIKIT.EQ.32) RATE = .0758
IF(KIT.EQ.33) RATE = .06~5
YY(LCITY) = YYH(LCITY)
DO 1203 JS=I,IEX
1203 YY(LCITY) = YYILCITY) . 11.0 + RATE)
GO TO 1207
120~ YYCLCITY) =
DO 1206 LLJ
1206 YY(LCITY) =
GO TO 1207
1205 CONTINUE
YY(LCITY) = YYHILCITY).YRON
1207 YFAC = YY(LCITY) /YYHILCITY)
XNBAR(LCITY) = XNBAH(LCITY).(YY(LCITY)-XXILCITY,20»/
1 CYYHILCITY)-XXHILCITY,20»
~MS
Rr1S
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
Rr>1S
RMS
RMS
RMS
RMS
Rr1S
RMS
RMS
RMS
RMS
RMS
YYH(LCITYI
= It1EX
IYGROWILCITY-53) + 1.0) . YY(LCITY)
72
722
723
72~
725
726
727
728
729
730
731
732
733
73~
735
736
737
738
739
7~0
7 ~ 1
7~2
7~3
7~~
7~5
7~6
7~7
7Lf8
7Lf9
750
751
752
753
75~
755
756
757
758
-------
SN 0669
ISN 0670
ISN 0671
ISN 0672
ISN 0673
ISN 067lf
ISN 0675
ISN 0676
ISN 0677
ISN 0678
ISN 0679
ISN 0680
ISN 0681
ISN 0682
ISN 0683
ISN 068lf
ISN 0685
\J1
'"
ISN 0686
ISN 0688
ISN 0689
ISN 0690
ISN 0691
ISN 0692
ISN 0693
ISN 0691f
ISN 0695
ISN 0696
ISN 0697
ISN 0698
ISN 0700
310
1 2 1 2
C
C
C
C
C
"(LC TY) = "H(LCITY).Y~AC
POP(LCITY ) = POPH(LCITY) .(I.+POPGR(LCITY»
TP(LCITY) = TPHILCITY) . YFAC
APV(LCITY) = APVH(LCITY) . YFAC
TO(LCITY) = TOH(LCITY) . YFAC
TA(LCITY) = TAH(LCITY) . YFAC
G(LCITy) = GH(LCITY).T(LCITY)!TH(LCITY)
C(LCITY) = CH(LCITY)~YFAC
XNT(LCITY) = XNBAR(LCITY)+SUMN(LCITY,20)
U(LCITY) = .Olf
XL(LCITY) = XNT(LCITY)!(I.-U(LCITY»
QM(LCITY) = QMH(LCITY)-VFAC
QBAR(LCITY) = QBARH(LCITY).(YY(LCITY)-XX(LCITY,20»!
1 (YYH(LCITY)-XXH(LCITY,20»
QC(LCITY) = QCH(LCITY)-C(LCITY)/CH(LCITY)
QT(LCITY) = QC(LCITY)+QBAR(LCITY)+QM(LCITY)
CONTINUE
FORMAT(FIO.3,60X,I5,I31
DATA TABLES 1970-1780
RMS
RMS
RM5
RMS
RMc;
RMS
Rr-1S
RMS
RMS
RMS
RM5
RMS
RMS
RMS
RMS
Rr-1S
RMS
RMS
RMS
RMS
Rr-1S
Rr-1S
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
EMP -NOW INRMS
Q-NOW AQCR')RMS
RMS
RMS
RMS
RMS
RMS
IF(ISBUG2.EQ.0) GO TO If 00
DO 1350 JRKK = I,INC
JJRK = ICHO(JRKK)
JRK = KCITY(JJRK)
WRITE(LWRIT,1315)
WRITE(LWRITtl3\ 1) (IAQCR(JRK,KRKl ,KRK=l ,20) ,KYY
1311 FORMAT(lH ,IIf,19AIf,'YEAR = ',IIf,'= NOW (T)',/I
WRITE(LWRIT,1312)
1312 FORMAT(IH ,'IND VA 1967 VA NOW VA (T-l)
1V-67 INV NOW WB 1967 WB NOW PI NOW
1315 FORMAT (lHl)
DO 13lfa KRK = 1,20
IF (XXH(JRK,KRK).EQ.o.O) GO TO 13lfO
~RITE(LWRIT,1310) KRK,XXH(JRK,KRKI,XX(JRK,KRK),XXT(JRK,KRK),
1 SUMNH(JRK,KRK),SUMN(JRK,KRK),XIH(JRK,KRK),XI(JRK,KRK),
EMP(M)
Q-67
759
760
761
762
71>3
761f
765
71>6
767
768
769
770
771
772
773
771f
775
776
777
778
779
780
781
782
783
78q
785
786
7P.7
788
789
790
791
792
793
791f
795
796
-------
lSN 0701
ISN 0702
ISN 0703
ISN 0705
ISN 0706
ISN 0707
ISN 0708
ISN 0709
ISN 0710
ISN 0711
ISN 0712
ISN 0713
ISN 071Q
0'
o
ISN 0715
ISN 0716
2W3H JR( ,K~(),WB(JRK,KR(),P JRK,K~K),
3 EOLDH(JRK,KRK) ,EOlD(JRK,KRK,5) ,JRK
1320 FORMATIIH ,10Fl0.1,2X,13)
13"10 CONTINUE
IF(YYH(JRK).EQ.O.O) GO TO 1350
WRITE(LWRIT,1323)
1323 FORMAT(IH ,II
WRITE(LWRIT,1321)
1321 FORMAT(IH,' INC 1967 INC NOW EMP(C)-67 EMP«()NOW REV 1967
IV NOW EXp 19A7 EXP NOW (ONSMP-67 (ONSMP NOW AQ(R')
WRITEILWRIT,137.0) YYH(JRK) ,YY(JRKI ,XNBAH(JRK) ,XNBAR(JRK),
1 TH(JRK) ,T(JRKI ,GHIJRK) ,G(Jr~KI ,(H(JRKI ,(JRK) ,JRK
WRITE(LWRIT,13231
WRITE(LWRIT,13221
1322 FORMAT(IH ,0 EMP(T)-67 EMp(TINOW LABOR-67 LAROR NOW QM 1967
1M NOW Q( 1967 Q( NOW QBAR 1967 QBAR NOW AQCR')
WRITE(LWRIT,13?01 XNTHIJRK),XNT(JRK) ,XLH(JRKI,XL(JRK) ,QMH(JRK),
lQM(JRK) ,QCH(JRK) ,Q(JRI
-------
MCN", NPT(2) ,NO;)~1 ,LPNT1, (LPNT2( I),
2 tlPRNI tlPRN2 dPRN3 tISBUG
3,ICTAB,IGOV
= ,NOPE 1),~ UMC IT
,ISBUG2
RMS
Rr-,S
RM5
RMS
RMS
C
C
e...........~O.Q..Qo~e..o.oo...........................................oRMS
e RMS
ISN 0717 00 7183 I = 1,2 RMS
ISN 0718 7183 READILREAD,71B'i) (T 1 TLE (J, I ), J=I,20) RMS
ISN 0719 718'4 FORMAT(ZOA'i) Rr1S
C RHS
e RMS
e READ IN FIXED DELTe FOR 3 1 2 ( A) REPORT RMS
C Rr1S
ISN 0720 IF(IA312.EQ.O) GO TO 7175 RMS
ISN 0722 REAOILREAo,7212) BEN RMS
ISN 0723 7212 FORMAT(AF5.0) RMS
ISN 072«4 ADJUST = BENEF.(BENIJYEAR-3)/BEN(8) RMS
ISN 0725 READ(LREAD,717'i) QGIVEN RMS
ISN 0726 717'i FORMAT(7FI0.'i) RMS
ISN 0727 7175 CONTINUE RMS
~ lSN 072A KPRNl = lPRNl RMS
ISN 0729 KPRN2 = IPRN2 Rr1S
ISN 0730 KPRN3 = lPRN3 RMS
ISN 0731 KSBUG = lSBUG RMS
lSN 0732 IF(INPT(l )+INPT(2).EQ.O) STOP RMS
ISN 073«4 NOPEN = NOPE + NOPEI RMS
lSN 0735 NUMCIT = NUMCIT + 1 RMS
ISN 0736 IHELP = 0 RMS
ISN 0737 700 CONTINUE RMS
ISN 0738 NUMCIT = NUMCIT - 1 RMS
ISN 0739 NEwCIT = 1 RMS
ISN 07'iO NBEFOR = 0 RMS
ISN 07'il NBEFOI = 1 RMS
ISN 07'i2 NBEF03 = 0 RMS
ISN 07'i3 NBEF02 = 1 RMS
ISN 07'i'4 IPRNl = KPRNI RMS
ISN 07«45 IPRN2 = KPRN2 Rt.'S
835
836
837
838
839
8'iO
8 ~ 1
8'42
8~3
8'i'4
8~5
8'i6
8L17
8L18
8Lf9
850
851
852
853
85'i
855
856
857
858
859
860
861
867.
863
86'i
865
866
867
868
869
870
871
872
-------
ISN 07.6
ISN 07Q7
ISN 07&+8
ISN 07&+9
ISN 0750
(SN 0751
ISN 0752
ISN 075Q
ISN 0755
ISN 0757
ISN 0758
ISN 0759
ISN 0760
ISN 0761
ISN 0762
ISN 0763
(SN 076Q
ISN 0765
ISN 0766
ISN 0767
0" I SN 0768
N°ISN 0769
ISN 0770
ISN 0771
ISN 0772
ISN 0773
ISN 077tf
ISN 0775
ISN 0776
ISN 0777
ISN 077B
ISN 0779
ISN 0780
ISN 0781
ISN 0787.
ISN 0783
PR~3 = K~{1/3
ISBlJG = KSBUG
NOPEl = NOPEN-NOPE
NOK = NOPE
(HELL = 0
705 IHELP = IHELP + 1
IF(IHELP.LT.3) GO TO 706
IHELP = 1
IF (NOPEN.EQ.l)
LSTRAT = 8
IKSW = I
DO 7999 K = I,MCNT
[K = MC[TY = XX192,20) + XXII,7.0)
PII92,20> = PII92,20> + PI(I,?O>
XII9Z,ZO> = XI(9Z,ZO> + XI(I,ZO)
XK(9Z,ZO) = XK(92,ZO> + XK(I ,ZO)
SUMNI92,ZO) = SUMN(97.,20) + SUMNII,20)
XN6AR(92) = XNBAR(9Z> + XNBAR(I)
YY(92) = YY192> + YY(I)
(92) = (9Z) + (I)
XNT(92) = XNTI92> + XNT(I>
XL(92) = XL(9Z> + XL(I)
QTI9Z) = QTI9Z) + QTII)
QM(92) = QMI9Z) + QMII)
Q8AR(92) = Q8AR(9Z) + QBARII)
QC(92) = QC(9Z) + QC(I)
G(92) = G(92) + G(I>
T(92) = T(92) + T(I)
TP(92) = TP(9Z) + TP(I)
TO(92) = TO(92) + TO(I)
TA(92) = TA(9Z) + TA(I)
POP(92) = POP(92) + POP(I)
APV(92) = APV(92) + APVII>
RMS
RMS
RMS
RMS
Rr-1S
R r1 c;
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
Rt-1S
RMS
RMS
RMS
RMS
RMS
RMS
R~1S
R~'S
Rr.,s
RMS
RMS
RMS
RMS
RMS
RMS
RMS
GO TO 70tf
C
C
873
871.f
875
876
877
87A
879
8AO
881
887.
883
88Lf
885
886
887
8AB
889
890
891
892
893
89q
895
896
897
898
899
900
901
902
903
90Lf
905
906
907
908
909
910
-------
ISN 0784
ISN 078S
ISN 0786
ISN 0788
ISN 0790
ISN 0791
ISN 0793
ISN 079~
0'
W
ISN 0795
ISN 0796
I5N 0797
ISN 0798
ISN 0799
ISN 0800
C
7999 CONTINUE
IKSW = 0.0
704 IF(NUMCIT.LT.l.AND.IDYEAR.NE.1967)
IF(NUMCIT.EQ.O) GO TO 'i00
706 INTT = INPT(IHELP}
IF(INTT.EQ.OI GO TO 705
GO TO (~81 ,10) .INTT
C
C
C
6001
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
RMS
Rt-15
RMS
RMS
RMS
RMS
RMS
RM5
RHS
RHS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
TWO DIGIT INDUSTRY INFORMATION .. TWO ALTERNATIVES... RMS
(I) DETAILED SIC INFORMATION ONE SMSA AT A TIME (LPNTI = lRMS
(21 ONE TWO DIGIT SIC ( AT A TI~E ) FOR ALL CITIES (LPNT1=RMS
WITHIN EITHER CATEGORY. THE ECONOMY CAN BE AFfECTED IN ONE Of TRMS
DEPENDING ON THE OAT A GIVEN. RM5
(1) COST FIGURES ARE GIVEN DIRECTLY (LPNT2=1) RMS
(7.) PRICE CHANGES FOR FU~L ARE GIVEN (LPNT2=2) RMS
NOX= NO. OF CITIES FOR WHICH TWO DIGIT SIC COST DATA IS GIVEN RMS
(LPNTI=I) RMS
RMS
RM5
RMS
RMS
RMS
RMS
RMs
RMS
RMS
RMS
RM5
Rt-15
GO TO 13S5
FORMAT(II0.SFIO.2)
10 CONTINUE
720 KOUNT = 0
IHELL = IHELL + I
LPNTT = LPNT2(IHELL)
LsTRAT = LPNTT + 5
C
SS DO 70 1=1,91
9 1
912
913
91 'i
91S
916
917
918
919
97.0
971
922
923.
92~
97.S
92b
927
92~
929
930
931
932
933
93~
93S
936
937
938
939
9'i0
9''11
942
943
94'i
94S
946
947
948
-------
SN 0801 DNR( ) = 1. ( RMS 949
ISN 0802 PON(I) = 0.0 RMS 950
ISN 0803 PPID(I) = 0.0 RMS 951
ISN 0801t PXID(I) = 0.0 RMS 952
ISN 0805 PXKD(J) = 0.0 Rr1S 953
ISN 0806 PXD(J) = 0.0 Rr1S 9S,+
ISN 0807 EO(I,5) = 0.0 RMC; 9S5
ISN 0808 PIDII) = 0.0 RMS 956
ISN 0809 '(181)(1) = 0.0 RMS 957
ISN 0810 XIDII) = 0.0 RMS 958
ISN 081 1 XKD(I) = 0.0 RMS 959
ISN 0812 PDN ( I) = 0.0 RMS 960
ISN OB13 70 X D ( I) = 0.0 RMS 961
C RMS 962
C IF LPNTI = 1 ° ... COST DATA fOR EACH CITY IS PRECEEDED BY THE RMS 963
C ACTUAL CITY CODE NUMBER AND THE NO. OF IN D. IN THAT CITY FOR RMS 96'+
C WHICH COST DATA IS GIVEN. RMS 965
C RMS 966
C IF LPNTI = 2 COST DATA FOR EACH I ND. IS PRE(EEDED BY THE SIC CORMS 967
( NUNBER AND THE NO. OF (ITIES FOR WHICH THAT IND. EXISTS AND HARMS 968
0" ( COST DATA. RMS 969
,j:>.o ( RMS 970
ISN 081'+ ,.,CNT = 1 RMS 971
ISN 0815 I I = 1 RMS 972
ISN 0816 NOTE = I RMS 973
ISN 0817 MAX = 0 RMS 97'+
ISN 0818 I F (LPNTT .NE. 3) GO TO 5'+ RMS 975
C RMS 976
( RMS 977
( GENERAL FEEDBACK SEGMENT RMS 978
C OPTION NUMBER EIGHT RMS 979
( RMS 980
ISN 0820 5992 FORMATI2(7FID.If,/) ,5FI0.'+,23X,13) RMS 981
ISN 0821 KOUNT = KOUNT + 1 RMS 982
ISN 0822 NOI = 0 RMS 983
ISN 0823 IFEED = 1 RMS 98,+
ISN 082'+ I K = ,., C I T Y ( 1 ) RMS 985
C I K = MMCITY(KOUNT) RMS 986
-------
S' 0825
ISN 0826
ISN 0827
ISN 0829
ISN 0830
ISN 0831
ISN 0832
ISN 083~
ISN 0836
ISN 0837
ISN 0838
ISN 0839
ISN 0840
ISN 08~1
ISN 0842
ISN 0843
ISN 0844
ISN 0845
ISN 08'16
ISN 08'17
'" I SN 08'18
IJ1
ISN 08'19
ISN 0850
ISN 0851
ISN 0852
ISN 0853
ISN 085'1
ISN 0855
ISN 0857
ISN 0858
ISN 0859
ISN 0860
LCITY = KC '.YII
-------
ISN 0862
ISN 0863
ISN 086'+
ISN 0865
ISN 0866
ISN 0867
ISN 086R
ISN 0869
ISN 0870
ISN 0871
ISN 0872
ISN 087.3
ISN 087'+
ISN 0875
ISN 0876
ISN 0877
ISN 0878
ISN 0880
ISN 0881
0' ISN 0882
0' ISN 0883
ISN 088'+
ISN 0885
ISN 0886
ISN 0887
ISN 0889
ISN 0890
ISN 0891
ISN 0892
ISN 0893
ISN 089Lf
ISN 0895
ISN 0896
ISN 0897
ISN 0899
REA) LREA), LflG) JCITY,NOI
GO TO 7153
7151 NOI = b
MIN = NOI
READILREAO,7152) HHI,JCITY
LCITY = KCITYIJCITY)
MCITYIII) = JCITY
REAO(LREAD,715Z1 HHA,JCITy
FUIIJCITYpJYEAR-31=FUIIJCITY,JYEAR-3)+HHIC7)
FUCJCITY,JYEAR-3)=FUIJCITy,JYEAR-3)+HHAC7)
DELTCCLCITY) = Oon
7152 FORMATI7FIO.Lf,3X,I3)
GO TO 725
7153 CONTINUE
C
RMS 1025
RMS 1026
RMS 1027
RMS 1028
RMS 107.9
.RM5 1030
RM5 1031
RMS 1032
RM5 1033
RMS 103'+
RMS 1035
RMS 1036
RMS 1037
Rt~S 1038
RMS 1039
R r1 5 1 0 'I 0
RMS 10'11
RMS 1042
RM5 10Li3
RMS 10LfLf
RMS 10'15
RMS 10Lfb
RMS 10Lf7
RMS 10LiA
RMS 10Lf9
RMS 10S0
RMS 1051
RMS 1052
RMS 1053
RMS 105Lf
RMS 1055
RMS 1056
RMS 1057
RMS 1058
RMS 1059
RMS 1060
Rt~S 1061
RMS 1062
LCfTY = KCITYIJCITY)
INDI = NOI -20 + 120/NOI)
IF(LPNT1.EQ.2) GO TO 1
HCITYCII) = JCITY
MIN = NOI
GO TO 2
1 MCNT = JCITY
MIN = JCITY
MAX = JKEEPIINDl)
c
2 ITT = 5
IFILPNTT.NE.I) GO TO 725
ITT = 1
READILREAD,730) DELTCILCIrY)
DELTCILCITY) = DELTCILCITY) I
730 FORMATII0X,FIO.2)
725 CONTINUE
QTILCITY>
c
IH = 0
J = 0
DO 556 JJ=I,MIN
I F I I A 3120 EQ. 1) GO
READ ILREAD,60011
TO 715'1
I K
-------
S~ 090)
ISN 0901
ISN 0902
ISN 0903
ISN 090£+
ISN 0905
ISN 0907
ISN 0908
ISN (1909
ISN 0910
ISN 0911
ISN 0912
ISN 0913
ISN 0915
ISN 0916
ISN 0917
ISN 0918
0" ISN 0919
~ ISN 0920
ISN 0921
ISN 0922
ISN 0923
ISN 0925
ISN 0921,
ISN 0927
ISN 0928
ISN 0929
ISN 0930
ISN 0931
ISN 0932
ISN 093£+
C
C
HOLDI ) = 0.0
HOLD(2) = 0.0
HOLD!) = 0.0
HOLD(L!) = 0.0
HOLD(5) = 0.0
NO MORE FUEL CALCULATION FOR HOLDIL!)
READ(LREAD,6001) IK,(HOLDIIS), I5=I,ITT)
IF(LPNT1.EQ.I) GO TO :3
IUN = KCITYIIK)
LCITY = IUN
MCITYIJ) = IK
GO TO £+0
3 INDI = IK - 20 +(20/IKI
NEXT = JKEEPIIND1)
IF (MAX.LT.NEXTI MAX = NEXT
IUN = INDI
INDIJI = INDI
£+0 GO TO (£+1,£+1, LPNTT
HS 1063
RMS 106£+
RMS 1065
RMS 1066
RMS 1067
RMS 1068
RMS 1069
R t-1 5 1 0 7 0
Ri1S 1071
RMS lcJ72
RMS 107)
RMS 107£+
RM5 1075
RMC; 1076
RMS 1077
RMS 1078
RMS 1079
RMS 10RO
RMS 1081
RM'5 1082
RMS lOR)
RMS 108£+
RMS 1085
RMS 10Rb
RMS IOR7
Rt-1S 1088
RMS 1089
RMS 1090
RMS 1091
RI1S 1092
RMS 1093
RMS 109'1
RM'5 1095
RMS 1096
RMS 1097
RMc; 109R
RMS 1099
RMS 110'0
C
715'+
IK = ISPINIJJ)
INDI = IK - 20 + (20/IK)
HOLD(1) = HHAIJJ)
HOLD(Z) = HHIIJJ)
IUN = INDI
IF ( I NDC I (lC I T y, I NO 1 ) . EQ. 0)
J = J + 1
INDIJI = INDI
GO TO 'II
IH = IH + 1
GO TO 556
GO TO 7155
7155
C
C
C
LPNT2 = 1 ... COST FIGURES GIVEN DIRECTLY
C
£+1 CONTINUE
LSM = iiI
P = I.
IFIIGOV.EQ.O) GO
HERE = HOLO(I)
TO 36£+1
-------
ISN 09:5
ISN 0936
ISN 0937
ISN 0938
ISN 0939
ISN 091i0
ISN 09lil
ISN 091t2
ISN 091i3
ISN 091i1f
ISN 091f5
ISN 091i6
ISN 091t7
ISN 09ltA
ISN
ISN
ISN
ISN
g;ISN
ISN
ISN
ISN
ISN
ISN
ISN
091f9
0950
0951
0952
0953
0951f
0955
0957
0958
0960
0962
C
C
C
C
C
C
C
C
C
C
"fA-iA = .OLDIZ)
HOLDll) = P . HERE
HOLDIZ) = P . HAHA
36lfl CONTINUE
PPIDIIUN) = -HOLD II) - DELTCILCITY> . EOLDILCITY,INDl,S) . P
QPOLD = PIILCITYdNDI)
QPNEW =QPOLD + PPIDIIUN)
QPAR = QINVIINDI,1 ).IXIHILCITY,INDI )-IOIINDI).XKHILCITY,INOI»)
QIOLD )=QPAR+QINVIINDl,Z).IQPOLO-PIHILCITy,INDI»+QINVIINDI,3.
1 XKILCITY,INDl)
QINEW =QPAR+QINVIINDl,Z).IQPNEW-PIHILCITy,INDI »+QINVIINDl,).
1 XKILCITV,INDI) -HOLD(2)
PXIDIIUN) = QINEW-QIOLD
PXKOIJUN) = PXIDIIUN)/II.+DIINDl»
K = LCITY
PXDIIUN) =-XXILCITY,INDl)+XXILCITY,I~DI).IPXKDIIUN)+XKIK,INDI»
1 /XKILCITY,INDl)
EDIIUN,S) = ABIIUN).PXDIIUN)
GO TO S56
It WRITE ILWRIT,S)
S FORMAT I1Hl,'OOPS-TRIED TO ACCESS FUEL SECTION AT RMS00858',IIII)
556 CONTINUE
MIN = MIN - IH
IF(MIN.EQ.0.AND.IA312.EQ.l) GO TO 5555
LSM = 556
IF(LPNTl.EQ.l) NOI = MIN
IF(LPNTl.EQ.Z) JCITY = MIN
IFILPNTT .EQ. 1) GO TO 17
..... EXTINCT .....
LPNTZ = Z ... CHANGE IN FUEL PRICES GIVEN
FOURTH RECURSIVE BLOCK - FUEL DATA
TWO DIGIT SIC - 12 EQUATIONS
ENDOGENOUS VARIABLES ARE...
EOLD - QUATITY OF EACH TYPE OF FUEL USED BY A GIVEN IND.
CTOT67 - TOTAL COST OF FUEL 1967
XX - VALUE ADDED 1967
R~S 1 )1
RMS 1102
RMS 1103
RMS I11J1t
RMS 1105
RMS I1n6
Rf-1S 1107
RMS 1108
RMS 1109
RMS 1110
Rf-1S 1111
RMS 1112
RMS 1113
Rf-1S I111t
RMS 1115
RMS 1116
RMo; 1117
RMS 1118
RMS 1119
RMS 1120
RMS 1121
RMS 1122
RMS 1123
RMS 1121f
RMS 1125
RMS 1126
RMS 11'27
RMS 1128
RMS 1129
RMS 1130
RMS 1131
RMS 1132
RMS 1133
RMS 113'+
RMS 1135
RMS 1136
RHS 1137
RHS 1138
-------
ISN 096'"
ISN 0965
ISN 0966
ISN 0967
ISN 096R
ISN 0969
ISN 0971
0"
'-.0 ISN 0977.
ISN 0973
ISN 097Li
ISN 0975
ISN 097b
ISN 0977
ISN 0979
ISN 0980
ISN 0981
ISN 0982
ISN 098,..
ISN 0985
ISN 0986
ISN 0987
ISN 0988
ISN 0989
c
c
C
C
C
C
C
C
C
C
C
C
C
C
Zz - TO"AL QUA" "Y OF Fl EL CO~ SU '1ED
XI,XK, AND PI
EXOGENOUS VARIARLES...
POLO - PRICE OF FUEL 1967
PNEW - PROJECTED PRICE OF FUEL FOR 1967
WITH POLUTInN COSTS
GAMA - FUEL EXPONENTS - CnNSTANTS FOR 1967
XKT - 1963
A - THE COEF. OF THE FUEL FUNCTION
ALIT - THE COEF. FROM THE FUNCTION Z=A.V
D,BETA,PICOEF, AND XX/XK
EXTINCT
.....
....~
WRJTECLWRIT,6)
6 fORMATCIHl ,'OOPS-TRIED TO ACCESS FUEL SECTION AT RMS00917',/////)
C
17 CONTINUE
M = 0
DO 22 IT = I,MIN
IFCLPNTI.EQ.2) GO TO 18
INDI = INDI IT)
IUN = INDI
GO TO 19
18 IK = MCITYCIT)
LCITY = KCITYIJK)
IUN = LCJTY
19 IF(LPNTT.EQ.I.OR.LPNTT.EQ.3) GO TO 21
IVRITECLWRIT,7.3)
23 FORMATCIHI ,'OOPS TRIED TO ACCESS FUEL SECTION AT RMS00998'.////)
21 I = LCITY
IFCWII.TNDI ).SUMNCI,20).EQ.0.0) GO TO 39
PDNC IUN) = SETAC JNDI ).PXDC IUN)fWI I, INDI)
DNRII) = DNRII) -IPDNCIUN)/SUMNII,20»
39 CONTINUE
37 FORHATCIH ,2I3,FI2.I,FI2.,..)
PIDCI) = PIDCI) + PPIOIIUN)
XIDCI) = XIDCJ) + PXIDIIUN)
RMS 39
RMS 11'+0
RMS 11'+1
RM5 11"'2
RMS 11"'3
RMS II"''''
RMS II"'S
RMS 11"'6
RMS 11"'7
RMS 11"'8
RMS 11'+9
RMS 1150
RMS 1151
Rt1S 1152
Rr1S 1153
RMS 115'+
RMS 1155
RMS 1156
RMS 1157
RMS 1158
RMS 1159
RMS I 160
RMS 1161
RMS 1162
RMS 1163
RMS Ilb,+
RMS 1165
RMS 1166
RMS 1167
RMS 1168
RMS 1169
RMS 1170
RMS 1171
RMS 1172
RMS 1173
RMS 117'+
RMS 1175
RMS 1176
-------
ISN 0990
ISN 0991
ISN 0992
ISN 0993
ISN 099'+
ISN 0995
ISN 0996
ISN 0997
ISN 0998
ISN 0999
-.J
o
ISN 1000
ISN 1001
ISN 1002
ISN 1003
ISN 100'+
JSN 1005
ISN 1007
ISN 1009
ISN 1011
XKDCI) = XK)()
XOCI) = XDCI) +
XOICI) = XOCI)
wBOq) = WBOCI)
02(I) = WSOCI)
22 CONTINUE
NOTE = 1
GO TO 5555
+ PXKDCIUN)
PXDC IUN)
R~ 5
RMS
Rf-1S
RMS
RMS
RM5
RMS
RMS
RMS
RM5
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RM5
RMS
RMS
RM5
RMS
GIVEN - PROFIT AFFECRMS
GIVE - INCOME AFFECTRMS
RMS
RMS
Rt1S
RMS
RM5
RMS
RM5
RMS
RMS
RMS
R/1S
RMS
+ BETACIND1).PXOCIUN)
C
C
(
(
C
C
C
(
C
C
INPT(2) = 1
ALL CITy INFO. FOR ALL IND. REQUESTED
li81 (aNT I NUE
NOTE = NOPE
C
C
C
(
C
C
C
C
C
NOTE»>AROVE READ INDICATES THE TYPE OF CONTROL (OsT DATA GIVEN
(\) NOPE> TYPES OF POLLUTION IMPACT TO 8E EXAMINED FOR EACH CITY
= 2»INCOME CHANGE GIVEN
= 'fo.ElECTRIC COST INCREASE
= S..ELECTRIC COST INCREASE
= »>TAX (HANGE GIVEN
(3) NUMBER OF CITIES FOR WHICH DATA IS GIVEN - MCNT
IF (NOPE-5) 102,102,li82
'f82 ~RITECLWRIT,'f83) NOPE,CLPARAC 11, 1=I,NOPE)
'f83 FORMATCIHI ,'BAD CONTROL COST PARAMETER CARO',lil'f)
STOP
10200 103 I=I,NOPE
IFCLPARACI).GT.5.0R.LPARACI).LT.l) GO TO '+82
IFC1.EQ. NOPE) GO TO 103
IF (LPARACI+l) .EQ. LPARA(I» GO TO 'fa2
103 (ONTINUE
1 177
1 1 7 A
1 179
1 1 A 0
1 \ 8 1
1 1 A 2
1 183
l1A'f
1 1 AS
1 1 A 6
1 1 A 7
1 1 A A
1 189
1 190
1 1 9 1
1 192
1 193
1 19,+
1 195
1 196
1 197
1 1 98
1 199
1200
1 2n 1
\202
1203
12nli
\205
1206
1207
\2118
1209
1 2 \ 0
1 21 1
1212
1 213
1 2 1 '4'
-------
C ~MS 1 2 5
C ....- MAKE EASILY ADJUSTABLE READ RMS 1216
C RMS 1 217
C 10~ REAOILREAO,105) HCITYII) .(COSTIJ,I),J=1 ,5) RMS I 218
C RMS 1 219
ISN 1012 IKSW = 0 RMS 1220
C RMS 122 1
C ZERO LOOp REMOVED FROM HERE RMS 1222
C RMS 17.7.3
ISN 1013 I F I I A 312. EQ. 1) GO TO 7100 RMS 122~
C RMS 1275
C RMS 1226
C COST INPUT FOR 1970-1980 SIMULATION RMS 1227
C RMS 1228
C RMS 1229
ISN 1015 DO 10Li I=I.MCNT RMS 1230
ISN 1016 SCIN = 0.0 Rt~S 123 1
ISN 101 7 SRIN = 0.0 RMS 1232
ISN 1 0 1 8 SC = 0.0 RM5 1233
I5N 1 n 1 9 SR = 0.0 RMS 123~
-J ISN 1020 109 READILREAD,105) RIN,R.CIN.CC,MCITY(I) RMS 1235
...... IS N 1021 RISUM = RIN + R + CIN + CC RMS 1236
ISN 1022 IF(RISUM.EQ.O) GO TO 1 1 0 Rt-lS 1237
ISN 102~ GO TO 1 1 1 RMS 1238
ISN 1025 1 1 0 NUMCIT = NUMCIT - 1 RMS 1239
ISN 102ft NEWCIT = 1 RMS 12110
ISN 1027 IF(NUMCIT.LT.l) GO TO 1355 RMS 1211 1
ISN 1029 GO TO 109 RMS 12112
ISN 1030 1 1 1 CONTINUE Rr-lS 12113
ISN 1031 K I = MCITY(I) RMS 12'1Lf
ISN 1032 L = KCITY(KI) RMS 12115
ISN 1033 IFIL.EQ.O) GO TO 10Lf RMS 12LJ6
C RMS 12'17
C USED FOR COST STRAT NO. 3 ONLY RMS 12LfR
C RMS 12Lf9
C COST (Lf , I) = COSTILJ.I) -.5 RMS 1250
C COST (3. I) = 30. Rr-lS 1251
ISN 1035 C 05 T I 1 . I) = CC-(R-.b) RMS 1 2 S"i'
-------
ISN 1036
ISN 1037
ISN 103A
ISN 1039
ISN 101i0
ISN 10lil
ISN 10'12
~ ISN 101i3
ISN 10li'l
ISN 10'15
ISN 10'16
ISN 10'17
ISN 10'18
ISN 10'19
ISN 1050
ISN 1051
ISN 1053
ISN 1051i
ISN 1055
ISN 1056
ISN 1057
(OFIII) = (IN-IRIN'.6)
(OSTI'I,I) = R'.6+RIN'.6
(OSTI3,1) = 0.0
DELT(II) = (OSTI'I,I)/QTIL)
RMS 1253
RMS 125'1
RMS 12S5
RMS 1256
RMS 1257
RMS 125A
RMS 12S9
RM5 1260
RMS 1261
RMS 1267.
RMS 1263
RMS 126'1
RMS 1265
RMS 1266
RMS 1267
RMS 126A
RMS 1269
RMS 1270
RMS 1271
RMS 1272
RMS 1273
Rt-1S 127'1
RMS 1275
RM5 1276
RMS 1277
RMS 1278
RMS 1279
RM5 1280
RMS 1281
RMS 12R2
RMS 1283
RM5 128'1
RMS 1285
RMS 1286
RMS 1287
RMS 12B8
RM5 1289
RMS 1290
(
(
C
U5ED ONLy FOR COST STRAT. NO.2 AND 3
(05T(I,I) = (05Tll,I)-(67.3.7-C05T(I,I)/808.067 )
C
(
C
C
(
C
(
USED FOR (OST STRAT NO.3 ONLY
C 05 T ( 1 , I) = (05 T I 1 , I) '.5
U5ED ONLy FOR C05T STRAT. NO.2 AND 3
(OFlll) = COFlll)-1623.7.COFII [)/Sna.067
C
C
(
(
U5ED FOR (OST STRAT NO.3 ONLY
(OFI(I) = (OF III I ".5
105 FORMAT(20X,2FI0.1/20X,2FI0.1,33X,I3,'X)
10'1 (ONTINUE .
GO TO 5555
(
(
C
7100
INPUT FOR 312(A) REPORT
DO 710'1 1=I,MCNT
READ(LREAD,7101) (IN,51N,RIR,R,MCITY(I)
7101 FORMATI'IFI0.5,33X,I3)
7111 CONTINUE
KI = r"cITYlII
L = K(ITY(KII
IF(L.EQ.O) GO TO 710'1
C05Tll,l) = SIN-(l.+AAA)'.'1 +R
DELTCII) = .60-IICIN-II.+AIA»+15IN.II.+AAA»)/QTIL)
COFIII) = CIN-II.+AIA)..'O +RIR
FUIIKI,JYEAR-3)=FUIIKI,JYEAR-3)+COFIII)
FU(KI,JYEAR-31=FUIKI,JYEAR-3)+COSTII,II+DELTC(I)'QTIL)
-------
ISN }( 58
ISN 1060
ISN 1061
ISN 1062
ISN 1063
ISN 106'1
ISN 1065
,-.]
vi
ISN 1066
ISN 1067
ISN 1068
ISN 1069
ISN 1070
ISN 1072
ISN 107'1
ISN 1075
ISN 1076
F(IGOV.EQ.o) GO.O 7 01 ~.,S 1291
HERE = (05T(I,I> RMS 1292
HAHA = caFI (I) RMS 1293
C RMS 129't
(.....~.....&...~....o...o............o.................................RMS 1295
C RMS 1296
DELTCII) = P .DELTCII) RMS 1297
(05T(I,I> = P . HERE RMS 1298
COFIII> = P . HAHA RMS 1299
C RMS 1300
C RMS 1301
C[[[RMS 1302
C RMS 1303
710't CONTINUE RMS 130'1
C ..&.. EXTINCT ..... RMS 1305
C RMS 1306
C RMS 1307
C COST I NPUT FOR 1967 RUNS RMS 1308
C RMS 1309
C ....0 EXTINCT ..... RMS 1310
C RMS 1311
C RMS 1312
5555 CONTINUE RMS 1313
LSM = 5555 RMS 131't
DO 5000 III = I,NOTE RMS 1315
II = 1 RMS 1316
IF(INTT.EQ.2.AND.MIN.EQ.0.AND.IA312.EQ.1> GO TO 1000 RMS 1317
IFIINTT.EQ.2) GO TO 100 RMS 1318
LSTRAT = LPARA(III) RMS 1319
II = LSTRAT RMS 1320
GO TO (1000,2000,3000dOOO,2000),LSTRAT RMS 1371
RMS.1322
RMS 1323
. RMS 132't
RMS 1325
R t~ c; 1 3 2 6
RM5 1327
RMS 1328
C
C
C
C
C
C
C
ALL CITY INFORMATION ALL IND......EFFECT OF DIRECT REDUCTION
OF INDUSTRY'S PROFITS BY AIR PLOUTION CONTROL EQUIPMENT COSTS
-------
ISN 1077
ISN 1078
ISN 1079
lSN 1080
ISN 1082
ISN 1083
ISN 1085
ISN 1086
lSN 1087
ISN 10BA
ISN 1089
ISN 1090
lSN 1091
15N 1092
ISN 1093
ISN 109Lf
-.J
,p..
ISN 1096
ISN 1097
ISN 1098
ISN 1099
ISN 1100
ISN 1101
ISN 1102
ISN 1103
ISN 110Lf
ISN 1105
ISN 1106
ISN 1107
~MS 1329
RMS 1330
RMS 1331
RMS 1332
00 IDOl K=I,MCNT RMS 1333
KI = MCITY(K) RMS 133Lf
I = KCITY(KI) RMS 1335
IF (1.f.Q.O) GO TO IDOl RMS 1336
ONR(I) = 100 RMS 1337
IF(INTT.EQ.2.ANO.MIN.EQ.0.AND.IA312.EQ.I) GO TO 7192 RMS 1338
IFILSTRAT-Lf) 311,300,300 RHS 1339
300 COSTlll,K) = DELTCIK) . QM(I) RMS 13LfO
LSH = 300 RMS 13Lfl
COFI(K) = 0.0 RMS 131.42
GO TO 311 Rt..,S 13Lf3
7192 COST(II,K) = HHA(7) RMS 13LfLf
LSM = 7192 RMS 13Lf5
COFl
-------
ISN IIOR
ISN 1109
ISN 1110
ISN 1111
ISN 1112
-J ISN
U1 ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
1 1 1 3
1 1 1 'f
1 1 1 5
1 1 1 ~
1 1 1 7
1 1 1 R
1 1 1 9
1 120
1 1 2 I
1127.
1 1 2 'f
1 125
1 126
1 127
1 128
ISN 1129
ISN 1130
1001
C
C
C
C
C
C
C
C
C
C
C
C
C
X)I I ) = XD )
~8DII) = BETA(20).XDII)
D2II) = WBDII)
CONTINUE
FIRST SIMULTANEOUS SYSTEM
Iy,G,T,
WHERE THE ENDOGENOUS VAR. ARE...
YY ~ TOT. REGIONAL INCOME
T - GOVERNMENT REVENUE
C - CONSUMPTION 119671
THE EXOGENOUS VAR. ARE ,..
Xx - VALUE ADOED IACTUAL VALUE
G - GOVERNMENT EXPENDITURE
C
C
100 LCHECK = 0
LVAR = 'f
MCASES = MCNT.'i
LVAR = 3
MCASES = MCNT.LVAR
120 MM = 0
KPIST = 'i
M = 0
DO 1002 K=I,MCNT
KI = MCITYIK)
I = KCITYIKI)
IJCIT = I
IF (I.EQeO) GO TO 1002
YIM+I) = YYII)
YIM+2) = GIll
YIM+) = CII)
XIMM+l) = XXII,20) + XDIIJ)
XIMM+2) = TII)
XIMM+) = FUDGE(I,I)
X(MM+3) =0.0
XIMr.,+'f) =0.0
C
AND C)
R~S 1367
RMS 1368
RMS 1369
RM5 1370
RMS 1371
RMS 1372
RM5 1373
RMS 137'f
RMS 1375
RM5 137/)
RM5 1377
RMS 137B
RM5 1379
RMS 13BO
RMS 1381
RMS 1382
RM5 1383
RMS 138'f
RMS 1385
RMS 1386
RMS 1387
RMS 1388
RM5 1389
Rr.,s 1390
RM5 1391
RMS 1392
RMS 1393
RMS 139'f
RM5 1395
RMS 1396
RMS 1397
RMS 1398
RMS 1399
RM5 1'i00
RMS 1 'f1J 1
RMS I 'f(12
RMS 1 'i03
RMS l'iO'i
AND + DELTA)
-------
ISN 1131
ISN 1132
ISN 1133
ISN 113'+
ISN 1135
ISN 1136
ISN 1137
ISN 113R
ISN 1139
ISN 11'+0
ISN 1 1'+ 1
ISN 11'+2
ISN 11'+3
ISN 11'+'+
ISN 11'+5
ISN 11'+7
~ ISN 11'+~
ISN l1Lf9
ISN 1150
ISN 1151
ISN 1152
ISN 1153
ISN 115'+
ISN 1155
ISN 1156
ISN 1157
ISN 1158
ISN 1159
ISN 1160
ISN 1161
ISN 1162
ISN 1163
ISN 1165
ISN 1166
C
XC"1M+1) = FU)GE 1,2)
M = M + 3
MM = MM + '+
1002 CONTINUE
LSM = 1002
R~ 5 11 05
RMS 1'+06
RMS 1 Lf07
RMS lLf08
RMS 1 Lf09
RM5 1'+10
RMS 1'+11
RMS 1'+12
RM5 1'+13
RMS 1'+'1 Lf
RM5 1'+15
RMS lLf16
RMS 1'+17
RMS 1'+18
RMS 1'+19
RMS I Lf20
RMS lLf21
RM5 lLf22
RMS 1'+23
RM5 1'+2Lf
RMS 1'+25
RMS 1'+26
RMS 1 Lf27
RMS 1 '+2R
RMS 1 '+29
RMS lLf30
Rf151'+.11
RMS 1'+32
RMS 1 Lf33
RMS 1 '+.'3Lf
RMS 1 '+35
RMS 1 '+36
RMS 1 Lf37
RMS 1'+3R
RMS 1 Lf39
RMS ILf'+O
RMS 1 Lf'+ 1
R M S 1 If If ';f
C
MARK = 1
C
1 0 1 '+
951
8015
C
I 0 1 5
IF CICTAB) ROI5,BOI5,tOl'+
DO 951 IM=I,LVAR
ARABCIM) = SCCICIM)
CALL SIMUL
M = 0
IF(LCHECK)1003,1003,,+000
1003 DO 1001f K=I,MCNT
KI = MCITYCK)
I = KCITYCKI)
IF CI.EQ.O) GO TO 100'+
YYDCI) = Y(M+l)
GD(I) = YCM+2)
(D(I) = YCM+3)
XD1(J) = 0.0
M = M + 3
100'+ CONTINUE
LCHECK = 1
NBEFOR = NBEFOI
GO TO 120
'+000 M = 0
LSM = '+000
NEWCIT = 0
NBEFOI = 0
DO 1007 K= 1 ,11CNT
KI = r1CITY(K)
I = KC I T Y C K I )
IF (I.EQ.O) GO TO 1007
YYO(I) = YYOCI) - YCM+l)
GOCI) = GOCI) - YCM+2)
-------
ISIJ 1167
ISN 1168
ISN 1169
ISN 1170
ISN 1171
ISN 1172
ISN 1173
ISN 117ti
ISN 1175
ISN 1176
ISN 1177
ISN 1178
ISN 1179
ISN 1180
ISN 1181
ISN 1182
ISN 1183
ISN 118ti
ISN 1185
ISN 1186
-J
:-JISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
1 187
1 18 A
1 189
I 190
1 192
1 193
1 1 9 ti
1 195
1 196
1 197
1 19 A
1 199
1200
1201
1202
1203
()II) = CD I) - YIM...3)
M = M + LVAR
APVDII) = 0.0
TPDII) = .05608.APV(I)oAPVD(I)
TODII) = oOI223-YYD(I)
TADII) = .Oti71oYYDII)
TDI!) = TPDII)+TOD(I)+TAD(I)
QCDII) = .Oti718 . (DII)
QMOII) = .1991S . XD(I)
QBARDII) = .079ti7oIYYDII)-XDII»
QTOII) = QCO(I)+QMD(I)+QBARD(I)
XD 1 I I) = Y Y D I I )
IF ILSTRAT-ti) 320,330,330
320 DNBARII) = .I03239oIYYDII)-XDIJ»)
Dill) = ONBAfnI)
GO TO 1007
330 DNBARII) = .103239 ~ (YYD(I)-XD(I) - DELTC(K) . QBARD(I»
DI(I) = DNBARII)
1007 CONTINUE
(
~MS II ti3
RMS 1 tititi
RMS 11.J1.J5
RMS II.Jti6
RMS I titi7
RMS I titiA
RMS lti"'9
RMS 1"'50
RMS II.\SI
RMS 1 ti52
RMS I 1.J53
RMS 1 tiSI.J
RMS 1 1.J55
RMS 1 1.J56
RMS 1 tiS7
RMS I,..S8
RMS 1 I.JS9
RMS 11.\(,0
RMS 1"'61
RMS I ti62
RMS 11.J63
RMS 1I.J6I.J
RMS 1 1.J65
RMS 1 1.J66
RMS 1"'67
RMS 1 ti68
RMS 1 1.J69
RMS 1 1.J70
RMS 1"'71
RMS 11.J72
RMS lti73
RMS 1 1.J7,..
RMS 1 1.J75
RMS 1'+76
RHS 1 '+77
RMS 1I.J78
RMS 1,..79
RMS 1,..80
LSM = 1007
(
DO 1012 K=I,MCNT
KI = MCITYlK)
I = KCITYIKI)
IF II.EQ.O) GO TO 1012
SUMNDII) = D2(1)/WII,7.0)
XNDII) = DNRII). SUMNDII) + Dl(l)
XLDII) = 1.03921.\ 0 XNQ(!)
XLL = XLII) + XLDII)
UDII) =-IXNOII)/XLL)
Dill) = 0.0
D 2 I I) = O. 0
1012 CONTINUE
73I.J8 FORMATI3FI003,FI0.,..,33X,I3)
IF(IPRN3)351,3I.J8,351
351 CONT 1 NUE
GO TO 131.J3,3I.J,..,3I.JI.J,3I.J3,3I.J3,3I.\Q,3"'''',3,..I.J),LSTRAT
-------
ISN 120~
ISN 1205
ISN 1206
ISN 17.07
ISN 1208
ISN 1209
ISN 1210
ISN 1212
ISN 1213
ISN 121~
JSN 1215
JSN 1216
ISN 1217
ISN 121A
ISN 1220
JSN 1222
ISN 122~
JSN 1225
JSN 1226
ISN 1227
ISN 1228
-J ISN 1229
00 ISN 1230
ISN 1231
JSN 1232
JSN 1233
JSN 123~
JSN 1235
ISN 1236
ISN 1237
JSN 123A
ISN 1239
ISN 12~0
ISN 12~1
ISN 12~2
ISN 12~3
ISN 12~~
ISN 12~5
3~3 :- (IPRNl) 3qA,3~8,3~q
3~~ DO 3~5 K = I,MCNT
K I = 11 C I T Y I K )
1 = KCITYIKII
CALL PRINI III
3~5 CONTINUE
3~8 IF (NOPEN.EQ.I) GO TO 5001
ISTOP = 0
33~8 CONTINUE
ISTOP = ISTOP + 1
DO 1013 K = I,MCNT
KI = MCITY(K)
1 = KCITY(KI)
IFIISTOP.EQ.I) M = 1
IFIISTOP.EQ.2) M = 92
IF 11.EQ.0) GO TO 1013
PIN(MI = PIN(M) + PIDII)
XIN(M ) = XINIM) + XIOII)
XKN(M) = XKNIM) + XKD(I)
XDNIM) = XDN(M) + XD(I)
YYN(M) = YYN(M) + YYD(I)
GNIM) = GN(M) + GD(J)
TN(M) = TN(M) + TD(J)
TPN(M) = TPN(M) + TPD(I)
TON(H) = TON(H) + TOO(I)
TAN(M) = TANIM) + TAD(I)
OHM) = CN(M) + CD( I)
DNBARNIM) = DNAARN(M) + ONBAR (I)
XNN(M) = XNN(M) + XND(I)
XLN(M) = XLN(M) + XLD(I)
UNIM) = UN(M) + UD(J)
SUMNN(M) = SUMNN(M) + SUMND(J)
QBARN(M) = QBARNIM) + QBARDIJ)
~TN(M) = QTN(M) + OTO(I)
QCN(M) = OCN(M) + QCD(J)
QMN(M) = QMN(M) + QMD(I)
1013 CONTJNUE
LSM = 1013
RMS 1~81
RMS 1 ~82
RHS I'H\3
RM5 I~A~
RMS 1 'fAS
RM<; 1 ~R6
RMS 1 ~A7
RMS I11B8
RMS l11R9
RMS 1 ~90
RMS 1~91
RMS 11192
RMS 1'+93
RMS 1 '+9~
RMS 1 ~95
RMS 1 ~9b
RM5 1 ~97
RMS 1 ~98
RM5 1 ~99
RHS 15()0
RMS 1501
RI1S 1502
RHO:; 1503
RM5 1 SOli
RMC; 1505
RMS 1506
RMS 1507
RHS 1508
RMS 1509
RMS 1510
RHS 1511
RMS 1512
RHS 1513
RMS 151~
RMS ISIS
RM5 1511>
RM5 1517
RMS 1518
-------
ISN 12'+6
JSN 12'+7
ISN 12'i9
ISN 1250
ISN 1251
ISN 1252
ISN 1253
ISN 125'i
ISN 1255
ISN 1256
ISN 1257
ISN 1258
ISN 1259
ISN 1260
-..)ISN
-.D ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
1261
1262
1263
126'i
1266
1267
126A
1269
1270
127 t
ISN 1272
ISN 1273
ISN 127'i
ISN 1275
ISN 1276
GO TO 133,+8,350) ,ISTOP
350 IF(IPRN3.EQ.O) GO TO 5001
IF (IPRN2) 5001,5001,3,+6
3'i6 IF (LSTRAT-5) 500l,3'f7,5001
3£+7 IKSW = 2
DO 13'+7 K = I,MCNT
Kl = MClTyIK)
I = KCITYIKI)
CALL PRINtll)
13£+7 CONTINUE
lKSW = 0
GO TO 5001
C
2000 CONTINUE
LSM = 2000
C
C
C
C
C
R1S
RMS
RM5
RMS
RMS
RMS
Rt"lS
RMS
RM5
RMS
RMS
RMS
RM5
RMS
Rt"lS
RMS
RMS
Rr1S
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
I 366816.).CBENEF/I0000.)RMS
RMS
RMS
RMS
RMC;
Rr15
RMS
RM5
RMS
CITY INFORMATION C ALL-IND)
EFFECT
ON ECONOMY THROUGH CHANGE IN REGIONAL INCOME
DO 2001 K=I,t-'1CNT
KI = MCITY(K)
I = KCITYCKI)
IF CI.EQ.O) GO TO 2001
IF CLSTRAT-5) 2009,2008,2009
2008 COSTIII,K)=-DELTCCK)oQCII)-DELTRIJYEAR-3).YY(I)
FUI IKI ,JYEAR-3)=FlJI CKI ,JYEAR-3)+DFUI (JYEAR-3).YY( I)
FUIKI,JyEAR-3)=FUIKI,JYEAR-3)+DFUA(JYEAR-3).YYCI)
GO TO 377
2009 COSTCIl,K) =-911.3 . IIYYCI) - XXCI,20»
1 -DELTRCJYEAR-3)eYYII)
377 YRCI) = IYYCl)+COSTCII,K»/YYCI)
TRCI) = 1.0
DNRII) = 1.0
2001 CONTINuF
GO TO 3002
c
c
CITY INFORMATION - ALL IND
1519
1520
t 521
1522
1523
157.£+
1525
1526
1527
1528
1529
1530
1531
1532
1533
153q
1535
1536
1537
153R
1539
15'iO
15'1 1
15'12
15'13
15'1'1
15'15
15'16-
15'17
15'1A
15'19
1550
1551
1552
15C,3
155'i
1555
1556
-------
C EFFECT RMS 1557
C ON ECONOMY THROUGH INCREASED TAXATION RMS 1558
C RMS 1559
ISN 1277 3000 CONTINUE RMS 1560
ISN 1278 LSM = 3000 RMS 156 1
ISN 1279 DO 3001 K=I,MCNT RMS 1562
ISN 1280 K I = MCITYIK) ~ RMS 1563
ISN 12 B 1 I = KCITYIKI) RMS 156~
ISN 1282 I F I I . EQ. 0) GO TO 3001 RMS 1565
ISN 17.8~ I F IIAPROP.EQ.I) GO TO 5050 RMS 156[,
ISN 1286 TRII) = I TlI) - COSTI I I ,K) )/TI I) RMS 1567
ISN 1287 Y R I I) = 1 eO RMS 1568
ISN 1288 DIIIRII) = 1.0 RMS 1569
ISN 1289 GO TO 3001 RMS 1570
ISN 1290 5050 TPDII) = . 05608-APVD I I ) -APV I I) RMS 157 1
ISN 1291 T R I I) = IT ( I ) + TPD I I ) ) I T I I ) RM'5 1572
ISN 1292 Y R I I) = 1.0 R t1S 1573
ISN 1293 DNRII) = 1. 0 RMS 157~
ISN 129~ 3001 CONTINUE RMS 1575
C RMS 1576
C THIRD SET OF SIMULTANEOUS EQUATIONS WHERE RMS 1577
C Y , G, AND C fiRE ENDOGENOUS R t.., '5 1578
ex> C T AND X ARE EXOGENOUS RMS 1579
o RMS 1580
C
ISN 1295 3002 MARK = 3 RMS 1581
ISN 129[, LSM = 3002 RMS 1582
ISN 1297 LVAR = 3 RMS 1583
ISN 129A MCASES = MCNT - LVAR RMS 15R~
ISN 1299 LCHECK = 0 RMS 15AS
ISN 1300 2006 CONTINUE RI1S 1586
ISN 1301 LSM = 2006 RMS 1587
ISN 1302 KPIST = 6 RMS 1588
ISN 1303 MM = 0 RMS 1589
ISN 130Lf t-1 = 0 RMS 1590
ISN 1305 DO 2002 K=I,MCNT RMS 1591
ISN 1306 K I = MCITYIK) RMS 1592
ISN 1307 1 = KCITYIKI) RMS 1593
ISN 1308 IF I Ie EQ. 0) GO TO 2002 RMS 159q
-------
ISN 1310
ISN 1311
ISN 1312
ISN 1313
ISN 1311f
ISN 1315
.ISN 1316
ISN 1317
ISN 1318
ISN 1319
ISN 1320
ISN 1321
ISN 1322
ISN 1323
ISN 1321f
ISN 1325
00 ISN 1326
I-' ISN 1327
ISN 1328
ISN 1329
ISN 1330
ISN 1331
[SN 1332
ISN 1333
ISN 1335
ISN 1336
ISN 1337
ISN 133A
ISN 1339
ISN 13'+0
JSN 13'+1
ISN 13'+2
ISN 13'+3
C
C
2002
C
Y(M+l) = YY( )
Y(M+2) = G(I)
Y(M + 3) = (I)
X(MM+l) = TR(I)
X(MM~2) = YR(I)
X(MM+3)= XX(I .20)
X(Mt.,+,+) = T(I)
X(MM+5) = FUDGE(I,I)
X(MM~5) = 000
X(MM+6) = FUOG(I.Z)
X(MM+6) = 0.0
M = M + 3
MM = MH + 6
CONTINUE
LSM = 200Z
~1S 1595
RMS 1596
RH5 1597
RMS 1598
RMS 1599
RMS 1600
RHS 16nl
RMS 1602
Rt-1S 1603
R 1-1 5 1 6 0 '+
RMS 1605
RMS 1606
RHS 1607
RMS 160R
RMS 1609
RHS 1610
RHS 1611
RMS 1612
RMS 1613
RMS 161'+
RMS 1615
RHS 1616
RHS 1617
RMS 1618
RMS 1619
RMS 1620
RMS 1621
RMS 1622
RMS 1623
RMS 1621f
RMS 1625
RMS 1626
RMS 1627
RMS 167.8
RMS 1629
RMS 1630
RMS 1631
RMS 1632
IF (ICTAB) 8019,8019,1018
1018 DO 952 [H=I,LVAR
952 ARAB(IM) = SCC1(IM)
C
8019 CALL SIMUL
LSH = 8019
C
1019
2003
ZOO If
1'1 = 0
IF(LCHECK)Z003,Z003,2005
DO 200'-+ K= 1 ,MCNT
KI = MCITY(K)
[ = KC[TY(KI)
IF (I.EQ.O) GO 10 2001f
YYO([) = Y(M+l)
GO(I) = Y(M+2)
(0(1) = Y(M+3)
TR(I) = 1.0
YR([) = 1.0
M = 1'1+3
CONTINUE
LSM = 200'f
LCHECK = 1
-------
ISN 13"f"f
ISN 13"f5
ISN 13"f6
ISN 13'+7
ISN 13"f8
ISN 13'+9
ISN 1350
ISN 1351
ISN 1352
ISN 1353
ISN 1355
ISN 1356
ISN 1357
ISN 135R
ISN 1359
ISN 1360
ISN 1361
ISN 1362
ISN 1363
00 ISN 1361.f
,N ISN 1365
I ISN 136b
ISN 1367
ISN 1368
ISN 1369
ISN 1370
ISN 1372
ISN 1373
ISN 1371.f
ISN 1376
ISN 1377
ISN 1379
ISN 1380
C
2005
2007
c
5001
5000
C
C
C
5003
NBEFOR = NBEF03
GO TO 2006
M = 0
NEWCIT = 0
L5M = Z005
NBEF03 = 0
DO 2007 K= 1, MCNT
KI = MCITYCI()
I = KCITYCKI)
IF (IoEQ.O) GO TO 2007
XD(I) = .31.f0'160(YYD(I)-Y(M+l»
M=t1+3
WBDCI) = ALPHA(ZO) .XO(I)
PID(I) = XDCI) .(lo-BETA(20»
QPOLD = PI(I,ZO)
QPNEW = PI(Ip?O)+PID(I)
QPAR =QINV(ZO,1 ).CXIH(I,20)-(D(ZO)-XKHCI,ZO»)
QIOLD =QPAR+QINV(ZO,Z)Q(QPOLO-PIHCI,20) )+QINVC20,3).XKCI,ZO)
QINEW =QPAR+QINV(20,Z).(QPNEW-PIH(J.ZO) )+QINV(ZO.3).XK(I.zn)
XID(I) = QINEW-QIOLD
XKDCI) = XIO(I)/(I.+DCZO»
D2(I) = WBD(I)
CONTINUE
L5M = Z007
GO TO '+000
IF (INTT.EQ.l) NOK = NOK - 1
CONTINUE
CONTINUE
IF (NOK.EQ.O) IHELP = IHELP + 1
NOK = 1
IF(INPT(Z).EQ.Z) NOPEI = NOPEI - 1
ID = NOPEI + INPT(I)
IF(ID.EQ.O) GO TO 700
~1S 1633
RMS 163'1
RMS 1635
RMS 1636
RM5 16.17
RI1S 163A
RMS 1639
RMS 16'+0
RMS 16'11
RMS 16'+Z
RMS 16'13
RMS 16'1,+
R 11 S 1 6 Ii 5
RMS 16'+6
RMS 16'17
RMS 16'18
RMS 16'19
RMS 1650
R r... S 1 6 5 1
RMS 1652
RMS 1653
RMS 165l.f
RMS 16S5
RMS 1656
RMS 1657
RMS 1658
Rf1S 1659
RMS 1660
RMS 1661
RMS 1667.
RMS 1663
RMS t 66Lf
RMS 1665
RMS 1666
RMS 1667
RMS 1668
RMS 1669
RMS 1670
-------
ISN 1382
ISN 1383
ISN 138'i
ISN 1385
ISN 1386
ISN 1387
ISN 138A
I 00
IW
INTT = INPT(2)
GO TO 10
95 I = 1000000
IS = ICHO(I)
STOP
600~ FORMAT(IHl,' INDUSTRY-',I3,' DOES NOT EXIST FOR CITY-',I3)
END
RMS 1671
RMS 1672
RMS 1673
RMS 167~
RMS 1675
RMS 1676
RMS 1677
.~
-------
LEVEL 20.1
ISN 0002
ISN 0003
ISN ODOlf
ISN 0005
I5N 0006
I5N 0007
ISN OOOB
00
~
(MAY 71)
05/360
FORTRAN H
COMPILER OPTIONS - NAME= MAJN,OPT=OO,LJNECNT=LfO,SIZE=OOOOK,
SOURCE,BCO,NOLJST,NODECK,LOAD,MAP,NOEOIT,NOID,NOXREF
SUBROUTINE PHELP{XX,XD,XIP)
IFIXX) 10dO,B
8 XJP = {XD/XX)$100.
RETURN
10 XIP = 0.0
RETURN
END
RMS 1678
RMS 1679
RMS 1680
RMC; 16AI
RMS 1687.
RMS 16B3
RMS 16ALf
DATE
72.22'
-------
LEVEL 20.1 (MAY 71)
05/360
FORTRAN H
COMPILER OPTIONS - NAME= MAIN,OPT=OO,LINfCNT=~O,SIZE=OOOOK,
SOURCE,BCO,NOLIST,NODECK,LOAD,MAP,NOEDIT,NOID,NOXREF
SUBROUTINE PRINl (JCtTY)
ISN 0002
ISN 0003
ISN OOOtt
ISN OOOS
IsN 0006
ISN 0007
ISN 0008
ISN 0009
IsN (1010
ISN 0011
ISN 0012
ISN 0013
ISN 001L1
ISN 001S
ISN 0016
ISN 0017
ISN 001A
~ ISN 0019
ISN 0020
ISN 0021
ISN 0022
ISN 0023
ISN OOZLI
ISN 0025
ISN 0026
ISN 0027
IsN OOZR
ISN 0029
ISN 0030
c
c
C
C
C
C
REAL KJNDX,LJNDX
DIMENSION (92),
DIMENSION DNBAR(92)
DIMENSION DNRARN(92),
DIMENSION G(92),
DIMENSION IAQCR(92,20),
DIMENSION LARELCI0,20),
DIMENSION PIC92,20),
DIMENSION PPID(91),
DIMENSION PXKD(91),
DIMENSION QBARN(92),
DIMENSION QCN(92),
DIMENSION QMN(9Z),
DIMENSION QTN(92),
DIMENSION SUMNN(9Z),
DIMENSION TN(92),
DIMENSION UN(92),
DIMENSION XDN(92),
DIMENSION XIN(9Z),
DIMENSION INPT(2)
DIMENSION XKN(92),
DIMENSION XLN(92),
DIMENSION XNN(92),
DIMENSION YY(92),
DIMENSION TITLE(20,Z)
DIMENSION PERUN(92)
DIMENSION TAf9Z) ,TAN(9/.1 ,TAD(92) ,TO(92) ,TON(92) ,TOD(92) ,TP(92),
1 TPN(92) ,TPD(92)
DIMENSION OTO(92) ,OTA(9Z) ,OTP(92) ,OAPV(9Z) ,OPOP(92) ,POPGR(91)
CD(92) ,
CN(92)
EOLD(9l,20,S) ,
GO(92) ,
INDCI(92,20),
MCITY(91),
PID(92) ,
PXD(91),
QF3AR(92),
QC(92),
Q ~, ( 9 2) ,
QT(92) ,
SUM ~J ( 9 2 , 2 0 ) ,
T ( 92) ,
U (921 ,
WN(92) ,
XI(92,20),
XK(92,20),
EO(91,5)
GN(92)
KCITY(100)
PDtH91 )
PIN(92)
PXID(91)
QBARD(92)
QCD(92)
QMD(92)
QTO(92)
SUt.,ND(92)
TD(92)
UD(92)
XD(92)
XID(92)
XKO(92)
XL(92) ,
XNB,\R(92),
XNT(92) ,
YYDI92> ,
XLD(92)
XND(92)
XX(92,20)
YYN(92)
RMS 168S
RMc; 16A6
RMS 16A7
RMS 1688
RMS 16B9
RMS 1690
RMS 1691
Rr.,s 1692
RMS 1693
RMS 169L1
RMS 1695
RMS 1696
RMS 1697
RMS 1698
RM5 1699
RMS 1700
Rt1S 1701
R t1 5 1 7 0 2
RMS 1703
RMS 170q
RMS 1705
RMS 1706
Rt1S 1707
R t-1 S 1 7 0 R
RMS 1709
RMS 1710
RMS 1711
RMS 1717-
RMS 1713
R r1 S 1 7 1 LI
RMS 1715
RMS 1716
RMS 1717
RMS 171A
Rt-1S 1719
RMS 1720
DA"E
72.22
-------
,--'-~---'-
ISN
ex> ISN
O'ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
ISN
SN 0031
ISN 0032
ISN 0033
ISN 003'i
ISN 0035
ISN 0036
ISN 0037
ISN 0038
ISN 0039
ISN 00'i0
ISN 00'i1
00'i2
00'+3
00'i'i
00'i5
00'i6
00'i8
00'+9
0050
0051
0052
0053
005'i
ISN 0055
ISN 0056
ISN 0057
ISN a05A
ISN 0059
ISN 0060
COMMON /A"AXI T;','PD,'O,TO),TA,'AD
COMMON/CPRIN/PXD,PXID,PXKD,PPID,PDN,EOLD,ED
COMMON/BUG/ISBUG,LREAD,LWRIT
COMMON /APRIN/PI,PID,XI,XID,XK,XKD,XX,XD,C,CD,yy,YYO,SUMN,SUMNO,
1 XNBAR,DNBAR,XNT,XND,U,UD,XL,XLD,G,GD,T,TD,QT
2 ,QTD,QC,QCD,QM,QMD,QBAR,QBARD
COMMON / ZEROT I PIN,XIN,XKN,XDN,YYN,GN.TN,CN,ONBARN,XNN,XLN,WN,
1 UN,SUMNN,QBARN,QMN,QCN.QTN
2 .TPN,TON,TAN
COMMON/NPRINI LABEL,MCITY,KCITY,
COMMON IPRITTI TITLE.AENEF
COMMON / IPRIN IINTT,LPNTT ,JYEAR,IAQCR,IJCIT,MIN,IA312
C~r"'MON/EVERY/OWB(92) ,OW(9Z) ,OX(92) ,01 (92) ,OY(92) ,OT(92) ,OG(92),
I OC(92),OXBAR(92),OQC(92),OQM(92),OQB(92), GRFAC,XINC,VK,
2 YGROW.IZS
3 ,OTP,OTA,OTO,OPOP,OAPV,POPGR
COMMONI WNEW~I WAGECO(20,3)
COMMON/ALET/PERXX.PERW.PERWB,PERSN,PERP,PERXI,PERXK,PERYY,PERNB,
1 PERT,PERG,PERC,PERQM.PERQB,PERQC,PERTP,PERTO.PERTA
EQUIVALENCE (INDCI(I),XX(I»
I = I C I T Y
WRITE(LWRIT,IOO)
LPYEAR = 1969 + JYEAR
IF(ICITY.EQ.92) GO TO 760
IF(IKSW-l) 1,750,750
1 WRITE (LWRIT,II0) LSTRAT,(LABEL(LSTRAT,J), J = 1,20)
~RITE(LWRtT,120) (IAQCR(I,J),J = 1,20) ,LPYEAR
ViRITE(LWRIT,130)
WR J TF: (LWR t T, 1 'in)
WRITE(LWRtT,150)
WRITE(LWRIT,160)
c
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
LSTRAT,INPT,IKSW,IHELP,IFEEDRMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
Rr..,s
RMS
RHS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
CALL OLETS(I)
PERXN = PERSN + PERNS
PERQT = PERQC + PERQM + PERQB
PERXL = PERXN I .96
CALL PHELp(PI(I,20),PIO(I) ,PIP)
CALL .PHELp ( X I ( I ,20) ,X T 0 ( I ) ,X J P)
1721
1722
1723
172'i
177.5
1726
1727
1728
1729
1730
173 I
1732
1733
173'i
1735
17 36
1737
1738
1739
17'i0
17'+ 1
17'i2
17'i3
17'i'i
17'+5
17'i6
17'i7
17'iA
17'i9
1750
17; 1
1752
1753
175'i
1755
1756
1757
1758
-------
S~ 006
ISN 0062
ISN 0063
ISN 006'+
ISN 0065
ISN 0066
ISN 0067
ISN 0068
ISN 0069
ISN 0070
ISN 0071
ISN 0072
ISN 0073
ISN 007'+
ISN 0075
ISN 0076
ISN 0077
ISN 0078
ISN 0079
00 ISN 0080
--.J ISN 0081
ISN 0082
ISN 0083
ISN 008'+
ISN on85
ISN 0086
ISN 0087
ISN 0088
ISN 0089
ISN 0090
ISN 0091
ISN 0092
ISN 0093
ISN 009'f
ISN 0095
ISN 0096
IFIU(I» 58,58,56
56 UP = UDII) lUll)
PERUX = '+.0
GO TO 59
58 PERUX = '+.0
UP = 0.0
59 CONTINUE
UX = UI1>.100.
UXD = lJD(I).100.
ViRITEILWRIT,190)
WRITE(LWRIT,170)
WRITE(LWRIT,180)
WRITE(LWRIT,200)
WRITE(LWRIT,210)
WRITE(LWRIT,160)
WRITEILWRIT,2'+0)
iI,R I TE (UVR IT dbO)
WRITE(LWRIT,260)
WRITE(LWRIT,160)
CALL
CALL
CALL
CALL
CALL
CALL
CALL
CALL
CALL
CALL
CALL
CALL
CALL
CALL
CALL
CALL
CALL
C
C
PHELp X((I,20),XKDII),XKP)
PHELpIXXII,20) ,XDII) ,XXP)
PHELp (C(I),CDII),CP)
PHELpIYY(I),YYDIII,YYPI
PHELPISUMN( I ,20) ,SUMNOI I) ,SUMNP)
PHELp(XNBAR(I),DNBARII),PNBAR)
PHELp(XNTII) ,XNDII),XNP)
PHELpIXLII) ,XLD(I) ,XLP)
PHELPITII) ,TDII),TPQ)
PHELPITPII),TPDII),TPP)
PHELpITO(I),TODII),TOP)
PHELp(TAII),TADII) ,TAP)
PHELp(QCII) ,QCD(I) ,QcP)
PHELPIG(I),GD(I) ,GP)
PHELp(QM(I),QMDIII,QMPI
PHELp(QBARIII,QBARD(I),QBARP)
PHELp(QTII),QTD(II,QTP)
R~S 759
RMS 1760
RMS 1761
RMS 17b2
RMS 1763
RMS 17bq
RMS 1765
RMS 1766
RMS 1767
RMS 17b8
RMS 1769
RMS 1770
RMS 1771
RMS 1772
RMS 1773
RMS 177,+
RMS 1775
RMS 1776
RMS 1777
RMS 1778
RMS 1779
RMS 1780
RMS 1781
RMS 1782
RMS 1783
RMS 178q
RMS 1785
RMS 1786
RMS 1787
RMC; 17AA
RMS 1789
RMS 1790
RMS 1791
RMS 1792
RMS 179)
RMS 179q
RMS 1795
RMS 1796
PERXX,XX( ICITY,lO) ,XDI ICITY> ,XXP
PERP ,PI I ICITY,lQ) ,PIOI ICITY) ,PIP
PERXI,XIIICITY,20),XIDIICITY),XIP
PERXK,XKI ICtTy,lO) ,XKO( ICITY> ,XKP
PERSN,SUMN( ICITY,20) ,SlJMND( ICITY> ,SUMNP
PERNB,XNBAR( ICITY> ,DNBAR( ICITY> ,PNBAR
-------
ISN 0097
ISN 0098
ISN 0099
ISN 0100
ISN 0101
ISN 0102
ISN 0103
ISN 010'1
ISN 0105
ISN 0106
ISN 0107
ISN 0108
ISN 0109
ISN 0110
ISN 0111
ISN 0112
ISN 0113
ISN 011'1
ISN 0115
ISN 0116
i~ISN 0117
ISN 0118
ISN 0119
ISN 0120
ISN 0121
ISN 0122
ISN 0123
ISN 012'1
ISN 0125
ISN 0126
ISN 0127
ISN 0128
C
WRITE(LWR "1300) PERYY,YY(IC "Y),YYD( C "Y>,YY:)
~J R I T E ( L W R IT, 2 9 0) PER C , C ( I C I T Y ) ,C D I I C I T Y ) , C P
~RITEILWRIT,310) PERXN,XNTIICITY),XNDIICITY),XNP
WRITEILWRIT,320) PERUX,UX,UXD,UP
WRITEILWRIT,330) PERXL,XLI ICITY) ,XLD( ICITY) ,XLP
IVRITE(L'JIIRITt3'10) G(ICITY>, GO(ICITYJ, GP
~jRITE(LWRITt350) PERT,T( ICITY> ,TDI ICITY> ,TPQ
VlRITE(LWRIT,351) PERTP,TP( ICITY) ,TPDI ICITY) ,TPP
WRITEILWRIT,3S7.) PERTO,TO(ICITY) ,TOO(ICITY),TOP
WRITE(LWRIT,3S3) PERTA,TAIICITY) ,TADIICITY),TAP
WRITE (LWRIT,35S)
WRITFILWRITt360) PERQT,QTI ICITY) ,QTDI ICITY) ,QTP
WR.ITEILWRIT,220) PERQM,QM( tCITY) ,QMD( ICTTY) ,QMP
WRITEILWRIT,270) PERQ8,QBAR( ICITY) ,QBARO( ICITYI ,QBARP
WRITF.:(LWRIT,370) PERQC,QCI ICITY) ,QCD( ICITY) ,QCP
100 FORMAT"IIHl)
110 FORMATIIIIIH ,'CONTROL STRATEGY NO. '11,' - ',20A'I,/)
120 FORMAT( IH ,'AQCR ',13,IX,19A'I,3X,' FOR ',1'1)
130 FORMAT(IH ,7'1X,'WITHOUT',7X,' WITH',7X,'NET',AX,'PERCENT')
1'10 FORMATIIH ,7'1X,'CONTROL',7X,'IT-I)',oX,'CHANGE',6X,'CHANGE')
ISO FORMAT(IH ,'X,'MANUFACTURING INDUSTRIES')
160 FORMAT(IH )
170 FORMAT(IH ,6X,'PROFIT IMILLIONS)','8X,IX,3IFIO.3,2X),FI0.')
180 FORMATIIH ,6X,'INVESTMENT (MILLIONS)',''X,IX,3IFI0.3,2X),
1801FI0.')
190 FORMATIIH ,6X,'VALUE ADDED (MILLIONS)',''IX,3IFIO,3,2X),FIO.')
200 FORMAT(IH ,6X,'CAPITAL STOCK IMILLIONS)','2X,3(FI0.3,2X),FI0.'I)
210 FORMATIIH ,6X,'EMPLOYMENT (1000 S)','6X,IX,3IFI0.3,7.X) ,FIO.')
7.20 FORMAT(IH ,6X,'ELECTRICITY USED BY MANUFACTURING INDUSTRIES 110 M
lKWH)',IIX,3(FI0.3,2X) ,FI0.')
2'10 FORMAT(IH ,'X,'OTHER INDUSTRIES')
260 FORMATIIH ,6X,'EMPLOYMENT (1000 S)','6X,IX,2IFI0.2,2X),FI0.3,2X,
1 F 1 0 . 'I )
270 FORMATIIH ,6X,'ELECTRICITY USED BY OTHER INDUSTRIES (10 M KWH)',
1 19X,3(FI0.3,2X),FI0.'I) .
290 FORMATIIH ,'X,'REGIONAL CONSUMPTION (MILLIONS)',3'1X,IX,FI2.3,2X,
29012(FI0.3,2x),FI0.')
300 FORMATIIH ,'X,'TOTAL PERSONAL INCOME FOR THE REGION (MILLIONS)'
~MS 1797
RMS 1798
RMS 1799
RMS 1800
RMS 1801
RMS 1802
RMS 1803
RMS 180Q
RMS 180S
RMS 1806
RMS 1807
RMS 180B
RI-1S 1809
RMS 1810
RMS 1811
RMS 1812
RMS 1813
RMS 181'1
RMS 181S
RMS 1816
RMS 1817
RMS 1818
RMS 1819
RMS 1820
RMS 1821
RMS 1822
RMS 1823
RMS 182'1
RMS 1825
RM5 1826
RMS 1827
RMS 187.8
RMS 1829
RMS 1830
RMS 1831
RMS 1832
RMS 1833
RMS 183'1
-------
ISN 0129
300 .9X,F 2.3,2X,2(;" 0.3,2X),F O.j) R'1S
310 FORMAT(IH ,~X,'TOTAL REGIONAL EMPLOYMENT (1000 S)',33X,IX,3(FI0.3,RMS
31012X) ,FIO.q) RMS
320 FORMAT(IH ,qX,'REGIONAl UNEMPLOYMENT (PERCENT)',36X,IX,2(FI0.3,2X)RMS
3201,FI0.q,2X,FIO.q) RMS
330 FORMAT(IH ,~X,'TOTAL LABOR FORCE (1000 S)',qlX,IX,3IFlO.3,2X), RMS
3301FIO.q) RMS
3QO FORMAT(IH ,qX,'GOVERNMENT EXPENDITURE FOR THE REGION (MILLIONS)', RMS
3QOlI9X,lX,3(FI0.3,2X) ,FI0.Q) RMS
350 FORMAT(IH ,qX,'GOVERNMENT REVENUE FROM THE REGION (MILLIONS)',22X,RMS
35011X,3IFI0.3,2X) ,FI0.q) RMS
351 FORMAT(IH ,qX,'GOVERNMENT REVENUE FROM PROPERTY TAXES (MILLIONS)',RMS
11~X,lX,3(FI0.3,2X) ,FI0.Q) RMS
352 FORMAT(IH ,qX,'GOVERNMENT REVENUE OTHER THAN PROPERTY TAX IMILLIONRM5
lX)',IQX,IX,3IFI0.3,2X) ,FIO.q) RMS
353 FORMATIIH ,QX,'INTRAGOVERNMENT AID TO THE REGION (MILLIONS)',23X, RMS
1 lX,3(FI0.3,2X),FI0.q) RMS
355 FORMAT (/,IH ,qX,'ELECTRIC POWER DEMAND' ,I) RMS
360 FORMATIIH ,6X,'TOTAL ELECTRIC CONSUMPTION FOR THE REGION liD M KWRMS
3601S)',13X,3(FI0.3,2X),FIO.Q) RMS
370 FORMATIIH ,6X,'REsIDENTIAL CONSUMPTION IN THE REGION (10 M KWH)',RMS
370117X,)IFI0.3,2X),FI0.4) RMS
380 FORMAT(IH ,'SIC DETAIL FOR MANUFACTURING INDUSTRIES ACTIVE IN THISRMS
3801 AQCR'/) RMS
qOO FORMATIIH ,'SIC DETAIL: NO DETAIL AVAILABLE FOR THIS AQCR') RMs
Q 51 FORMA T (1 H ,'[[[ RM S
ISN 0130
ISN 0131
ISN 0132
IsN 0133
ISN 0134
IsN 0135
ISN 0136
ISN 0137
IsN 013A
~ ISN 0139
ISN 01qO
ISN 01ql
ISN 01q2
ISN 01Q3
1 . . . . . . . 0 . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . RMS
q52 FORMAT I1H ,'...... -VALUE ADDED- . -INRMS
IVESTMENT- . .CAPITAL STOCK- .' RMS
q53 FORMA T (1 H ,'. SIC................................................ RMS
ISN 01QQ
ISN 0145
1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .' ) RMS
45Li FORMAT 1 I H ,'.COOE. NO CONTROL. NET . PERCENT . NO CONTROL. RMS
1 NET . PERCENT . NO CONTROL. NET . PERCENT . ' ) RMS
LiS5 FORMAT (I H ,'...... (MILLIONS). CHANGE . CHANGE . (MILLIONS). RMS
1 CHANGE . CHANGE . IMIlLIONS). CHANGE . CHANGE . ' ) RM5
Q56 FORMAT 1 1 H t'...... -PROFIT- . -NO. ORMS
1 F EMPLOyEES- . -ELECTRIC POWER CONSUMPTION- .' ) RMS
Q57 FORMAT 1 1 H ,'...... (MILLIONS). CHANGE . CHANGE . llano 5) . RMS
ISN 0lq6
ISN 01Q7
ISN 01Q8
1835
1836
IB37
1838
18.19
IBLiO
I B Q 1
IB~2
IBQ3
1 8q 4
IB~S
18'+/)
ISQ7
18qa
1849
IB50
185 1
18S2
IB53
185Q
lass
18S6
1857
1858
18S9
1860
1861
1862
1863
-------
1 CHA'GE . CHNG:: . I 0 ~ (WH). CHA'GE . CHANG:: .'
IFIINTT.EQ.2.AND.MIN.EQ.O,AND.IA312.EQ.I) GO TO ~61
IFIINTT.EQ.2.AND.MIN.EQ.O) GO TO ~60
IFIIHElP.EQ.2) GO TO 500
RETURN
~60 WRITE ILWRIT,160)
WRITE IlWRIT,lbO)
WRITE (L~RIT,~OO)
RETURN
~61 WRITEILWRIT,~62)
~62 FORMATIIIIH ,'NO DETAilED SIC AVAILABLE FOR
1 COSTS USED INSTEAD'I' IMpACT TRANSFERED TO
RETURN
500 WRITE ILWRIT,100)
WRITE ILWRIT,160)
WRITE IlWRIT,~SI)
WRITE ILWRIT,Q52)
WRITE (LWRIT,~53)
WRITE ILWRIT,q5~)
WRITE ILWRIT,QS5)
WRITE IlWRIT,~51)
530 FORMATIIH ,IH.,13,2H .,)(IX,F9.3,2H .,F9.3,2H .,F9.~,2H .))
DO 600 JND = 1,19
IF IINDCIII,JND)) 510,600,510
510 INDI = JNO+20-II/JNO)
XXP =IPXDIJND) IXXII,JND)).100.
IF IXI (I.JND)) 515,515,520
515 XIP = 0.0
GO TO 525
520 XIP = IPXIDIJNO) IXI(I,JNO)).100.
525 IF IXKII,JND)) 527,527,535
527 XKP = 0.0
GO TO 5~O
535 XKP = (PXKO(JND) IXKII,JND)).100.
5~0 WRITE (LWRIT,530) INDI ,XXI I ,JNO) .PXDIJND)
IPXIDIJND) ,XIP,XKII,JNO),PXKDIJND) ,XKP
600 CONTINUE
WRITEILWRIT,100)
ISN 01Lf9
ISN 0151
ISN 0153
ISN 0155
ISN 0156
ISN 0157
ISN 0158
ISN 0159
ISN 0160
ISN 0161
ISN 0162
ISN 0163
ISN 016~
ISN 0165
ISN 0166
ISN 0167
ISN 0168
ISN 0169
ISN 0170
ISN 0171
~ ISN 0172
ISN 0173
ISN 017~
ISN 0175
ISN 0176
ISN 0177
ISN 0178
ISN 0179
ISN 0180
ISN 0181
I . ISN 0182
ISN 0183
ISN 018~
ISN 0185
ISN 0186
R~ 5
RMS
RMS
RMS
RMS
RMS
RM'5
RMS
RMS
RMS
THIS AQCR - TOTAL MANURMS
REGIONAL ANALYSIS') RM5
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
,XXP,XIII,JND),
1873
187Lf
1875
1876
1877
1878
1879
1880
1881
1882
1883
188~
1885
1806 .
1887
1888
1889
1890
189 I
1892
1893
189~
1895
1896
1897
1898
1899
1900
1901
1902
1903
190~
1905
1906
1907
190B
1909
1910
-------
WRITEILWRIT,160)
WRITEILWRIT,,+SII
WRJTEILWRJT,IiS6)
WRJTEILWRJr,~53)
WRITEILWRIT,,+S,+I
WRITEILWRIT,l.fS71
WRITEILWRIT,,+Sl )
DO 700 JND = 1,19
IF IINDCIII,JNDII 610,700,610
610 INDI = JND~20-II/JNDI
IF IPIII,JNDI) 615,615,620
615 PIP = 0.0
GO TO 625
620 PIP = IPPIDIJND) IPIII,JND).100.
625 IF IEOLDII,JND,5») 635,630,635
630 PCT = 0.0
GO TO 6,+0
635 PCT = IEOIJND,5)/EOLOII,JND,5».100.
6'+0 IF ISUMNII,JND» 6,+5,6,+5,650
6'+5 SUMNP = 0.0
GO TO 655
650 SUMNP = IpDNIJND) ISUMNII,JND» -100.
655 WRITE ILwRIT,S30) INDI,PIII,JND),PPIDIJND),PIP,SUMNII,JND),
1 PDNIJND) ,SUMNP,EOLDI I ,JND,S) ,EDIJND,S) ,PCT
700 CONTINUE
IF IINTT.NE.2.0R.LPNTT.NE.2) RETURN
WRITEILWRIT,701 )
701 FORMATIIH1,'OOPS-TRIED TO ACCESS FUEL PRINTOUT IN
RETURN
750 GO TO 1760,770),IKSW
760 BENT = SENEF/I000.
WRITEILWRIT,911) IITITLEIIMI,JMI),IMI=1,20),JMI=I,2)
WRITE ILWRIT,900)
GO TO 790
770 WRITE ILWRIT,910)
790 WRITE ILWRITtl20) IIAQCRII,J),J=I,20) ,LPYEAR
WRITEILWRIT,130)
I"iRITEILWRIT,I'+OI
ISN 0187
ISN 0188
ISN 0189
ISN 0190
ISN 0191
ISN 0192
ISN 0193
ISN 019ii
ISN 0195
ISN 0196
ISN 0197
ISN 0198
lSN 0199
ISN 0200
15N 0201
ISN 0202
ISN 0203
ISN 0201.1
ISN 0205
ISN 0206
-.D ISN 07.07
~.tSN 0208
ISN 0209
ISN 0210
ISN 0211
ISN 0213
ISN 021'+
ISN [1215
ISN 0216
ISN 0217
ISN 021A
ISN 0219
ISN 0220
ISN 0221
ISN 0222
ISN 0223
ISN 022'+
HS
RMS
RMS
RM5
RMS
RMS
RMS
RMS
RMS
Rr1S
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RI1S
R r1S
RMS
RMS
RMS
RMS
RMS
PRINl,RMSOI758')RMS
RMS
RMS
RMS
RMC)
RMS
RMS
RMS
RMS
RMS
RMS
191 1
1912
1913
1911f
1 915
1 9 1 6
191 7
1918
1 9 1 9
1920
192 1
1922
1923
1921f
1925
192h
1927
192R
1929
1930
1931
1932
1933
193'+
1935
1936
1937
1938
1939
19'+0
19'+ 1
191.12
19l.f3
19'+t.t
19145
19146
19'+7
19l.fA
-------
ISN 0225
ISN 0226
ISN 0227
ISN 0228
ISN 0229
ISN 0230
ISN 0231
ISN 0232
ISN 0233
ISN 023Lf
ISN 0235
ISN 0236
ISN 0237
ISN 07.38
ISN 0239
ISN 02'iO
ISN 02Lfl
ISN 02"12
ISN 02Lf3
-.D ISN 02Lf'i
N ISN 02"15
ISN n2'i6
ISN 02"17
ISN 02'i8
ISN 02'i9
ISN 0250
ISN 0251
ISN 0252
ISN 0253
ISN 025'i
ISN 0255
ISN 0256
ISN 0257
ISN 0258
WR T::(LWRI", 50)
W R I T E ( L ';1 R IT, 1 6 n )
~MS 19'i9
RMS 1950
RMS 1951
RMS 1952
RMS 1953
RMS 195'i
RMS 1955
RMS 1956
RMS 19S7
RMS 1958
RMS 1959
RMS 1960
RMS 1961
RMS 1962
RMS 1963
RMS 196"1
RMS 1965
RMS 1966
RMS 1967
RMS 1968
RMS 1969
RMS 1970
RMS 1971
RMS }972
RMS 1973
RMS 197'i
RMS 1975
RMS 1976
RMS 1977
RMS 1978
RMS 1979
RMc; 1980
RMS 1981
RMS 1982
RMS 1983
RMS 198'i
RMS 1985
RMS }986
C
C
CALL OLETS ( I )
PERXN = PERSN + PERNS
PERQT = PERQC + PERQM + PERQB
PER XL = PfRXN / .96
CALL PHELp(PIII,20),PINII),PIP)
CALL PHELp(XIII,20) ,XINII),XIP)
CALL PHELPIXKII,lO) ,XKNII),XKP)
CALL PHELPIXXII,20),XDNII),XXP)
CALL PHELp ICII) ,CNII),CP)
C A'L L PH E L P ( Y Y I I ) , Y Y N I I ) ,Y Y P )
CALL PHELp I SUMN ( I ,20) ,SUMNN I I ) ,SUMNP)
CALL PHELpIXNBARI I),DNBARNII) ,PNBAR)
CALL PHELPIXNTII),XNN(I) ,XNP) .
CALL PHELp(XLII) ,XLNII) ,XLP)
CALL PHELp(GII),GNII) ,GP)
CALL PHELPIT(I),TNII) ,TPQ)
CALL PHELPITPII),TPNII),TPP)
CALL PHELpITAII) ,TANII) ,TAP)
CALL PHELp(TOII),TONII) ,TOP)
CALL PHELPIQCII),QCN(T) ,OCP)
CALL PHELPIQMII),QMNII),QMP)
CALL PHELplQBARI I) ,QBARNI I),QBARP)
CALL PHfLPIQTII) ,QTNII),QTP)
c
c
85'i IF (UII» 858,858,856
856 UP = IUNII) / UI I» . Ion.
GO T0860
858 UP = 0.0
860 CONTINUE
PERUN(I) = .O'i . II. + IUP/l00.»
PERUX = 'i.O
UX = UII)-100.
UXQ = UNII).100.
-------
ISN 0259
ISN 0260
ISN 0261
ISN 0262
ISN 0263
ISN 026Lf
ISN 0265
ISN 0266
ISN 0267
ISN 0268
ISN 0269
ISN 0270
ISN 0271
ISN 0272
ISN 0273
ISN' 027Lf
ISN 0275
:~ ISN 0276
ISN 0277
ISN 0278
ISN 0279
ISN 0280
ISN 0281
ISN 0282
ISN 0283
ISN 028Lf
ISN 0285
ISN 0286
ISN 0287
ISN 0288
. ' )
ITONS)'
)
C
C
C
WRITEILWRIT,190)
WRITEILWRIT,170)
WRITEILWRIT,IS0)
WRITECLWRIT,200)
WRITEILWRIT,210)
WRITEILWRIT,160)
ViRITEILWRIT,2LfO)
WRITECLWRIT,160)
WRITECLWRIT,260) PERNB,XNBARIICITY),DNBARNIICITY),PNBAR
WRITECLWRIT,160)
WRITFCLWRIT,300) PERYY,YYIICITY),YYNIICITY),YYP
\'IRJTEIL'NRIT,290) PERC,CI ICITY) ,CNI ICITY) ,CP
WRITECLWRITdI0) PERXN,XNTI ICITY) ,XNNI ICITY) ,XNP
WRITECLWRIT,320) PERUX,UX,UXD,UP
WRITEILWRIT,330) PERXL,XLI ICITY) ,XLNI ICITY) ,XLP
WRITECLWRIT,3LfO) GIICITY), GNIICITY), GP
WRITECLWRIT,3S0) PERT,TI ICITY> ,TNC ICITY) ,TPQ
WRITECLWRIT,351) PERTP,TP( ICITY) ,TPNI JCITY) ,TPP
WRITEILWRIT,352) PERTO,TOI JCITY) ,TONI ICITY) ,TOP
WRITE(LWRITd5:\) PERTA,TAI ICITY) ,TANI ICITY) ,TAP
WRITE ILWRIT,3SS)
I',RITEILWRITd6IJ) PERQT,QTI ICITY) ,QTNC ICITY) ,QTP
WRITECLWRIT,220) PERQM,QMIICITY),QMNC ICITY) ,QMP
WRITECLWRIT,270) PER(~B,QP,ARC ICITY> ,QBARNI ICITY> ,QBARP
WRITECLWRIT,370) PERQC,QCI ICITY) ,QCNI ICITY) ,QCP
PERXX,XXI ICITY,ZO) ,XDNC ICITY> ,XXP
PERP ,PI C ICITY,20) ,PINC ICITY) ,PIP
PERXI,XICICITY,20) ,XINCICITY) ,XIP
PERXK,XKC ICITY,20) ,XKNC ICITY) ,XKP
PERSN,SUMNI ICITY,lO) ,SUMNNI ICITY) ,SUMNP
C
C
911 FORMATIIH ,112C20X,20ALf,/»
900 FORMATCIII,lH ,'TOTAL NET EFFECT
1 IN THIS RUN',/)
910 FORMATIIII,lH ,
1000 FORMAT CIH ,'......
ITITY OF COAL- .
1005 FORMATIIH ,'......C$ PER
lCHANGE . CHANGE .IPER
1010 FORMAT CIH ,'......
OF ALL CONTROL STRATEGIES
, EFFECT OF 3121A) CONTROL COSTS',/)
-COST OF COAL-
-COST OF OIL-
TON). CHANGE. CHANGE.
M GAL). CHANGE. CHANGE .'
-QUANTITY OF OIL-
.
~ \ 5
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
Rr~s
RMS
RMS
RMS
RM5
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
Rr-1S
RMS
PURSUEDRMS
RMS
RMS
-QUANRMS
RMS
. RMS
RMS
-CORMS
1987
1988
19A9
1990
1991
1992
1993
199'1
1995
1996
1997
199R
1999
2000
2001
20n2
7.003
20rJ'I
2005
2006
2007
7.00A
7.009
2010
201 1
7.012
2013
201'1
7.015
2016
2017
7.01R
2019
2020
7.021
2022
7.023
202/;>
-------
5" OF GAS- . -QUANTITY OF GAS- . ' ) R'S 2025
ISN 0289 1015 FORMAT(IH ,'...... (1000 GAL). CHANGE . CHANGE .(PER ClI.FT). RMS 2026
lCHANGE . CHANGE .(lOOOCU.FT). CHANGE . CHANGE . ' ) RMS 2027
ISN 0290 1020 FORMAT (1 H t'...... -COST OF ELECTRICITY- . -CORMS 7.028
1ST OF COKE- . -QUANTITY OF COKE- . ' ) RMS 7.029
ISN 0291 1025 FOR i1 A T ( 1 H ,0...... (PER M KWH). CHANGE . CHANGE (PER TON) . RMS 2030
lCHANGE CHANGE . (TONS) . CHANGE . CHANGE . ' ) RMS 7.0:q
ISN 0292 RETURN RMS 2032
ISN 0293 END RMS 2033
-.D
~
-------
LEV:- ZO.1
lSN OOOZ
ISN 0003
lSN OOOtf
ISN 0005
ISN 0006
ISN 0007
ISN OOOB
ISN 0009
ISN 0010
ISN 0012
ISN OOltf
ISN 0015
ISN 0016
ISN 0017
-.D ISN 0019
\.}1 ISN 0021
ISN 0022
ISN 0023
ISN OOZtf
ISN 0026
ISN 0027
ISN 0029
ISN 0030
ISN 0032
ISN 0033
ISN 003tf
ISN 0035
ISN 0036
ISN 0037
MAY 711
OS/360
fOR. HN H
COMPILER OPTIONS - NAM:= MA N,OP "=OO,LINECNT=qO,5 ZE=OOOOK,
SOURCE,BCD,NOLIST,NOOECK,LOAD,MAP,NOEDIT,NOID,NOXREF
SUBROUTINE OLETS(I)
REA L K J N D X , L J ~JO X
DIMENSION GRFAC(100,20,3),XINC(60,3) ,VK(20),YGROW(3B)
DIMENSION OTP(9Z1 ,OTA(9Z) ,OTO(92) ,OPOP(92) ,OAPV(92)
C n f111 0 N / EVE R Y / 0 W B ( 9 2 ) ,0 'N ( 9 Z ) ,OX ( 9 2 ) ,0 I ( 9 2 ) , 0 Y ( 9 2 ) , 0 T ( 9 2 ) ,0 G ( 9 2 ) ,
1 OC(92) ,OXBAR(9Z1 ,OQC(92) ,OQM(921 ,OQB(97.), GRFAC,XINC,VK,
2 YGROW,IZS
3 ,OTP,OTA,OTO,OPOP,OAPV,POPGR
COMMON/ WNEWW/ WAGECO(ZO,3)
COMMON /IPRIN/ INTT,LPNTT,JYEAR
COMMON/ALET/PERXX,PERW,PERWB,PERSN,PERP,PERXI,PERXK,pERYV,PERNB,
I PERT,PERG,PERC,PERQM,PERQB,PERQC,PERTP,PERTO,PERTA
IF (I.EQ.9?) GO TO Z500
IF (IZS.EQ.OI GO TO 2600
200 LCITY = I
lEX = JyEAR+3
IKX = JVEAR + 2
IF(OX(II.EQ.O.OI GO TO lZ00
IF(GRFAC(LCITV,20 ,I).EQ.O.O) GO TO 1200
GFACI = GRFAC(I,20,1)
GFACZ = GRFAC(I ,20,2)
GFACD = GFACZ - GFACI
IF(GFACI.GT.O.O) AFI = ALOGIO(GFACI )/3.
RAFI = 10.cH>AFI
IF (GFACl.GT.O.O.AND.GFACZ.GT.O.O) AFZ = (ALOGI0(GFACZ)-
I ALOG10(GFAC11 )/5.
RAFZ = 10.-QAFZ
IF (GRFAC(LCITV,20,31.GT.0.0.AND.GFAC2.GT.0.OI
1 AF3=(ALOGIO(GRFAC(LCITV,20,3) )-ALOGI0(GFACZ»/5.
RAF3 = 10...AF3
GO TO (1120,1130,1130,1130,1130,1130,11 qO, 11 qO, 11 qO, 11 qO, 11 qO) ,
1 JVEAR
1120 VFAC = RAFl~.3.
GO TO 1150
1130 KJNDX = JYEAR-l
LJNDX = 0
RMS Z03Q
RMS Z035
RHS 2036
RMS Z037
Rf1S 2038
RMS Z039
RMS ZO'iO
R r1 S Z 0 'i 1
RMS ZOqz
RMS ZO'i3
RMS ZO'iq
RMS ZO'i5
RMS ZO'i6
RMS ZO'i7
RMS ZO'iA
RMS 20'i9
RMS Z050
RMS 2051
RMS205Z
RMS 2053
RMS ZOSq
RMS 2055
RM5 2056
RMS 2057
R r1 5 7 0 5 8
Rr1S 2059
RMS 2060
RMS 7061
RMS 2062
RMS 2063
RMS 206ti
RMS 2065
RMS 2066
RMS 2067
RMS 2068
RMS 2069
)A T ~
72.22
-------
S~ 0038
ISN 0039
ISN OO~O
ISN OOql
ISN 00Q2
ISN 00Q3
ISN OO~Lf
ISN no~s
ISN OOLi6
ISN 00~7
ISN 00~8
ISN 00Q9
ISN 0050
ISN 0051
ISN 0052
ISN 0053
JSN 005~
ISN 0055
JSN 0056
ISN 0058
ISN 0060
-.D ISN 0061
0'
ISN 0063
ISN 0065
JSN 0067
ISN 0068
ISN 0069
ISN 0070
JSN 0071
ISN 0072
ISN 0073
ISN 007~
ISN 0075
ISN 0076
ISN 0077
ISN 0078
ISN 0079
JSN 0080
VFAC = (RAF..3.).(RAF2..KJ~DX)
GO TO 1150
11~0 KJNDX = S
LJNDX = JYEAR-b
VFAC = (RAFl.~3.).(RAF2..KJNDX).(RAF3..LJNDX)
1150 PERXX = OX(I).VFAC
PER'N = OW(I)
DO 1161 JS = I,IEX
1161 PERW = PFRW~I.018
PERWB = OWB(I)GVFAC
2000 FORMATIIH ,I~,5Fl0.2)
PERP = PERXX-PERWB
PERXI = 0I(I)oVFAC
PERXK = PERXX/VK(20)
10 FORMAT I1H ,I~,JFIO.3)
PERSN = PERWB/PERW
1770 FORMAT(IHl,QI5)
1200 CONTINUE
IFILCITY.GT.5J) GO TO 120~
IF(XINCILCITY,l).NE.O.O) GO TO 1205
RATE = 0.0
IFII.FQ.22) RATE=.0762
IF II.EQ.31) RATE=.0758
IFII.EQ.Jl) RATE=.06~5
PERYY = OY(I)
PERYY = PERYY-Il.0+RATE)..IEX
GO TO 1207
120'+ PERYY=Oyll)
PERYY = PERYY-(YGROWII-5J)+1.0).~IEX
GO TO 1207
lZ05 CONTINUE
YRAGI = XINCILCITY,I)/OYII)
YRAG2 = XINCILCITY,2)/XINC(LCITY,I)
YRAG3 = XINC(LCITY,3)/XINC(LCITY,Z)
YRAMI = ALOGI0lYRAGI )/3.
YRALI = 10.$~YRAMI
YRAM2 = ALOGIO(YRAG2)/5.
YRAL2 = 10.-eYRAM2
R r~ 5 Z 0 7 0
RMS Z071
RMS 7072
RMS 2073
RMS 207q
Rr~~ 2075
RMC; 2076
RMS 2077
Rr~s 2078
RMS 2079
RMS ZOAO
RI-1S 2081
RMS 2082
RMS 2083
RMS 7.0Aq
RMS 2085
RMS 2086
RMS ZOA7
RMS 2088
R r~ 5 7. 0 8 9
Rr~5 2090
RMS 2091
RMS 2092
RMS 2093
RMS 7.09~
RMS 2095
RMS 2096
Rr~s 2097
RMS 2098
RMS 2099
Rt1S 7.100
RMS 2101
RMS 2102
RMS 210J
RMS 210q
RMS 2105
RMS 2106
RMS 210'/
-------
IS' (08
ISN 0082
ISN 0083
ISN 008Lf
ISN 0086
ISN 0087
ISN 0089
ISN 0090
ISN 0091
ISN 0092
ISN 0093
ISN 009Lf
ISN 0095
ISN 0096
JSN 0097
ISN 009R
ISN 0099
ISN 0100
ISN 0101
ISN 0102
..!) ISN 0103
-J ISN 010'+
ISN 0105
ISN 0106
ISN 0107
ISN 0108
ISN 0109
ISN 0110
ISN 0111
ISN 0112
ISN 0113
ISN 011Lf
ISN 0115
ISN 0116
ISN 0117
ISN 0118
ISN 0119
ISN 0120
lZ06
1207
121 0
17.1 2
YRAM3 = ALOG 0(YRAG3)/5.
YRAL3 = 10.o~YRAM3
PERYY = OY(I)t>YRAL1..3
IF (KJNDX.EQ.o.) GO TO 1206
PERYY = PERYYeYRALZ..KJNDX
IF (LJNOX.EQ.O.) GO TO 1206
PERYY = PERYY ~YRAL3~eLJNDX
PERYY = PERYY-l.1LfLf
YFAC = pERYY/OY(I)
PERNB = OXBAR(I)~(PERYY-PERXX)/(OY(I)-OX(I»
PERT = OT(I)<:>YFAC
PERTP = 0TP(J)oYFAC
PERTO = OTO(I) . YFAC
PERTA = OTA(I) . YFAC
PERG = OG(I).YFAC
PFRC = OC(I).YFAC
PERQM = OQM(I)oVFAC
PERQB = OQB(I).(PERYY-PERXX)/(OY(I)-OX(I»
PERQC = OQC(I)oYfAC
CONTINUE
FORMAT(Fl0.3,60X,15,I3)
ZERXX = ZERXX + PERXX
ZERW = ZERW + PERW
ZfRNB = ZERWB + PERWB
ZfRSN = ZERSN + PERSN
ZERP = ZERP + PERP
ZERXI = ZERXI + PERXI
ZERXK = ZERXK + PERXK
ZERYY = ZERYY + PERYY
ZERNB = ZERNB + PERNS
ZERT = ZERT + PERT
ZERTP = ZERTP + PERTP
ZERTO = ZERTO + PERTO
ZERTA = ZERTA + PERTA
ZERG = ZERG + PERG
ZERC = ZERC + PERC
ZERQM = ZERQM + PERQM
ZERQB = ZERQB + PERQB
R~S 2 08
RMS 7.109
RMS 2110
RMS 2111
RMS 2112
RMS 7.l13
RMS 211Lf
RMS 2115
RMS 2116
RMS 2117
RMS 2118
RMS 2119
RMS 2120
RHS 2121
RMS 2122
Rt1S 2123
RHS 7.12Lf
RMS 2125
R~1S 2126
Rt-1S 2127
Rt-1S 2128
RMS 2129
RMS 2130
RMS 2131
RMS 7.132
RMS 2133
RMS 213Lf
Rt-1S 2135
Rt-152136
RMS 2137
RMS 2138
RMS 2139
RMS 21LfO
RMS 21Lfl
RMS 21Lf2
RMS 7.1Lf3
RMS 21Lfq
RMS 21Lf5
-------
ISN 012
ISN 0122
TSN 0123
ISN 01lLf
ISN 0125
ISN 0126
ISN 0127
ISN 0128
ISN 0129
ISN 0130
ISN 0131
I
ISN 0132
ISN 0133
ISN 013Lf
ISN 0135
ISN 0136
ISN 0137
ISN 0138
ISN 0139
ISN 01LfO
-..0 ISN 01Lfl
00 ISN 01Lf2
ISN 01Lf3
ISN 01LfLf
ISN 01Lf5
ISN 01Lf6
ISN 01Lf7
ISN 01Lf8
ISN 01Lf9
ISN 0150
ISN 0151
ISN 0152
ISN 0153
ISN 015Lf
ISN 0155
ISN 0156
ISN 0157
ISN 0158
2500
2600
ZE~QC = ZERQC + ~ERQC
RETURN
PERXX = ZERXX
PfRW = ZERW
PERI'IB = ZERWB
PfRSN = ZERSN
PERP = ZERP
PERXI = ZERXI
PERXK = ZERXK
PERYY = ZERYY
PERNB = ZERNB
PERT = ZERT
PERTP = ZERTP
PERTO = ZERTO
PERTA = ZERTA
PERG = ZERG
PERC = ZERC
PER Q M = Z E R Q 11
PERqB = ZERQB
PERQC = ZERQC
ZERXX = 0.0
ZERW =0.0
ZER\"vB=O.O
ZERSN=O.O
ZERP =0.0
ZERXI=O.O
ZERXK=O.O
ZERYY=O.O
ZERNB=O.O
ZF:RT =0.0
ZERTA = 0,,0
ZERTP = 0.0
ZERTO = 0.0
ZERG =0.0
ZERC =0,,0
ZER(~M=OoO
ZERQB=O.O
ZERQC=O.O
H1S 2 . 6
RMS 21'17
RMS 21 LfB
RMS 21Lf9
RMS 2150
R'''S 2151
RMS 2152
RMS 21S3
RMS 215Lf
RMS 2155
RMS 2156
RMS 2157
RMS 2158
RMS 2159
RMS 2160
RMS 2161
RMS 2162
RMS 21A3
RHS 216Lf
RMS 2165
RMS 2166
RMS 2167
RMS 2168
RMS 2169
RMS 217r)
RMS 2171
RMS 2172
RMS 2173
RMS 217Lf
RMS 2175
RMS 2176
RMS 2177
RMS 2178
RMS 2179
RMS 2180
RMS 2181
RMS 2182
RMS 2183
-------
ISN 0159
ISN 0161
ISN 0162
ISN 0163
ISN 0161f
ISN 0165
-.D
-.D
2700
F (IZS.NE.O
IZS = 1
GO TO 200
CONTINUE
RETURN
END
RETURN
~MS 2181f
RMS 2185
RMS 2186
RMS 2187
RMS 2188
RMS 2189
-------
LEV~- 20.
ISN 0002
ISN 0003
ISN OOoq
ISN 0005
ISN 0006
ISN 0007
ISN 0008
ISN 0009
ISN 0010
ISN 0011
ISN 0012
ISN 0013
I5N 00lQ
ISN 0015
ISN 0016
I
I~
o
o
caMP
MAY 7 )
05/360
FORTRAN H
_~R OP' .ONS - NA' ~= ~A N,O:>.'=Oo,LINf.CNT=QO,SIZE=OOOOK,
50URCE,BCO,NOLIST,NODECK,LOAD,MAP,NOEDIT,NOID,NOXREF
SUBROUTINE VALUE
COMMON ISIMI M1 ,M2,MCA5ES,Y( 1000) ,XI 1000)
1 ,MARK,LVAR
GO TO(100,200,300,QOO,500) , MARK
100 CALL VALUEI
GO TO 1000
200 CALL VALUE2
~O TO 1000
300 CALL VALUE3
GO TO 1000
'100 CALL VALUELf
GO TO 1000
500 CALL VALUES
1000 CONTINUE
RETURN
END
RMS 2199
RMS 2200
RMS 2201
RMS 2202
RMS 2203
RMS 7.20'1
RHS 2205
RMS 2206
RMS 2207
RMS 2208
RMS 2209
RMS 2210
RMS 2211
RMS 2212
RMS 2213
RMS 221'1
)A" ~
72.22
":>
-------
LEVE- 20.1
(HAY 7 )
05/360
FORTHII H
COMP . - ~ ~ OP > IONS - NAM::= "A N, OPT=OO, - NEC1\I'>=1 0, S . ZE=OOOOK,
SOURCE,BCO,NOLI5T,NODECK,LOAD,MAP,NOEOIT,NOIO,NOXREF
SUBROUTINE VALUEI
ISN 0002
15N 0003
ISN OOO~
ISN 0005
ISN 0006
ISN 0007
ISN 0008
ISN 0009
ISN 0010
15N 001]
ISN 0012
ISN 0013
......
o
......
c
c
COMMON 151MI
1 ,MARK,LVAR
COMMON IFATI
COMMON I BUG
Ml ,MZ,MCASES,Y( 1000) .X( 1000)
F( 1000) ,LSTORE,KPIST
I ISBUG,LREAO,LWRIT
J = LSTORE + 1
DO 1 I=Ml,MZ,LVAR
F(I) = -VI]) + .~2836.XIJ) +
FII+I) = -Y(I+l) + 23.~227+
1 + X(J+Z)
FII+2) =-Y(I+2) + 137.9922 + .bI3256.Y(I)
1 + X(J+3)
1 CONTINUE
RETURN
END
2.83038eY(I+l) + .9~812.Y(I+2)
.9 q 215 .'X ( J + 1 )
RMS 2215
RMS 2216
RMS 2217
RMS 7.218
RMS 2219
RMS 2220
RMS 2221
RMS 7.222
RMS 2223
RMS 7.22q
RMS 227.5
RMS 2226
RMS 2227
RMS 222A
RMS 7.229
RMS 2230
RMS 2231
OA" ::
72.27
'>
-------
1-
LEVEL 20.1
IMAY 71)
05/360
FORTRAN H
COMPILER OPTIONS - NA'::= MAl' ,OP "=00,LINECNT=,+0,5IZE=0000K,
SOURCE,BCD,NOLIST,NODECK,LOAD,MAP,NOEDIT,NOID,NOXREF
SUBROUTINE VALUE2
ISN 0002
ISN 0003
ISN 000'+
ISN 0005
ISN 0006
ISN 0007
ISN 0008
ISN 0009
ISN 0010
I5N 0011
ISN 0012
ISN 0013
ISN 001'+
ISN 0015
I-'
o
N
C
c
COMMON ISIMI
1 ,MARK,LVAR
COMr1QN IFATI
COMMON I BUG
Ml ,M2,MCASES,YI 1000) ,x( 1000)
FI 1000) ,L5TORE,KPIST
I ISBUG,LREAD,LWRIT
J = LSTORE
DO 1 I=MI,M2,LVAR
FII) =-yII) + XIJ+,+)8YII+,+) + XIJ+I)
FII+l) = -YII+I)-13.95875 + 1.0392'+5-YII)+ 361.37'+-YII+3)
1 + X (J+ 7 )
FII+2) =-YII+2) + 1.8Q02 + ( .8377-XIJ+3» +(.1116oYII+3»
1 + X(J+6)
FII+3) =-YII+3) + «(YII+!)-Y(I)/Y(J+l)
F(I+,+) =-YII+q) + IX(J+2)/Y(I+2»
1 CONTINUF.:
RETURN
END
RMS 2232
RMS 2233
RMS 22314
RMS 2235
RMS 2236
RMS 2237
RMS 2238
RMS 2239
RMS 22'+0
RMS 22'+1
RMS 22'+2
Rr1S 22'+3
RMS 7.2'+,+
RMS 22'+5
RMS 22'+6
RMS 22'+7
RMS 22'+8
RMS 22'+9
RMS 2250
DATE
72.22
. )
-------
LEVE- 20.
(MAY 71)
05/360
fORTRAN H
(OMP I - ~R 0:>" ONS - NA'1E= MA N, 0 :>'=00, - N :CNT=qO, 5 I ZE=OOOoK,
SOURCE,BCD,NOLIST,NODECK,LOAD,MAP,NOEDIT,NOID,NOXREf
SUBROUTINE VALUE3
15N 0002
ISN 0003
15N ODOLf
ISN 0005
ISN 0006
ISN 0007
ISN 0008
ISN 0009
ISN 0010
ISN 0011
ISN 0012
ISN 0013
..-
o
w
(
COHMON 151MI
1 ,MARI(,LVAR
COHMON IFATI
COMMON I BUG
Ml,M2,HCASES,YI 1000) ,XI 1000)
FI 1000) ,LSTORE
I 15BUG,LREAD,LWRIT
c
c
J = LSTORE
00 1 I=Ml,M2,LVAR
F(I) = -VII) + .qZ836.X(J+3)+ 2.83038.YII+1) +
FII+I) = -YII+1) + 23.Lf227+ .9q215.XII+3).X(J+1)
1 + XIJ+q)
F(I+2) =-Y!I+2) +
1 + X(J+5)
1 CONTINUE
RETURN
END
.9Q812.Y(I+2)
137.9922 + .613256.Y(I) .XIJ+2)
RMS 2251
RMS 2252
RMS 2253
RMS 7.2Sq
RMS 2255
RMS 2256
RMS 2257
RMS 2258
RMS 2259
RM5 2260
RMS 2261
RMS 2262
RMS 2263
RMS 226Lf
RMS 221>5
RMS 2266
RMS 2267
RMS 2268
DATE
72.221
':>
-------
LEV~- 20.1
ISN 0002
IsN 0003
IsN OOO~
IsN [JOOS
IsN 0006
IsN 0007
rSN 0008
IsN 0009
ISN DOlO
ISN nOll
ISN 0012
IsN 0013
IsN OOI~
ISN 0015
ISN 0016
IsN 0017
IsN 0018
ISN 0019
:; ISN 0020
,p. ISN 0022
IsN 0023
ISN 002Lf
ISN 0025
IsN 0026
ISN 0027
ISN 0028
ISN on29
ISN 0030
ISN 0031
ISN 0032
ISN 0033
ISN 003~
ISN 0035
ISN 0036
I MA Y 7 )
05/360
FORTRAN H
COMPI_~R OPT ONs - NAM~= MAIN,OPT=Oo,.IIJECN"=LfO,sIZE=O:OOK,
sOURCE,BCD,NOLIsT,NODECK,LOAD,MAP,NOEDIT,NOID,NOXREF
SUBROUTINE VALUELf
DIMENSION A(Lf)
C
COMMON 151MI MI ,M2,MCASES,YI 1000) ,X( 1000) ,MARK,LVAR,MKEEP
COt1MON IFATI F(1000),LsTORE
MI., = LSTORE + 1
MT = 5 - MKEEP
MT1 = 2~MT
KD = MKEEp - 1
KD1 = KD02 + 1
DO 20 I=Ml,f'12,LVAR
Ii = I-MT
MR = MM - MT1
SUM = 0.0
B = ALOGIYIII+8)) - ALOGIX(MR+ILf))
00 lO KK=l,KD
K = KK - 1
SUM = SUM + XIM~1+K).YII+K)
AIKK) = 1.0
IFIXIMM+9)~XIMM+K).EQ.0.0) GO TO 9
AIKK) = IXIMM+KD)/XIMM+K)).(XIMM+K+MKEEP)/XIMM+KDI))
9 CONTINUE
B = B - XIMM+K+MKEEP) - ALOGIAIKK))
10 FII+K) = -YII+K) + AIKK).YIII+Lf)
F(II+6) = -YIII+6) + SUM + X(MM+KD).yll+KD)
MM = MM - MT1
FIII+Lf) =-YIII+Lf) + EXPIB)
FIII+5) = -YIII+5) + IXIMM+10)+YIII+10))/11.0+XIMM+11))
FIII+7) = -YIII+7) + YIII+S).XIMM+16)
FIII+8) =-Y(II+8)+XIMM+15).YIII+7)
FIII+9) = -YIII+9) + XIMM+12).YIII+7)
FIII+I0) = -YIII+10) + XIMM+17)+ XIMM+13).Y(II+9)
20 CONTINUE
RETURN
END
RMS 2269
RMS 2270
RMS 2271
RMS 2272
RMS 2273
RMS 227Lf
RMS 2275
RMS 7276
RMS 2277
RMS 227B
RMS 2279
RMS 2280
RMS 2281
RMS 2282
RMS 2283
RMS 228Lf
RMS 22BS
RMS 2286
RMS 2287
RM5 2288
RM5 2289
RMS 2290
RMS 2291
Rt1S 2292
RMS 2293
RMS 229~
RMS 2295
RMS 2296
RMS 2297
RMS 2298
RMS 2299
RMS 2300
Rt1S 23n 1
RMS 2302
RMS 2303
DAT~
72.22
-------
COMPILER OP"O~S - NA'1E= 'A N,OP"=QO,LINECNT=LfO,SIZE=(OOOK,
SOURCE,BCD,NOLIST,NODECK,LOAO,MAP,NOEOIT,NOIO,NOXREF
SUBROUTINE VALUES
COMMON ISIMI MI ,M2,MCASES,Y< 1000) ,X( 10001
1 ,MARK,LVAR
COMMON IFATI F(1000),LSTORE,KPIST
COMMON I BUG I ISBUG,LREAO,LWRIT
J = LSTORE
00 I I = MI,M2,LVAR
FII) = -YIII+.056084XIJ+I)
F(I+l) = -Y(I+l)-21.Lf3+.01223.XIJ+Z)
F(I+2) = -Y( 1+2)+317.S'f-815.72'fa(Y(I)+Y(I+1 II/Y(I+3)-
1 .0272$X(J+21/XIJ+31+1.0SS~(YIII+Y( 1+1) I/X(J+3)
F<'I+3) = -Y(I+31+YII+I )+Y(I+2)+Y(I)
1 CONTINUF.
RETURN
END
LEV~- 20.1
ISN 0002
ISN 0003
ISN 000..
ISN 0005
ISN 0006
ISN 0007
ISN OOOS
ISN 0009
ISN DOlO
ISN 0011
ISN 0012
ISN 0013
ISN' 001..
......
o
\.J1
(MAY 71)
05/360
FORTRAN ~
RMS 230'f
RMS 2305
RMS 2306
RMS 2307
RMS 230S
RMS 2309
RMS 2310
RMS 2311
RMS 2312
RMS 2313
RMS 231'f
RMS 2315
RMS 2316
RMS 2317
RMS 2318
)A" ~
72.22
" ,)
-------
LEVE - 2C.
ISN 0002
lSN 0003
ISN 0001f
ISN 0005
ISN 0006
ISN 0007
ISN 0008
ISN 0009
ISN (lOIO
ISN 0011
ISN 0012
ISN 0013
I-'
o
0'
(MAY 71)
05/360
FOR'RA' H
COMPILER OPTONS - NAME= MAIII,OPT=OO,LNECNT=IfO,5 ZE=OOOOK,
SOURCE,BCD,NOLIST,NODECK,LOAD,MAP,NOEDIT,NOID,NOXREF
SUBROUTINE DEBUG(y,MCNT,LVAR)
DIMENSION Y(lOOO)
VJRITE(6,21 )
MCASES = LVAR~MCNT
K = 0
DO 10 I=I,MCASES,LVAR
K = K + LVAR
10 WRITE(6,20) (Y(J), J=I,K)
20 FORMAT(!H ,10Fl1.1)
21 FORfHTIIHl)
RETURN
END
RMS 2319
RMS 2320
RMS 2321
.RMS 2322
RMS 2323
RMS 2321f
RMS 2325
RMS 2326
RI1S 2327
RMS 2328
RMS 2329
RMS 2330
)A"[
72.2.
-------
LEV~- 20.
ISN 0002
ISN 0003
ISN 0001./
lSN ODDS
ISN 0006
IsN 0007
IsN 0008
ISN 0009
IsN 0010
IsN 0011
ISN 0012
IsN 0013
ISN 0011./
ISN 0015
...... ISN 0016
1:3 ISN 0017
ISN 0018
IsN 0019
IsN 0020
lsN 0021
IsN 0022
IsN 0023
ISN 0021./
IsN 0025
ISN 0026
ISN 0027
IsN 0028
ISN 0029
ISN 0031
(MAY 71)
05/360
FORTRAN ~
COMPILER O:>"O~s - NAME= 'A N,O:>'=OO,- NECNT=l.fO,SIZE=OOOOK,
SOURCE,BCD,NOLIST,NODECK,LOAD,MAP.NOEDIT,NOID,NOXREF
SUBROUTINE 51MUL
c
C
C
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
RMS
IX = tCTAB + LCHECK + NBEFOR RMS
IF(IX-3) 69,60,69 RMS
60 IF(NEWCIT) 71,70,71 RMS
71 I = IJCIT RMS
WRITE(LPR,I20) (IAQCR(I,J), J=I,20),LPYEAR RMS
120 FORMAT(lHl,111125X,'AQCR ',13,IX,t9Al.f,SX,'FOR ',Ii.f112SX, RMS
1 'CONVERGENCE PROCESS FOR ESTIMATED VALUES WITHOUT POLLUTION CONTRRMS
10LS') RMS
70 CONTINUE RMS
WRITE(LPR,50) MARK RMS
IF(MARK.EQ.2) GO TO 900 RMS
50 FORMATIIH ,III./X,'...............~.................................RMS
MACRO MODEL SIMULATION PROGRAM I - MAINLINE
DIMENSION DEL(1000)
DIMENSION IAQCRI92,20)
COMMON ISIMI Ml ,M2,MCASES,YI 1000) ,X( 1000)
1 ,MARK,LVAR
COMMON IFATI FlIOOO) ,LSTORE,KPJST,ARAB( 10) ,LCHECK
COMMON ICATI A(20,20)
COMMON ISENSI ICTAB,NBEFOR,NEWCIT
COMMON I IPRIN IINTT,LPNTT ,JYEAR,IAQCR,IJCIT
COMMON I BUG I IsBUG,KRD,LPR
LSTORE = 0
KREP = 0
106 FORMAT(II./)
MCNT = MCASES/LVAR
U1N = LVAR
IMAX = 1./
LPYEAR = JYEAR + 1969
M 1 = l'
H2 = MINOILMN,MCASES)
IT = 0
C
2331
2332
2333
2331./
2335
2336
2337
233A
2339
23'10
23'11
231./2
23l.f3
23'11./
231./5
231./6
23 '+7
23l.f8
23'19
2350
23S1
2352
2353
235'f
2355
2356
2357
2358
2359
2360
2361
2362
2363
236'f
2365
2366
DATE
72.22
-------
ISN 0032
ISN 0033
ISN 003Lf
ISN 0035
ISN 0036
ISN 0037
ISN 0038
ISN 0039
ISN DOLfO
ISN 00'11
ISN OO'i2
ISN 00'13
....... IS N
~ ISN
ISN
ISN
ISN
OOLfLf
00'15
00'16
00'17
00'18
ISN OOLf9
ISN 0050
ISN 0051
ISN 0052
ISN 00S3
ISN OOSLf
ISN 0055
ISN 00S6
ISN 0057
ISN 0058
ISN 0059
ISN 0060
.
I[[[RMS 2367
2'/T5.'.'.TI20.'.'/T5.'.'.25X,'SIMULTANEOUS SYSTEM NUMBER - '.11, RMS 2368
3 TI20.'.'/T5,'.'.TI20.'.') RMS 2369
RMS 2370
RMS 2371
RMS 2372
RMS 2373
RMS 237'1
RMS 2375
RMS 2376
RMS 2377
RMS 237a
RMS 2379
RMS 2380 .
RMS 2381
RMS 2382
RMS 2383
Rf-1S 238'4
RMS 2385
RMS 2386
RMS 2387
RMS 2388
RMS 2389
RMS 2390
RMS 2391
RMS 2392
RMS 2393
RMS 239'1
RMS 2395
RMS 2396
RMS 2397
RMS 2398
RMS 2399
RMS 2'100
RMS 2'101
RMS 2Lf02
RMS 2Lf03
RMS 2'iOLf
(
2 WRITEILPR.9) (ARAB(J),J=t,LVAR)
WRITE(LPR,IO) (YIJ), J=MI,M2)
9 FORMATIIH ,T5,'O',5X,'ENDOGENOUS VAR.
I TI20,'.')
10 FORMATIIH ,T5,'.',aX,'ACTUAL VALUE ',3IF9.2,8X),TI20".')
II FORMATIIH ,TS,'.',8X,'ITERATION-',I2,2X,3IF9,2,F8,S),TI20,'.')
GO TO 69
',3(SX,ALf,'
CONVI),
(
900 CONTINUE
WRITEILPR,99)(ARAB(J),J=I,LVAR)
WRITEILPR,91D) IY(J), J=Ml,M2)
99 FORMATIIH ,TS,'~',5X.'ENOOGENOUS VAR.
1 TI20,'.')
910 fORHATIIH ,Ts,'.',aX.'ACTUAL VALUE ',SIF9.2,8X),T120,'.')
911 FORMATIIH .T5.'.',8X,'ITERATION-',I2,2X,S(F9.2,F8.S),TI20,'.')
',S(SX.ALf,'
CONVI),
C
5,690,690
C
69 IF IIT-IHAX)
S IT = IT & 1
3 CALL GRADN
CALL VERT
CALL VALUE
COMPUTE NEW VALUES AND TEST FOR CONVERGENCE
7 IAS = IT
IT = IMAX
NN = 0
00'1 I = Hl,H2
NN = NN + 1
YY = YII)
Y I I) = o. 0
KK = 0
-------
..... ISN
o ISN
...0
ISN
ISN
ISN
ISN
ISN
0082
0083
0085
0086
0087
0088
0089
YD 5 = ABSIYI )/YY RMS 2tf05
YOIS = ABSIYOIS-I.) RMS 2tf06
DELlI) = YDIS RMS 2~07
IFIYDIS.LT.O.009)GO TO tf2 RMS 2tf08
IT = IAS RMS 2tf09
~2 CONTINUE RMS 2~10
~ CONTINUE RMS 2~ll
IFIIX.NE.3) GO TO 69 RMS 7.~12
688 CONTINUE RMS 2~13
IFIMARK.EQ.2) GO TO 673 RMS 2~1~
'IIRITEILPR,11) IAS,IYIJ),DELIJ), J=MI,M2) RMS 7.Lf15
GO TO 69 RMS 2tfl6
673 ~IjRITE(LPR,911) IAS,(Y(J) ,DEL(J), J=Ml,M2) RMS 2tf17
689 GO TO 69 RMS 2~IA
690 IF(IX.NE.3) GO TO 692 RMS 2~19
691 ~'iRITE ILPR,10S) IAS RMS 2'120
1 0 5 FOR 11 A T I I H ,T 5 , ' . , ,T 1 2 0 , ' . , , I T 5 t ' . , , ~ X , ' I T ERA T ION 5 T 0 CON V ERG E = , , R M 5 2 ~ 7. 1
1 12tTI20t'o't/~X,'o~""'."'o~"'o""""""""""".""",oRMS 2tf22
2""..o~'.......o.......a......'"",""""o""",""""""'RMS 2~23
Lf"') RMS 2Lf2'1
692 CONTINUE RMS 2'125
IF 1M2 .EQ. MCASES) RETURN RMS 2'126
MI = HI + LMN RMS 2Lf27
LSTORE = LSTORE + KPIST RMS 2'128
H2 = MINOIM2+LMN,MCASES) RMS 2~29
GO TO I RMS 2~30
END RMS 2Lf31
ISN 0061
ISN 0062
ISN 0063
ISN 006~
ISN 0066
ISN 0067
ISN 0068
ISN 0069
ISN 0071
ISN 0072
[SN 007'1
[SN 0075
ISN 0076
ISN 0077
[SN 0078
[SN 0080
ISN 0081
-------
LEVEL 20.1 (MAY 7 )
OS/3b(
:-OR"RAN H
COMP L~~ O;)TIONS - NAM~= MA IJ,O;)'=OO,LINECN"=LfO,S ZE=OOOOK,
SOURCE,BCO,NOLIST,NODECK,LOAD,MAP,NOEDIT,NOIO,NOXREF
SUBROUTINE GRAON
ISN 0002
ISN 0003
ISN ODOLf
ISN 0005
ISN 0006
ISN 0007
ISN 0008
ISN 0009
ISN 0010
ISN 0011
ISN 0012
ISN 0013
ISN 001Lf
ISN 0015
ISN 0016
:::: ISN 0017
o ISN 0018
ISN 0019
ISN 0020
ISN 0021
ISN 0022
ISN 0023
ISN 002/f
ISN 0025
ISN 0026
ISN 0027
c
C
C
MACRO SIMULATION I - GRADN
OIMENSION 0(1000)
COMMON ISIMI Ml ,M2 ,MCASES,Y( 1000) ,X (1000)
I ,MARK,LVAR
COMMON IFATI F(IOOO),LSTORE
COMMON ICATI A(20,20)
COMMON I BUG I ISBUG,KRD,LPR
106 FORMAT (IH ,Il.f)
NN' = 0
00 1 1 = M 1 ,M2
NN = NN + 1
XD = .0000IeY(I)
IF (XD-.OoOOI) If, If, 5
/f XD = .001
5 YII) =y(I) & XD
CALL VALUE
DO 2 J = M 1 ,M2
2 DIJ) = F(J)
Y(I) = Y(I) - 2.0.XO
CALL VALUE
KK = 0
D 0 3 J = 1'1 1 ,M 2
KK = KK + 1
3 A(KK,NN) = (DIJ)-F(J»/(XO.2.)
1 Y(I) = Y(J) & XD
RETURN
END
RMS 2Lf32
RMS 2/f33
RMS 2/f3/f
RMS 2/f35
RMS 2/f36
RMS 2/f37
RMS 2/f38
RMS 2/f39
RMS 2/f/f0
RMS 2/f'H
RMS 2/f/f2
RMS 7./fl.f3
RMS 2/f1f/f
RMS 2/fLf5
RMS 2/fLf6
RMS 2'1/f7
RMS 2l.f/f8
Rt'1S 2'1lf9
RMS 2'150
RMS 7.'151
RMS 2l.f52
RMS 2Lf53
RMS 2'1SLf
RMS 2/f5S
RMS 2/fS6
RMS 2/fS7
RMS 2Lf58
RMS 2/fS9
RMS 2/f60
RMS 2/f61
OA"E
72.22
-------
LEV~- 20.
(MAY 71)
05/360
FORTRAN H
COMPILER OPTIONS - NAM~= MAIN,O;)'=OO,- NECNT='fO,5IZE=0000K,
SOURCE,BCD,NOLIST,NODECK,LOAD,MAP,NOEDIT,NOID,NOXREF
SUBROUTINE VERT
ISN 0002
ISN 0003
ISN OOO'f
ISN 0005
ISN 0006
ISN 0007
ISN 0008
ISN 0009
ISN 0010
ISN 0011
ISN 0012
ISN 0013
ISN 001Li
ISN 0015
...... ISN 0016
...... ISN 0017
......
ISN 0018
ISN 0019
C
C
C
MACRO SIMULATION I - VERT
COMMON ISIMI HI ,M2,MCASES,Y( 1000) ,X( 1000)
1 ,MARK,LVAR
COMMON ICATI A(20,20)
COMMON I BUG I ISBUG,KRD,LPR
7 00 1 K=I,LVAR
R = 1.0/A(K,K)
A(K,KI = leD
DO 2 J=I,LVAR
2 A(K,J) = RoA(K,J)
00 1 1=I,LVAR
IF (K-I) 3dt3
3 AIK = A(I,K)
A(I,KI = 0.0
00'+ J=I,LVAR
'f A(I,J) = A(I,J) - AIK.A(K,J)
1 CONTINUE
RETURN
END
RMS 2'f62
RMS 2Li63
RMS 2'f6'f
RMS 2'f6S
RMS 2'f66
RMS 7.'f67
RMS 2'+68
RMS 2'f-69
R '1 5 2 'f 7 0
RMS 2'f71
RM5 7.'f72
RMS 2'f73
RMS 2'f7'f
RMS 2'f75
RHS 2'f76
RM5 2'f77
RMS 2'f78
RMS 2'f79
RMS 2'f8Q
RMS 2'fRI
RMS 2'f82
RMS 2'fA3
OA"~
72.22
'OJ
-------
.
--------- ---
3.0
PROGRAM lOA (INPUT-OUTPUT ANALYSIS)
3. 1
Introduc tion
Program lOA contains a lOO sector United States input-output
flow table, and is designed to conduct the following functions:
First of all, this input-output table can be aggre-
gated into a new table according to the numbe r of
sectors desired by the user. The maximum number of
the sectors after aggregation must not exceed 49.
Then the program will estimate a new input-output
coefficient table and the corres ponding Leontief
inverse matrix.
Take changes in the final demand by sector as the
input, and calculate the changes in the gross
domestic products by sector (in terms of value of
shipment).
1£ change in the total consumption of the United
States is given, the program will estimate
changes in consumtpion by sector acccording to
average propensity to consume by sector stored
in the program and calculate changes in the gross
domestic products by sector (in terms of value of
shipments).
A detailed description of input preparations and the program is pro-
vided in the following sections.
112
-------
I-~----~---"""---"-
3.2 Job Deck Setup
The ProgranJ lOA (Input-Output A nalys is) is a 535 statement
FORTRAN program which uses approximately 184K of core and
30 seconds on the IBM 360/75.
Includes a set of data ~
and six sets of
parameter cards
/*
Data Input Deck
/ / GO S YSIN DD>.'<
Includes the f Program lOA
mainline and / / (Input-Output Analysis)
t h r e e s ub - ~
routines ((
(535 card(3)
Program Deck
/ / Job
Job Control Language
(JCL) Cards
113
-------
3.3
Data Input Deck
This deck includes seven sets of parameter and input data cards,
a total of 1,391 cards':' arranged in the following manner.
1.
Original Input-Output Flow Table: 1,296 cards
Flow table is arranged by reading the first 10 columns of
the table, then the next 10 columns, etc.
Cols. 1-70
Cols. 73-75
Cols. 78 -80
2.
Dimension Card:
1 OF7. 0, each card reads 10 fields, that is,
first 10 columns of first row of the input-
output table and then first 10 columns of
second row of the table, and so on.
13, identification of the tenth of the column
of the 1-0 flow table; for example, 10th,
20th, 30th, etc.
13, number of the row of the 1-0 flow table,
i.e., 1,2,3, etc.
1 card
This card specifies the matrix dimensions of the original
and aggregated 1-0 tables.
Cols. 1-4
Cols. 5-8
115, total number of columns in the original
table
108, total number of rows in the original
ta ble
The original 1-0 table included a detailed final demand sector,
aggregation of rows or columns and some misdevious sectors.
':'Number of cards will be
cified by the user. At present,
sector input-output table.
changed according to the sectoring s pe-
the data input deck is prepared for a 42
114
-------
Cols. 9-12
14, total number of columns in the new 1-0
table. At present, it includes 42 sectors
but column 43 is consumption, therefore,
it is 43.
Cols.13-16
14, total number of rows in the new 1-0
table. Again, row 43 is value-added by
sector, therefore, it is 43.
3.
Sector Specification Card(s}: 3 cards (for 42 sector table)
Since the new 1-0 table is an aggregation of the original 1-0
table, each sector in the new table contains one or more than
one sector(s) from the original table. These cards specify
the number of sectors which will be aggregated into each
"new" sector in the form of format (2014). For example, if
the first sector of the new table aggregates three sectors and
the second sector contains two sectors from the old table,
punch 3 on column 4 and 2 on column 8 and so on.
4.
Sector Ordering Cards: 6 cards (for 42 se ctor table)
Sectors in the original table must be lined up according to the
order of the sector specification cards in the form of format
(2014). For example, if sectors 1, 4, 7 of the original table
will be aggregated into the first sector of the new table, punch
1 on column 4, 4 on column 8, 7 on column 12 and so on.
5.
Label Cards: 42 cards, one card for each new sector (42
sectors)
Cols. 1-28
Label of the sector
6.
Change of United States Consumption Expenditures:
1 card
Cols. 1-6
F6. 1 Benefit from. cleaner air in terms of
change in total consumption expenditures by
($ million). It will be distributed by average
propensity to consurne by sector.
115
-------
7.
Change in Final Demand by Sector: 42 cards, one card for
each sector in the new I -0 table
Cols. 1-10
FIO.2 Change in final demand by sector,
if it is a decrease of the demand, add a
negative sign ($ million).
3.4
Description of the Program Deck
The program consists of a mainline and three subroutines.
The
. cards are sequenced along the right hand edge in columns 73 -80 with
the prefix lOA.
The program begins by reading the input-output table (l00 sectors)
from the cards.
The program then immediately transfers control to
subroutine AGGRE.
The subroutine immediately reads in the first
parameter card specifying the number of rows and columns in the
original 1-0 table and the number of rows and columns which are to be
in the aggregate matrix.
Following this, cards are read in which tell
AGGRE how many original sectors are to be aggregated in each of the
new sectors and which of the original sectors are to be aggregated into
which new sectors.
AGGRE then proceeds to do the actual aggregation
of old sectors into new sectors.
When the aggregate sectors have been
computed, control. is returned to the main program.
The main program now reads and prints the labels for the sec-
tors of the new 1-0 table calculated by AGGRE.
116
-------
The program next tabulates and prints the new 1-0 flow table.
Next, the new 1-0 coefficients table is generated and printed.
The pro-
gram now calls subroutine MINVR which calculates the Leontief inverse
matr ix of the 1-0 table.
When control is returned from MINVR, the
Leontief inverse matrix and its determinant are printed.
The program
now reads the estimated averc.ge consumption by sector, calculates and
prints each sectorls share of the increased consumption.
Next, a
series of cards, which contain the change in final demand by sector,
are read.
Changes in dome stic products by sector due to the air pol-
lution control are now calculated and the resultant cost and benefit to
the structure of the national economy is printed, and the program term-
inates.
The flow chart of Program lOA is given in Figure 3.1 and the pro-
gram listing is provided in Section 3. 5.
117
-------
FIGURE 3. 1
Sta rt
- - ---
Read labels
for new sec-
tors matrix
Print the
new I-a
flow table
Generate the
new I-a
coefficient
table
Print the
new I-a
coefficient
table
118
Aggregate
original I-a
. sectors into
I new sectors
,i
-------
':'
FIGURE 3. 1 (continued)
Calculate the
Leontieff inverse
of the aggregate
1-0 table.
-.-- -
Print the
Leontieff
ead change in
consumption
and final
dema nd data
Calculate the change in
the gross domestic
product by sector
Print the effects of
air quality control
on the structure of
the national eco-
nom~
~
119
-------
3.5
Frog,ram lOA (Input-Output Analysis)
':"
120
-------
i--
~
~
['
L
IJ
C
L
I:
I.,
I~
c.
!->
c'
"
*******************************************************0****
*
_.
~.
PROGRAM IDA
PREPARED RY CONSAD RESEARCH CORPORATION
APRIL 1971
*
*
*
*
*
THIS PROGRAM I
-------
......
N
N
?q TO SACCO MAN F~CTUR:S
30 RROAD fiND NA~ROW FAA~ CS, YARN AN) "HR :1\0 M -_S
31 ~Isr. T[XTILE GOJDS AND FLOO~ COVERINGS
37 APPAREL
33 MISC. FAnRICATED TEXTILE PRODUCTS
3/t LLJr.H~ER "NO \>IOOD PRODUCTS, EXC. CONTAINERS
35 wonnE~ CO~TAINERS
36 HOUSEHOLD FUR~lITURE
,7 OTHER FURNITURE AND FIXTlJRES
3H rAPFR '\"JD ALLIFr) PRODUCTS, EXC. CONTAINERS
3q PAPERROARO CONTAINERS AND ROXES
40 PRINTING AND PUBLISHING
41 CHEMICALS AND SELECTED CHEMICAL PRODUCTS
42 PLASTICS AND SYNTHETIC t-1!1TERI ALS
41 DRIJGS, CLEANING, f\ND TOILET PREPARATIONS
44 PAPHS AND hLLIEJ) PRODUCTS
45 PFT~OLEU~ REFI~I~G A~D RFLATED INDUSTRIES
46 RU~RFR AND MlSCELLtNFOUS PLf\STICS PRODUCTS
47 LEATHER TANNING AND INDUSTRIAL LEATHER PRODUCTS
48 FOOTWEAR AND OTHER LEATHER PRODUCTS
4q GLL\SS :'\ND GLASS PRODIICTS
50S T (] N E AND C LAY P q 0 D UC T S
51 P RIM £l R Y I R 0 t-J AND S TEE L M i\ N U F ACT U R I N G
52 COPPER t4ANUFI\CTLJR ING
53 AllJ~INU~ '<1ANUFACTlJR ING
54 nTH~R NONFERROUS ~-1ET"LS MANUFACTURING
55 METAL CONTAINERS
56 HEATING, PLUMRI~Jr, AND STPUCTURAL t-1ETAL PRODUCTS
57 STAMPINGS, SCRfW MACHINE PR00UCTS AND BOLTS
5 q 0 T Ii E R F i\ f3 R I C " TED M F TAL PRO 0 UC T S
5q ENGINES nND TURBINES
60 FA~M MACHINERY AND fOUIPMENT
61 CONSTRUCTION, ~INING, AND OIL FIELD MACHINERY
6 2 ~H T F R I A L S I H\ N D L I N G 1"1/\ CHI N E R Y AND ~ QUI P MEN T
6 3 ", E T 1\ lIrHl R KIN G t1 A CHI N E R Y i\ N D E 0 U 1 P ,-, E N T
64 SPECIAL INDUSTRY MACHINERY I\ND EQUIPMENT
65 GENERAL INDUSTRIAL MACHINERY AND EQUIPMENT
66 MACHINF SHOP P~O~UCTS
67 n,FFICE, CO'-1PlJTIN~ AND ACCOUNTING r~ACHINFS
6R SERVICE INDUSTRY ~ACHINES
69 ELFCTRIC INDUSTRI AL EOUIP~~ENT AND APPARATUS
70 HClIJSEHOLD f\PPL II1.NCFS
71 [LFCTRIC LIGHTIN~ !INn WIRING EQUIPMENT
72 r~AnIO, TV AI\1O C(1"1~lJi\JICATInN EOUIDME~JT
71 E L crT Q. I C C CJ "11-' fJ I'J ~ \! r SAN 0 1\ C C E S S rJ R. I F S
7 Lt r..q c; C. f I F. C T!~ rr. A L MAC H I N F P. Y, E 0 tJ T. P M E f\j T II NOS U P P LIE S
HA
lOA
IDA
IDA
lOA
lOA
10/\
IDA
lOA
Ion
IDA
IDA
IDA
lOA
lOA
IDA
IDA
lOA
lOA
lOA
IDA
lOA
IDA
IDA
IDA
10/\
10/\
IDA
IDA
IDA
10/\
lOA
lOA
IDA
1 0.1\
IDA
lOA
IDA
lOA
!(lA
lOA
10/\
IDA
IDA
IDA
lOA
;' f\,
r; ~!~ "'~I
'I):,! 1
f'
I
,I, '. ,
"
05q
060
061
062
063
064
065
066
067
068
069
070
071
072
073
074
075
076
077
078
07Q
080
ORl
082
083
084
085
OR6
087
ORB
089
090
091
092
093
09/t
095
096
097
098
OQg
100
101
102
103
10/+
".;'
,I, c--
-------
~ .....
'N
:W
76 AI ~CR.AF' AN) PAR.' S
77 n.'Hf~ '~ANSDO~TA'ION EQUP~EN.'
78 SCIENTIFIC AND CONTR~LLING INSTRUMENTS
79 OPTICAL, GPHTHALMIC AND PHOTOGRAPHIC EQU!PMfNT
f30 I.1ISCELlAI\JEGUS ,"1 l\\JIIF/\CTtW ING
Pol TR!\NC:;PrJRTATION I\"JD ~-:AREHOUSING
R2 COMMIJNTC:\T!ONS, F:XC. RADD AND TV BROADCASTING
83 RA~IO "NO TV BROADCASTING
84 E L U: T R I C lJT I L IT.I F. S
85 r;l\S lJT!LITIFS
86 WATER ANO SI\NITA~Y SERVICES
87 WHOLESALE AND RETAIL TRADE
88 FPJANCE AND If\JSU~ II.NCE
R 9 R F. l\ L EST 1\ T F A \I f) R. F: I\! TI\ L
90 HOTELS, PFRSO~Al AND REPAI~ SERVICES EXC. AUTO
Yl RlISrNF:SS SfRVICfS
92 RESEl\~CH A~() DEVELOPMENT
93 AUTOMOHILE REPAIR AND SERVICES
94 AMUSEMENTS.
95 MFnJCII.L, ~DUCl\TInNAL SERVICES, AND NONPROFIT ORG.
96 FEn~RAL GOVERNME~T ENTEQoRISFS
97 ST~TE ~NO LOCAL GOVERNM~NT ENTERPRISES
9A GROSS rMDORTS OF GOODS AND SERVICES
99 BUSII\JESS TPAVEL, ENTERTAINMPH, t\f'JD GIFTS
100 OFFICE SUPPLI=S .
; >1 ..
'~T,¥'":
"J.
ICA .06
OA 107
IDA 108
101\ 109
IDA 110
IDA 111
IDA 112
IDA 113
IDA 114
lOA 115
101\ 116
[OA 117
IOI\ 118
lOA 119
IDA 120
10.11. 121
IDA 122
TiJA 123
101\ 124
IDA 125
InA 126
IDA 127
lOA 123
IDA 12g
IDA 130
IDA 131
IfJA 016
lOA 017
!iJA 01S
IDA 019
IDA 020
InA 021
IDA 022
InA 023
IDA 024
lOA 025
IDA 026
lOA 132
lOA 133
c.}
IDA 134
lOA 135
IDA 1~6
IilA 1.37
IIlJ\ 1,R
JO" 139
FJJ\ 140
:1 .~ ,,:; "
() I r-' t: t~ S I rJ ~J A!{ ( 1 5 ° )
DP"ENSIO~ Cf1LSU(lO)
DIMENSIn~ RC50,50),C(22,31),SUM(31)
DIMENSION LABfLCIO,50)
01MFNSrnN STORF(SO,3J,P(SO,50),SI50,SO) ,VECTOR(SOJ
COM~O~ ~K( 103,150)
C (H'~ :''10 N N ARK, I T EST, I R 0 ( 50 ) , I R 0 vI, S T E P
CO ;".;-10 N !\ IN ( 3000 ) , A ( 1 20 , 1 50 ) ,N C , N R , I C , I R
COM~ON KRO,LP~,KP
F () U I V AL FN C E (r", N,r R ) , ( S ( 1,1 ) ,A C 1 ,1 ) ) ,( 1\ K ( lJ ,VE C TOR ( 1 » ,
1 C K K ( 1 . 1 ) , p ( 1 , 1 ) ) , (1\ ( 51 ,51 ) , S T OR E ( 1 , 1 ) )
KqD = 5
LPR = 6
KP = 7
STFP = 0.0'
nn 10 1=1,50
10 IROCI) = I
IQrr:=o
\1IJLT=lO
Ifln')=()
-------
) 0 70 M = 1 , 12
I(K=:)
[F (M-12) 21,22,22
22 fv1UI.T=')
I h I)D= 4
JH[=540
21 0'1 25 J=1,108
I l. 0= ( J - 1 ) * ~~U L T + 1
IH[=JLO+IAf)f)
25 P.EIiJ IKRD,7!)3) (AIN(I),I=ILO,IHU
703 FOP.~I\T 110F7.1)
1)0 h5 1=1,10
65 CO L SUI T) = :). 0
00 611'1'/\01=1,106
00 60 M40=1,~ULT
KK=KK+1
60 COLSU(MAO) = COLSU(MAO)+AINIKK)
61 CONTI~UF
\oJ RITE (L P R , BOO) (C (1 L SUI J J) , J J = 1 , MU L T )
ROO F 0 R~'i l\ T (1 X ,1. 0 r- 10. U
f) 0 30 I t1 = 1 , i'W l T
IREC=IRCC+l
.1=0
00 '50 I=P1,.1HI,MULT
.1=.1+1
50 A(Ir.:EC,J) = AIN(I)
30 CO~H I "JUF
20 CONTINUE
~
N
~
AGGREGATES I/O MATRIX
Cl\LL :\GGPF.
IR1=IR - 1
IR2 = r~ + 1
I~3 = I~ + 2
ICI = IC - 1
WRITE (6,4)
4 FORI1;ATI 1H1,?5X,'THF. INDUSTRIAL SECTORS REFEREO
1 GENE~AL'/26X,'FLOW TABLE AND THE INPUT/OUTPUT
226X,'AS FOLL!)\-JS...'/II
1)0 3 J=1,Jp.l
REA[)15,l) (LARELlI,J), 1=1,7)
3 wRITEI6,n ILl\RFL(I,.1), 1=1,7),.1
1 FOR\1.i\T (7f}J~)
2 F n R ~-1 h T I 11~ ,30 X ,7 A It, 3 X, ,- SEe T OR ' , I 3 )
i~ HIT ElL f) K , 1 1 6) r R 1
: ~ {~ ;:: I ~ l'": ',:
r: );' ~
,/i.l~;\: T~ ~~
- ,. ,-. ,-
~ '
) ~ ,-.-, ';.:
:': ,; , r: I;' 1\,1
; (I
TO BY
TARL E
f: '; n .~"
OA
:OA
IDA,
t O.t:J..
10A
lOA
101\
lOA
I (l A
[0/1
lOA
IDA
lOA
lOA
lOA
101\
IDA
IDA
IDA
IflA
lOA
lOA
III A
InA
IDA
lOA
IDA
IDA
IDA
IDA
lOA
IOI\
lOA
IflA
IDA
IDA
IDA
NIJM~ER IN THEIGA
AQ.E DEFINED'/IOA
lOA
IDA
rrJA
InA
lOA
r Il"
InA
42
.43
lit 4
145
146
147
14.'3
149
1'50
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
1 73
174
175
176
177
178
179
180
1 B 1
1R2
1~3
18ft
18'>
186
1r.7
'J
n
t f:J ~:::
-------
r-
'RF-J"S nUT THF ~:\-I G:NE~AL =_OW TABLE
M = ()
501 00 601 I=I,IC
0'1 602 J=I.IR
."1 = M + 1
602 g ( I , J) = II. IN ( r~ )
601 CflNTJNUE
TABULATES THE COLUMN SUMS
00 603 J = l,IC
B(J,IQ21 = 0.0
DO 603 T = 1, I R 1
603 O(J,IR2) = (3(J,IR21
RIIC,IR) = 0.0
+ B(J,I)
I R OW = I R 2
ITEST = IC
STEP = 0
CALL PRINT«(3}
....
'N
'1J1
GENERATES THE AGGREGATE I/O TABLE
WRITE(LPR,655) JRl
655 F1R~AT(IHl,35X,'INPUT/OUTPUT
~o 606 N=I,IC .
CON = B(N,IR) + B(N,IRZ)
B(N,IR3) = 0.0
IF(CON)6,606,6.
6 DO 607 I=l,Iq
~ ( N, I) = R I ~J, I ) I C 0\1
627 BIN,IIU) = F~IN,rP'3) .. B(N,I)
607 CDNf!I~\J~
R ( N, I !< 2) = '1 ( N , r R ? ) I C J N
606 CONTINUE
TABLE (',I2,' - SECTORS}'//)
T R n'A = I R 3
fTEST = IC
STfP= 1
CALL PRfNT (R)
WRITE(6,6666) IR1,IR1,IRl,IRl
6666 FnR~t\T(/145X, 'NOTE..~OW-',12,'+l = VALUE ADDED (PERCENTAGE)'
1 151X,'R{rtJ-' .17,'+-2 = INTE~r~tDrI\TF P~Of)UCT (PERCFNTI\GF, SUM
./ : ':"i :f~
) I~ ':! I~;: I [ ~i,:
~ !~.:;.. ~\ " 'l
,; r!
OA .89
f(A .90
IDA 191
IDA lq2
10/\ 193
IDA 194
IDA 195
InA 196
IDA 197
IDA 198
lOA 199
IOL\ 200
10/\ 201
lOA 202
IDA 203
IDA 204
lOA 205
lOA 206
InA 207
IDA 208
fOA 209
InA 210
10t\ 211
IDA 212
10/\ 213
IDA 214
lOA 215
10/\ 216
IDA 217
lOA 218
IDA 219
IDA 220
10/\ 221
1011. 222
10/\ 273
IDA 224
fOil. 225
rOA 226
lOA 227
IDA 22A
IDA 229
IDA 230
10/\ 23L
!(J/\ 232
InA 233
OF RIC/\ 23't
I; I;;;
-------
OA
IDA
101\
111"
IDA
lEONTIEF INVERSE MATRIX OF THE INPUT/OUTPUT TAIOA
)'11> IDA
IDA
lOA
IDA
IDA
IDA
lOA
IDA
IDA
IDA
lOA
IDA
IDA
lOA
H1A
IOA
TOA
lOA
IDA
IDA
lOA
Infl
lOA
lOA
IDA
lOA
OOllARS',5XIOA
lOA
IrJA
lOA
lOA
lOA
InA
10"
IDA
InA
lOA
Inl\
lOA
lOA
THE LEONTIF-F INVERSE
WPITE(LPR,742) JRl
742 FfJRM"HIHl,?IX,'THE
ISLE 1',I3,' SECTORS
on 7701 1=1, IC 1
on 7702 J = 1,ICI
7702 fHI,J) = -[3II,J)
7701 H( [,I) = 1.0 ~ 8(1,1)
MA = IC 1
DETA = 0.0
CALL ~INVR (A,MA,DETA)
WRITE (LDR,773) DETA
773 FORMAT(IH ,'DFTERMINA~T=',FI2.10)
I ROW = r C 1
ITEST = ICt
STFP = 1.0
CALL PRINT(R)
-
N
0'
MOOIFIEO CONSUMPTION VECTOR
CON = ESTIMATED MAqGINAL INCREASE IN CONSUMPTION
~(IC,2) = PERCENTAGE CONSUMPTION OF EACH SECTOR
AK( [) = fACH SECTORS SHARE OF THE INCREASED CONSUMPTION
111 FOR!v1ATIF6.U
REAO(KRO,111) (ON
\-4 p, IT E ( 6, 6it 1) CON
647 FORMAT(lHl,18X,' CONSUMPTION BY SFCTOR'/25X,
I'BE~EFIT - CHANGE IN CONSUMPTION'
l'=',F6.1,' MILLION DOLLARS'1135X,~PERCENT',7X,'MIlLION
2'SECTOR OEFI"JITION')
00 643 I=l,IRl
AK(I) = CON~:I3( IC,I)
643 WRITE(6,644) I,BIIC.I),AKII),(LABEl(J,I), J=I,7)
644 FOR~ATIIRX,'SECTnR-',I2tFI4.6,8XtFI0.2,8X,7A4)
MULTIPLIES THE LEONTIEF INVERSE BY
1. MrlOIFIED CONSUMDTION VECTOR
2. CONTROL COST VECTOR
MN = IRI
STEP = ?
TROW = I R
JTfST"" IF
236
237
238
239
240
241
242
~43
244
245
2it6
247
248
249
250
251
252
253
254
255
256
251
258
259
2&0
261
262
263
264
265
266
267
268
269
270
271
272
213
274
21'>
276
2 1..1
218
279
280
281
. >
-------
OA
IDA
lOA
lOA
IDA
POLUTION ClNTROL COST - CHANGE IN FINAL OEMANDIOA
IDA
InA
lOA
IDA
rr1A
IDA
IDA
lOA
lOA
IDA
IDA
lOA
lOA
lOA
lOA
I f) A .
CONTROL. EFFECT nF COST ON STRUCTUREIGA
roll.
FJA
IDA
IDA
AIR POLUTION CONTROL. BENEFICIAL EFFECT ON STRUCTIOA
ECONO~Y'II) lOA
IDA
IDA
IDA
InA
In/\
IDA
IDA
lOA
roll.
IDA
JO/\
IHl,30X,' AIR POLUTION CONTROL. NET FFFECIOA
NATIONAL ECONOMY'//) 10/\
IDA
InA
STRUC lOA
I [J 1\
MA RK = 1
GO Tn 444
1000 MARK = ~ARK + 1
WI~ITE(6,777)
777 FOP'-1AT{IHl,25X,'AIQ.
1'111
READ(5,1!) (VECTOP.(J), J=l,MN)
71 FOR\1A T( 30X,' SFCTOR-' f 12 ,Fll.3 ,4X, 7A'i')
11 FOR~I\T(FI0.2,70X)
r)o 7777 J=I,MN
':IRITE(!),7l) .J,VfCTOR(J),(U\BEU I,J), 1=1,7)
7777 CONTINUE
444 CONTIr\IU'=
00 23 IK = I,MN
5 TOP, E ( I K , !"1 AR K) = 0.0
00 ? 3 J = 1 , ~<11\1
S e J, I K) -:: R ( J , I K) f,c VE C TOR ( J)
STORE(IK,f'1ARK) = STORE(IK,MARK) + SeJ,IK)
23 CONTINUE
IF eM!\RK-l) 32,31,32
32 \tJQ. JTE (6,33)
~31 FOH:>1I\T(IHl,?'OX,' AIR pnLUTION
-.J 1. or Nt\T TUNAL ECONOMY'I I)
GO TO 34
31 WRITE(6,35)
35 F n R i., A T ( 1 HI, ? 0 X , ,
1 U~ E nF NI\ n ONAL
DO 16 ICT = 1,MN
00 36 Ie P = 1, Mi~
36 P(ICT,ICr) = S(ICT,ICP)
34 CALL PRINT (S)
IF (MARK-2) 1000,,7, ?OOO
3 7 f'1 ARK = ,
DO 38 1= 1,f'.1N
STnR~(I,3) = STORE(I,l) + STORE(I,2)
DO 3 A J = 1 , ~-1N
38 S ( I , J) = P ( I , J) + S ( I , J )
\.J~ ITE (6, 330)
: 330 FCP~Ml\T(
IT O~ STQUCTURE OF
GO TO 3'-t
2000 CONTINUE
Rl FORI>1r\T( 1Hl,7.0X,'SUMf~"RY OF FFFECT OF AIR POLUTION. CONTROL ON
1 T U \~ F 0 F N A Tl n N 1\ LEe 0 N C):'.1 Y , I I) a x " ( I N MIL L r 0 N S ) , I I
.---;, ."
J,: ,.' \"0 ,')
:, "
,.:1.
~:;I
.c:. I~ :..; ,0
.~ 'il ~
2R3
2R4
285
286
287
288
289
290
291. .
292
293
?CJ4
795
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
32,
324
325
326
327
32R
"-,('
,
-------
23X,' ~E'\!EF.. ',3X,'
C:S., ,7X, 'EFFECT'I
:. A
HA
IDA
IDA
lOA
lOA
lOA
IDA
101\
IDA
lOA
IDA
101\
lOA
IDA
lOA
10/\
101\
10,"
IDA
lOA
101\
BEEN READ IN AND STORED lOA
FILE CARDS IDA
101\
lOA
10/\
IDA
IDA
101\
IDA
COLUMN lOA
lOA
lOA
InA
lOA
lOA
IDA
101\
10./\
101\
lOA
101\
101\
lOA
lOA
00 2222 K=1,3
7222 STORE(IR2,K) = 0.0
W~ITE(6,Rl)
o 0 2 ? 2 I = 1 , ~~ N
hi RITE ((), 4 l-t ) 1,( 5 TOR E ( 1 , MAR K) ,M ARK = 1 , 3 ) ,( LA BEL ( J, 1) , J = 1 , 7 )
DO 222 K=1,3
STORF.(IR2,KI = STORE(IR2,K) + STORE( I,K)
222 cnNTI"II.Ir:
WRITF.(6,4444) (STORE(IR2,~/\RK),MARK=1,3)
44 FfJRMAT(ZOX,'SECTOR-',I2,3F12.2,5X,7A4)
4444 FOR ~A A T ( ~ 4 X , ' TOT 1\ L ' , 3 F 1 2 .2 )
5 TOP
E"JD
SURROUTINE AGGRF-
DIIv1["ISION ~(150) ,JCOM( 100) ,KCm1C 150)
COMMON BK ( 108, 150)
C:OMf>10"! '1ARK, ITF-STdRO! 50), JRm.J, STEP
CO M,-"O N A I N ( "30(0) , I'd 1 20 , 150 ) ,NC, N R , I C , I R
COMMON KRD,LPR,KP
. ~PROGRAM ASSUMES MATRIX TO BE AGGREGATED HAS
: coBY COLUMNS ON 01 SK AND AL Sf") THI\T THE DEF INE
IN THE MAIN PROGRAM ARE I\S FOLLOWS
NC = TOTAL = OF COLUMNS IN INPUT MATRIX
NR = TOTAL = OF ROWS [N INPUT MATRIX
IC = TOT~L = OF C:OLU~NS IN OUTPUT MATRIX
lR = TOTAL = OF ~OWS IN OUTPUT MATRIX
** W~RNING ** THF LAST ROW IN THE OUTPUT MATRIX MUST BE
lOR (VALUE ADDED) IN THE INPUT TBL. THE LI\ST
IN THE OJTPUT MATRIX MUST BE COLUMN 108
(CDNSUMPT ION I IN THE INPUT MATRIX.
READ (K~D,7001 NC,NR,IC,IR
\~ RITE (L P R , 70 0) N C , N R, I C, I R
00 5 [ = 1, I C
Of) I) J = 1,.\JR
i3K(I,J) = 0.0
5 CONTI NIJF
700 F[)R~1A T (2014)
IF (IC-IR) 11,10,10
10 I~AX=IC
t\JM fiX = NC
Gn fl:J l?
11 Irv1AX=IR
t.j~ /'.
f.'i,j
330
331
332
~3~
334
335
336
337
33R
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
3?1
372
373
37/.,
375
'C>
-------
12 ~EAD K~D,700} (JCm~(I ),1= .d~AX
READ «RD,700) (KUli"1n ),I=1,NMAX
K=O
[")(1 '30 I=l,Ir.
J J = J C 0:-1 ( I )
DO 35 J=1,JJ
I<=K+1
IRFC=KCflM(K)
DO 16 L=1,N:t
36 R K ( I, L) = BK ( I ,l ) {,A ( I K EC , L )
35 r:ONT I r~lJF
30 CONTINUE
ITOP = IR*IC
fFJ 41 l=l,ITOP
41 AIN(I 1 = 0.0
M=O
DO 40 1= 1, IC
K=Q
I) D ItO J = 1 ,I I{
JJ=JCOM(J)
M=M+1
1)0 43 l=t,JJ
~ K=K+l
-.J::> IF(K-10'".}) 48,47,47
It7 WRITE(6,491 K
49 FOR:~,~T(lHl,'E~ROR...T~IED TO ADD ROW-',I3,' AND THERE
1 ' / 1 X, 'r~ I: A N IN G F UL ROW SIN THE I N I TI A L I /0 TAB L E ' / / )
STOP
48 CONTJNUE
K K =1«( OM ( K )
4, l\INCMI = ATNCi'O&!3K(I,KK)
40 C (1 NT I 1\)t.J r
OUTPUT MATRIX ~RRANGED IN AIN BY COLUMNS
RETtJRN Tn t-1AIN PD.OGPAM TD WRITE AGGREGATED MATRIX
RETURN
F ~J f)
SURP.OUTiNE PRINT(S)
DIMENSIr'N B( ')0,50)
Cowmr\! SK( lOR, 150)
C (1 WI!J f\J rJ, ARK, IT EST, I R (1 ( 50 ) , T ROW, S T E P
CCJWHV~ AIf\J(3000),A( 120,1501 ,NC,NR,IC,IR
CO~~ON KRo,Lr~,KP
700 FnR\1AT CIH .n:,9C5X,'COL',I3)
7 ()? FOR '~A T C 1 X , ' R r; '".J' , I ., , 2 X , 1 0 F 1 1 . 0 )
7 0 <; FOP f ~ 1\ T C) x , I R 01001' ,r 3 , 1. X ,t 0 r 1 1 . 4 )
7 0 ~ F n R '~f\ T (1 x , I ~ () ',oJ I , r 3 , 2 x , 1 0 F 1 1 . 7)
;'
:;::: 1-"
, "
! ,', 1.
~- ~ c,: j
( A 377
:nA 378
IDA 379
lOA 380
lOA 3R1
IDA 382
IDA 383
10/\ ,84
lOA 385
lOA 386
IDA 387
IelA 388
lOA 389
IDA 390
18A 391
lOA 392
IDA 393
lOA 394
IDA 395
lOA 396
IDA 397
[0 A 398
lOA 399
IDA 400
lOA 401
ARE ONLY 108IOA 402
101\ 403
IDA 404
10/\ 405
lOA 406
lOA 407
IO/\ 408
IDA 409
lOA 410
lOA 411
IDA 412
lOA 413
10/\ 414
lOA 415
lOA 416
lOA 417
IO/\ 418
IDA 419
lOA 420
rrJ/\ If 21
10/\ 422
f~I';: 0:;'1
-------
I < = 1
MARK = (ITEST&8)/9
on R02 J=1 ,~ARK
IJ = J>:'9
IF([J - JTEST)~05.904,804
R04 IJ = ITEST
805 \4RITF. (U'R,700) (IRO(I)d=IKdJ)
IFCSTEP-l.0)333.332.334
333 on 801 v=l,IROW
801 vi R I T F: (L P R ,702) I PO C M) ,( fH I ,101 ) d = I K,I J )
GO TO 99
334 00 822 M=1,tROW
32? WR ITE(lPR,705) IRO(M), (B( I,M), I=IK,IJ)
GO TO 99
132 00 811 M=l,IROW
1;/ R I T E C L P R , 703) I RO C M) , C B ( t . M ) , I = I K, I J )
i 811 CONTI NUE
99 CONTINUF
IK = tJ & 1
\~RITE CLPR.701)
I R02 CQNTI NUF
. 701 FOP'1AT( lHl)
RETIJRN
F. f\I 0
5URROUTINE MTNV~ CA,MA,DETA)
f) I MFN S ION A ( 50,5 0) , I R ( 50) . I C C 50 )
COMMON nK(10B,150)
COM~ON "1ARK, ITf=ST.IRO( 50), IROW,STEP
C 0 w~o N III N U 0 00) , [) C 120 , 150 ) ,NC . NR , I D ,I K
COMMON KRD,lPR,KP
on 1 I = 1 . '<\ A
IR(I)=O
1 Ie ( [ ) =0
DETA = 1.0
TOl = . a00001
5=0.0
R = M/\
::'BEGIN*
2 1=0
J=O
TFST=O.O
[)~l 1- 5 K = 1 , MA
IFCIR(K»25,39,25
19 CONTINlJJ:
DO ;U.. l= 1, /-1~
IF(ICCL»?4,3R,24
11,°-
, J ,; J~
(A 424
IDA 425
InA '.26
10/\ 427
Ir)A 428
IDA 429
LOA 430
IDA 431
InA 432
lOA 433
rOA 434
lOA 435
IDA 436
lOA 437
IDA 438
lOA 439
lOA 440
lOA 441
IDA 442
IDA 443
IDA 't It 4
IDA 445
IDA 446
lOA 447
InA 't48
lOti 449
lOA 450
lOA 451
IDA 452
IDA 453
IDA '..54
InA 455
IDA 456
10/\ 457
IDA 4.58
rnA It 59
IDA 460
lOA 461
rnA 462
IDA 463
IDA 464
IDA 465
IDA 466
IDA 467
lOA 468
IOA 469
r 'i /. 1(,
~
w
o
-------
X=AI3S(lUK,L) )
.:CX-TES" ~(t,37,37
37 CONT I NU:::
I=K
J=L
TEST=X
24 CONTINUE
25 CnNTINUF
THIS IS THE TEST FOR A DEPENDENT TOLERANCE
I F (TO L ) 73, 74,74
7, TOL=A8S{TOL*TEST)
74 CO NT I NU t
,~
'-~ $
*END"~ SUBMXS $$
.~
......
w
......
PI V= A C I , J )
DETI\=PIV~:DETA
IF CAR S C P I V ) - T rJ L ) 1 7, 17, 35
35 CONTINUE
IR(I)=J
IC(J)=I
nrV=l.O/PIV
ACI,J)=PIV
rJU5K= 1, ~"'A
IFCK-J)36,5,36
36 AC I,K) = I\CI,K)*PIV
5 C (] NT I NIJ E
D09K=1,MA
TFCK-I )33,9,33
33 CONTI"JUF.
PIV1=ACKtJ)
6 DnRL= 1 t r4A
IFCL-J)31t8,.,1
31 ACK,l) = ACKtL) - PIV1*A(I,U
A CONTI NIJF
9 CONTI NUE
OOllK=1.MI\
IF C K- I ) 32, 11 ,32
32 ACKtJ) = -PIV*ACK,J)
11 CONT I I-,JUF
$=5':1.0
IFC<;-P,)2t12,12
12 :)(,16J=1,MA
1<=1((1)
e'L ['~. i: r ;.
~
( A
lOA
IDA
lOA
lOA
lOA
lOA
lOA
. lOA
TrJA
lOA
IfJA
IrJ/\
IDA
IDA
101.\
lOA
10A
lOA
lOA
lOA
J 0/\ .
lOA
lOA
10/\
lOA
InA
lOA
lOA
JrlA
101.\
lOA
lOA
lOA
IflA
IDA
10A
lOA
lOA
lOA
lOA
lOA
lOA
IDA
lOA
lOA
:1 f~1
47.
472
473
474
475
476
477 .
478
479.
4AO
481
4H2
483
484
485
486
487
488
489
490
4(11
492
491
494
495
496
497
If 9R
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
51'5
'516
-------
=CK-)41, .6,4.
4. CONT ;-.,Jl~
OETA=-[)FTA
~F"n It l r, 1 . M A
Tr:i'~D=A(K,L)
A( K,L )=/\( I ,L)
1 4 A ( I, L ) = T F to1 P
D015L=1,MA
TEIv1P=ACL,IvI)
!\CL,M)=!l.CL,I)
15 AC L, t )=TEMP
ICC''')=K
IR(K)=~
16 CONTINUF
LE~VE WITH A WARNING OF ERROR
17 If\JVFR = S
RETURN
END
IJA
:C A
YOA
InA
! t1 II
In"
IDA
IDA
IDA
101\
IDA
IDA
IDA
IDA
lOA
. 101\
IDA
lOA
.....
'w .
N
,
5 .8
5 .9
520
5Z1
522
523
524
~25
526 .'
51.7
,28
529
530 .
531
532
533
534
535
-------
o
4.0 PROGRAM FBE (Feedback Effects)
4. 1 Introduction
Program FBE contains a regional market share matrix devel-
oped in Appendix D of Volume I of CONSAD's May 15 Report.
The purp,ose
of the program is to take the output from Program lOA and then distribute
the interregional feedback to each AQCR.
These interregional feed-
backs are measured in terms of changes in value-added by two-digit
SIC detail for manufacturing indus tries.
Output from Program FBE
also punches a deck of cards which will be part of input deck to the
Program RMS to be used for simulation option 8 as described before.
Detailed specifications of the program preparation is given in
the following sections.
133
-------
4.2
Job Deck Set- Up
The Program FBE (Feedback Effect) is a 109 statement FORTRAN
program which uses approximately 60K of core.
'-
/,..
'>'
Includes three sets
of input data cards
(248 cards)
Data Input Deck
/ / GO .SYSIN DD>:~
109 cards
Program F BE
/ / Job
Job Control Language (JCL)
Cards
134
/
-------
4.3
Data Input Deck
This deck includes three sets of input data cards, a total of 248
cards, arranged in the following manner.
1.
The ratio of value -added /value of shipment by sector for
42 sectors under consideration (6 cards)
Format: 8FIO.4
Output from Program lOA, 42 cards, containing the cost and
benefit effects of air quality control measured by changes in
value of shipments by sector (one card/sector).
2.
Gols. 1-10
Gols. 11-20
F 1 0.4, benefit effects
F 1 0.4, cost effects
The regional market share matrix, 200 cards, contains the
share of national value-added by industry for each of the
100 AQCRs under consideration.
3.
First 100 cards contain market share of first 11 two -digit
SIGs by AQGR (SIC 20-30) (one card for each AQGR).
Gols.
Gols.
Gols.
Gols.
Gols.
Gols.
Gols.
Gols.
Gols.
Gols.
Gols.
Gols.
1-5
6-10
11-15
16-20
21-25
26-30
31-35
36-40
41-45
46-50
51-55
56-60
15, AQGR code
F5.4, share for industry 20
F5.4, share for industry 21
F5.4, share for industry 22
F5.4, share for industry 23
F5.4, share for industry 24
F5.4, share for industry 25
F 5.4, share for industry 26
F5.4, share for industry 27
F5.4, share for industry 28
F5.4, share for industry 29
FS.4, share for industry 30
The second 100 cards contain market share of second nine
two-digit SIGs by AQCR (SIC 31-39) (one card for each AQCR)
135
-------
Cols.
Cols.
Cols.
Cols.
Cols.
Cols.
Cols.
Cols.
Cols.
Cols.
1-5
6-10
11-15
16-20
21-25
26-30
31-35
36-40
41-45
46-50
15, AQCR code
F5.4, share for industry 31
F 5.4, share for industry 32
F5.4, share for industry 33
F5.4, share for industry 34
F5.4, share for industry 35
F5.4, share for industry 36
F5.4, share for industry 37
F5.4, share for industry 38
F5.4, share for industry 39
4.4
Description of the Program Deck
Program FBE is a FORTRAN program consisting of 109 state-
ment cards.
The cards are sequenced along the right hand edge
(columns 73 -80) with the prefix FBE.
The program begins by reading in data set one (the ratio of value-
added/value of shipment).
Immediately thereafter the cards produced
by lOA (cost and benefit effects of air quality control) are read.
The
program then multiplies the cost and benefit for each sector by its re-
s.pective ratio of value-added/shipment in order to convert the value of
shipment by sector from Program lOA to value -added.
The next section of the program aggregates the 42 cost and bene-
fit figures calculated above into 20 two-digit SICs and one "all other'!
industry classi£ic~tion.
The program now.punches and writes output
~ards containing (1) the benefit and cost to "all other" industries (1 card),
a.nd (2) the. cost and benefit to the 20 two-digit SICs (20 cards).
136
-------
L-
!
Next the program reads in the regional market share matrix
(data set 3).
The program processes one AQCR at a time for the next
set of calculations.
For each AQCR, the cost and benefit for each two-
digit SIC is multiplied by the share for that AQCR/SIC combination in
order to determine the cost and/or benefit effect associated with in-
dustry in that AQCR.
The net effect is calculated as the sum of the
cost and benefit figures.
When all 100 AQCRs have been processed, the program punches
card and writes output in three steps.
First, the cost effect by AQCR
is punched and written, seven two-digit SICs per card (or line) fol-
lowed by the AQCR code.
Thus, three cards are punched for each
,
AQCR.
The format used is (7FIO.4, 3X, 13).
(300 cards)
Second, a deck, identical in size and format, containing the bene-
fit effect is punched and the output for benefit written.
Finall y , the
deck containing the net figures is punched and written, and the pro-
gram terminates.
137
-------
Start
Read ratios of
value-added/
value of ship-
ment
Read cost,
benefit as ap-
plied to value
of shipment
Multiply ratios >:'
cost & benefit.
figures to convert
to value -added
Aggregate 42 sec-
tors to 20 sectors
and one "all other"
1
138
FLOW CHART OF
PROGRAM FBE
~
-------
1
Write resul-
tant cost and
benefit fig-
ures
Read distribu-
tion matrix
Calculate cost &
benefit share for
each AQCR/SIC
combination
Punch resul-
tant shares
Write
results
End
139
FLOW CHART OF
PROGRAM FBE
(continued)
-------
4.5
Program F BE (Feedback Effect)
."
140
-------
~: 1.: ::~ ':~ ::= ::: '::: ~:' :::~ =:' ;:= =:: ::;:: ;:: :::: :~ :-;:~: ,~ :;-= :::: =:: ::= ;:-: ~:: :;: ~:': ::-: :::: :::' ::~::: )::::= * f.: ::-:: :::: ~: ::: :::: ::: I.:: ~:~ =:; .::: ::: ;::: ;;-: ~:::: ::::::: ':: ::: 1,~ :!.: =-:: ::: 1.:::::
,}:
::;c
p r:> r. r, D c.. t" !- '\ F:
P~FPARf~ py rrNSAD DEsr~RCH COPP~R~TTn~
f\P:'!IL 107]
::<
,}:
,>:
...
THIS pr , 2 0 ), I) i S T R, 1\ ( 1 0 () , ? 0 )
i) 1'1 T " C 'J S T, G F \! f': / ? (' ':: n . (', t C ~: 0 . r; /
I) I . 1 r t'l S I U (\! D r S T L1 C I t ~ 'J , '2 (I )
nATA v / '2 * 0.0 /
1.0 F:l[)MAT(RF1.0.I-d
~::L\[1 (S,lC) 0(1), r = 1,4?)
J ,l ? (' I = 1, It ?
"'E I\() I 1),1')) Y(I), l(T)
Y ( I) = Y IT) ':: X ( T )
lIT) = llf) ':' XCI}
.... ? 0 C IJ 'H I r,: l J !=
t nn ~o I = 1,q
.V(J) = V(i) + YITI
l~ VCZ} = V(?) + 7( I)
cnST(l) = Z(10)
lE:-Jf(l) = Y(lO)
CnST(?) = lCll)
'1F\J[(?) = YI1I)
Ci)STL~) = .1) ,:' ZC]?}
,3ENECl) = .1) ::' Y(12)
cnST(4) = .s * 7(1;)
:3 F ~.! [= I 4) = .!) ':' Y ( I? )
in l-t [i T = 1:3 , ? 1
J = I - q
CflSTIJI = l(TI
!3c~IFIJI = Y(f)
1+0 CIl;n f\H!F
nn 1)0 I = 2?,71)
enST( lit) = criST I 14)
~r:\Jr {14} = '\t\Jt.:( J/~)
') c; Ul ~H I :--1 ' J r:
j:l A0 T = ?f,,7'i
J = r - 11
C ") ,; T ( .J) = 7 If)
.q: ~ II ( _I) ." V ( n
+ Z ( T )
+ Y ( T )
~:
::'
:~
':<
~:
y.:
):::
:::
::::
=p, E 1)01
FR F 002
FKF 003
Ff,[ 0 C'/~
FRF 00')
FBF 006
Fn[ n07
FRE COP,
FilE 009
FBF. 010
FRF OIl
FRE 012
FRE 013
FflE I) 1I~
FBE 015
FRE 016
FH[ 017
FP,F 0]3
F Ii F. O]q
F [)-E 020
FRE 011
FRE 022
F[)E 023
F 11 F 074
FP,[ 02')
FP. E 0'2",
FBE n27
FBE 02'1
FBE 029
FF\C 030
FBE 011
FRF 032
FGE () ~3
FB F () 3 It
FRE ()3')
FfH: O~6
FRf.: 017
Ff~ F 018
FBE 03')
F8 E 040 'J
FRF 041
fi3 E () 1+ 2
FRF (t/+ 1
F f\ F n 't It.
FRE 0 't ')
FR~ 04f-,
F P r~ r)I/7
-------
..
FR E 04!3
F1E 049
FI3E 05.J
Ff3E ')51
FR F 05~
FE'E 051
FGE 054
FRF 0')5
FRr: 01)6.
FgE 057
FP [ 0')8
FB E C59
FBE 06')
FRE 061
FBI: 062
F-B<: () 6 '.
FR F (A4
FP,E ')6'3
FI3 E () 6 /)
FRf 0/)7
F~E ~6g
FB E 069
FB E 070
FRE 071
Ff1E 'J7?
Ff3E 073
H\f 074
FI3 E 075
FBE 076
F p, E 077
FB E: 078
FR E 079
FI3E 0'3')
FRF: OHI
FRF OR?
FR E 0Bl
Ffl E 1)P'4
Ff1!: 08')
I=R [: Or.,s
FP,E l)R7
F r, E ()B8
F8!: 08Q
F~; E 0(1)
F[3f 0'11
rRE O'l?
F r, r:: ')T~
Ff\f- ~1 f.) 't
I) C C (1 'I)T T ~JlII=
CDS 'I~) = 7.(20) -I- 1(30)
3[i'JElli) = YIZ ) I ') T <; T ,-. "J, ( I , J ), J = Y, 1 '5 ). I
\',10 r T F 16. 11. Z) I ,) 1ST '1. G I I , .1 ), .J = 1. A, ? C), I
'..J PIT [ ( 7 , ]/i 7) ( i) I S T '{!1 ( I , J ), .J = 1 (, , 2 ') ), I
1 5 r. UJ \!T.r . i: ,r
iY\ 1 f) f) T = 1, 1 ,., ('\
. J :", T T:- I '-" 11, ()) : j 1ST P 1\ I I , 1 ). (r) 1ST ~ I' ( I , .J ), J = '3, 8), I
, J '{ I T r: ( 7 , 1 4 ~j) I') 1ST P '\ ( T . 1 ), (r;] S T I' 1\ ( J "J ), J = .~, ~~ ), I
vi" t T r I f I, II, r:) I:) r S T '} 1\ ( J , .J ), .J = (', I C', ), 1
'tJl, [ T I I 1 , 1/, :J) (:l 1ST i> f\ IT, .J ), .J = '), 1 ') ), I
./ P, 1 T r I h, 1/, (') (In s 1 '.' 1\ ( J , .J ), .J = 1 (, , ? 0 ), I
i" r 1 r (r ,) 1 I, -, \ () J ',1 rl{ I ,J \ .J ..,. 1 f I ? (' ) I
,~ il E /oj E I J )
':' [(1STI J)
+ OI STPfl,{ T ,J)
-------
160 CGI\;T 1 "11 IE
f1n 1. 7 C I = I,] (;i")
.~,... I T f- ( h, 1 4 (') r'! C) T ~ C ( I, 1 ), (') 1ST D ( ( I , J ), J = 3, R ), I
'<1 q I T r ( 7 , 1 4 C;) D I <, T R C ( I , 1 ), ( r) I <) r ~ r. ( I , J ), J = " 13), I
H r::. I T F ( h, 1 l, 0) (f) r C; T i~ C ( I , J ), J = C), 1 S ), T
I..J Q I T ~ ( 7 , J l-t 5 ) ( :) T S T <:' r ( T , J ), J = C), 1 'j ), I
,~r) r T f" ( '" 14?) () I S T i} r. ( I , J ), J = 16,70), I
;-JRIlr(r,1/.7) ('iIST'J((I,J), J = 16,?C), I
1 7'J r: (I \I T I ". I I r
140 FORM/IT( 1'-j ,7FJO.4, ~x,I3)
145 FOP'1!', T( 7J:ln .4, 3X, 1'3)
I/-t? F (I P ~i t. T ( 1 H ,5 F 10 . It , ? 3 X , T "3 )
147 r:fJK:'v'-\T(5F1G./i,2 n, rl)
<;Tnr
FNf)
FBF 0g5
FB E I) 96
Ff3E n97
FI1F 09~
F[-\E 099
FRF: lOr)
FRF 1111
r-rJC 102
r- f, ..
FBE 1.01
FRF: 104
FRE 105
F~E 10h
F~F 1.07
FBE lOR
FI3 E 1.09
.....
,j:>.
w
-------
APPENDIX
EPA CONTROL COST ESTIMATES
..
144
-------
The control costs used in this report are the preliminary estimates
corresponding to the control implied by Clean Air Amendment of 1970,
provided by EPA for this study.
These estimates are termed as EPA
cost estimates in the report.
However, the EPA cost estimates are
given in the form of national estimates by industry category.
The con-
trol costs by AQCR were then estimated by the production capacity of
the corresponding industry and use of electricity in each AQCR in pro-
portion to the national estimates of control costs.
The EPA cost estimates reflect the air pollution control of sta-
tionary sources for five pollutants, namely, particulates, oxides of
sulfur, carbon monoxide, hydrocarbon, and oxides of nitrogen.
The
stationary sources covered in these estimates are: >:<
Industrial Proce s s
Grain Milling and Handling
Kraft Pulp
Nitric Acid
Sulfuric Acid
Phosphate
Petroleum Refining & Storage
Asphalt Batching
Cement
Gray Iron Foundries
Iron and Stee 1
Primary Copper
Primary Lead
Primary Zinc
Primary Aluminum
Secondary Non-ferrous Metals
Stationary Combustion
Electric Power
Other Stational Combustion
Solid Waste
SIC2042
SIC2611
SIC28l9
SIC287l
SIC29ll
SIC295l
SIC324l
SIC332l
SIC3323
SIC333l
SIC3332
SIC3333
SIC3334
SIC334l
>,'< See Table 3.3, Volume I of this report.
145
-------
1- -
Regional cost estimates are, in general, calculated by the pro-
duction capacity of each four-digit SIC industrial categories in each
AQCR in proportion to national production capacity of the corre spond-
ing industries.
The measurement of regional production capacity is
based on the value-added definitions of GBE.
Thus:
N
C = C
ijt jt
V..t
-2:L
N-
V.
Jt
w he r e C.. is control cost of industry j in year t in ith AQCR
1Jt
N
Cjt is control cost of industry j in year tin the nation
V. .t is the value-added of industry j in year tin ith AQCR
1J
VN is the value-added of industry j in year t in the nation
jt
Ideally, such cost estimates should be based on four-digit SIC
industry categories.
However, in some occasions, a three-digit or
two-digit production capacity estimate (that is Vijt/V~) has been used,
because, in some AQCRs, 4- or 3-digit information is not available
to the public. ~<
~< There is a disclosure problelll of making available any indus-
trial data in the census publication unless at least 4 or more plants
are located in the same geographic unit in order to protect each indi-
vidual firm in the defined area.
146
-------
Regional control costs of steam electricity were estimated by
total usage of electricity in the region rather than production capacity
Besides industrial usage. the residential use of electricity was also
estimated in proportion to population.
. Other stationary combusion and solid waste costs were allocated
to each AQCR according to manufacturing production capacity to the
nation.
..
. 147
------- |