FINAL REPORT
VOLUME II: REVISED USER'S GUIDE
Prepared for:
Environmental Protection Agency
Office of Air Programs
Research Triangle Park
North Carolina
January 7, 1972
-------
FINA L REPORT
VOLUME II: REVISED USER'S GUIDE
THE OAP REGIONAL ECON011IC MODEL UTILIZA TION --
PHASE I
Prepared for:
Environmental Protection Agency
Office of Air Programs
Research Triangle Park
North Carolina
Prepared by:
CONSAD Research Corporation
121 North Highland A venue
Pittsburgh, Pennsylvania 15206
January 7, 1972
-------
TABLE 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
Simulation) .,. .
.......
9
2. 1. Int roduction. . . . . . . . . . . . . . . . . . . . 9
2. 2. Job Deck Set - Up. . . . . . . . . . . . . . . . . . 10
2. 3. Data Input Deck. . . . . . . . . . . . . . . . ... 11
2.4. Description of the Program Deck. . . . . . . . . . 32
2.5. Program RMS (Regional Model Simulation) . . . 38
3. O. PROGRAM lOA (INPUT-OUTPUT ANALYSIS).
. . . .107
3. 1. Introduction. . . . . . . . . . . . . .107
3.2. Job Deck Setup. . . . . . . . . . . . . . . . . .108
3.3. Data Input Deck. . . . . . . . . . . . . . . . . .109
3.4. Description of the Program Deck. . . . . . . . . .111
3.5. Program lOA (Input-Output Analysis). . .115
4. O. PROGRAM FBE (Feedback Effects) . .. .'. . .
. . . .
. .116
4. 1. Int roduction. . . . . . . . . . . . . . . . . . . . 128
\
4.2. Job Deck Set-Up. . . . . . . . . . . . . . .129
4.3. Data Input Deck. . . . . . . . . . . . . . . . . .130
4.4. Description of the Program Deck. . . . . . . . . .131
4.5. Program FBE (Feedback Effect). . . . . . . . . . .135
APPENDIX: EPA Control Cost Estimates. . . . . . . . . . . .139
ii
-------
"-~---- "---- --- --~" -
1. 0 USER'S GUIDE: A GENERAL DESCRIPTION
1. 1 Introduction
This is a revised user's guide of the OAP Economic Model System
to facilitate an efficient assessment of control strategies. ~:~ As des-
cribed and demonstrated in Volume 1, air pollution control strategies
are first transfQrmed into a set of appropriate model inputs. ,;<>:' The
computer program accepts the se inputs as "exogenous shocks II to the
model system and generates as output, from the model system, meas-
urements of the "changesll in the key economic variables in the AQCR IS
under study.
The computer simulation program of the OAP Economic Model
System includes three programs, namely, Program RMS (Regional
Model Simulation), Program lOA (Input-Output Analysis), and Program
FBE (Feedback Effects).
As indicated in Volume 1, the input-output
model and the interregional feedback effects were introduced in order
to measure the external market of an AQCR.
The main part of the
simulation program, consequently, is the Program RMS -- the simu-
lation tool for the Regional Economic Model.
~, CONSAD Research Corporation, An Economic Model System for
The Assessment of Effects of Air Pollution Abatement, Volume II,
User's Guide to the OAP Economic Model System, May 15, 1971.
>:,* A Iso, see CONSAD, op. cit., Volume 1. Development and
Demonstration Phase.
1
-------
Program RMS has been revised so that the growth rates by
AQCRs could be included to provide economic projection of future years.
In strategy simulation, cumulative effects over time of any given
strategy were also included.
At the same time, a new set of control
costs which are the preliminary estimates corresponding to the control
implied by Clean Air Amendment of 1970 was provided by EPA.
The se
estimates are termed as EPA cost estimates in the rest of the report.
,
.,
2
-------
1. 2 A General Description
PROGRAM RMS
The Regional Economic Model developed in Volume I cons is ts
of five major modules of equations, namely:
a.
manufacturing industry equations
income and regional economy equations
b.
c.
other industry equations
d.
labor market equations
fuel demand submodel equations
e.
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 I'simula-
tion options" which are related to a number of control inputs cor res-
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 te rms of increase in the regional
disposable income by AQCR, resulting from changes in property values,
savings in health expenditures, etc. from cleane r 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 increase to industry.
This option takes the price increases for the electric power which is
, being passed to industries.
Option 5: Electricity price increase to the residents.
This option takes price increas es for the elect:dc powe r indus try which
are being passed to the residents as a reduction in their disposable
income.
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: Fuel price change to two-digit SIC.
This option takes price changes in fuel oil, natural gas, coke and coal
as control input.
Option 8: Feedback Effects.
This option takes the output from Program FBE to measu re the inte r-
regional feedback effects.
Each option
consists of modules A, B, C, D, except 7 which
includes fuel demand s ubmodel (module E).
Howeve r, depending upon
whether aggregate or two-digit SIC manufacturing sectoring and other
4
-------
simulation specification are available, ten blocks of equations, each
related to one module, were included in this Program RMS.
A block
flow chart for the simulation options available in Program RMS is
shown in Figure 1. 1, and Table 1. 1 shows the relation between-simu-
lation option module of the Regional Economic Model and equation
blocks included in the block flow chart.
For example, in simulation option 1 the industry control costs
are read in fron~ EPA cost estimates control input data.
This leads to a
I
reduction in both profits and investment of the manufacturing industry
thus affecting the entire manufacturing industry module, and in turn
I
I.
I
I.
regional income, consumption, government sector of module Band
other industry of module C and finally labor market of module D. If other
simulation options were tr ied, the sequence will refe l' to corres ponding
blocks in the Program and tlnettl effects will be summarized in the pro-
g ram output.
Usually each strategy simulation will involve more than
one simulation option.
An example of the simulation 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 lOO sector
i
,
,
':
u.S. input-output table to a desired num.ber of sectors. Currently,
5
-------
(j'
Set 3 E
Fuel submodel
(2-digit SIC)
6
Set 2 A
Manufacturing
industries
(2-digit SIC)
I
I
I
I
I
I
I
I
I
I
I
I
I
I
- - - -- -
Set 1 A 0,
Manufacturing
industries
(2 -digit SIC)
Set 4 A
Manufacturing
industries
(2-digit SIC)
- - - - --
(Two-digit SIC detail)
o
= simulation options
FIGURE 1. 1
BLOCK FLOW CHART OF THE
PROGRAM RMS
Set 5 A
Manufacturing
industries
Set 6 B
Income deter-
mination
Set 7 B
Income deter-
mination
Set lOA
Manufac tur ing
indus tries
Set 7 C
Other industrie
Set 8 D
Labor market
-------
Ta:>.e ..1.
Rela :ions between Simuation 0:> :ion anc
~quation B .oc ,s in ::>rogram R vrs
-J
Simulation 1 2 3 4 5 6 7 8
Option of Industry Benefit Local Electricity Electricity Cont rol Fuel price Feedback
. Program Control by Taxa- price in - price in- Co s t by change to effects
EquatIon RMS Cost AQCR tion creas e to crease to 2-digit 2-digit
Blocks in
Program RMS indus try res identials SIC SIC
Block No. 1 A
Block No. 2 A
Block No. 3 E
Block No.4 A
Block No. 5 A A
Block No. ,6 B B B B B
Block No. 7 C C C C C C C C
Block No. 8 D D D D D D D D
Block No. 9 A A A
Block No. 1'0 B B B
Module:
A
B
C
D
E
Manufacturing Industry
Income and Regional Economy
Other Industry
Labor Market
Fuel Demand Submodel
-------
--:------,'-~~-;:-~.o.:..-',;---;-....~'-~" ...', .,',
a 42- sector 1-0 model has been used (see Appendix D. Volume I of the
May 15 Rep'ort).
The program then provide s a new 1-0 coefficient
matrix and the corresponding Leontief inverse matrix with new sector-
ing.
Finally, the program takes changes in final demand, eithe-r change
by sector or change by average propensity to consume by sector and pro-
duces the change in gross domestic product by sector.
Thus, the change s in the final demand, both as a re sult of costs
and benefits of air pollution control, will cause the change s in the
gross 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
transfers the value of shipment into value added by sector and then
distributes the feedback effect to each AQCR with the regional market
share matrix.
The output from Program FEE 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 pro-gram
with a 2100 statement program which includes nine subroutines to solve
both linear and non-linea r simultaneous equation sys tem. For the non-
linea r s ys tem, an iterative algorithm- -an extended Newton I s method>:'
has been introduced.
Since the program has been designed for a broad range of
policy 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 prog ram prepai'ation is given in
the following sections.
>.'( See M. K. Evans and L. R. Klein, The Wharton Econometric
Forecas ting Model, Unive rs ity of Penns ylvania, 1967, Cha pte r IV.
>:0:' 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.
Program
Cards j
FBE Deck
Benefit
Card
One Set for
Each Year
Requested
One Set for Each Run
Constant
Modifiers
Years para-
meter Card
Prog. RMS /
(Reg. Mode
Simulation)
Program
Deck
II Job
Job control
lang. (JCL)
10
-------
2.3 Data Input Tape
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 r e
are two ways of reading through the tape, depending on whether .the model
is to be used as before, for the 1967 simulation, or if the updated 1970-
1980 time series version is decided instead.
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
Cols. 5-10
.IS, two-digit SIC code
Cols. 11-20
FI0.4, labor share from production
function by industry
Cols. 21-30
Cols. 31-40
Cols. 41-50
F 1 0.4 capital share by industry
F 1 0.4, depreciation rate by industry
F 1 0.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 numbe r (blank for the 92nd
card)
19A4, name of AQ~R as s ociated with the code
in cols. 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 2 014, relative codes for AQCRs 61-80':<
Card 5 2014, relative code s for AQCRs 81-100':<
4.
Fuel Submode 1 Coefficients: 58 ca rds
(a) Share coefficients from fuel equation by two -digit SIC
(19 cards)
Cols. 9 -1 0
Cols. 11-20
Cols. 21-30
Cols. 31-40
Cols. 41-50
Co1s. 51-60
12, code of the two-digit SIC for which this
data applies
FI0.5, coefficient of coal
F 10. 5, coefficie nt of coke
F 1 O. 5, coefficient of oil
FlO.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 coefficients 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 applie s
FIO.5, constant determined in the invest-
ment equation
FIO.5, constant multipJier from the fuel
equation
FIO.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
FI0.6, coefficient
AQCR Economic Data: 455 cards. five cards for each AQCR
for which data was available
First Card:
Cols. 1-10
Cols. 11-20
Cols. 74-76
SecoJ;ld Card
Cols. 1-10
Cols. 11-20
F 1 O. I, labor force in each AQCR
F 10. 3, unemployment rate in each AQCR
13, AQCR code
FIO. I, 1967 employment in manufacturing
industries in each AQCR
FIO. I, 1967 wage bill for manufacturing
industrie s
13
-------
Gols. 21-30
Gols. 41-50
Gols. 74-76
Third Card:
Gols. 1-10
Gols. 11-20
Gols. 21-30
Cols. 74-76
Fourth Card:
Gols. 1- 10
Cols. 74-76
Fifth Card:
Gols. 1-10
Gols. 11-20
Gols. 21-30
Gols. 74-76
FIO.l, 1967 value-added in manufacturing
industries
FIO.l, 1967 investment in manufacturing
industries
13, AQCR code
FIO.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
F10.1, 1963 value-added for manufacturing
industrie s
13, AQCR code
14
-------
'-
6.
Electric Consumption (10 million KWH) by AQCR: 91 cards
Cols. 1- 10
FlO.O, total electric consumption in each
AQCR
Cols. 11-20
FlO.O,. residential electric power consump-
tions
Cols. 21-30
FlO.O, electric power consumed by manu-
facturing industries
Cols. 31-40
F 1 O. 0, electric power consumed in other
industries
Cols. 74-76
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
Cols. 1-10 FI0.l, number of employees
Cols. 11-20 FI0. 1,. wage bill
Cols. 21-30 F 1 O. 1, value -added
Cols. 74-76 13, AQCR code
Cols. 77-78 12, two-digit SIC code
A 119'sll card follows the 1963 data.
(b) Data for 1967
Cols. 1-10 FlO.l, number of employees
Cols. 11-20 F 1 O. 1, wage bill
Cols. 21-30 F 1 O. 1, value - added
Cols. 41-50 F 10. 1, investment
15
-------
Co1s. 74-76
Co1s. 77-78
13, AQCR code
12, two-digit SIC code
A 19's" card follows the 1967 data.
Detailed Fuel Information by Two-Digit SIC by AQCR--
There are two cards for each AQCR/two-digit SIC combina-
tion (2, 058 cards)
8.
Card 1
Co1s. 1-10
Co1s. 21-30
Cols. 31-40
Cols. 41-50
Co1s. 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
FlO.O, total cost of fuel consumption (all
type s)
F 1 O. 0, quantity of coal
F 1 O. 0, cost of coal
F 1 O. 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
FlO.O, quantity of oil
FlO.O, cost of oil
F 1 O. 0, quantity of natural gas
F 10.0, cost of natural gas
FlO.l, quantity of electricity
FIO.O, cost of electricity
16
-------
Cols. 74-76
13, AQCR code
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
applies
7F10.4, feedback/distribution to SICs 8-14
13, code of the AQCR for which this data
applies
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 0.3, feedback distr ibution
17
-------
Gr owth rate s for 1970: 104 cards, lwo 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
10.
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
11.
Growth rate s for 1975:
104 cards, same format as above.
12.
Growth rate s for 1980:
104 cards, same format as above.
13.
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
F 1 O. 0, Total per sonal income for 1970
13, AQCR code
F 1 O. 0, total per sonal income for 1975
13, AQCR code
FI0.0, 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
.OBE growth factors were available
Cols. 11-20
IlO, AQCR code
Cols. 31-40
--
F10.2, total personal 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 Input Parameter Cards
1.
The Years Parameter Card, for the time series run (with
constant modified switch)
Gols. 1-4 Total number of years included
Cols. 6-16 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.
Col. 80
1 constant modifiers added to the input
card deck afte l' this
Blank - no constant modifiers; next cards
are for testing convergence of city codes
2.
The Constant Modifiers are calculated from the difference
between the actual and estimated values without control in
1967 and are used to adjust the constant terms in the simul-
taneous systems in the time series convergence analysis.
Cols. 1-10
F 1 O. 2, regional adjustment constant for
the government expense equation
Gols. 11-20
F 1 0.2, regional adjustment constant for
the consumption equation
Gols. 21-30
F 10. 2, regional adjustment constant for the
wage equation, not currently in use
Cols. 31-40
F 1 O. 2, regional adjustment constant for
the labor force equation, not currently in use
Cols. 77-78
12, AQCR number
There are 91 cards in this set, one for each AQCR.
20
-------
3.
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.
4.
The Benefit Card
Cols. 1-10
F10.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)
Cols. 11-20
FIO.2, estimated national benefit less
compact of cost due to government assis-
tance. If there is not government assis-
tance, this number should equal the one
in the previous ten columns (in millions)
Cols. 21-30
F 1 0.2, normalizer: if benefit = $10 billion,
it equals 2,000. If benefit = $15 billion, it
equals 3, 000.
Cols. 31-40
FIO.2, percentage reduction in manufac-
turing costs resulting from government as'-
sistance
5.
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.
21
-------
6.
Input Parameter Card
A 11 of the following cards are nece s sary for each year
requested.
6. a. 1
Card Column
Variable Name
Des cription
1-4(14 )
INPT (1)
Information for aggregate
manufacturing industries
des ired.
(O or 1) - control input
data will be supplied.
to) - No
(1) - Yes
5-8(14)
NOPE
Number of simulation
options to be purs ued
for agg regate manufac-
turing industries (0, 1,
2, 3, 4, 5).
9-12
+ (NOPE - 1)
* 4 [NOPE ~,
14]
LPARA(I),
1= 1, NOPE
Type of simulation options
to be pursued in this run.
Specify "NOPE I( options.
o - none
1 - Industry control cost
2 - Benefit by AQCR
3 - Change in local taxa-
tion
4 - Electricity-price in-
creas e to indus try
. 0
5 - Electricity-price in-
crease to residential
us e rs
N through
N + 3
MCNT
If INPT(l) = 1, MCNT
specifies the number of
cities for which cost
changes are given or if
INPT( 1) = 0, MCNT °is
specified as 1.
22
-------
Card Column
(N + 3) through
(N-+ 7)
[14]
(n + 8) through
(N + 11)
[14]
Variable Name
INPT (2)
NOPEl
(N + 12) through LPNT
(N + 1 5 )
[14]
(Ii + 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
information to be included
(0 or 2) and simulated
o - not wanted
2 - des ired
Number of simulation op-
tions to be pursued for
detailed two -dig it SIC
data (O~ 1, 2, 3)
= 1
Types of detailed two-
digit SIC simulation op-
tions to be purs ued,
specify I'NOPE 1" options.
o - none
1 - control cos t by 2 dig it
SIC. This is s imu-
lation option number
6.
2 - impact of fuel price
'cha ng e
3 - national feedback/
dis tribution affect.
This indicates simu-
lation option numbe r
8 is to be pursued.
The number of cities fo r
which this control card
specification applies.
-------
~----
6.a.2
Card Column
L through
(L+3 )
[14]
(L+4) through
(L+7)
~4]
(L+8) through
(L+ 11)
[14]
(L+12)
(L+ 15 )
[14]
(L+ 16)
(L+19)
14
(L+20)
(L+23 )
[14]
Variable Name
1PRN1
1PR N2
1RPN3
1SBUG
1SBUG2
1CTAB
24
Description
The next five variable s
control the printo_\lt
options.
=0 312(A) summary table
requested but delete
printout of options 1,
4 & 5.
=1 Print results of all
strategies
=0 Do not print 312(A)
table (i. e., summary
table of 1 st five strate-
gies
=0 overrides all preceed-
ing reque sts and prints
out only the net sum-
mary tables for all op-
tions requested
= 1 Turns on debug print-
out patche s
=0 Skips option
=1 Detailed data tables
per year are given
=0 Skips option
=1 Convergence tables for
simultaneous systems
design
=0 Skips option
-------
Card Column
Variable Name
Description
(L+24 )
(L+27)
[14]
IGOU
= 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
Des cr iption
4
I
.Aggregate control cost for manufac-
turing, etc.
8
2
Two simulation options are to be pur-
.sued in this strategy simulation for
aggregate part of the model system.
12
4
Option 4 is chosen.
16
5
Option 5 is also chosen.
19-20
.1
Control input data for 91.AQCR'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 purs ued.
39-40 .
91
Two -digit SIC data available for 91
AQCR's.
-------
AN EXAMPLE OF PARAMETER CARD SPECIFICATION
Card Column Punch Description
41-48 B LA NK The IPRNSl, IPRNZ' printout options
are bypas sed
49-52 0 Final summary table requested
53-56 BLANK 1',';) debugging switches needed
57-60 B LA NK Avoide update tables
61-64 BLANK Delete convergence tables
N 65-68 1 Impact of cost with government
-...] assistance requested
-------
6.b.
Title Cards. On the next two cards, specify the heading
for the final summary table. Example:
Card 1
Cols. 1-80
Five Year Extended Implementation
Card 2
Cols. 1-80
Without Government Assistance
6. 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
6. d.
One Blank Card Necessary Here:
The initial parameter card required for each year makes it pos-
sible to reque st 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 reque sted 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 user
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 five) 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-dig:Ct 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
FIO.5, investment cost - stationary combus-
tion
Cols. 11-20
Cols. 21-30
Cols. 71-73
Cols. 74-76
FIO.5, annual cost - stationary combustion
FIO.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 FIO.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
Co Is. 31-40 FI0.4, investment cost for SIC 29
Cols. 41-50 FIO.4, investment cost for SIC 32
30
-------
Cols. 51-60
FI0.4, investment for SIC 33
Cols. 61-70
FI0.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 number six)
Cols. 1-10
Cols. 11- 20
Cols. 21-30
Cols. 31-40
Cols. 41-50
Cols. 51-60
Cols. 61-70
FI0.4, annual cost for SIC 20
FIO.4, annual cost for SIC 26
F 1 0.4, annual cost for SIC 28
FI0.4, annual cost for SIC 29
FIO.4, annual cost for SIC 32
FI0.4, annual cost for SIC 33
FI0.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 (IGOU =0 or 1), 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 Des cription of the Program Deck
The prog ram RMS cons ists of a MAINLINE and nine subroutines,
PRIN1, VALUE, VALUE1, VALUE2, VALUE3, VALUE4, SIMJ)"L,
GRADN and VERT, a total of 1700 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
cards.
The program parameters are then read from the input para-
meter card(s).
The paramete rs specify to the prog ram how many
and which of the various simulation options have been chosen by the
user.
The program now reads frorn 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,
VALUE2, 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 equations needed to implement the desired simulation option have
been solved.
The program now calls subroutine PRINl which prints
out the original values, the values after iteration and the percent
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 showing the total effect of the combined results of the simu-
lation options and the program terminates.
If all specified options
have not been pursued, the program 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
-------
MACRO FLOWCHART
OF RMS
Read stan-
dard input
data
eadpara-
meter input
cards
No
Yes
ead air pollu-
tion control cost
data for AQCR's
re uested
8
G)
Yes
Go to (2,3,4)
based on in-
put param-
eters
8
3
Cost affects
profit direct-
ly or through
increase in
ost of elec.
Cost affects
income
Cost of con-
trol affects
local taxation
34
-------
11
No
No
Go to 9, 10,11
based on inpu
par ameter s
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 industrie s
specified above
Solve system of
equations using
new costs
-------
Print out
detailed SIC
information
11
36
-------
olve system of
equations for new
values from
changed costs
Print out
re suIts
7
Print out net
results for all
simulation
options
37
-------
2.5 Program RMS (Regional Model Simulation)
38
-------
F ITRTR A '
0001
0002
0003
000'+
0005
0006
0007
0008
0009
0010
0011
0012
0013
001'+
0015
0016
0017
0018
0019
002U
0021
0022
0023
002'+
0025
0026
. .
G~EL
(
(
(
(
(
(
(
(
(
(
(
(
, )
MA IJ
DA"E = 72000
10/1.j1.j/0
[[[
. .
. A SIMULATION OF ECONOMIC EFFECTS OF AIR POLLUTION
ON AIR QUALITY CONTROL REGIONS
PROGRAr1 RMs
PREPARED BY (ONsAD RESFAR(H CORPORATION
NOVEMBER 1971
.
.
.
.
.
(ONTROL .
.
[[[
REAL KJNDX,LJNDX
DIMENSION A(20),
DIMENSION ALPHA(20),
DIMENSION BEN(8)
DIMENSION C(92),
DIMENSION COFI(91),
DIMENSION CTOT67(91,20),
DIMENSION D2(91),
DIMENSION DNBARN(92),
DIMENSION EOLDHI9I,20)
DIMENSION EOLD(9I,20,S),
DIMENSION GD(92),
DIMENSION IAQCR(92,20),
DIMENSION INPT(2),
DIMENSION KEE~IS),
DIMENSION LPARA(IO),
DIMENSION PDN(91) f
DIMENSION PICONIZO),
DIMENSION PNEWD(91 ,5),
DIMENSION PXD(91),
DIMENSION PZD(91),
DIMENSION QBARN(92),
DIMENSION QCN(92),
DIMENSION QMN(92),
DIMENSION QTN(92),
DIMENSION QGIVEN(7)
AB(91),
8(91,20),
(D(92) .
(OsTI91,5).
0(20) .
DELTC(91),
DNR(91),
G(92).
GN(92) .
IND(20) ,
J/(EEPI20) ,
LABEL! 10,20),
LPNT21),
PI(92,20),
PID(92),
POLD(91 ,20,5),
P X I D 1 91) ,
QBAR 1 92),
QC(92),
QM(92),
QT(92),
RNETI9I) ,
ALlT(20)
BETA(20)
CN(92)
(TOTD7(91)
D I 191 )
DNBAR(92)
EDI91,5)
GAMAIZO,S)
HOLDIS)
IND(II92,20)
KCITYII00)
LKEEPI3,5)
MCITYI91>
PICOEF(20)
PINI9Z)
Pp1D(91)
PXKD(91)
QAARD(92)
QCDI9Z)
QMD(92)
QTD(92)
sIJMN(92,20)
:)AG ~ JOO
.
.
~'1s00001
~'S00002
RMS00003
RMSOOOO't
RMsOOOOS
RMs00006
RMs00007
RMS00008
RMs00009
RMSOOOI0
RMSOOO 11
RMS00012
RMS00013
RMSOOOll.j
RMSOOOIS
RMS00016
RMS00017
RMS00018
RMS00019
RMS00020
RMS00021
RMS00022
RMS00023
RMS0002'i
RMS0002S
RMS00026
-------
---- ---- ----- -
------------ ----
FORTRAN I V G LEVEL ZO MA N DATE = 7Z006 10/QQ/08 ) AG ~ 000
C RMS00039
0027 DIMENSION SUNND(92), SUMNN(92), T(92) RMSOOO'O
0028 DIMENSION TO(9Z), TN(92), TR(9I) RMSOOO'+l
0029 DIMENSION U (92). UD(92), UN(92) RMSOOO'+Z
0030 DIMENSION VK(ZO). W(91 ,7.0), WB(91,ZO) RMSOOO'+3
0031 DIMENSION WBD(91), WD(91), WN(92) RMSOOO'+'+
003Z DII-IE"NSION WT(91 ,20), XD(92), XON(92) fH1S000'+S
0033 DIMENSION XD1(91), XI(92,20), XID(92) RMSOOO'+6
003'+ DIMENSION XIN(9Z), XK(92,20), XKD(92) RMSOOO'+7
0035 DIMENSION XKN(92), XKT(91,20), XL(92) RMSOOO'+8
0036 DIMENSION XLD(92), XLN(92), XND(9Z) RMSOOO'+9
0037 DII-IFNSION XNBAR(92), XNN(92), XNT(9Z) RMSOOOSO
OU38 DIr~ENSION XX(92,20), XX63(91 ,20), Y R ( 9 I) RMSOOOSI
OU39 DIMENSION YY(92), YYD(92), YYN(92) RMS00052
00'+0 DIMENSION YGROW(38) RMSOOOS3
00'+1 DIMENSION ZERO(16S6), ZZ(91,20) RMSOOOS'i
e RMSOOOSS
e RMS00056
OO'+Z DII-IENSION PERPIN(92). PERTN(92), PERSUM(92) RMS00057
00'+3 DIMF.:NSION PERXIN(92). PERCN(92), PERQBR(92) RMS00058
00'+'+ DIMENSION PERXKN(92), PERONB (92), PERQTN(92) RI'1S000S9
00'+5 DIMENSION PERXDN(92). PERXNN(9Z), PERQCN(92) RMS00060
00,+6 DIMENSION PERYYN(92), PERXLN(92), PERQMN(92) RMS00061
00,+7 DIMENSION PERGN(9i), PERUN(92) RMS00062
e RMS00063
00,+8 DIMENSION PER-RC ( 156'i) RMS0006'i
C RMS0006S
00'i9 DIMENSION PMED(ZO,S) ,KGOR(ZO,S) RMS00066
e RMS00067
OOSO DIMENSION XLH(92) ,UH(92) ,SUMNH(92,20) ,WBH(9Z'20) ,XXH(92,20) RMS00068
0051 DIMENSION XIH(92,20) ,YYH(92) ,TH(92) ,GH(92) ,CH(92) ,XNTH(92) RMS00069
00S2 DII-IENSION XNBAH(92) ,QTH(92) ,QMH(92) ,QCH(92) ,QBARH(92) ,XINC(60.3) RMS00070
00S3 DIMENSION GRFAC(100,20,3) ,IYEAR(Il) ,ICHO(91),XXT(92,20) RMS00071
e RMS00072
e RMS00073
OOS'+ DIr1ENSION SCel (3) ,SCC2(S) RMS0007'i
e RMSOUQ7S
e RMS00076
-------
---- -------- ----
FO~"~AN V G _::VEL
)055
c
c
0056
0057
0058
C
0059
C
0060
C
0061
C
0062
0063
C
006tf
C
0065
0066
C
0067
C
0068
C
0069
0070
0071
C
C
C
C
20
MAIN
)A'E = 72006
:)AGE 000
1011 't/0
~ISOO)77
~'S00078
RMS00079
RMS00080
RMS00081
RMS00082
RMSOOOB3
RMSOOOBLf
RMSOOQ85
RMS00086
RMS00087
RMS00088
RMS00089
RMS00090
RMS00091
RMS00092
RMSOU093
RMS0009Lf
RMS00095
RMS00096
RMS00097
RMS00098
RMS00099
RMSOOI00
RMSOOI0l
RMSOOI02
RMSOOI03
RMSOOI0'+
RMSOOI05
RMSOOI06
RMSOOI07
COMMON IPRITTI TITLE,BENGOV RMSOOI0B
COMMON/NPRINI LABEL,MCITY,KCITY, LSTRAT,INPT,IKSW,IHELp,IFEEORMSOOI09
COMMON I IPRIN IINTT,LPNTT ,JYEAR,IAQCR,IJCIT,MIN,IA312 RMSOOII0
RMSOOlll
RMS00112
RMSOOl13
RMSOO 11 Lf
DIMENSION :'lDGE 9,'+)
DIMENSION HHI(7),HHA(7)
DIMENSION ISPIN(71
DIMENSION TITLE(20,21
COMMON I ZEROT I PIN,XIN,XKN,XON,YYN,GN,TN,CN.ONRARN,XNN,XLN,WN,
1 UN,SUMNN,QBARN,QMN,QCN,QTN
COMMON IBLOXKll
I
2
PERPIN,PERTN,PERSUM.PERXIN,pERCN,PERQBR,
PERXKN,PERDNB,PERQTN,PERXDN,PERXNN,PERQCN,
PERYYN,PERXLN,PERQMN,PERGN,PERUN
COMMON IFATI F( 10001 ,LSTORE,KPIST,ARAB( 10) ,LCHECK
COMMON ISENSI ICTAB,NBEFOR,NEWCIT
COMMON ICATI Z(20,201
COMMON ISIMI Ml,M2,MCASES,Y( 10001 ,X( 1000) ,MARK,LVAR,MKEEP
COMMON I BUG I ISBUG,LREAO,LWRIT
COMMON IBUGGI IPRN1,{PRN2,IPRN3
COMMON I CPRIN I PXD, PXID, PXKD, PPIO, CTOTD7' pDN, CTOTb7
I ,EOLD,ED,POL~,PNEWD
COMMON IAPRIN/PI,PtD,XI,XID,XK,XKD,XX,XD,C,CD,YY,YYD,SUMN,SUMNO,
1 XNBAR,DNBAR,XNT,XND,U,UD,XL,XLD.G,GD,T,TD,QT
2 ,QTD,QC,QCD,QM,QMD,Q8AR,QBARD
-------
FOR"RM V G LEVEL 20 MAIN OAT:: = 72)06
0)72 EQUIVALENCE I ZERO I 1 I ,P I NIl I I
C
0073 EQUIVALENCE IPERRCll I,PERPIN(I II
C
0071f EQUIVALENCE (I NDC I ( 1 I , XX ( 1 ) I
C
C
C
1(/11/08
:JAG:: [(0'
C
C
C
BEWARE ... SIC AT A TIME WILL BE DIFF.
RMSOO . S
RMSOO .6
RMS00117
RMS00118
RMSOOl19
RMS00120
RMS00121
RMS00122
RMS00123
RMS00121f
RMS0012S
RMS00126
RMS00127
RMS00128
RMS00129
RMS00130
RMS00131
RMS00132
RMS00133
RMS00131f
RMS0013S
RMS00136
RMSOU137
RMS00138
RMSOU139
RMSOOIIfO
RMSOOIIfI
RMSOOllf2
RMS001"'3
RMS0011f'i
RMSOOllf5
RMSOOllf6
RMS001".7
RMSOOllf8
RMSOOllf9
RMS00150
RMS00151
RMS00152
0075
0076
0077
DATA LKEEP/ 3,l,I,'i,3,2,5,"',3,O,5,,,.,O,O,51
DATA SCC1/' Y',' G', 'C'I
DATA SCC2/' N-T',' L',' W',' U','SUMN'I
C
C
0078
DATA ISPIN/20,26,28.29,32,33,,,.0/
C
0079
0080
0081
DATA JKEEP/20.2/
DATA PMED,KGOR/I00.0.0,100.01
DATA FUDGE/36'i-0.01
C
0082
0083
008'i
0085
0086
0087
0088
ADJUST = 2000.
NOX = -I
LLLLL = 0
LLKKK = 0
LSTRAT = 1
LPNTT = 0
INTT = 0
0089
0090
0091
0092
0093
009'i
0095
0096
0097
00 106 I = Itl656
106 ZERO(I) = 0.0
DO 63'i0 1=1 tl56'"
63'i0 PERRCII) = 1.0
ISBUG = 1
LREAD = 5
LWRIT = 6
AlA = 0./897.
AAA = 960./111f5.
-------
FOR"RAN V G LEVEL 20 ~ A : I
C
0)98 IIOCS = 1
0099 IIOCS = 0
C
0100 ISBUG2 = 1
0101 ISHUG2 = 0
C
0102 I A 3 I 2 = 0
0103 IA312 = 1
C
010~ LTAPE = 9
0105 IDYEAR = 1967
0106 JYEAR = -2
0107 IDYEAR = 1971
0108 IFEED = 0
0109 q05 FORMAT(20A~)
01 10 ~10 fORMAT(20I~J
01 1 I ~qO FORMAT( IOX.5FI0.QJ
0112 DO ~OQ I = 1.9 1
01 13 DO. ~O~ J = 1,20
01 1 q XXCI,JJ = 0.0
0115 XXHCI.JJ = 0.0
01 16 DO 61 K=I.5
01 1 7 6 1 EOLD(I.J.KJ = 0.0
01 18 ~OQ CONTINUE
C
C
C GENERAL INPUT VARIABLES
C
C
01 19
0120
0121
0122
0123
012Q
0125
0126
)ATE = 72 )06
10/~~/08
PAG ~ )00
(IAQCRC I ,JJ .J=I.Z0J
RMS00153
RtiS001S£f
RMS001S5
RMS001S6
RMS001S7
RMS001S8
RMS001S9
RMS00160
RMS00161
RMS00162
RMS00163
RMSOOI6~
RMS0016S
RMS00166
RMS00167
RMS00168
RMS00169
RMSOOl70
RMS00171
RMS00172
RMS00173
RMS0017Q
RMS0017S
RMS00176
RMS00177
RMS00178
RMS00179
RMS00180
RMS00181
RMS00182
RMS00183
RMS0018Q
RMS0018S
RMS00186
RMS00187
RMS00188
RMSOo189
RMS00190
DO ~30 I = I.Z0
READ CLTAPE.qqOJ ALPHACI),BETACI),O(I) ,PICOEF(I),VKC1)
If (ISBUGJ Q30,Q30.qZO
~20 WRITE(LwRtT.~~O) ALPHACIJ.SETACI),D(I).PICOEF(IJ.VK(I)
Q30 CONTINUE
DO ~60 I = 1.92
READ(LTAPE.~61 J
QS8 fORMAT C IH1)
-------
F mITR A N--rv G L"£'I E L
0127
0128
0129
0130
0131
0132
0133
013&f
0135
0136
0137
0138
0139
01&f0
C
o 1 &I 1
01&f2
01&f3
01&1&f
01&15
01&16
01&f7
01&f8
01"19
0150
o 1 5 1
0152
0153
0151.f
0155
0156
0157
20
OAT:: = 72006
10/1.f1.f/08
'AGE 000
M A . ~
&159 FORMAT(IHl,I/IH
lUOY',/)
1.f60 CONTINUE
&f61 FORMAT (1&f,19A&f)
WRITE (lWRIT,1.f59)
KITTY = 0
00 &fbLf I = 1,91
WRITE (lWRIT,&f62) ( IAQCR(I,J),J=I,20)
KITTY = KITTY + 1
IF (KITTY-31) 1.f61.f,&f63,Lf63
1.f62 FORMAT(tH ,23X,13,2X,19A&f)
&f63 KITTY = 0
WRITE IlWRIT,1.f5B}
&f6&f CONTINUE
WRITE (lWRIT,&f58)
,'AQCR S :oOR WHICH OA.A WAS AVA _A3lE :OOR THIS
READIlTAPE,&fl0) (KCITYII) ,1=1,100)
IF(LLKKK.EQ.l) GO TO 120
LLKKK = 1
READIlTAPE,1.f90) I(GAMA(IKIT,JKIT), JKIT = 1,5), IKIT = 1,19)
&f90 FORMATltOX,5FI0.5)
00 Lf92 KIT = 1,20
REAO(lTAPE,&f91) PICON(KIT), AIKIT), AlIT(KIT)
1.f91 FORMAT( tOX,3FI0.5)
&192 Corn I NUE
DO 799 KIT = 1,19
799 READ(LTAPE,798) AB(KIT)
7 9 8 FOR r-1 A T ( lOX, F 1 0 . 6 )
IF( IOYEAR.EQ.1967) Go TO 9&10&f
C
C
C
C
1970-1980 INPUT SEGMENT
DO &f50 I = 1,91
REAO(lTAPE,&I&f5) XlHI I) ,UHf I} ,SUMNH( I ,20) ,WBH( I ,20) ,XXH(I ,20),
1 XIH( 1,20) ,YYHI I) ,THI I) ,GH( I) ,CHI J)
IF (ISBUG) &f&l7,1.f1.f7,I.fLf&f
LfI.fLf WRITEILWRIT,I.fLf6)XLH(I) ,UHII) ,SUMNH(I,20),WBH(I,20),XXHII,20),
STRMS0019
RMS00192
RMS00193
RMSOOI9&f
RMS00195
RMS00196
RMS00197
RMS00198
RMS00199
RMS00200
RMS00201
RMS00202
RMS00203
RMS00201.f
RMS0020S
RMS00206
RMS00207
RMS00208
RMS00209
RMS00210
RMS00211
RMS00212
RMS00213
RMS0021&f
RMS0021S
RMS00216
RMS00217
RMS00218
RMS00219
RMS00220
RMS00221
RMS00222
RMS00223
RMS00221.f
RMS00225
RMS0022b
RMS00227
RMS00228
-------
FOR'RAN
o 58
0159
0160
016.1
0162
0163
0161i
0165
0166
0167
0160
0169
0170
0171
0172
0173
0171i
0175
0176
0177
0178
0179
0180
0181
0182
0183
0181i
0185
0186
0187
'V G LEVEL
20
DA ';: = 72006
M A . ~
10/~1.f/08
'AGE 000
1 XIH(I.2r),YYH(I).T~( ).G~( ),CH ~'1S)0229
. .
Ii Ii S FOR M A ., F 1 0 . , F 1 0 ." 3 , / . 3 :- 1 O. ,1 J X ,:' O. 1 . / , 3 :- 1 ( . 0 , / . :- 1 0 . 0 , / , 8 0 X ) R ' 5 0 0 2 3 0
1.f1.f6 FORMATIIH ,FI0.1,IX'F10.3.1X,Ii(Fln~1 ,IX) ,I.fIFIO.O.IX)) RMS00231
'1'17 XNTHII) = (1.0-UH(l}).XLH(I) RMS00232
XNRAH I I) = XNTH ( I ) -SUMNH I I ,20) RMS00233
'150 CONTINUE RMS0023/.f
00 /.f70 I = 1,91 RMS00235
READILTAPE,1.f66) QTH(I) ,QCHIII,QMH(I) ,QBARHII) RMS00236
WRITE(LWRIT,466) QTH(I},QCH(I),QMHII} ,Q8ARHtI} RMS00237
FORMATIIH ,2UIS) RMS00238
COIHINlJE RMS00239
FORMATI~F10.0) RMS0021.f0
REAOILTAPE,5991.f) XKIT,KIND RMS0021.fl
FORMAT (FI0.1,66X,12) RMS0021.f2
IFIKIND-99) 5993,5996,5996 RMS0021.f3
5996 READILTAPE,5997) XKIT.XKIT1,XKIT2,XKIT3,ICT,KIND RMS0021.f1.f
5997 FORMATI.1FI0d dOX,F10el ,23Xtl3d2) RMS0021.f5
IF IKIND-99) 5998,S999,S999 RMS0021.f6
5998 KIND = KIND-20+120/KINO) RMS0021.f7
1CITY = KCITYIICT) RMS0021.f8
IF IICITY.EQ.O) GO TO 5996 RMS0021.f9
'/1 B H I I C ! T Y , K ! NO) = X KIT 1 R M 5 0 0 2 5 0
SUMNHI!CITY,KINO) = XKIT RMS00251
XXHI!CITY,KIND) = XKIT2 RMSOU252
XIHI!CITY,K!ND) = XKIT3 RMS00253
GO TO 5996 RMS002511
C RMS00255
C[[[RMS00256
'165
1.f7(J
'166
5993
599'1
C
C
C
C
9'10'1
RMS00257
RMS00258
RMS00259
RMS00260
CONTINUE RMS00261
DO 9'150 1=1,91 . RMS00262
REAO(LTAPE,91.fI.fS) XLII),UII),SUMN(I.20) ,WBII,20I,XXII,20), RMS00263
IX! I I ,20) ,YY( I) ,TI I) ,G( I I ,C( I), XNTI I), WTI I ,20) ,XX63( I .20)RMS002611
IF I ISBUG) 9'1'17,9'1'17,9444 RMS00265
-------
FOH~AN
0188
0189
0190
0191
0192
0193
019'+
0195
0196
0197
0198
0199
0200
0201
0202
0203
0201f
0205
0206
0207
0208
0209
0210
021 I
0212
0213
021'1
0215
0216
0217
0218
0219
0220
V G LEVEL
20
MA N
DA' ~ = 72006
lXI(I.201.YYI.TI .G(II,C :1,
9lflf7 PI(I,201 = XXII,~0)-WB(I.20)
XK(I,20) = XXII,201/VKI201
WII,20) = WB(J.20)/SUMNII,201
WT(I.201 = WT(I,201/XNTII)
XNT(I) = 11.0-U(III.XLlII
XNBAR(J) = XNT(I)-SUMNII,20)
9'+50 CONTINUE
. ,
9lfqS FORMAT(FIO.l,FI0.3,I,FI0.1 ,FI0.l,FIO.I,10X,Fln.I,I,
13FI0.0,I, FI0.0,I,FIo.I,FIO.I,FI0.1,10X,FIO.11
9qlf6 FORMATI IH .FIU.l,FI0.3,I,FI0.I,FI0~1 ,FIO.l rl0X,FI0.l,l,
13FIO.0,1. FIO.O,I.FI0.I.FI0.I.FIO.1 ,IOX.FIO.I)
C
9'+65
9lf70
1f67
C
9993
9995
9996
9998
WT( 1,201 ,XX631 I ,201 ~MS00267
RMS00268
RMS00269
RMS00270
RMS00271
RMS00272
RMS00273
RMS00271f
RMS00275
RMS00276
RMS00277
RMS00278
RMS00279
RMS00280
RMS00281
RMS00282
RMS00283
RMS0028&f
RMS00285
RMS00286
RMSOU287
RMS00288
RMS00289
RMS00290
RMS00291
RMS00292
RMS00293
RMS0029Q
RMS00295
RMS00296
RMS00297
RMS00298
RMS00299
RMS00300
RMS00301
RMS00302
RMS00303
RMS0030Q
x ~ T (
I ,
Do 9Q70 1=1,91
READ(LTAPE, 1f661 QTI I I ,QCIII,QMIII,QBAR(II
IF(ISBUG) 9lf70,9Q70,9'i65
WRITE ILWRlT,'+67) QT(II ,QC(I) ,QMI II .QAARIII
COrH I NUE
FORMAT(IH ,lfFI0.0)
READ(LTAPE,999QI XKIT, XKITI, XKIT2,
IFIKIND - 991 9995,9996.9996
KIND = KIND - 20 + 120 I KIND)
ICITY = KCITY(ICTI
IF II CITy .EQ. 01 GO TO 9993
WTIICITy,KINO) = XKITI / XKIT
XX631ICITY,KINDI = XKIT2
XKTIICITY,KIND) = XKIT2 I VKIKIND)
GO TO 9993
READILTAPE,9997) XKIT, XK,ITl, XKIT2.
IFIKIND - 99) 9998,5999,5999
KIND = KIND - 20 + 120 / KINDI
ICITY = KCITYIICTI
IF II CITy .EQ. 01 GO TO 9996
WIICITY,KINDI = XKITI I XKIT
WBIICITy.KIND) = XKITI
SUMNIICITy,KIND) = XKIT
PIIICITy,KINDI = XKIT2 - XKITI
ICT, KIND
XKIT3, IcT, KIND
10/qlf/08
)AG: 0)0
-------
FORTRAN I V G - ~VEL
0221
0222
0223
022li
0225
0226
RMS003CS
RMS00306
RMS00307
RMS00308
RMS00309
RMS00310
RMS00311
RMS00312
RMS00313
RMS0031li
,KIND,IXII),YII),1=3,5)RMS00315
RMS00316
RMS00317
RMS00318
RMS00319
RMS00320
RMS00321
RMS00322
RMS00323
RMS0032li
RMS0032S
RMS00326
RMS00327
RMS00328
RMS00329
RMS00330
RMS00331
RMS00332
RMS00333
RMS00331.f
RMS00335
RMS00336
RMS00337
RMS00338
RMS00339
RMS003liO
RMS003'il
RMS003'i2
20
M A I II
DATE = 72006
= X(ITZ I VK K 'D)
= XKIT2
= XKIT3
XKIICI"y,K ND)
XXIICITy,KINO)
XI I ICITy,KIND)
GO TO 9991,
FORMATI3FI0.l,1.f3X,I3,IZ)
FORMATI3FI0.1,10X,FI0.1,23X,I3,I2)
999'i
9997
C
C...........
0227
0228
0229
0230
0231
0232
0233
023"+
OZ35
0236
0237
0238
0239
02,,+0
02"+1
02'12
02'i3
02"+1.f
02liS
02,,+6
02'17
02,,+8
02li9
0250
0251
0252
0253
025"+
02S5
C
5999
6002
6103
610"+
6000
6003
6006
2301
2302
CONTINUE
READILTAPE,6000) CT,IXII) ,Y(I),I=I,2) ,ICT
IF (KINO-99) 6002,6006,6006
KIND = KIND - 20 + 120 I KIND)
ICITY = KCITYIICT)
CTOT67IICITY,KINO) = CT
006003 1=1,5
EOLOIICITy,KIND,I) = XII)
IF IGAMAIKIND,I).NE.O.O) GO TO 6103
POLDIICITy,KINO,I) = 0.0
GO TO 6003
IF (X 1 I ) .y 1 I ) . NE. 0.0) GO TO 610'1
POLOIICITy,KIND,I) = 0.0
GO TO 6003
POLDIICITY,KINO,I) = YIII/XII)
FORMATIFI0.0,10X,"+FI0.0,13X,I3,I2,1,'FI0.0,Fln.1,Fld.O)
PMEDIKIND,I) = PMEDIKIND,II+POLD(ICITy,KINO,I)
KGO~IKIND,I) = KGORIKINO,II+l
CONTINUE
GO TO 5999
CONTINUE
DO 2302 JNK = 1,19
007.301 KGO = l,S
IFIKGORIJNK,KGO).EQ.O) GO TO 2301
PMEDIJNK,KGO) = PMEOIJNK,KGO)/KGOR(JNK,KGO)
C 0 ~JT I N U E
CONTINUE
DO 2350 ICK = 1,91
00 23"+5 JNK = 1,19
1( II /.f/08
PAG~ :(0'>
-------
,-- ..
I
FOR'RAN IV G LEVEL
20
MA I
0256
0257
0258
0259
0260
0261
0262
0263
026~
0265
0266
0267
0268
n269
0270
0271
0272
0273
027~
0275
0276
0277
0278
0279
0280
0281
0282
0283
028'i
0285
0286
0287
0288
0289
0290
0291
0292
KG = JKEEP(JrK
00 23~0 KNK = 1,5
KGO = LKEEP(KG,KNK)
IF (KGO) 23~5,23~5,2310
IF (POLDltCK,JNK,KGO» 2320,2315,2320
POLDIICK,JNK,KGO) = PMED(JNK,KGO)
IF IEOLO(ICK,JNK,KGO» 23'i0,2325,23~0
EOLDltCK,JNK,KGO) = .001
CONTINUE
CONTINUE
CONTINUE
PERUN(92) = .O'i
00 20 I I = 1,7
READ ILTAPE,202) K,IK
K = K - 20+(20/K)
201 JKEEP (K) = IK
202 FORMAT (213)
IfILLLLL.EQ.l) GO TO 38
LLLLL = 1
DO 5991 ICT = 1,91
I = ICT
PERUN(I) = .O~
IfIIDYEAR.EQ.1967) GO TO 5991
00 5990 MJ= 1 tl9
5 9 9 0 E 0 L D H ( I C T , M J) = E OLD ( I C T , 1-1 J , 5 )
EOLDHIICT,20) = QTH(ICT)
REAOILTAPE,S992) IB(I,J), J = 1,19)
5991 COfH I NUE
DO 298 I = 1,91
READ (LTAPE,299) RNETII)
298 CONTINUE
00 ~Dl 1=1,8
READ(LTAPE,'i05)ILABELII,J), J = 1,20)
tfOI CONTINUE
299 FORMATI20X,FIU.3)
38 CONTINUE
LPNTI = 0
2310
2315
2320
232!:>
23tfO
231f5
235U
c
OA"E = 72Jn6
lC/.1/08
)AGE : Olr
RMS003'i3
RMS0031 tf
RMS003tfS
RMS003tf6
RMS003tf7
RMS003tf8
RMSOU3tf9
RMS00350
RMS00351
RMS00352
RMS003S3
RMS003Stf
RMSOD3SS
RMS00356
RMS00357
RMS003S8
RMS00359
RMS00360
RMSOD361
RMS00362
RMS00363
RMS0036'i
RMS00365
RMS00366
RMS00367
RMS00368
RMS00369
RMS0037D
RMSDD371
RMSD0372
RMSD0373
RMSD037'i
RMSOD375
RMSD0376
RMSD0377
RMS00378
RMSD0379
RMS0038D
-------
FUKI~AN IV G LEVEL
C
C
0293
029£i
0295
0296
0297
0298
0299
0300
0301
0302
0303
030£i
0305
0306
0307
0308
0309
0310
0311
0312
0313
031£i
0315
0316
0317
7010
7015
lU2U
1050
1060
C
0318
0319
0320
0321
0322
0323
753£i
7533
C
C
C
C
£0
MAIN
DATE = 72006
10/4 '/08
;)AG~ (01
GROW. .
~AC .O~ CALCULATIONS FO~
970-1980 ONLy
IFCIDYEAR.EQ.1967) GO TO £i00
7009 READCLREAD,1070) INY,CIYEAReJRK) ,JRK=I ,INy)
WRITEIL~RIT,1070) INY,CIYEARIJRK) ,JRK=I,INY)
IMII = IYEARCINY) - 1968 - 110(5
1070 fORMATCl£i,IX,1215)
IF I1NY.EG.0) STOP
DO 1020 JRK = 1,5£i
READ CLTAPE,7010) CGRFACCJRK,KRKd) ,KRK=I ,20)
FORMAT Cl~F5.0,/,6F5.0)
FORMATC I~F5.3,/,6F5.3)
CONTINUE
DO 1025 JRK = 1,5£i
DO 1025 KRK = 1,20
GRFACIJRK,KRK,I) = GRFACeJRK,KRK,I)/100.
1025 CONTINUE
00 1030 JRK = 1 ,5~
READILTAPE,7015) IGRFACIJRK,KRK,2),KRK=I,20)
1030 CONTINUE
00 10~0 JRK = I,S~
READCLTAPE,7015) IGRFACIJRK ,KRK,)) ,KRK=1 ,20)
10'10 COtHINUE
DO 1060 JRK = 1 ,S~
READILTAPFd050) IXINCeJRK,KRK) ,KRK=1 ,3)
FO~MATII0X,FI0~0,/,10X,FI0.0,1,10X.FI0.0)
CONTINUE
DO 7S3~ J=I,91
READILREAD,7533) (FUDGE(J.I) , 1=1,'1)
IFIIS8UG.NE.I) GO TO 753'1
WRITEILWRIT,7533) (FUDGEeJ.I) , 1=1,'1)
CONTINUE
FORMAT(~FI0.2)
FOR AQCRS NOT INCLUDED IN OBE REPORT ICODE GREATER THAN S'I )
~'S0038
R'S00382
RMS00383
RMS0038'f
RMS00385
RMS00386
RMS00387
RMS00388
RMS00389
RMS00390
RMS00391
RMS00392
RMS00393
RMS0039'f
RMS00395
RMS00396
RMS00397
RMS00398
RMS00399
RMSOO'fOO
RMSOO'lOI
RMSOO'l02
RMSOO'103
RMSOO'lO'f
RMSOO'l05
RMSOO~06
RMSOO'f07
RMSOO~08
RMSOO~09
RMSOO'l10
RMSOO'lII
RMSOO'l12
RMSOO'f13
RMSOO'lI'1
RMSOO'fIS
RMSOO~16
RMSOO'll7
RMSOO'lIB
-------
F~AN~ G L~v~L'
C
C
C
C
C
032Lf
0325
0326
1051
131 U
C
C
0327
0328
0329
0330
0331
0332
0333
033Lf
0335
0336
1052
1053
105Lf
C
0337
0338
0339
03LfO
ZtJ
MA IJ
OAT::: = 72006
REAoILTAPE,1051) YGROW
FORMATI30X,FI0.2,LfOX)
FORMATIIH ,13,12F9.1,I5)
00 105Lf JRK = 1,20
00 105Lf KRK = 1,3
00 1053 KC = 55,100
LCITY = KCITY(KC)
IF(LCITY.EQ.OI GO TO
GRFACILCITY,JRK,KRK)
GO TO 1053
GRFACILCITY,JRK,KRK)
CONTINUE
CONTINUE
1052
= GRFACI32,JRK,KRKI
= 0.0
00 1055 IJ=I,38
WRITEILWRIT,1070) IJ
1055 YGROWIIJ) = YGROWIIJI1100.
REAoILREAD,lOUOI INC,(ICHOIJRKI,JRK=I,INC)
C
03Lfl
03Lf2
031.f3
7211
C
C
C
031.fLf
031.f5
03Lf6
03Lf7
031.f8
7225
C
IFIIA312.EQ.OI GO TO 7213
READ(LREAD,7211) B~NGOV,BENEF,THIS1,P
FORMATII.fF10.2)
...
OVERWRITES FEEDBACK MATRIX IF BENEFIT NOT EQ~ $10 BILLION
IF(BENGOV.EQ.I0000.1 GO TO 7213
00 7225 ICT = 1,91
READ(LREAD,59921IBIICT,JI,J=I,19)
VVRITE(LWRIT,59921 (B( ICT,J) ,J=l, 19)
CONTINUE
10/1. /08
:JAGE 001
~MS)OI.f 9
R~SOOI2:
RMSOOLf21
RMSOOl.f22
RMSOOLf23
RMSOOl.f21.f
RMSOOl.f25
RMSOOl.f26
RMSOOl.f27
RMSOOl.f28
RMSOOl.f29
RMSOOl.f30
RMSOOl.f31
RMSOOl.f32
RMSOOl.f33
RMSOOl.f31.f
RMSOOl.f35
RMSOOl.f36
RMSOOLf37
RMSOOl.f38
RMSOOl.f39
RMSOOLfLfO
RMSOOLfl.fl
RMSOOl.f1.f2
RMSOOl.f1.f3
RMSOOI.fI.fLf
RMSOUI.fLf5
RMSOOl.f1.f6
RMSOOl.f1.f7
RMSOOLfI.fB
..RMSOOI.fLf9
RMSOOl.f50
RMSOOLf51
RMSOOLfS2
RMSOOLf53
RMSOOLf51.f
RMSOOLf5S
RMSOOl.f56
-------
F ORTR A ~
03,+9
0350
0351
0352
0353
035'1
0355
0356
0357
0358
0359 '
0360
0361
0362
0363
036'1
0365
0366
0367
0368
0369
0370
0371
0372
0373
037'1
0375
\ G LEVEL
7213
C
C
C
C
C
C
20
)AE = 72006
10/II/Of
MA N
CONTINUE
ELIMINATES REQUESTS FOR INFO. ABOUT CITIES FOR WHICH NO DATA IS
AVAILABLE (EXAMPLE - NO DATA FOR CITY NO. 131
NIL = INC
DO 1063 1=1 tlNC
106'1 KIT = ICHOIII
lAB = KCITYIKIT>
IF(IABI 1063.1061.1063
1061 NIL = NIL - 1
DO 1062 NIPP = I.NIL
ICHOINIPPI = ICHO(NIpp + II
1062 CONTINUE
WRITEILWRIT,10651 KIT
GO TO 106'1
1063 CONTINUE
INC = NIL
FOr-
-------
Fa ~'.RAN
0376
0377
0378
0379
0380
0381
0382
0383
038'1
0385
0386
0387
0388
0389
0390
0391
0392
0393
039'1
0395
0396
0397
0398
0399
O'+OU
0'+01
0'+02
0'+03
0'+0'+
0'105
'V G LEVEL
(
(
(
(
6106
1 1 1 5
(
C
C
C
20
DATE = 72006
MAIN
~ERXKNll) = I. + IXKNII) / XKII.2)) . PERXKNII)
PERXDN(I) = II. + IXONII)/XXII,20») . PERXDNII)
PERYYNII) = 11. + (YYNII)/yYII») . PERYYNII)
PERGNII) = II. + IGNII)/GII») . PERGNII)
PERTNII) = 11.+ ITNII) / TII») . PERTNII)
PERCNI}) = 11. + I(NI!)/ (II») . PERCNII)
PERONBII) = II. + IONBARNI}) / XNBARI I») . PERDNBI})
PERQBRI}) = (1. + (QBARNII) / QBARII») . PERQBRI})
PERQ(NII) = (1. + IQ(NII) / Q(II»I . PERQCN(I)
PERQMN( I) = (1. + IQMNII) / QMII»I . PERQMNII)
PERSUM(I) = (1. + (SUMNN( I) / SUMNI!,20») . PERSUM(I)
DO 6106 1=1.1656
ZEROII) = 0.0
(ONTINUE
IF IJYKNT.EQ.INY) Go TO 7009
JYKNT = JYKNT + 1
JYEAR = IYEARIJYKNT)
lEX = JyEAR + 3
... TO GENERATE OCS SET TEST TO 8 - 1977 - NOW )978
...
IF(JYEAR.EQ.IMII) GO TO 95
KYy = JYEAR+1969
XX(92,20) = 0.0
PI192,20) = 0.0
XI(92,20) = 0.0
XK192,201 = 0.0
SUMNI92,20) = 0.0
XNBAR(92) = 0.0
YY(92) = 0.0
XNT(92) = 0.0
(192) = 0.0
XL(92) = 0.0
10/'+'1/08
PAG~ 001!
RMSOO,+9S
RMSOOl.f96
RMSOO'+97
RMSOO,+98
RMSOO'+99
RMSOOSOO
RMSOOSOI
RMS00502
RMS00503
RMS0050'+
RMSOOSOS
RMSOOS06
RMSOOS07
RMS00508
RMSOOS09
RMSOOSI0
RMSOOSll
RMSOOS12
RMS00513
RMSOOSll.f
RMSOOS15
RMS00516
RMSOOS17
RMSOOS18
RMS00519
RMSOOS20
RMSOOS21
RMS00522
RMSOOS23
RMSOOS2'1
RMsoaS2S
RMS00526
RMSOOS27
RMS00528
RMS00529
RMSOOS30
RMSOOS31
RMSOOS32
-------
FOR'1AN I V G - ~VEL
20
MA N
DA"E = 72006
1(/~~/08
('f06
0'f07
0'f08
0'f09
O'fl0
O'f 1 1
0'f12
Q'I921 = 0.0
QMI921 = 0.0
QBARI921 = 0.0
QCI921 = 0.0
GI921 = 0.0
TI921 = 0.0
GO TO (1120,1130,1130,1130,1130,1130,II'fO,II'fO,II'f0,11~0,11~01,
1 JYEAR
1120 KJNDX = 0
LJNDX = 0
GO TO 1150
1130 ,KJ~JDX = JYEAR-l
L J tW X = 0
GO TO 1150
11'f0 KJNDX = 5
LJNDX = JYEAR-6
1150 DO 1210 MPCITY= I,INC
KIT = ICHOIMPCITYI
LCITY = KCITY(KITI
DO 1200 JGE = 1,20
IFIXXHILCITY,JGEI.EQ.O.OIGO TO 1200
IFIGRFAC(LCITy,JGE,II.EQ.O.O) GO TO 1200
W(LCITY.JGE) = WBHILCITY,JGE)/SUMNH(LCITY,JGE)
116U GFACI = GRFAC(LCITY,JGE,I)
GFAC2 = GRFAC(LCITY,JGE,21
GfACD = GFACZ'- GFACI
XXILCITY,JGEI = XXHILCITY,JGEI.(GFAC1+KJNDX/5..GFACD +
1 (LJNDX/5..(GRFAC«(CITY,JGE,3)- GFACZ»)
VFAC = XXILCITY,JGE)/XXH(LCITY,JGlI
DO 1161 JS = I,IEX
1161 WILCITY.JGEI = W(LCITY,JGE) . 1.018
"BILCITy,JGEI = WBH(LCITY,JGEI.VFAC
IfISUMNH(LCITY,JGEI.EQ.O.O.OR.W(LCITY,JGE).EQ.O.o)
1 WRITE(LWRIT,17701 KIT,JGE
1770 fORMATltHI,'fI51
SUMNILCITY,JGE) = WB(LCITy,JGEI / W(LCITY,JGEI
PIILCITY,JGEI = XXILCITY,JGEI-WB(LCITy,JGE)
XIILCITy,JGE) = XIH(LCITY,JGEI.VFAC
O'f13
O'fl~
0~15
0~16
0~17
O~I8
0'+1 9
0~20
0~21
0~22
0~23
O'f2'+
0~2S
0'f26
0'f27
O~28
O'f29
O~30
O'f31
O'f32
0'f33
0'f3'f
0~35
0'f36
0~37
0"+38
0"+39
0,,+,,+0
:tAGE 00I~
~'S00533
RMS00531
RMS00535
RMS0053b
RMS00537
RMS00538
RMS00539
RMS005"+0
RMS005~1
RMSOOS"+2
RMS005'f3
RMS005'f'l
RMS005'1S
RMS005'1b
RMSOOS'I7
RMS005'18
RMS005'+9
RMS00550
RMS00551
RMS00552
RMSOU553
RMS0055'+
RMS00555
RMS005Sb
RMS00557
RMS00558
RMS00559
RMSOOSbO
RMSOOSbl
RMSOOS62
RMSOOS63
RMSOOSb~
RMSOOSb5
RMSOOS66
RMSOOSb7
RMS00568
RMS00569
RMSOU570
-------
20
DA":: = 72006
IC/jl/08
FOR .RAN I V G - ::VEL
O. . 1
0'iq2
O'i'i3
Oqq'i
O'iqS
C
C
Oqq6
0'i'i7
O'iQ8
0'i'i9
O'i50
0'i51
0'i52
OQ53
O'iS'i
0'i55
0'i56
0'iS7
0'i58
0'i59
0'i60
0'i61
0'i62
0'i63
0'i6'i
0'i65
0'i66
0'+67
0'i68
0'+69
0'i70
0'+71
0'+72
0'i73
MAl'
XKILCITY,JGE) = XXI -C.Y,JG::)/VKIJG::
EOLDILCITY,JGE,5) = EOLOHILCITY,JGE) - VFAC
1200 CONTINUE
IFILCITY.GT.S3) GO TO 120'i
IFIXINCILCITY,I).NE.O.O) GO TO 1205
YYILCITY) = 0.0
GO TO 1210
RATE = 0.0
IFIKIT.EQ.23) RATE = .0762
IFIKIT.EQ.32) RATE = .07S8
IFIKIT.EQ.33) RATE = .06'i5
YYILCITY) = YYHILCITy)
00 1203 J5=ltlEX
1203 YYILCITY) = YY(LCITY) - 11.0 + RATE)
GO TO 1207
120'+ YY(LCITY) =
DO 1206 LLJ
1206 YYILCITY) =
GO TO 1207
1205 CONTINUE
YY(LCITY) = XINC(LCITY,I)+KJNDX/5~(XINC(LCITY,2)-XINC(LCITY,I))+
1 LJNDX/S.(XINCILCITY,3)-XINCILCITY,2))
YY(LCITY) = YYILCITY) -1.1'i'i
1207 YFAC = yYILCITY) /YYH(LCITY)
XNBAR(LCITY) = XNBAH(LCITY)~(YY(LCITY)-XXILCITY,20))/
1 (YYHILCITY)~XXHILCITY,20))
T(LCITY) = TH(LCITY).YFAC
G(LCITY) = GH(LCIT~).TILCITY)/THILCITY)
C(LCITY) = CHILCITY).YFAC
XNTILCITY) = XNBAR(LCITY)+SUMNILCITY,20)
U(LCITYJ = .0'+
XL(LCITY) = XNTILCITY)/II.-U(LCITY))
QMILCITY) = QMH(LCITY).VFAC
Q~ARILCITY) = QBARH(LCITY)-(YY(LCiTY)-XX(LCITY,2n))/
1 IYYHILCITY)-XXHILCITY,20))
QCILCITY) = QCHILCITY).CILCITY)/CHILCITY)
QT(LCITY) = QCILCITY)+QBARILCITY)+QM(LCITY)
1210 CUNTINUF.
YYH(LCITYI
= ItlEX
(YGROVJ(LCITY-53)
+ 1.0)
- YY(LCITy)
)AGE 001~
~'SOOS7
RMSOOS72
RMSOUS73
RMSOOS7'i
RMS0057S
RMSOOS76
RMSOOS77
RMSOOS78
RMSOOS79
RMSOOS80
RMSOOS81
RMS00582
RMS00583
RMS0058'+
RMSOOS8S
RMSOOS86
RMS00587
RMSOOS88
RMS00589
RMSOOS90
RMSOOS91
RMS00592
RMSOOS93
RMSOOS9'i
RMSOOS9S
RMS00596
RMS00597
RMSOOS98
RMS00599
RMSOObOO
RMS00601
RMS00602
RMS00603
RMS0060'i
RMS0060S
RMS00606
RMS00607
RMSOOb08
-------
FOR . ~AN IV G _EVEL
DATE = 72(06
1014 'f/OB
20
MA IJ
O. 7'+
0'+75
212 FORMA"(F O.3,60X,I5t13
GO TO (1220tl230tl235tl235,235tl2J5tl2J5,1Zlo, 2'+0, Z'+Otl2/i0,
1 12/i0),JYEAR
1220 DO 1227 MPCITY= 1,INC
KIT = ICHOCMPCITY)
~RITE(6,1070) KIT
LCITY = KCITY(KIT)
DO 1225 JGE = 1,20
IFCXXHCLCITY,JGE).EQ.O.O) GO TO 1225
XXT(LCITY,JGE) = XXH(LCITy,JGE)+2./3~.(XXCLCITY,JGE)-XXH(LCITY,
1 JG E) )
1225 CONTINUE
1227 CONTINUE
GO TO 1260
123U KJNOX = 0
LJNDX = 0
GO TO 12S0
1235 KJNDX = JYEAR-2
LJNDX = 0
GO TO 1250
1 2 If 0 I< J rw X = 5
LJNDX = JYEAR-7
12S0 DO 1257 MPCITY= I,INC
KIT = ICHOCMPCITY)
LCITY = KCITY(KIT)
DO 1255 JGE ='1,20
IF (XXHCLCITY,JGE).EQ.O.O) GO TO 1255
XXTCLCITY,JGE) = XXH(LCITY,JGE).(GRFAC(LCITY,JGE,I)+
I(KJNDX/5.CGRFAC(LCITY,JGE,2)-GRFAC(LCITY,JG~'I»)
1 +CLJNDX/S.CGRFAC(LCITY,JGE,3)-GRFAC(LCITY,JGE'2»»
1255 CONTINUE
1 2 S 7 C 0 ~n I N U E
1260 DO 1262 MPCITY= I,
. KIT = ICHOCMPCITY)
LCITY = KCITY(KIT)
DO 1261 JGE = 1,20
IF CXXHILCITY,JGE).EQ.O.O) GO TO 1261
XKTILCITY,JGE) = XXT(LCITy,JGE)/VK(JGE)
0/i76
0'+77
0/i78
Oq79
0/i80
0'+81
Oq82
0'+83
0'+8q
OQ85
OQB6
0/i87
0'+88
0/iS9
OQ90
OQ91
OQ92
OQ93
O/i9Q
OQ95
OQ96
0'+97
OQ98
O/i99
0500
0501
0502
0503
050Q
0505
0506
0507
INC
:»AGE 001/
~'SOOb09
R~S:Ob a
RMS00611
RMS0061Z
RMSOQ613
RMS0061/i
RMS00615
RMSOOb16
RMSOOb17
RMS00618
RMS00619
RMSOOb20
RMS00621
RMSOOb22
RMS00623
RMS0062'+
RMS00625
RMSOOb26
RMSOIJ627
RMS00628
RMSOOb29
RMS00630
RMS00631
RMS00632
RMS00633
RMS0063'f
RMSOOb35
RMS00636
RMS00637
RMS00638
RMS00639
RMS006'fO
RMS006'f1
RMS006'f2
RMS006Q3
RMS006'f'f
RMSOOb'f5
RMS006'f6
-------
FO~"RAN
0508
0509
0510
051 1
0512
0513
051~
0515
0516
0517
0518
0519
0520
0521
0522
0523
052~
0525
0526
0527
0528
0529
0530
0531
0532
0533
053~
0535
0536
0537
V G LEVEL
20
MA . ,
(OHINUE
(ONTINUE
1261
1262
(
(
(......
( REMOVE AFTER THIS RUN
(
659't
6595
6596
WRITE(6.6591i)
WRITE(6,6595) PERRC
WRITE(6,6591i)
VIRITE(6,6596) XX
FORMAT(IHl)
FORMAT(lH ,20F5.2)
FORMAT(lH ,10FI0.2)
C..........
C
C
635'1
FORMAT(lH tiS)
00 6355 HI = I.IN(
IK = ICHO(MI)
I = KCITY(II()
PI(I,20) = PI(I,20) . PERPIN(I)
XI(I,20) = XI(I.20) Ot PERXIN(I)
XK(I,20) = XK(I.20) . PERXKN(I)
XX(I,20) = XX(I.20) . PERXON(I)
(I) = C(I) .-PERCN(I)
YY(I) = YY( I) . PERYYN(I)
SU~1N(I,20) = SUMN(j'20). PERSUM(I)
XNBAR(I) = XNBAR(I) . PERONB(I)
XNT(I) = SUMN(I.20) + XNBAR(!)
UtI) = PERUN(I)
G(I) = G(I) . PERGN(I)
T(I) = T(I) . PERTN(I)
QC(I) = QC( I) . PERQCN(])
QH(t) = QM( I) . PERQHN(I)
QBAR«() = QBAR(I) . PERQBR(I)
XL(I) = XNT{J) 1 (1. - PERUN{J»
QT«() = QC(I) + QM(I) + Q8ARII)
DATE = 72006
1011 1 108
'AGE (OlL
RMS006C7
RMS006Lff
RMS0061i9
RMS00650
RMS00651
RMS00652
RMS00653
RMSOU6S't
RMS00655
RMS00656
RMS00657
RMS00658
RMS00659
RMS00660
RMS00661
RMS00662
RMS00663
RMS00661i
RMS00665
RMS00666
RMSOQb67
RMS00668
RMS00669
RMS00670
RMS00671
RMS00672
RMSOOb73
RMS0067't
RMS00675
RMS00676
RMS00677
RMS00678
RMSOU679
RMS00680
RMS00681
RMS00682
RMS00683
RMS0068Lf
-------
FUKTKAN IV G LEVEL
20
MAIN
)A"E = 72006
10/1 'U08
:)AG~ (01
R'S(06~5
R~s00686
RMS00687
RMs00688
RMS00689
RMS00690
RMS00691
RMS00692
RMS00693
Rr'1S0069't
RMS00695
RMS00696
RMS00697
RMS00698
EMP -NOW INRMS00699
Q-NOW AQCR')RMS00700
RMS00701
RMS00702
RMS00703
RMS0070't
RMS00705
RHS00706
RMS00707
RMS007U8
RMS00709
RMS00710
RMS00711
RMS00712
RMS00713
RERMS0071't
RMS0071S
RMS00716
RMS00717
RMS00718
RMS00719
QRMS00720
RMS00721
RMS00722
0538
6355
C
C
C
C
C
CONTINUF:
DATA TABLES 1970-1780
0539
05't0
05'tl
05'i2
05'i3
05't'i
05't5
05'i6
05'i7
0552
0553
055't
0555
0556
0557
0558
IFIISBUG2.EQ.0) GO TO 'tOO
00 1350 JRKK = 1,INC
JJRK = ICHO(JRKK)
JRK = KCITY(JJRK)
WRITE(LHRIT,1315)
VI R I T E ( L 'J R IT, 1 3 1 1) (I A Q C R I J R K , K R K ) ,K R K = 1 , 2 0 ) ,K Y Y
1311 FORMAT(IH ,1't,19A't,'YEAR = ',I't,'= NOW (T)',/)
WRITEILWRIT,1312)
1312 FORMAT(IH ,'IND VA 1967 VA NOW VA IT-I)
IV-h7 INV NO~ WB 1967 WB NOW PI NOW
1315 FORMAT (IHl)
00 13'i0 KRK = 1,20
IF (XXH(JRK,KRK).EQ.O.O) GO TO 13'+0
WRITE(LWRIT,1310) KRK,XXH(JRK,KRK),XX(JRK,KRK),XXT(JRK,KRK),
1 SUt-1NH(JRK,KRK) ,SUMN(JRK,KRK) ,XIH(JRK,KRK) ,XI (JRK,KRK),
2WBH(JRK ,KRK) '~IB(JRK,KRK) ,pI (JRK,KRK).
3 EOLDH(JRK,KRK) ,EOLD(JRK,KRK,5),JRK
132U FURMAT(IH ,10FI0.l,2X,I3)
13't0 CONTINUE
I F ( Y Y H ( J R K ) . E Q.. 0 . 0) GOT 0 1 3 5 0
WRITE(LWRIT,1323)
1323 FORMAT(IH ,I)
WRITE(LWRIT,1321)
1321 FORMAT(IH,' INC 1967 INC NOW EMP(C)-67 EMp(C)NOW REV 1967
IV NO~ EXP 1967 EXp NOW CONSMP-67 CONSMP NOW AQCR')
WRITE(UVRIT,1320) YYH(JRK) ,YY(JRK) ,XNBAH(JRK) ,XNRAR(JRK).
1 TH(JRKI ,T(JRKI ,GH(JRK) ,GIJRKI ,CHCJRK) ,CIJRKI ,JRK
WRITE(L~RITtI323)
WRITECLWRIT,1322)
1322 FORMATCIH " EMPCT)-67 EMP(T)NOW LABOR-67 LABOR NOW QM 1967
1M NO~ QC 1967 QC NOW QBAR 1967 QBAR NOw AQCR') I
I
WRITECLWRIT,1320) XNTHCJRK) ,XNTIJRK) ,XLHIJRK) ,XLlJRK) ,QMHCJRK),
EMpCM)
Q-67
05'i8
05'+9
0550
0551
0559
0560
0561
0562
0563
-------
F O""'RTR A ~
056'i
0565
0566
0567
0568
0569
0570
0571
0572
I. G LEVI:.L
zo
MA 'J
DATE = 72006
10/li'/OB
3AGE : 02;
JRI<) , QAAR 1 JRK) , JRK
lQMIJRK) ,QCHIJ~() ,QCIJRK) ,QBARH
1350 CONTINUE
C
C
C
C
C
RMS00723
RMS00721
RMS00725
RMS00726
RMS00727
RMS00728
GO TO 1355 RMS00729
'i00 CONTINUE RMS00730
C RMS00731
C RMS00732
C RMS00733
C[[[RMS0073'i
C
C
C
C
C
C
C
C
C
C
MAKE OPTIONAL - ONLy DATA TABLES GIVEN
IF ISBUG2= 1 AS IT IS HERE
GENERAL PARAMETER CARD
RMS0073S
RMS00736
RMS00737
RMS00738
RMS00739
RMSOQ7'i0
RMS007'il
RMSOU7'i2
Rr~S007'i3
RMS007'i'i
RMS007'iS
RMSOU7'i6
RMS007'i7
RMS007'i8
ISBUG2=1 DATA TABLES FOR 1970'5 REQUESTED
ICTAB = 1 CONVERGENCE PROCESS FOR NEWTON RAPSON REQUESTED
IGOV = 1 FOR GOVERNMENT ASSISTANCE
REAOCLREAD,'il0) INPT(1),NOPE,CLPARACI), I=I,NOPE),
1 MCNT,INPT(2),NOPEl,LPNT1,ILPNT2II), I=I,NOPEI ),NUMCIT
2 ,IPRNI dPRN2 ,IPRN3 ,rSBUG dSBUG2
3 dCTAB,IGOV
c......................................~................................RMS007'i9
C
7183
718'+
C
C
C
C
00 7183 1= 1,2
REAOCLREAO,71BIi)
FOf~MATC20A'+)
RMS00750
RMS007S1
RMS00752
RMS00753
RMS007S'i
RMS007S5
RM500756
RMS00757
RMS00758
RMS00759
RMS00760
ITITLECJ,I), J=I,20)
-------
FO n ~AN
0573
057Lf
0575
0576
0577
0578
0579
0580
0581
0582
0583
n58Lf
0585
0586
0587
0588
0589
0590
0591
0592
0593
059~
0595
0596
0597
0598
0599
0600
0601
0602
0603
060~
0605
0606
0607
0608
0609
0610
V G _:::VEL
20
MA N
ADJUST = fEII:::;-.(B:' (JY:A~-3)fBEN(8
REA)I_R::AD,717Lf) QG VE~
717Lf FORMAT(7FIO.Lf)
7175 CONTINUE
KPRNI = IPRNI
KPRN2 = IPRN2
KPRN3 = IPRN3
KSBUG = 15BUG
IFIINPTII)+INPTI2).EQ.O) STOP
NOPEN = NOPE + NOPEl
NUMCIT = NUKCIT + I
IHELP = 0
700 CONTINUE
NUMCIT = NUMCIT - 1
NEiiCIT = 1
NBEFOR = 0
NBEFOI = 1
NBEF03 = 0
NBEF02 = 1
IPRNI = KPRNI
IPRN2 = KPRN2
IPRN3 = KPRN3
ISSUG = KSBUG
NOPEI = NO PEN-NOPE
NOK = NOPE
IHELL = 0
705 IHELP = IHELP + 1
IF(IHELP.LT.3) GO TO 706
IHELP = 1
IF (NOPEN.EQ.1)
LSTRAT = 8
IKSW = 1
DO 7999 K = I,MCNT
IK = MCITY(K)
I = KC1TYIIK)
CALL PRIN1 (I)
XXI92,20) = XX(92,20) + XX(I,20)
PI(92,20) = PI192,20) + PIII,20)
GO TO 70Lf
)A': = 72006
:»AG::: (02l
10/' 'f/oa
R'1S0J76
R'1S00762
RMS00763
RMS0076'f
RMS00765
RMS00766
RMS00767
RMS00768
RMS00769
RMS00770
RMS00771
RMS00772
RMS00773
RMS0077Lf
RMS00775
RMS00776
RMS00777
RMS00778
RMS00779
RMS00780
RMS00781
RMS00782
RMS00783
RI'1S0078'i
RMS00785
RMS00786
RMson787
RMS00788
RMS00789
RMS00790
RMS00791
RMS00792
RMS00793
RMS0079'i
RMS00795
RMS00796
RMS00797
RMS00798
-------
FO nRAN
06.1
06 2
0613
061'i
0615
0616
0617
0618
0619
0620
0621
0622
0623
062'i
0625
0626
0627
0628
0629
0630
0631
0632
0633
063'i
0635
0636
0637
0638
0639
06'iO
06'i1
06'f2
06'f3
06'f'f
06'f5
.V G LEVEL
(
7987
7988
(
(
7999
20
MAN
DA"::=72006
X 92,20 = XI(92,20) + XI( ,20)
X(192,20) = XK(92,ZO) + XK(I,20)
SUMNI92,20) = SUMNI92,20) .SUMN(I,20)
XNRARI9Z) = XNBAR(92) + XNBAR(I)
YY(92) = YYI9Z) + YYII)
(192) = (19Z) . (II)
XNT(92) = XNT(92) + XNTII)
XL(92) = XL(92) . XL(I)
QT(92) = OT(92) + OT(J)
QM(92) = QM(92) + QM(I)
QBAR(92) = QBARI9Z) . QBARII)
QC(92) = Q(192) + Q(J)
G(921 = G(92) + G(I)
T(92) = T(92) + T(J)
PERPIN(J) = (1. + (PIN(I)/PI(I,20») . PERPIN(I)
PERXIN(I) = (1. . (XIN(I)/XIII,20») . PERXIN(I)
PERXKN(I) = II. + (XKNII) / XK(I,zO») . PERXKN(I)
PERSUM(I) = (1. + (SUMNN(I) / SUMN(J,ZO») . pERSUM(I)
PERXDN(I) = (1. . (XON(J)/XXIJ,20») . PERXDN(I)
PERYYNII) = (1. + (YYN(J)/YY(I») . PERYYN(I)
PERGNII) = (1. + (GN(I)/GII») . PERGN(I)
PERTNII) = (1.+ (TNIJ) / T(I») . PERTNII)
PERCNII) = (1. . (CN(I)I (II») . PER(NIJ)
PERDNB(I) = (1. + IDNBARN(I) / XNBARI I») . PERDNB(I)
PERQBRII) = II'. . (QBARN(I) / QBAR(I») . PERQBR(I)
PERQ(NII) = 11. . (QCN(I) / Q(II») . PERQ(NII)
IF(QMII» 7988,7988,]987
CONTINUE
PE~QMNII) = II. . (QMN(I) / QM(I») . PERQMN(J)
GO TO 7999
PERQMN(I) = 1.0
CONTINUE
IKSW = 0.0
. .
70'f IFINUMCJT.LT.l.AND.IOYEAR.NE.1967)
IFINUMCIT.EQ.O) GO TO 'fOO
GO TO 1355
10/' 'f/08
:JAG:: 0)2:
R'S00799
RMS00800
RMS00801
RMS00802
RMS00803
RMS0080'f
RMS00805
RMS00806
RMS00807
RMS00808
RMS00809
RMS00810
RMS00811
RMS00812
RMS00813
RMS0081'f
RMSOOB15
RMS00816
RMS00817
RMS00818
RMS00819
RMS00820
RMS00821
RMS00822
RMS00823
RMSOn82'f
RMS00825
RMS00826
RMS00827
RMS00828
RMS00829
RMS00830
RMS00831
RMS00832
RMS00833
RMS0083'1
RMS00835
RMS00836
-------
FOR'RAN I V G LEVEL
061 6
06..7
06..8
06..9
C
C
C
6001
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
0650
0651
0652
0653
065"
20
MAIN
DA'~ = 72006
10/""/08
)AG:: :)2
706 INTT = INP'(HELP)
IFIINTT.EQ.O) GO TO 705
GO TO ("8I,10),INTT
RMSOOl37
RMS00838
RMs00839
RMS008"0
RMSOua&fl
RMSOOaQ2
RMSOoaQ3
RMSOOaQq
RMSOOS..5
RMSOOS&f6
RMSOOSq7
RMSOOSqa
TWO DIGIT INDUSTRy INFORMATION .. TWO ALTERNATIVES... RMS008Q9
(I) DETAILED SIC INFORMATION ONE SMSA AT A TIME (LPNTI = lRMSOOS50
(2) ONE TWO DIGIT SIC ( AT A TIME) FOR ALL CITIES (LPNTl=RMsOOS51
WITHIN EITHER CATEGORY, THE ECONOMY CAN BE AFFECTED IN ONE OF TRMSOOS52
DEPENDING ON THE DATA 6IVEN. RMs00853
II) COST FIGURES ARE GIVEN DIRECTLY (LPNT2=1) RMsOOS5Q
(21 PRICE CHANGES FOR FUEL ARE GIVEN (LPNT2=2) RMsOQ85S
NOX= NO. OF CITIES FOR WHICH TWO DIGIT SIC COST DATA IS GIVEN RMs00856
(LPNT1=1) RMS008S7
RMs00858
RMSOOa59
RMSOOS60
RMsOOS61
RMsOOS62
RMSOOS63
RMS0086Q
RMS0086S
RMS00866
RMSOOS67
RMS00868
RMSOOS69
RMs00870
RMS00871
RMS00872
RMSOUS73
RMSOOS7&f
FORMAT(II0,5FI0.2)
10 CONTINUE
720 KOUNT = 0
IHELL = IHELL + 1
LPNTT = LPNT2(IHELL)
LSTRAT = LPNTT + 5
C
0655
0656
0657
0658
0659
0660
5S 00 70 1=1,91
ONR(I) = 1.0
PDN(I) = 0.0
PPIO(I) = 0.0
PXIDII) = 0.0
PXKDII) = 0.0
-------
FORTRAN
0661
0662
0663
0661i
0665
0666
0667
0668
0669
0670
0671
0672
0673
067/f
0675
0676
0677
0678
0679
0680
0681
0682
0683
.V G _~VEL
C
C
C
C
C
C
C
C
C
C
C
C
C
C
5992
C
2(
MA II
DA"E = 72006
011 I /oa
'X)I ) = 0.0
EDII,S) = 0.0
PIOII) = 0.0
WBDII) = 0.0
XIDII) = 0.0
XKDII) = 0.0
PDNII) = 0.0
70 XDII) = 0.0
IF LPNTI = I ... COST DATA FOR EACH CITY IS PRECEEDED BY THE
ACTUAL CITY CODE NUMBER AND THE NO. OF IND. IN THAT CITY FOR
WHICH COST DATA Is GIVEN.
IF LPNTI = 2 COST DATA FOR EACH IND. IS PRECEEDED BY THE SIC
NUMBER AND THE NO. OF CITIES FoR wHICH THAT IND. EXISTS AND
COST DI\TA.
MCNT = 1
I I = I
NOTE = 1
MAX = a
IF ILPNTT
.NE. 3) GO TO 51f
GENERAL FEEDBACK SEGMENT
OPTION NUMBER EIGHT
FORMATI217FI0.lit/);5F10.If.,23X,I3)
KOUNT = KOUNT + I
NOI = 0
IFEED = 1
IK = MCITYll)
IK = MMCITYIKOUNT)
LCITY = KCITYIIK)
00 S6 I = 1 tl9
IFIINDCIILCITYtI) .EQ. 0) GO TO 56
NOI = NoI + 1
INDINOI) = I
"AGE ou~'t
RMS:(875
RMSH876
RMS00877
RMsooa78
RMS00879
RMS00880
RMSOOB81
RMSOOBB2
RMS00883
RMSOOa81i
RMsouaa5
RMS00886
RMsooa87
CORMS00888
HARMS008B9
RMS00890
RHSOOB91
RHsooa92
RMS00893
RMSOOB91i
RMsooa95
RMSOOB96
RMSOOB97
RMsooa98
RMS00699
RMS00900
RMS00901
RMS00902
RMS00903
RMS0090/f
RMS0090S
RMS00906
RHS00907
RMS0090a
RMS00909
RMS00910
RMS00911
RMS00912
-------
FO~"RAN
068'
0685
0686
0687
0688
0689
0690
0691
0692
0693
069"+
0695
0696
0697
0698
0699
0700
0701
0702
0703
070"+
0705
0706
0707
0708
0709
0710
071 1
0712
0713
071"+
0715
0716
V G LEVEL
20
H A I I
DA":: = 7l)06
o fLIifI 0 8
lTAGE-uu25
56 CONTIIJU~ .
IFIIGOV.EG.I.AND.IJO .1'1::.0) GO"O 3659
IFINOI.EQ.O.OR.ADJUST.EQ.O) GO TO "+3
3659 CONTINUE
MIN = NOI
MCTTYIII) = IK
00 1f2 I = l,tHN
INDI = INDII)
IUN = INDI
PXDIIUN) = BILCITY,IND1) . ADJUST
EDIIUN,5) = ABIIUN).PXDIIUN)
PPIDIIUN) = 11.- BETAIIND1».PXDIIUN)
PXIDIIUN) = PICOEFI INDI ).ppiDI IUN)
1f2 PXKOIIUN) = PXIDIIUN)/Il. + DIIND1»
GO TO 17
1f3 IFf ED = 0
IFIIPRN3.EQ.0> GO TO 5003
WRITE ILWRIT,930)
930 FORMAT 1//,IH ,'INTERREGIONAL FEEDBACK EFFECT NOT AVAILABLE fOR
lIS AQCfP)
GO TO 5003
/ THISI
C
C
C
3121A) COST DATA - INPUT FOR DETAILED SIC
51f IFIIA312.EQ.I) GO TO 7151
READILREAD, 1.f(0) JCITY,NOI
GO TO 7153
7151 NOI = 6
MIN = NOI
READILREAD,7152) HHI,JCITY
LCTTY = KCITYIJCITY)
MCITYIII) = JCITY
READILREAD,7152) HHA,JCITY
DELTCILCITY) = 0.0
7152 FORMATI7FI0.I.f,3X,I3)
GO TO 725
7153 CONTINUE
C
RMS00913
RMS009 . 'f
"RMS00915
RMS00916
RMS00917
RMS00918
RMS00919
RMS00920
RMS00921
RHS00922
RMS00923
RMS00921f
RMS00925
RMS00926
RMS00927
RMS00928
RMS00929
RMS00930
THRHS00931
RMS00932
RMS00933
RMS0093'f
RMS00935
RMS00936
RMS00937
RMS00938
RMS00939
RMS009lfO
RMS0091.fl
RMS0091.f2
RHS0091.f3
RMS0091.f1f
RMS0091.f5
RMS0091.f6
RMS009'f7
RMS0091f8
RMS009'i9
RMS00950
-------
FOR'RAN
07 7
0718
0719
0720
0721
0722
0723
072'1
0725
072b
0727
0728
0729
0730
0731
0732
0733
073'1
0735
0736
0737
0738
0739
07'10
07'11
07'12
07'13
07'1'1
07'15
07'16
07'17
07"8
07'19
0750
0751
V G LEVEL
Zo
MA N
.CITY = KC"Y(JCITY)
NDI = NO I -20 + (ZO/NOI)
IFILPNTl.EQ.Z) GO TO 1
MC I T Y ( I I) = JC I TY
MIN = NOI
GO TO 2
1 MCNT = JCITY
MIN = JCITY
MAX = JKEEP(IND1)
C
2 ITT = 5
IFILPNTT.NE.l) GO TO 725
1 T T = 1
REAO(LREAD,730) DELTC(LCITY)
DELTC(LCITY) = DELTC(LCITY) I
730 FORMAT(10X,fl0.2)
725 CONTINUE
C
C
IH = 0
J = 0
DO 5Sb JJ=I,MIN
If(IA312.EQ.l) GO TO 715"
READ (LREAO,6001) IK
HOlO( 1) = 0.0
HOLO(2) = 0.0
HOLD(3) = 0.0'
HOlO(") = 0.0
HOlO(S) = 0.0
HOLO(") = .2-POLD(LCITY,IK,q)
REAO(lREAD,6001) IK,(HOLO(IS),
IFILPNTl.EQ.l) GO TO 3
IUN = KCITY(IK)
LCITY = IUN
11 C I T Y ( J) = I K
GO TO "0
3 INO] = IK - 20 +(20/IK)
NEXT = JKEEP( IN01)
IF (MAX.LT.NEXT) MAX = NEXT
DATE = 72006
10/q'U08
::lAG ~ 002(,
QT(LCITY)
R~S)095
R'1S00952
RMS00953
RMS0095Q
RMS00955
RMS0095b
RMS00957
RMS0095B
RMS00959
RMS00960
RMS00961
RMS009b2
RMS00963
RMS0096Q
RMSOQ9bS
RMS00966
RM500967
RMS00968
RMS00969
RMS00970
RMS00971
RMS00972
RMS00973
RMS0097Q
RMS00975
RMS00976
RMS00977
RMS00978
RMS00979
RMS00980
RMS00981
RMS00982
RMS00983
RMS0098"
RMS00985
RMS00986
RMS00987
RMS00988
15=1 dTT)
-------
C
C
C
0766
C
0767
0768
0769
0770
0771
0772
0773
077'1
0775
0776
0777
0778
0779
0780
0781
0782
0783
FORTRA~
0752
0753
075"1
0755
0756
0757
0758
0759
0760
0761
0762
0763
076"1
0765
V G LEVEL
20
DA.~ = 72006
MA ' I
10/'lCf/08
UN = NDI
NDIJ = IN)
'10 GO TO I'll,,!), LPNTT
C
715"1
IK = ISPINIJJ)
INDI = IK - 20 + 120/IK)
HOLD 1 1) = HHAIJJ)
HOLDIZ) = HHIIJJ)
IUN = INDI . .
IFIINDCIILCITY,IND1).EQ.O)
J = J + 1
INf)IJ) = INDI
GO TO '11
IH = IH + 1
GO TO 556
GO TO 7155
7155
LPNT2
COST FIGURES GIVEN DIRECTLY
= 1 ...
Cfl CONTINUE
p = 1.
IFIIGOV.EQ.O) GO TO 36Cfl
HERE = HOLDIl)
HAHA = HOLD(2)
HOLDIl) = P . HERE
HOLDIZ) = P. HAHA
36'11 CONTINUE
PPTDIIUN)
PXIDIIUN)
PXKDIIUN)
K = LCITY
PXOIIUN) =-XXILCITY,INDI )+xXILCITY,IND1).IPXKD(IUN)+XKCK,INOl))
1 IXKILCITY.INDl)
EDIIUN,5) = ABIIUN).PXOCIUN)
GO TO 556
Cf DO 5 K=I.5
IFIPOLDILCITY.INDI ,K)~EQ.O.)HOLD(K) =0.0
5 PNEWDCIUN,K) = POLDllCITYtlNOl.K) + HOLD(K)
556 CONTINUE
= -HOLDll) - OELTC(LCITY) . EOLDClCITY,IN01.S) . P
= PICOEFtINDl).PPID(IUN) - HOLD(2)
= PXIOIIUN)!11.+DCIND1))
)AG~ ):2
RMS00989
RMS00990
RMS00991
RMS00992
RMS00993
RMS0099'1
RMS00995
RMS00996
RMS00997
RMS00998
RMS00999
RMSOI000
RMSOI001
RMSOI002
RMSOl003
RMSOI00Cf
RMSOI005
RMSOl006
RMSOI007
RMSOl008
RMSOI009
RMSOI0I0
RMSOI0l1
RMSOI012
RMSOI013
RMSOIOICf
RMSOI015
RMSOI016
RMSOl017
RMSOl018
RMSOI019
RMSOI020
RMSOI021
RMSOI022
RMSOI023
RMSOl021.f
RMSOJ025
RMSOI026
-------
FORT~AN IV G LEVEL
( 781
0785
0786
0787
0788
0789
0790
0791
0792
0793
079li
0795
0796
0797
0798
0799
20
)A":: = 72006
10/«i'i/08
MA N
llN = ~ I' -
FIMP.EQ.O.AND.IA32.::Q.l)
IFCLPNTI.EQ.l) NOI = MIN
IFILPNTI.EQ.2) JCITY = MIN
IFCLPNTT .EQ. 1) GO TO 17
GO TO 5555
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
LPNT2
= 2 ...
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
ZZ - TOTAL QUATITY OF FUEL CONSUMED
XI,XK, ANI) PI
EXOGENOUS VARIABLES...
POLO - PRICE OF FUEL 1967
PNEW - PROJECTED PRICE OF FUEL FOR 1967
WITH POLUTION COSTS
GAMA - FUEL EXPONENTS - CONSTANTS 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
ICHECK = 0
MKEEP = MAX + 2
DO 27 I=I,MKEEP
27 KEEPCI) = LKEEPCMAX,I)
6 MM=O
KPlsT = 18
M = 0
DO 12 I=I,MIN
IFILPNT1.EQ.2) GO TO 7
INDI = INDCI)
IUN = INDI
:tAG:: OJ2f
RMSOl( 27
RMsOI02~
RMsOI029
RMSOI030
RMsOI031
RMSOI032
RMSOI033
RMsOI03tf
RMsOI035
RMSOI036
RMsOI037
RMSOI038
RMsOI039
RMsOIO'iO
RMsOI0'f1
RMsOIOlf2
RMSOIO'i3
RMSOIO'f&J
RMsOIO'fS
RMsOI0'i6
RMSOI0'f7
RMsOIO'l8
RMsOIOlf9
RMSOI050
RMSOI051
RMsOI052
RMsOI053
RMSOI051f
RMsOI055
RMsOI0S6
RMSOI057
RMSOI0S8
RMSOI0S9
RMsOI060
RMSOI061
RMSOI062
RMSOI063
RMSOI061f
-------
FO ~"RAN I V G - ~VEL
20
MA N
)A'~ :: 72006
10/~£f/08
)AG~ :02~
0830
0831
0832
0833
083'i
0835
0836
MARK:: 't
LVAR :: 11 + (MKEEP-5)
MCASES :: MIN.LVAR
(ALL SIMUL
IFIICHECK.EQ.I) GO TO 17
ICHECK = 1
M :: 0
RMSOI065
RMSOI066
RMSOI067
RMSOI068
RMSOI069
RMSOI070
RMSOI071
RMSOI072
RMSOI073
RMSOI07't
RMSOI075
RMSOI076
RMSOI077
RMSOI078
RMSOI079
RMSOI0BO
RMSOI081
RMSOI082
RMSOI083
RMSOI08~
RMSOI085
RMSOI086
RMSOI0B7
RMSOI088
RMSOI089
RMSOI090
RMSOI091
RMSOI092
RMSOI093
RMSOI09'i
RMSOI09S
RMSOI096
RMSOI097
RMSOI098
RMSOI099
RMSOII00
RMSOI101
RMSOI102
08 )(
080
0802
0803
080'i
0805
0806
0807
0808
0809
U810
081 1
0812
0813
D811f
0815
0816
0817
0818
0819
0820
0821
0822
D823
D821f
0825
0826
0827
0828
0829
GO TO 8
( :: ~C "Y( )
LeTTY:: KCITYIIK)
IUN :: LCITY
8 00 9 K::l,MKEEP
KK :: KEEP(K)
X(MM+K) :: POLDILCITY,IND1,KK)
IFIICHECK.EQ.l) X(MM+K) :: pNEWD(IUN,KK)
9 X(MM+K+MKEEP) :: GAr1AIIND1,KK)
MM :: MM+2.IMKEEP-S)
X (MM+ 1 1) :: X K T I LC I T Y , I ND 1 )
XIMM+12) :: DIIND1)
XIMM+13) :: 1.- BETA(IND1)
XIMM+llf) :: PICOEF(IND1)
XIMM+15) :: AIIND1)
X(MM+16) = ALITIIND1)
XIMM+17) = XXILCITY,IND1)/XK(LCITY,IN01)
XIMM+18) :: PICON(IND1)
00 11 K=I.MKEEP
KK = KEEPIK)
11 Y(M+K) :: EOLOILCITY,INDI ,KK)
M :: M + MKEEP - 5
YIM+6) :: XKILCITY,IN(1)
YIM+7) :: CTOT67ILCITY,IND1)
YIM+8) :: XXILCITY,IND1)
YIM+9) = ZZILCITY,IND1)
YIM+I0)= PIILCITY,IN(1)
YIM+ll) :: XIILCITY,IND1)
M :: M + 11
12 MM :: MM + 18
7
C
-------
FO n ~AN
08: 7
OS38
0839
08liO
08lil
081i2
081i3
081i1i
D81i5
OSli6
081i7
081i8
081i9
0850
0851
0852
0853
0851i
0855
0856
OS57
0858
0859
0860
0861
0862
0863
0861i
086S
0866
0867
0868
0869
0870
0871
0872
0873
V G _~V~L
20
MA N
00 16 = ,MIN
IFILPNT1.EQ.2) GO TO 13
INOI = INDII)
IUN = INOI
GO TO Iii
13 IK = MCITYII)
LCITY = KCITYIIK)
IUN = LCITY
1"+ 00 15 K=I,MKEEP
15 EOIIUN,K) = YIM+K)
M = M + MKEEP - 5
PXKOIIUN) = YIM+6)
CTOT07IIUN) = YIM+7)
PXO(IUN) = YIM+8)
PZOIIUN) = YIM+9)
PPIDIIUN) = yIM+IO)
PXIDIIUN) = YIM+l1)
16 M = M + 11
GO TO 6
c
17 CONTINUE
M = 0
00 22 IT = I,MIN
IFILPNT1.EQ.2) GO TO 18
INOI = INDIITI
IUN = INDI
GO TO 19
18 IK = MCITYIITI
LCITY = KCITYIIK)
IUN = LCITY .
19 IFILPNTT.EQ.I.0R.LPNTT.EQ.3) GO TO 21
DO 23 K=I,MKEEP
23 EDIIUN,K) = EDIIUN,K) -YIM+K)
M = M + MKEEP - 5
PXKDIIUN) = PXI
-------
FOR' ~AN IV G LEVEL
)874
0875
0876
0877
0878
0879
0880
0881
0882
0883
088"
0885
0886
0887
0888
0889
0890
0891
0892
0893
0894
20
MAIN
)A':: =.72006
:)AG:: (03
P~IDI UN) = PPIDIIUN) - YIM+IO)
~XID IUN) = PXIDIIUN) - XIM+ll)
M = M + 11
21 I = LCITY
IFIWII,IND1).SUMNII,20).EQ.0.O) GO TO 39
PDNIIUN) = BETAIINDl).PXDIIUN)/W(I,IND1)
DNRII) = DNR(I) -(PDNIIUN)/SUMNII.20»)
39 CONTINUE
37 FORMATIIH ,ZI3,FI2.1,FIZ.,,)
PIOII) = PID(I) + PPIOIIUN)
XI0II) = XIOII) + PXIDIIUN)
XKOII) = XKDII) + PXKOIIUN)
XDII) = XOII) + PXDIIUN)
XD1(1) = XOII)
~BDII) = WBDII) + BETAIIND1).PXDIIUN)
D2II) = W8DII)
22 CONTINUE
NOTE = 1
GO TO 5555
C
C
C
C
C
C
C
C
C
C
INPT(2) = 1
ALL CITY INFO. FCR ALL IND. REQUESTED
"81 CONTINUE
NOTE = NOPE
C
C
C
C
C
C
C
10/44/08
RMS01141
RMSOll,,2
RMSO 11 43
RMS011"4
RMSOll'iS
RMS011'i6
RMS011'i7
RMS011'i8
RMSOll'i9
RMSOl150
RMS01151
RMS011S2
RMSOl1S3
RMS01154
RMSOllSS
RMSOl1S6
RMSOl157
RMS01158
RMS011S9
RMSOl160
RMSOl161
RMSOl162
RMS01163
RMS0116'i
RMS0116S
RMS01166
RMS01167
RMS01168
RMSOl169
RMS01170
RMS01171
RMS01172
RMS01173
RMSOl17'i
RMS0117S
GIVEN - PROFIT AFFECRMsOl176
GIVE - INCOME AFF~CTRMsOl177
RMSOl178
NOTE»>ABOVE READ INDICATES THE TYPE OF CONTROL COST DATA GIVEN
(1) NOPE> TYPES OF POLLUTION IMPACT TO BE EXAMINED FOR EACH CITY
= 2»INCOME CHANGE GIVEN
= 'i..ELECTRIC COST INCREASE
= S..ELECTRIC COST INCREASE
= 3»TAX CHANGE GIVEN
-------
FOR.RAN IV G LEVEL
20
MA . \
DAT~ = 72006
H/"/OB
PAG~ (03~
C
C
3) \ UM3ER O~ C .. ~S ~OR W'i CH )A..A IS G VEN - MO.
RMSOl179
RMSOIIBO
RMSOI1Bl
RMSOIIB2
RMS01183
RMSOI1B'+
RMSOI1B5
RMSOIIB6
RMSOIIB7
RMSOI1B8
RMS01189
RMS01190
RMS01191
RMS01192
RMS01193
RMSOI19'+
RMS01195
RMS01196
RMS01197
RMS01198
RMS01199
RMS01200
RMS01201
RMS01202
RMS01203
RMSOI20~
RMS0120S
RMS01206
RMS01207
RMS01208
RMS01209
RMS01210
RMS01211
RMS01212
RMS01213
RMSOI21'+
RMS01215
RMS01216
0895
0896
0897
0898
U899
0900
0901
0902
0903
IF INOPE-5) 102,102,,+82
'+82 WRITEIL~RIT,~83) NOPE,ILPARAII), 1=I,NOPE)
~83 FORMATIIH1,'BAD CONTROL COST PARAMETER CARD',,+I~)
STOP
102 DO 103 1=I,NOPE
IFILPARAII).GT.5.0R.LPARAII).LT.I) GO TO ~82
IFII.EQ. NOPE) GO TO 103
IF ILPARAII+I) .EQ. LPARAII» GO 10 '+82
103 CONTINUE
C
C ..~.. MAKE EASILY ADJUSTABLE READ
C
C 10~ READILREAD,10S) MCITYII),ICOSTIJ,I),J=I,5)
C
090~
IKSW = 0
C
C
C
ZERO LOOP REMOVED FROM HERE
0905
0906
IfIIDYEAR.EQ.1967) GO TO 9909
IFIIA312.EQ.I) GO TO 7100
C
C
C
C
C
COST INPUT fOR 1970-1980 SIMULATION
0907
0908
0909
0910
091 1
. 0912
0913
091~
0915
0916
0917
DO 10lf 1=I,MCNT
SCIN = 0.0
SRIN = 0.0
se = 0.0
SR = 0.0
109 READILREAD,10S) RIN,R,CIN,eC,MCITYII)
RISUM = RIN + R + CIN + CC
IfIRISUM.EQ.O) GO TO J 10
GO TO 111
110 NUMCIT = NUMCIT - I
NE\(!jCIT = 1
-------
FORTRAN
0918
0919
0920
0921
0922
0923
092tf
0925
0926
0927
0928
0929
0930
0931
0932
0933
2 )
)A"~ = 72006
)AGE 003:
V G LEVEL
M A I ~
10/'+'f/08
~(NUMCIT._". I GO '.0
GO TO 09
111 CONTINUE
KI = MCITY(II
L = KCITYlKII
IF(L.EQ.O) GO TO 10'+
RMS012 7
RM5012 ~
RM501219
RMS01220
RMS01221
RM501222
RM501223
RM50122'+
RMS01225
RMS01226
RMS01227
RM501228
RM501229
RM501230
RMS01231
RH501232
RMS01233
RM50123'f
RM501235
RM501236
RMS01237
RM501238
RMS01239
RMSOI2'+0
RMS012tfi
RMSOI2'+2
RMS012'f3
RMS012'+'+
RMS012'+S
RM5012'f6
RM5012'+7
RMS012'f8
RMS012'+9
RMS012S0
RM5012S1
RMS01252
RM5012S3
RMSOI2S'+
3SS
C
C
C
C
(
USED FOR COST 5TRAT NO.3 ONLY
C05T(tf,l) = (05T(,+,I) -.5
COST(3,}) = 30.
COsT(l,l) = (C-(R-.o)
(OFI(I) = CIN-(RIN-.6)
(05T(,+,I) = R-.6+RIN_.6
(OST(3tI) = 0.0
DELTC(I) = COsT(,+,I)/QT(L)
c
C
(
USED ONLY FOR COST STRAT. NO.2 AND J
COST(I,I) = COST(1 ,1)-(623.7-(05T(1,1 )/808.067 )
C
(
(
(
(
(
C
USED FOR COST STRAT NO.3 ONLY
COS T ( 1 , I) = COS T ( 1 , I I -.5
USED ONLY FO~ COST STRAT. NO.2 AND 3
COFI(I) = COFI(I)-t623~7.COFI(I)/808~067
(
(
(
C
USED FOR COST STRAT NO.3 ONLY
(OF 1 ( I) = (OF I ( I) ..5
105 FORMAT(20X,2FI0.l/20X,2F10.1,33X,I3.~X)
10'+ CONTINUE
GO TO 5SSS
C
C
C
INPUT FOR 312(A) REPORT
-------
FORT Ur
0931f
0935
0936
0937
0938
0939
09lfO
o9lfl
o91f2
o91f3
091f1f
091f5
091f6
o91f7
o91f8
o91f9
0950
0951
0952
0953
0951f
0955
0956
0957
0958
0959
V G LEVEL
7100
7101
7 1 1 1
C
20
)A'E = 72006
PAGE 003
MA N
10/'+'+/08
)0 7101f I=.,MCN'
R:ADILREAD,71: 1 CIN,SIr ,R,rc "Y I)
fORMATI3FI0.5,'+3X,13)
(OfHINUE
KI = MCITYII)
L = KCITYIKI)
IFIL.EQ.O) GO TO 710,+
COSTI! ,I) = SIN.(I.+A~A)..,+ +R
DELTCI}) = .60.IICIN.ll.+AIA»+(SIN.II.+AAA»)/QTIL)
COFIII) = CIN-II.+AIA)..,+O
IFIIGOV.EQ.O) GO TO 710'+
HERE = (OSTll.1)
HAHA = COFIII)
RMS0125S
RMS01256
RMS012S7
RMS012S8
RMS01259
RMS01260
RMS01261
RMS01262
RMS01263
RMSOI26'+
RMS01265
RMS01266
RMS01267
RMS01268
C
C......................................$................................RMS01269
RMS01270
DELTCII) = P . DELTC(I) RMS01271
COSTllt1) = P . HERE RMS01272
COF}II) = P . HAHA RMS01273
C RMS0127'+
( RMSo1275
C.................................................~~~~~~~~~~............RMS01276
C RMS01277
7101f CONTINUE RMS01278
GO TO 5555 RMS01279
RMS01280
RMS01281
RMS01282
RMS01283
RMSOI28'+
RMSU128S
RMS01286
RMS01287
RMS01288
RMS01289
RMSOi290
RMS01291
RMS01292
C
C
C
C
C
9909
COST INPUT FOR 1967 RUNS
CONTINUE
00 910'+ 1=I,MCNT
SCIN = 0.0
SRIN = 0.0
SC = 0.0
SR = 0.0
00 9109 J=I,S
READILREAD,9105)
RIN,R,CIN,CC,MCITY(I)
-------
FOR .RA~
)960
096
0962
0963
096'1
0965
0966
0967
0968
0969
0970
0971
0972
0973
097Lf
0975
0976
0977
0978
0979
0980
0981
0982
0983
V G LEVEL
9109
910Lf
9105
(
(
(
5555
(
(
(
(
(
(
(
(
(
(
(
1000
20
MA N
)A" ~ = 72006
10/LfLf/OS
SR[N = SRIN + CRIN/5)..6
SCIN = S(IN + C[N/5
SC = SC + CC/5
SR = SR + CR/5)..6
Kl = M(ITYC[)
L = KCITY(KII
[FCL.EQ.O) GO TO 910Lf
COSTCltll = SC - SR
COFIC[) = SCIN - SRIN
COSTCI.f.[1 = SR + SRIN
OELT(III = COST(Lf.I)/QTCL)
CONTINUE
FORMATC20X.2FI0.1/20X,2FI0.1.33X.13.LfX)
(ONTINUE
005000 It I = I.NOTE
I I = 1 . .. .
IFCINTT.EQ.2.AND.MIN.EQ.0.AND.IA312.EQ.l)
IF(INTT.EQ.21 GO TO 100
LSTRAT = LPARAIIIII
II = LSTRAT
GO TO 11000.2000.3000,1000.20001.LSTRAT
ALL (ITY I~FORMATION ALL IND..::::EFFECT OF DIRECT REDUCTION
OF INDUSTRY'S PROFITS BY AIR PLOUTION (ONTROL EQUIPMENT.COSTS
GO TO 1000
PID = CHANGE IN PROFIT
XIO = CHANGE IN INVESTMENT
XKO = CHANGE IN CAP. STOCK
XD = CHANGE IN VALUE ADDEO
WBD = CHANGE IN WAGE BILL
DO 1001 K=I.MCNT
KI = M(ITYIKI
I = K(ITYCKII
)AG~ ):3~
~IS01293
R'SOI29'
RMS01295
RMS01296
RMS01297
RMS0129S
RMS01299
RMS01300
RMS01301
RMS01302
RMS01303
RMS0130Lf
RMS01305
RMS01306
RMS01307
RMSOl308
RMS01309
RMS01310
RMS01311
RMS01312
RMS01313
RMS01311.f
RMS01315
RMS01316
RMS01317
RMS0131S
RMS01319
RMS01320
RMS01321
RMS01322
RMS01323
RMS01321.f
RMS01325
Rr1S0 1326
RMS01327
RMS01328
RMS01329
RMS01330
-------
FO ~.. ~AN
(198'
0985
0986
0987
0988
0989
0990
0991
0992
0993
099"1
0995
0996
0997
0998
0999
1000
1001
1002
1003
100"1
V G _~VEL
2(
)A"E = 72006
MAIN
/
~ (I.EQ.D) GO '0 001
DNRII) = 1.0.. .
IFIINTT.EQ.2.AND.MIN.EQ.0.AND.IA312.EQ.I)
IF (LSTRAT-Lf) 310,300,300
300 COST(II.K) = DELTCIK) . QMII)
COFIIK) = 0.0
GO TO 310
7192 COSTIII.K) = HHA(7)
COFI(K) = HHI(7)
LSTRAT = 1
IFIIGOV.EQ.O) GO TO 310
C
GO TO 7192
10/' Lf/OB
:JAG ~ 003:
RMS0133
RMS01:32
RMS01333
RMSOI33'f
RMS0133S
RMS01336
RMS01337
RMS01338
RMS01339
RMS013LfO
RMSOl3LfI
RMS013Lf2
C...~[[[RMSOI3Lf3
C RMS013LfLf
COSTIII.K) = p. COST(II,K) RMS013Lf5
COFIIK) = P . COFIIK) RMS013Lf6
C[[[RMSOI3Lf7
C
310 PID(I) =-COSTIII,K)
XID(I) = PICOEF(20).PIDII) -COFI(K)
XKD II) = XID(I) /ll.+D(20»
XDII) = -XX(I,20) + XXII,20) . «XKDII)+XK(I,20» / XK(I,20»
XDI(I) = XDII)
WBDII) = BETA(20).XDII)
02(1) = WBD
-------
FOR.~Ar
005
1006
1007
1008
1009
1 0 1 0
1 01 1
1 0 1 2
1013
101'i
1 0 1 5
1016
1017
1018
1019
1020
1021
1022
1023
102'i
1025
1026
1027
1028
1029
1030
103 1
1032
1033
1031.f
1035
1036
V G LEVEL
20
)AT~ = 72)06
M A . r
I O/LfI /08
~AG~ (03
C
RMS01369
RMS01370
RHS01371
RMS01372
RMS01373
RMS0137Lf
RMS01375
RMS01376
RMS01377
RMS0137B
RMS01379
RMS01380
RMS01381
RMS01382
RMS01383
RMS0138Lf
RMS01385
RMS01386
RMS01387
RMS01388
RMS01389
RMS01390
RMS01391
RMS01392
RMS01393
RMS0139Lf
RMS0139S
RMS01396
RMS01397
RMS01398
RMS01399
RMS01LfOO
RMSOILfOI
RMS011.f02
RMSOILf03
RMSOILfOLf
RMSOILfOS
RMSOILf06
100 LCHEC( = 0
LVAR = Lf
MCASES = MCNT.1.j
LVAR = 3
MCASES = MCNT.LVAR
120 MM = 0
KPIST = I.j
M = 0
DO 1002 K=I,MCNT
KI = MCITY(K)
I = KeITY(KI)
IJeIT = I
IF (I.EQ.O) GO TO 1002
Y(M+l) = YYII)
Y(M+2) = G(I)
YCM+3) = e(I)
X(MM+l) = XX(I,20) + XOI(I)
X(MM+2) = T(I)
X(MM+3) = FUDGE(I,I)
X(MM+Lf) = FUDGEII,2)
M=r1+3
MM = MM + Lf
1002 CONTINUE
e
C
C
MARK = I
e
101Lf
951
8015
C
101 S
IF (ICTAS) B01S,80rS,101Lf
00 951 IM=I,LVAR
ARAB(IM) = SeCl(IM)
CALL SIMUL
M = 0
IF(LCHECK)1003,1003,LfOOO
1003 00 100Lf K= 1 ,MCNT
KI = MCITY(K)
I = K e I T Y ( K I )
IF (I.EQ.O) GO TO 100Lf
-------
FOR "RAN I V G _~VEL 20 M A I ' )A"E = 72006 10/'i1f/08 ;)AG ~ (03
)37 YYD(I) = Y ( 1 + ) RMS01')7
038 GD ( I) = Y(1+2) RMS011f08
1039 CDII) = YIM+3) RMS011f09
10'i0 XDIII) = 0.0 RMS011f10
1 O'i 1 M = M + 3 RMS011f11
10'i2 100'i CONTINUE RMSOllf12
10'i3 LCHECK = 1 RMS011f13
10'i1f NBEFOR = NBEF01 RMS011f1'i
10'i5 GO TO 120 RMS01'ilS
10'i6 If 000 M = 0 RMS01'i16
101f7 NEwCIT = 0 RMSOllf17
10'i8 NBEFOI = 0 RMSOllf18
10'i9 DO 1007 K=l,MCNT RMSOllf19
1050 K 1 = MCITYIK) RMS011f20
105 1 I = KCITYIKJ) RMS011f21
1052 1 F II.EQ.O) GO TO 1007 RMSOllf22
1053 YYDII) = YYDII) .. yIM+1) RMS011f23
105'i G D I I) = GO I I) .. Y(M+2) RMSOllf21f
lOSS CD I I) = CD I I) - YIM+) RMS01'i25
1056 M = M + LVAR RMS011f26
1057 TDII) = .10o'i66.YYOII) RMS011f27
1058 QCOIJ) = .01f718 . CO I I) RMSOllf28
1059 QMO I I) = .t991S . X 0 I I) RMS011f29
1060 QBARDII) = .0791f7.IYYDII)-XDII» RMSOllf30
1061 QTDII) = QCOII)+QMOII)+QBARDII) RMS011f31
1062 IF ILSTRAT-'i) '320.330.330 RMSOl1f32
1063 320 DNBARII) = .103239.IYYD(I)-XO(I» Rt1S01'i33
1061f DIll) = ONBARII) RMS011f31f
1065 GO TO 1007 RMSOllf3S
1066 330 ONAARII) = .103239 . IYYDI I )-XDI I) - DELTCIK) . QBARD(I» RMS011f36
1067 D 1 I I) = DNBARII) RMS011f37
1068 1007 CONTINUE RMS011f38
C RMSOllf39
1069 DO 1012 K=I,MCNT RMS011f1f0
1070 K 1 = MCITY(K) RMS011f1f1
1 071 1 = KCITYIKI) RMS011f'i2
1072 1 F I I. EQ. 0) GO TO 1012 RMS01'i1f3
1073 SUMNDII) = 0211 IIW( I ,20) RMS011f1f1f
-------
fO~..RAN .V G LEVEL
20
MAIN
071
)75
1076
1077
1078
1079
1080
108 I
1082
1083
108'i
1085
1086
1087
1088
1089
1090
1091
1092
1093
109'i
1095
1096
IU97
1098
1099
1 100
1 101
1 102
I 103
110'i
I 105
I 106
I 107
I 108
1 109
I 1 I 0
I I I 1
OA"::=72006
lOll I 108
)AG;: :)3
XNO(I) = DNR(I) . 5UMND(I) + ) ( )
XLO(I) = 1.0392'i . XNO(I)
XLL = XL (I) + XLO(I)
UD(I) =-(XNO(I)/XLL)
01(J) = 0.0
02(1) = 0.0
1012 CONTINUE
73'i8 FORMAT(3fI0.3,fl0.'i,33X,13)
IFC IPRN3)3S1 t3'iBt351
351 CONTINUE.
GO TO (3'i3,3'i'i,3'i'i,3'i3,3'i3,3'i'i,3'i'i,3'i'i),L5TRAT
3'i3 IF (IPRNI) 3'iB,3'i8,3'i'i
31f'i 00 3'i5 K = I,MCNT
KI = HCITY(K)
I = KCITY(KI)
CALL PRINI (I)
3lf5 CONTINUE
3'i8 IF INOPEN.EQ.I) GO TO 5001
ISTOP = 0
33'i8 CONTINUE
I5TOP = ISTOP + I
00 1013 K = l,MCNT
KI = MCITYIK)
I = KCITYIKI)
IFCISTOP.EQ.l) M = I
IFCISTOP.EQ.2)- M = 92
IF (I.EQ.O) GO TO 1013
PIN(M) = PIN(M) + PIO(I)
XIN(M ) = XIN(M) + XIOII)
XKNCM) = XKN(M) + XKOII)
XDNIM) = XDNCM) + XDCI)
YYN(M) = YYN(M) + YYO(I)
GNIM) = GN(M) + GOCI)
TNCM) = TNIM) + TDCI)
CNIM) = CNIM) + CDCI)
DNBARNIM) = DNBARNCH) + DNBAR II)
XNN(M) = XNNIM) + XNDII)
XLNIM) = XLNCH) + XLDCI)
RMS01115
RMSOl'i'ib
RMSOl'ilf7
RMSOI'i'i8
RMSOI'i'i9
RMSOI'iSO
RMSOl'i51
RMSOl'i52
RMSOI'iS3
RMSOl'i51f
RMSOl'i5S
RMSOl'i56
RMSOl'i57
RMSOI'i58
RMSOI'i59
RMSOI'ibO
RMSOl'ibl
RMSOllf62
RMSOllfb3
RMSOl'i6'i
RMSOl'ibS
RMSOllfb6
RMSOl'i67
RMSOl'i68
RMSOl'i69
RMSOl'i70
RMSOI'i71
RMS01'i72
RMSOI'i73
RMSOlIf7'i
RMSOllf75
RMSOllf76
RMSOI'i77
RMSOI'i78
RMSOl'i79
RMSOI'i80
RMSOllfBl
RMSOllfB2
-------
tORT~AN I~ G LEVEL
20
MAIN
1 1 1 2
1 1 1 3
1 1 1 'i
1 1 15
1 1 1 6
1 117
1 1 1 8
1 1 1 9
1 120
1 1 21
1 122
1 123
112'i
1 125
1 126
1 1 27
1 128
1 129
1 130
UNIMI = UNIMI + UDII)
sUMNNIMI = sUMNNIM) + SUMND. .
QBARNIMI = QBARNIM) + Q8ARDII)
QTNIMI = QTNIMI + aT DIll
QCNIMI = QCNIMI + QCUII)
QMNIM) = QMNIM) + QMDII)
1013 CONTINUE
GO TO 133'i8,3S01 , IS TOP
3S0 IF (IPRN3.EQ.0) GO TO 5001
IF (IPRN21 SOOI,SOOI.3'i6
3'i6 IF (LSTRAT-S) SOOI t3'i7,5001
3'i7 IKSW = 2 .
00 13'i7 K = I,MCNT
KI = MCITYIKI
I = KCITY(KI)
CALL PRIN1( II
13'f7 CONTINUE
IKSW = 0
GO TO SOOI
1 1 3 1
C
2000
C
C
C
C
C
CONTINUE
DA.~ = 720n6
10/'f"+/08
;)AG ~ 00'f
R~SOI'f83
RMS01'f8'f
RMS01"t8S
RMS01"+86
RMS01'f87
RMsOl"t88
RMsOl"t89
RMsOl'i90
RMS01'f91
RMS01'f92
RMS01'f93
RMS01"t9"t
RMSOl'f9S
RMSOI'f96
RMS01'197
RMS01'f98
RMS01'f99
RMS01S00
RMS01S01
RMS01S02
RMsOlS03
RMsOISO'l
RMsOISOS
RMsOIS06
RMS01S07
RMS01508
RMS01S09
RMS01S10
RMS01511
RMs01S12
RMS01513
RMS0151'f
RMS0151S
1 366816.).IBENEF/I0000.)RMS01S16
RMS01517
RMS01S18
RMS01S19
RMS01S20
CITY INFORMATION ( ALL-INO)
EFFECT
ON ECONOMY THROUGH CHANGE IN REGIONAL INCOME
1 132
1 133
113'i
1 1 35
1 1 36
1 137
1 138
1 1 39
1 1 'i 0
1 1 'I 1
11'12
1 1 "t 3
00 2001 K= 1 ,MCNT
KI = MCITY(K)
I = KCITYIKII
IF (I.EQ.O) GO TO 2001
IF ILSTRAT-S) 2009,2008,2009
2008 COSTIII,K) = -DELTCIK).QC(II
GO TO 377
2009 COSTI II ,K) =-911.3 . (IYYII) - XXII ,201)
377 YR(I) = (YYII)+COST(II,K)I/YYII)
TRII) = 1.0
ONRII) = 1.0
2001 CONTINUE
-------
tORT~AN
11'+'+
1 1'+ 5
11,+6
1 1 '+ 7
1 1'+ 8
1 1'+ 9
1 150
1 1 5 1
1 152
1 153
1 15'+
1 155
1 156
1 157
1 158
1 159
1 160
1 1 6 1
1 162
1 1 63
116'+
1 165
1 166
1 167
1 168
1 169
1 170
1 1 7 1
v G _::VEL
C
C
C
C
C
3000
3001
C
C
C
C
C
3002
2006
20
~ A ~
)A';: = 72006
10/'+'+/08
;JAG:: OOif
GO TO 3002
~IS01521
RMSOl522
RMS01523
RMSOI52'+
RMS01525
RMSOIS26
RMSOl527
RMSOl528
RMSOIS29
RMS01S30
RMS01531
RMSOIS32
RMSOIS33
RMSOI53'+
RMS0153S
RMSOIS36
RMS01S37
RMSOIS38
RMS01S39
RMSOlSifO
RMSOIS'+1
RMS015if2
RMS01Sif3
RMSOlS'i'+
RMS01S'f5
RMS01S'+6
RMS015'i7
RMSOIS'i8
RMSOIS'i9
RMSOISSO
RMSOl551
RMS015S2
RMSOISS3
RMSOl5Sif
RMSOISSS
RMS01SS6
RMSOl5S7
RMSOlSSB
CITY INFORMATION - ALL IND
EFFECT
ON ECONOMy THROUGH INCREASED TAXATION
CONTINUE
DO 3001 K=I.MCNT
KI = MCITYIK)
1 = KCITYIKI)
IF CI.EQ.O) GO TO 3001
TRCI) = ITCI) - COSTCII.K»/TII)
YRCI) = 1.0
DNRII) = 1.0
CONTINUE
THIRD SET OF SIMULTANEOUS EQUATIONS WHERE
Y,G, AND C ARE ENDOGENOUS
T AND X ARE EXOGENOUS
MARK = 3
lVAR = 3
MCASES = MCNT . lVAR
LCHECK = 0
CONTINUE
KPIST = 6
,.IM = 0
M = 0
DO 2002 K=I,MCNT
KI = MCITY'CK)
1 = KCITY'IKI)
IF CI.EQ.O) GO TO 2002
Y' I M + I) = YY I I )
YIM+2) = GII)
YIM + 3) = (II)
X I 1-1 M + 1) = T R C I )
XIMM+2) = YRII)
XIMM+3)= XXII,20)
-------
FORTRAN
1 172
1 173
1 17'+
1 175
1 1 7 b
1 177
1 178
1 179
1 180
1 181
1 182
1 183
118'+
1 185
1 186
1 187
1 188
1 189
1190
1 1 91
1 192
1 193
119'+
1 195
1 1 9 b
1 197
1 198
1 199
1200
1201
1202
1203
120'+
1205
V G _~VEL
2002
C
20
)A"E = 72006
~ A N
10/'+'+/( f
;)AG~ CO'f
X(MM+') =
X(MM+5) = FUDGE(I,I)
X(MM+b) = FUDGE(I,2)
M = M+3
MM = MM + b
CONTINUE
~MS)1559
RMS01560
RMS015bl
RMS01562
RMS01563
RMSOI56'f
RMS01565
RMS01566
RMS01567
RMS01568
RMS01S69
RMS01570
RMS01571
RMS01S72
RMSOlS73
RMS0157'i
RMS01575
RMS01576
RMS01577
RMS01578
Rt1S01579
RMS01580
RMS01581
RMS01582
RMS01583
RMSOI58'+
RMS01S85
RMS01586
RMS01S87
RMS01588
RMS01589
RMS01590
RMS01591
RMS01592
RMS01593
RMSOI59'f
RMS01595
RMS01596
IF (ICTAB) 8019,8019,1018
1018 DO 952 IM=I,LVAR
952 ARAB(IM) = SCC1(IM)
. C
8019
C
1 0 1 9
2003
200'f
C
2005
CALL SIMUL
M = 0
IF(LCHECK)2003,2003,2005
DO 200"1 K=I,MCNT
Kl = MCITY(K)
I = KCITV(KI)
IF (I.EQ.O) GO TO 200'f
YYOII) = Y(M+1)
GO (I) = YIM+2)
CD (I) = YIM+3)
TR(J) = 1.0
YR(I) = 1.0
M = M+3
COtH 1 NUE
LCHECK = 1
NBEFOR = NBEF03
GO TO 2006
M = 0
NEwCIT = 0
NBEF03 = 0
DO 2007 K=I,MCNT
KI = MCITYIK)
1 = KCITYIKI)
IF (1.[0.01 GO TO 2007
XDII) = .3'f0'f6.(YYO(I)-YIM+I»)
-------
FOR'RA~
206
1207
1208
1209
1 210
1 2 1 1
121 2
1213
121 ~
1 215
1 216
121 7
1 2 1 8
121 9
1220
122 1
1222
1223
122 ~
1225
1226
1227
1228
IV G L~V~L
2007
C
5001
SOOO
C
C
C
S003
2 )
MA . '
DAT~ = 72)06
M = M + 3
WBD(I) = ALPHAI201 .XOII)
PIDII) = XD(I) .ll.-BETA(20»
XID(I) = PICOEF(20) . PIDCI)
XKDII) = XIOII)/(1.+DI20»
D2(I) = WBDII)
CONTINUE
GO TO '1000
IF IINTT.EQ.l) NOK = NOK - I
CONTINUE
CONTINUE
IF INOK.EQ.O) IHELP = IHELP + 1
NOK = 1
IFIINPTI21.EQ.2) NOPEI = NOPEI - 1
ID = NOPEI + INPTll)
IFIID.EQ.O) GO TO 700
INTT = INPT(2)
GO TO 10
95 I = 1000000
IS = ICHOII)
SToP
600'1 FORMATIIH1,' INDUSTRY-'d3,' DOES NOT EXIST FOR CITY-'d3)
END
/
10/'1'1/08
)AGE :(I.L
RMSOlS97
RMS01S98
RMS01599
RMS01600
RMS01601
RMS01602
RMS01603
RMS01601.f
RMS0160S
RMS01606
RMS01607
RMS01608
RMS01609
RMS016I0
RMS01611
RMS01612
RMS01613
RMSOI61'1
RMS016lS
RMS01616
RMS01617
RMS01618
RMS01619
RMS01620
RMS01621
RMS01622
RMS01623
-------
fOR"RAN I V G - ~V ~L
0)) 1
0002
0003
000'"
0005
0006
0007
0008
0009
2 )
P~ E - '
SUBROU .'E PHELP(XX.X).'ERX,X
c
c
If(XX) 10,10,S
8 XOEP = XX 1 PERX
XIP = ( XD/XX) . 100.
RETURN
10 XOEP = 0.0
XIP = O~O
RETU~N
END
',XD~'
OA.~=72006
1 Of' . 10S--
~AbE UU01
RMS01621
RM~01625
RMS01626
RMS01627
RMS01628
RMS01629
RMS01630
RMS01631
RMS01632
RMS01633
RMS01631.f
-------
I - ~-~- --
FO~. ~AN
I V G _~VEL 2: PRINI )A"~ = 71)06 1 0 I 'i'f I 0 8- "AGE OUU
0)) SUBROJT NE 'R ~ 1 (I(ITYI RMS01635
( RM.SO 1636
( RMS01637
( RMS01638
0002 DIMENSION (911, (D(92I, (N(92) RMS01639
0003 DIr1ENSION (TOTD7(91I, (TOT67(91 ,ZOI, DNBAR(92) RMSOI6'i0
000'+ D I 11 ENS ION DNBARN 1 92 I, EOLDI9J ,20,SI, EO(91,5) RMS016'i1
0005 DIMENSION GI921, GD(921, GNI921 RMSOI6'i2
0006 DIMENSION IA(~CR(92,lO), IND(I 192,101, K(ITYIIOOI RMS016'i3
0007 DIMENSION LABELlI0,20I, MCITYI91 I, PDN(911 RMSOI6'i'f
0008 DIMENSION PJ192,201, PJD(921, PIN(921 RMSOI6'+S
0009 DIMENSION PPID(91), PXD(91), PXID(911 RMSOI6'i6
0010 DIMENSION PXKD(911, QBAR(92), QBARD(921 RMSOI6'i7
001 1 DIMENSION GBARNI91l, QC(92), QCD(92) RMS0161t8
0012 D I 1'1 ENS I ON QCN(92), QI1(92), QMD(92) RMSOI6'i9
0013 DIMENSION QMN(92), QT(92), QTD(921 RMS016S0
001'+ DIMENSION QTN(92), SUMNI92.20), SUMIW 1 92) RMSOl6S1
0015 DIMENSION SUMNN(92), T(92), TD(92) RMS016S2
0016 DIr1ENSION TN(92), U 1 92), UDI921 RMS016S3
0017 DIMENSION UNI921, WN(92), XO(92) RMSOI65'i
0018 DIMENSION XDN(92I, X I (92,20), XIDI921 RMS016SS
0019 DIMENSION XIN(92), XK(92,201, XKD(921 RMS01656
0020 [) I r1 ENS I ON I NP T ( Z I RMS016S7
0021 DIMENSION XKN(9ZI, XL(911, XLD(92) RMS01658
0022 DIMENSION XLN(92), XNBARI921, XNDI921 RMS016S9
0023 DIMENSION XNNI921, XNT(92). XXI92,ZO) RMS01660
002'+ DIMfNSION YY(921, YYD(92) t YYNI921 RMS01661
0025 DIr-'ENSION TITLEIZO,21 RMS01661
0026 DIMENSION POLD(91 ,ZO,SI ,PNEWDI20,SI RMS01663
( RMS01661t
0027 DIMENSION PERPINI921, PERTN (92 I, PERSUMI921 RMS01665
0028 DIMENSION PERXINI91l, PER(N I 92 I, PERQBR(91) RMS01666
0029 DIMENSION PERXKN(92), PERDNfJ (92), PERQTNI9Z) RMS01667
0030 DIMENSION PERXDN(92), PERXNNI9ZI, PERQCN(92) RMSOl668
0031 DIMENSION PERYYN(91), PERXLN(92), PERQMN(92J RMSOlbb9
0032 0 I r-1 ENS I ON PERGN(92), PERUNI921 RMS01670
C RMSOlb71
( RMSOlb72
-------
FOR "RAN IV G _~VEL
0033
C
C
003Li
0035
0036
0037
0038
0039
OOLiO
00/i1
00/i2
OOLi3
OD/i/i
DO/i5
00/i6
00/i7
00/i8
00Li9
0050
0051
0052
C
0053
005Li
0055
0056
0057
0058
0059
0060
0061
RMS01673
RMS0167'f
RMS01675
RMS01676
RMS01677
COMMON 1 CPRIN 1 PXO, PXID, PXKO, PPIO, CTOT07, pDN, CTOT67 RMS01678
1 ,EOLO.EO,POLO,PNEWD RMS01679
COMMON/BUG/ISBUG,LREAD,LWRIT RMS01680
COMMON IAPRIN/PI,PID,XI,XIO,XK,XKD,XX.XO,C,CD,YY,YYO,SUMN,SUMNO, RMS01681
1 XNBAR,DNBAR,XNT,XND.u.l1D,XL,XLD,G,GD,T,TO,QT RMSOl682
2 ,QTO,QC,QCD,QM,QMO,Q8AR,QBARO RMS01683
COMMON 1 ZEROT 1 PIN,XIN,XKN,XDN,YYN,GN,TN,CN,ONRARN,XNN,XLN,WN, RMSOI68/i
1 UN'SUMNN,Q8ARN,QMN,QCN,QTN RMSOl685
COMMON/NPRINI LABEL,MCITY,KCITY, LSTRAT'INPT,IKSW,IHELP,IFEEDRMSOI686
COMMON IPRITTI TITLE,BENEF RMS01687
COMMON 1 IPRIN IINTT,LPNTT ,JYEAR.IAQCR,IJCIT,MIN,IA312 RMS01688
EQUIVALENCE I INOCI 11) ,XXIII) RMS01689
I = ICITY RMS01690
~RITEILWRIT,100) RMS01691
LPYEAR = 1969 + JYEAR RMSOI692
IFIIClTY.EQ.92) GO TO 760 RMS01693
IFI IKSW-l) 1,750,750 RM50169/i
1 WRITE (LIf:RIT,110) LSTRAT,(LABEL(LSTRAT,J), J = 1,20) RMS01695
VIRITEILWRITtl20) IIAQCRII,J),J:: 1,20) ,LpYEAR RMS01696
WRITEILWRIT,130) RMS01697
WRITEILWRIT,l'iO) RMS01698
WRITEILwRIT,150) RMS01699
WRITE(LWRIT,160) RMS01700
RMS01701
RMS01702
RMS01703
RMS0170'f
RMS01705
RMS017U6
RMS01707
RMS01708
RMS01709
RMSOI710
20
COMMO'
1
2
CALL
CALL
CALL
CALL
CALL
CALL
CALL
CALL
(ALL
PR '
OAT~ = 72006
:tAGE : OO~
lC/11/08
IBLOX(11
PERP .N, :tERTN,PERSUM,PERXIN,pERO, :t~RQ3~,
'PERXKN,PERONB,PERQTN,PERXON,PERXNN,PERQCN,
PERYYN,PERXLN,PERQMN.PER~N,PERUN
PHELPIPIII,20).PIDII),PERPINII) ,PIP,pERP)
PHELPIXI(I,20),XIDII) ,PERXINII),XIP,pERXI)
PHELPIXKII,20),XKDIII,PERXKNII),XKP,pERXK)
PHELPeXX(I,20),XO{I),pERXDN(I) .XXP,PERXX)
PHELP ICII) ,(0(1) ,PERCNII) ,CP,PERC)
PHELP(YYII),yYD(I).PERYYNII) ,YYP,PERYY)
PHELP(SUMNII.20),SUMNO(I) ,PERSUMII) ,SUMNp,PERSN)
PHELPIXNBARII),ONBARII),PERDNBII),PNBAR.pERNB)
PHELPIXNTIl),XNO(I),PERXNN(I),XND,PERXN)
-------
fORTRA~
0062
0063
0061.!
0065
0066
0067
0068
0069
0070
0071
0072
0073
0071.!
0075
0076
0077
0078
0079
0080
0081
0082
0083
0081.!
0085
0086
0087
0088
0089
0090
0091
0092
0093
0091.f
0095
0096
.V G LEVEL
20
PRIN1
DATE = 72006
CALL PHELPIXLII),XlDII),PERXLNII),XlP,PERXL)
CAll PHELPIGII),GDII),PERGNII),GP,PERG)
CAll PHELPITII),TOII),PERTNII),TP,PERT)
CAll PHElPIQCII) ,QCDI I) ,PERQCNII),QCP,PERQC)
CAll PHELPIQMII),QMOII),PERQMNII),QMP,PERQM)
CALL PHELPIQBARII) ,QBARD( I) ,PERQBRII),QBARP,PERQB)
PERXN = PERSN + PERNB
PERQT = PERQC + PERQM + PERQB
PER XL = PERXN / .96
C
C
C
IfIUII») 58,58,56
56 UP = UD(I) / UII)
PERUX = 1.f.0
GO TO 59
58 PE.RUX = 't.0
UP = 0.0
5 9 C 0 ~JT I N U E
UX = U(I)-100.
UXD = UDII)-100.
~RITElli\RIT,190)
\"iR I TE I LI'IR IT, 170)
WRITEILvwRIT,180)
Vv R I TEl L VI R IT, 2 0 0 )
Vi R I TEl L VI R IT, 2 1 0 )
WRITEILWRIT,160)
WRITEILWRIT,2I.fO)
'IJ R I TEl l VI R IT, 1 6 0 )
WRITEILWRIT,260) PERNB,XNBARI ICITY) ,DNBARI ICITYI ,PNBAR
WRITEILWRIT,160)
WRITEILWRIT,300) PERYY,YYIICITY),YYOIICITY),YyP
WRITEILWRIT,290) PERC,CI ICITY! ,COI ICITY) ,cP
WRITEILWRIT,310) PERXN,XNTIICfTY),XNDIICITY),XNP
WRITEILWRIT,320) PERUX,UX,UXD,UP
WRITEILWRIT,330) PERXL,XL(ICITY),XLDI ICITY) ,XlP
WRITEILWRITdI.fO) GIICITY), GD(ICITYI, GP
\~ R I TEl L vi R IT, 3 5 0) PER T t T ( I C I T Y ) ,T D ( I C I T Y ) ,T P
WRITE IL~RIT,355)
PERXX,XX( 1(ITY,20) ,XDI ICITY) ,XXP
PERp ,PIIICITY,20),PIO(ICITY),PIP
PERXI ,XI I 1(ITY,20) ,XIDI I(ITY) ,XIP
PERXK,XK( ICITY,20) ,XKD( ICITY) ,XKP
PERSN,SUMNI lCITY,7.0) ,SUMND( 1(ITY) ,SUMNP
10/'t"/08
PAGE 000
~'S01711
~MS01712
RMS01713
RMS011l't
RMS0171S
RMS01716
RMS01717
RMS01718
RMSOI719
RMS01720
RMSOl721
RMS01722
RMS01723
RMS0172't
RMS01725
RMS01726
RMSOl727
RMS01728
RMS01729
RMSOl730
RMSOl731
RMSOl732
RMS01733
RMSOI73't
RMSOl735
RMSOl736
RMSOl737
RMSOl738
RMSOl739
RMSOI7't0
RMSOl71.fl
RMSOI7't2
RMS0171f3
RMS017't1f
RMSOl71.fS
RMSOl71f6
RMS0171f7
RMSOl71f8
-------
FOR.~M
20
)A"E = 720n6
'AGE (00
V G _~VEL
PRINI
10/11/08
0097
0098
0099
0100
0101
0102
:>..........01)
0103
010'f
0105
0106
0107
0108
0109
01 10
01 1 1
01 12
01 13
011'f
01 15
0116
01 1 7
01 18
:"..........01 )
01 19
0120
0121
WRITEI _w~IT,360)" :>~RQT,QT( IC.Y) ,QTD( IC.Y) ,Q"":>
VI R I TEl L W R IT, 2 2 0 'p E R Q M , Q M ( I C I T Y ) ,Q M D ( I C I " . Y ) ,Q M :>
\'1 R I TEl L Vi R IT, 2 7 0) PER Q B , Q BAR ( I C I T Y ) ,Q BAR D 1 I C I T y ) ,Q BAR P
W R I TEl L ~': R IT, 3 7 0) PER Q C , Q C ( I C I T Y ) , Q C D ( I C I T Y ) ,Q C P
100 FORMATIIHI)
110 fORMATI///IH ,'CONTROL STRATEGY NO~ 'II,' - ',20A'f,/)
$
IEYOO'fI COMMA[[[
120 FORMATI IH ,'AQCR ',I3,IX,19A'f,3X,' FOR ',I'f) RMSOl755
130 FORMATIIH ,7'fX,'VJITHOUT',7X,' WITH',7X,'NET',8X,'PERCENTI) RMS01756
1'f0 fORMAT(lH ,7'fX,'CONTROL',7X,'(T-l)',6X,'CHANGE',6X,'CHANGE') RMS01757
150 fORMATI IH ,'fX,'MANUfACTURING INDUSTRIES') RMS01758
160 fORMATIIH ) RMS01759
170 fORMATIIH ,6X,'PROFIT IMILLIONS)'.'f8X,IX,3(FIO.3,2X),FI0.'t) RMS01760
18U FORMATIIH ,6X,'INVESTMENT (MILLIONS)','f'fX,IX'3IFI0.3,2X), RMS01761
1801FI0.'f) RMSOl762
190 FORMAT(IH ,6X,'VALUE ADDED (MILLIONS)','f'fX,3IFI0.3,2X),fI0.'f) RMS01763
20U FORMATIIH ,6X,'CAPITAL STOCK (MILLIONS)','f2X'3IFI0.3,2X),FI0.'f) RMSOI76'f
210 FORMATIIH ,6X,'EMPLOYMENT (1000 S)t,'t6X,IXt3(FI0~3,lX),FI0.'f) RMS01765
210 FORMAT(lH ,6X,'ELECTRICITY USED BY MANUFACTURING INDUSTRIES (10 M RMSOl766
lKWH)',11X,3IF10.3,lX),FI0.'f) RMS01767
2't0 FORMATIIH ,'fX,'OTHER INDUSTRIES') RMS01768
26U FORMAT(lH ,6X,'EMPLOYMENT (IUOO S),,'t6X,IX,2(FI0~2,2X)tfl0.3,2X, RMS01769
1 FI0.'f) RMS01770
270 FORMAT(IH ,6X,'ELECTRICITY USED BY OTHER INDUSTRIES (10 M KWH)', RMS01771
1 19X,3(FI0.3,2X),FI0.'f) RMS01772
290 fORMAT( IH ,'tX,'REGIONAL CONSUMPTION (MILLIONS)',3'fX,IX,FIZ.3,2X, RMS01773
29012IFI0.3,2X),FI0.'t) RMSOI77'f
300 FORMATIIH ,'tX,'TOTAL PERSONAL INCOME FOR THE REGION IMILLIONS)' RMS01775
$
RMS01719
RMSOI75)
RMS01751
RMS01752
RMS01753
RMSOI75't
IEYOO'f! COMMA[[[
300119x,FI2.3,2X,2IFI0.3,2X) ,FI0.~) RMS01776
310 FORMATIIH ,'tX,'TOTAL REGIONAL EMPLOYMENT (1000 S)',33X,lX,3IFl0.3,RMSOI777
31012X) ,FI0.'f) . RMS01778
-------
FO~"~AN
o 22
0123
012~
0125
0126
0127
0128
0129
0130
01 31
0132
0133
013~
0135
0136
0137
0138
0139
01'fO
o 1 'f 1
01'!2
01'!3
OI~'I
01'f5
01'16
V G _:VEL
PRINI
)A"E = 72006
)AG~ JOl
20
10/'fLf108
:X)ENDITURE ~O~ THE REG ON (MILLIONS)', RMS01783
RMS0178Lf
REVENUE FROM THE REGION (MILLIONS)',22X,RMSOI78S
RMS01786
RMS01787
KWRMS01788
RMS01789
KWH)',RMS01790
RMS01791
MANUFACTURING INDUSTRIES ACTIVE IN THISRMS01792
RMS01793
RMS0179~
'i 51 FORMA T (1 H ,'[[[ RMSO 1 795
. . . . .
3'10 ~ORMAT(IH ,IX,'GOVERNMENT
3'10 9X,IX,3(FI J.3,2X),~10.Lf)
3S0 FORMAT(IH ,'iX,'GOVERNMENT
3Sul1X,3(FI0.3,2X),FI0.'i)
3S5 fORMAT (/,IH ,'iX,'ELECTRIC PO~ER DEMAND' ,II
360 FORMAT(IH ,6X,'TOTAL ELECTRIC CONSUMPTION FOR THE REGION (10 M
3601S)',13X,3(FI0.J,2X),FI0.'i)
370 FORMAT(IH ,6X,'RESIDENTIAL CONSUMPTION IN THE REGION (10 M
37UI17X,3(FI0.3,lX),FI0.Lf)
380 FORMAT(\H ,'SIC DETAIL FOH
3801 AQCR'/)
'100 FORMAT(IH ,'SIC DETAIL: NO DETAIL AVAILABLE FOR.THIS AQCR').
'iS2 FORMA T (1 H
IVESTMENT-
'iS3 FORMAT (lH
.
-VALUE ADDED- .
-CAPITAL STOCK- .'
. . . " . ...
, '. SIC................................................ RMSO 1 799
RMS01796
-INRMSOI797
RMS01798
, . . . . . . .
1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . )
'is'! FORMAT (IH ,'.CODE. NO CONTROL. . NET PERCENT. NO CONTROL.
1 NET. PERCENT. NO CONTROL. NET. PERCENT.' )
'i5S FORMAT (IH ,'...... (MILLIONS). CHANGE CHANGE (MILLIONS).
1 CHANGE CHANGE (MILLIONS). CHANGE. CHANGE .' )
'iS6 FORMAT (IH ,'...... -PROFIT-
IF EMPLOYEES- . -ELECTRIC POWER CONSUMPTION-
'iS7 FORMAT (IH ,'...... (MILLIONS). CHANGE. CHANGE
1 CHANGE. CHANGE. (10 M KWH). CHANGE. CHANGE
. .
If(INTT.EQ.2.AND.MIN.EQ.0.AND.IA311.EQ.I) GO TO ~61
IF(INTT.EQ.2.AND.MIN.EQ.0) GO TO ~60
IF(IHELP.EQ.11 GO TO SOD
RETURN
'i60 WRITE (LWRIT,lbO)
WRITE (lWRIT,16o)
WRITE (LWRIT,'iOO)
RETURN
'i61 WRITE(L~RIT,'i62)
'ibZ FORMAT(//IH ,'NO DETAILED SIC AVAILABLE FOR
1 COSTS USED INSTEAD'I' IMPACT TRANSFERED TO
RETURN
RMS018UO
RMSOl801
RMS01802
RMSOl803
RMSOI80'i
-NO. ORMSOIBOS
RMS01806
. RMS01807
RMSOIB08
RMS01809
RMS01810
RMSOIBll
RMS01812
RMSOIB13
RMS01811t
RMS0181S
RMS01816
RMS01817
THIS AGCR - TOTAL MANURMS01818
-------
fOR'"RAN IV G LEVEL
PR IJ
)A'"~ = 71006
20
o j 7
o '18
01'19
0150
o 1 5 1
0152
0153
015'1
0155
0156
0157
0158
0159
0160
o 1 6 1
0162
0163
016'1
0165
0166
0167
0168
sou WRITE (LWR .,100)
W R I T E (L W R I '", 60)"
~RITE (LWRIT,1.f51)
~RITE (LWRIT,1.f51)
WRITE (LWRIT,1.f53)
WRITE (LWRIT,'51.f)
WRITE (LWRIT,'55)
~RITE (LWRIT,1.f51)
530 FORMAT(IH ,IH.,13,lH .,3(IX,f9.3,lH .,F9.3,lH
DO 600 JND = 1,19
IF (INDCI(I,JND» 510,600,510
510 INDI = JND+20-(I/JND)
XXP =(PXD(JNO) /XX(I ,JND) ).100.
IF (XI(!,JND» 515,515,520
515 XIP = 0.0
GO TO 525
520 XIP = (PXIO(JND) /XI(I,JND».100.
525 IF (XK(I,JND» 527,527,535
527 XKP = 0.0
GO TO 5'10
53S XKp = (PXKD(JND) /XK(I,JND».100.
5'10 WRITE (L'/'JRIT,530) INDI ,XXI I ,JND) ,PXD(JND)
IPXID(JND) ,XIP,XK(I,JND) ,PXKD(JND) ,XKp
bOO CONTINUE
WRITE(LWRIT,100)
i"w' Ft I T E ( L W R I T tl 60 )
Vi R I T E ( L W R IT, 'I 5 1 )
IVR! TE (LWR IT, '156)
WRITE(LWRIT,1.f53)
WRITE(LWRIT,'SI.f)
WRITE(LWRIT,'57)
WRITE(lWRIT,1.f51)
DO 700 JND = 1,19
IF (INDCI(I,JND» 610,700,610
610 INDI = JND+20-(I/JND)
IF (PI(J,JND» 615,615,620
615 PIP = 0.0
GO TO 625
0169
0170
017 1
0172
0173
017'1
0175
0176
0177
0178
0179
0180
0181
0182
0183
.,F9.I.f,2H ~»
,XXP,XI (I ,JND),
10/'1'1/08
PAGE 000
RMSOIB2
RMS01822
RMS01823
RMS018211
RMS01825
RMSOIB26
RMS01827
RMSOIB28
RMS01829
RMS01830
RMS01831
RMSOIB32
RMS01833
RMSOI83't
RMS01835
RMS01836
RMS01837
RMS01838
RMSOIB39
RMSOIB'tO
RMSOIBl.fl
RMS018112
RMS018113
RMSOIB'f'l
RMSOIB'fS
RMSOI8't6
RMS0181.f7
RMS018'18
RMSOIB'f9
RMSOIBSO
RMSOIBSI
RMS01852
RMS01853
RMSOIB5't
RMSOIB55
RMS01856
Rr1S01B57
RMSOIB58
-------
FORTRAN IV G LEVEL
20
PRINI
)A'"E = 7Z006
10/~~/08
:JAG:: (00
018~
0185
0186
0187
0188
0189
019U
0191
0192
0193
620 PIP = PPID(JND) /PI(I,JND)1HOJ. RMS01859
625 IF (EO-D(I,JND,5t) 635,630,635 RMS01860
630 PCT = 0.0 RMS01861
GO TO 6~0 RMS01862
635 PCT = IEDIJND,5)/EOLDII,JND,5».100. RMS01863
6'1U If (SUMN(I,JND» 6'15,6'15,650 RMS0186'1
6~S SUMNP = 0.0 RMS01865
GO TO 655 RMS01866
65U SUMNP = (PDN(JND) /SUMN(I,JND» *100: RMS01867
655 V'IRITE ILWRIT,530) INDI,PI(I,JND),PPIDIJND),PIp,SUMNII,JND), RMS01868
1 PDN(JND) ,SUMNP,EOLDI I ,JND,S) ,EDIJND,S) ,PCT RMS01869
700 CONTINUE RMS01870
IF (INTT.NE.Z.OR.LPNTT.NE.Z) RETURN RMS01871
WRITE (LWRIT,100) RMS01872
WRITE (LWRIT,160) RMS01873
WRITE (LWRIT,Lf51) RMS0187~
WRITE ILWRIT,1000) RMS0187S
WRITE (LWRIT,'1S3) RMS01876
~RITE (LWRIT,'1S'1) RMS01877
WRITE (LWRIT,1015) RMS01878
WRITE ILWRIT,'151) RMS01879
DO 712 JND = 1,19 RMS01880
If (INDCI I 1 ,JND» 701,712,701 RMS01881
701 INDI = JND+20-(I/JND) RMS01882
IF (POLD( I ,JND,I» 703,702,703 RMS01883
702 PCT1 = 0.0 RMSOI88~
GO TO 70~ RMSOl885
703 PCTI = (PNEWDIJNDtl>!POLDII,JNDtl».100. RMSOl886
70'1 IF (EOLDII,JND,l» 706,705,706 RMSOl887
705 PCT2 = 0.0 RMS01888
GO TO 707 RMSOl889
706 PCT2 = IEDIJND,l)/EOLD(I,JNDtl»)*100. RMS01890
707 IF IPOLDII,JND,J» 709,708,709 RMS01891
70~ PCT3 = 0.0 RMS01892
GO TO 710 RMS01893
709 PCT3 = IPNEWD(JND,3)/POLDII,JND,3».100. RMSOI89'1
710 WRITE(LWRIT,530) INOI ,POLD( I ,JND, 1) ,PNEWDIJND, 1) ,pCTl ,EOLD( I .JND, 1 )RMSOI895
I.EDIJND,I) ,PCTZ,POLD( I ,JND.3) ,PNniDIJND,J) ,PCT3 ~ RMS01896
019'1
0195
0196
0197
0198
0199
0200
0201
0202
0203
020Lf
0205
0206
0207
0208
0209
0210
021 1
0212
0213
021Lf
0215
0216
0217
0218
0219
-------
FORTRAN
0220
0221
0222
0223
022Lf
0225
0226
0227
0228
0229
0230
0231
0232
0233
023Lf
0235
0236
0237
0238
0239
02LfO
02 'i 1
OZLfZ
02Lf3
OZLfLf
OZLfS
OZLf6
02Lf7
02Lf8
02Lf9
02S0
0251
0252
0253
025Lf
025S
0256
V G LEVEL
)A.~ = 72006
10/LfCf/08
20
PR N
712 COfHINUE
V, RITE (U'iRITtlOOf
WRITE (LWRIT,160)
WRITE (LWRIT,Lf51)
WRITE (LWRIT,1010)
WRITE (LWRIT,Lf53)
WRITE (LWRIT,LfSLf)
WnITE ILWRIT,101S)
WkITE (LWRIT,Lf51)
DO 72Lf JNO = 1,19
IF IINDCI(I,JNO» 713,72Lf,713
713 INDI = JND+20-ll/JND)
IF IEOLDII,JND,3» 71S, 71Lf, 715
71Lf PCTI = 0.0
GO TO 716
715 PCT1 = (EO(JND,3)/EOLD(I,JND,3»-100:
716 IF IPOLD(I,JND,Lf» 718,717,718
717 PCT2 = 0.0
GO TO 719
718 PCT2 = IPNEWDIJND,Lf)/POLOII,JND,Lf» -100.
719 IF IEOLDII,JND,Lf» 721,720,721
720 PCT3 = 0.0
GO TO 722
721 PCT3 = IEDIJND,Lf)/EOLD(I,JND,Lf) )-100:
722 WRITE (LWRIT,530) INDI ,EOLDI I ,JND,3) ,ED(JNDt3) ,PCT1 ,PaLD( I ,JNO,Cf),
1 PNEWD(JND,Lf) ,PCT2,EOLD( I ,JND,Lf) ,EO(JND,/.f) ,PCT3
72Lf CONTINUE
WRITE (LWRIT,100)
WRITE (LWRIT,160)
WRITE (LWRIT,Lf51)
WRITE ILWRIT,1020)
WRITE (LWRIT,1.f53)
WRITE ILWRIT,Lf5Lf)
WRITE (L'i/RITtl02S)
WRITE (LNRIT,LfSl)
DO 736 JND = 1,19
IF (INDCIII,JND» 725,736,72S
725 INDI = JNO+20-11/JND)
)AGE 000,r
RMS01897
RMS01898
RMS01899
RMS01900
RMS01901
RMS01902
RMS01903
RMS0190/.f
RMS01905
RMS01906
RMS01907
RMS01908
RMS01909
RMS01910
RMS01911
RMS01912
RMS01913
RMS0191'i
RMS01915
RMS01916
RMS01917
RMS01918
RMS01919
RMS01920
RMS01921
RMS01922
RMS01923
RMS0192Lf
RMS01925
RMS01926
RMS01927
RMS01928
RMS01929
RMS01930
RMS01931
RMS01932
RMS01933
RMS0193'i
-------
FO ~'RA~
0257
0258
0259
0260
0261
0262
0263
026'1
0265
0266
0267
0268
0269
0270
0271
0272
0273
02711
0275
0276
0277
0278
0279
0280
0281
0282
0283
028Lj
0285
0286
0287
0288
0289
0290
0291
IV G LEVEL
20
DAT~ = 72006
PR N
10/£f'UOB
~ I::JO-DII,JND,S)
726 PCTI = 0.0
GO TO 728
727 PCTI = IPNEWOIJND,5)/POLDII,JND,S».100.
728 IF IPOLDII,JND,2» 730,729,730
72Y PCT2 = 0.0
GO TO 731
730 PCT2 = IPNEWOIJND,2)/POLDII,JND,2».100.
731 IF (EOLDII,JND,2» 733,732,733
732 PCT3 = 0.0
GO TO 73Lj
733 PCT3 = IEDIJND,2)/EOLDII,JND,2».100.
73'1 WRITE ILWRIT,S30) INDI,POLDI I ,JND,S) ,PNEWDIJND,S) ,PCTI ,POLDI I,JND,
1 2) ,PNEWDIJND,2) ,PCT2,EOLDll,JNO,2) ,EDIJND,7) ,PCT3
736 CONTINUE
RETURN
750 GO TO 1760,770) tIKSW
760 BENT:: BENEF/IOOO.
WRITEILWRIT,911} IITITLEIIMJ,JM1), 1M 1 :: 1,20), JMl :: t,2},BENT
WRITE ILWRIT,900)
GO TO 790
770 WRITE ILWRIT,910)
790 WRITE ILWRIT,120) IJAQCRII,J),J=l,20) ,LPyEAR
WRITE
-------
FO~T~AN IV G ~~VEL
0292
)293
029'f
0295
0296
0297
0298
0299
0300
0301
0302
0303
030'f
0305
0306
0307
0308
0309
0310
031 1
0312
0313
031'f
0315
0316
0317
0318
0319
0320
0321
0322
0323
032'+
0325
20
PRINI
)A"E = 72006
CALL PHELP(XL( I) .XLN( I) .PERX-~ (I) ,X_;).p~~XLJ
CALL PHELP(G(I).GN(I) ,PERGN(I ,GP,PERG)
CALL PHELP(T(I) ,TN(I),PERTN(I) ,TP,PERT)
CALL PHELP(QTII),QTN(I) ,PERQTN( I),QTP,PERQT)
CALL PHELP(QCII) ,QCNII),PERQCNI I) ,QCP,PERQC)
CALL PHELPIQM(I).QMNII),PERQMNI I) ,QMP,PERQM)
CALL PHELPIQBARI I) ,QSARNI I) .PERQBRI I) .QBARP.PERQB)
PERXN = PERSN + PERNB
PERQT = PERQC + PERQM + PERQB
PERXL = PERXN / .96
C
C
85~ IF CUll» 858,858,856
856 UP ~ IUN(I) / U(I» . 100.
GO T0860
858 UP = 0.0
860 CONTINUE
PERUN(I) = .O'i . (I. + (UP/IOO.»
PERUX = '+.0
UX = UII).IOO.
UXO = UN(I).IOO.
C
C
WRITEILWRIT,190)
~J R I T E ( L IV R IT, 1 70 )
'iiJ.(ITE(LWRIT,180)
'lJRITE(LWRIT,200)
WRITEILWRIT,210)
WRITE(LWRIT,160)
WRITEILWRIT,2'f0)
WRITEILWRIT,160)
W R I T E ( L Ifi R IT, 2 6 0 )
WRITECLWRIT,160)
WRITE(LWRIT,300)
'lJRITE(LWRIT,290)
W R I TEl L \'J R IT, 3 1 0 )
WRITEILWRIT,320)
rVRITEILWRIT,330)
PERXX,XX( ICITY,20) dON( ICITY) ,XXP
PERp ,PIIICITY,20),PIN(ICITY>,PIP
PERXI ,XI I ICITY,20) ,XIN( ICITY> ,XIP
PERXK,XK( ICITY,20) ,XKN( ICITY> ,XKP
PERSN,SUMN( ICITY,20) tSU~NN( ICITY) .SUMNP
PERNB,XNBAR( ICITY> ,DNBARNI ICITY> .PNBAR
PERyy,yy( ICITY> .YYN( ICITY) ,yyP
PERC,CC ICtTY) ,CNI lelTY) ,eP
PERXN,XNT< IeITY) ,XNN( lelTY) ,xNP
PERUx,UX,UXD,UP
PERXL,XLC ICITY) .XLN( ICITY) ,XLP
10/1 'f/08
;)AG~ 0011
~MS01973
RMSOI97'f
RMS01975
RMS01976
RMS01977
RMS01978
RMS01979
RMS01980
RMS01981
RMS01982
RMS01983
RMSOI98'f
RMS01985
RMS01986
RMS01987
RMS01988
RMS01989
RMS01990
RMS01991
RMS01992
RMS01993
RMSOI99'i
RMS01995
RMS01996
RMS01997
RMS01998
RMS01999
RMS02000
RMS02001
RMS02002
RMS02003
RMS0200'f
RMS02005
RMSOZ006
RMS02007
RMSOZ008
RMS02009
RMS02010
-------
fa n HN
V G _:VEL
20
PRIN1
DATE = 72006
10/1f1f/01
)AGE 001
C
0326
0327
0328
0329
0330
0331
C
C
0332
0333
033'"1
0335
0336
0337
0338
0339
03'"10
03'"11
03'i2
~RITEILWRIT,3'i0) GIICITY), G'I CITY), GP
VI R I TEl L W R IT, 3 5 0) . PER T , T 1 I C I . . Y ) , T NIl C I T Y ) , T P
WRITE ILWRIT,35S)
WRITE(LWRIT,360) PERQT,QTI ICITY) ,QTNI ICITY) ,QTP
W R I TEl L \'J R IT, 2 2 0) PER Q,-1 , Q M ( I C I T Y ) ,Q M NIl C I T Y ) ,Q M P
WI
-------
FORTRAN IV G LEVEL
0001
0002
0003
000'+
ODDS
0006
0007
0008
0009
20
ABC)
)AT~ = 72000
10/IPfI:a
)AGE 000
SUBROU" NE A~CD
COMMO~ /BUG/ ISBUG,-READ,LWRIT
COMMON /BUGG/ IPRN1,IPRN2,IPRN3
IPRNI = 1
IPRN2 = 1
IPRN3 = 1
ISBUG = 1
RETURN
END
~IS02J38
RMS02039
RMS020't0
RMS020'tl
RMS020't2
RMS020't3
RMS020't't
RMS020'fS
RMS020't6
-------
Fon~AN
000
0002
0003
ODOlf
0005
0006
0007
0008
0009
DOlO
001 1
0012
n013
001Lt
0015
V G .:VEL
DATE = 72)(6
10/" 108
)AGE : JO
20
VA_JE
SUBROU. .Nf VALU::
CO t1 M 0 N /5 I M / M 1 ,M"2 , M CAS E 5 , Y ( 1000) , X ( 1 000 )
1 ,MARK,LVAR
GO TO(100,200,300,LtOo,SOO) , MARK
100 CALL VALUEI
GO TO 1000
200 CALL VALUE2
GO TO 1000
300 CALL VALUE3
GO TO 1000
LfOO CALL VALUELf
GO TO 1000
500 CALL VALUES
100U CONTINUE
RETURN
END
RMS020' 7
RriS020't8
RMS020Lf9
RMS020S0
RMS020S1
RMS020S2
RMS020S3
RMS020SLt
RMS020SS
RMS02056
RMS020S7
RMS020S8
RMS020S9
RMS02060
RMS02061
RMS02062
-------
'FORTRAN
00: 1
0002
0003
ODOLf
0005
0006
0007
0008
0009
0010
0011
0012
V G _:VEL
20
VA_UE
DA"E = 72006
O/If'i/OB
5JBROU' I!E VALUE
c
COMMON 151MI
1 ,MAHK,LVAR
COr'1MON IF A T I
COMMON I BUG
Ml,M2,MCASE5,Y(1000),X(1000)
FI 1000) ,LSTORE,KPIST
I 15BUG,LREAD,LWRIT
C
J = L5TORE + 1
DO 1 I=Ml,M2,LVAR
FII) = -YII) + .Lf2B36-X(J) +
F(I+I) = -YII+l) + 23.Lf227+
1 + XIJ+2)
FII+2) =-Y(I+2) + 137.9922 + .6132S6.Y(I)
1 + X I J + 3 )
1 CONTINUE
RETUHN
END
2:83038.yll+l) + .91f812.Y(I+2)
.91f21S.XIJ+l)
Pl' G E""""U'C"O 1
RMS:2063
R'S02061f
RMS02065
RMS02066
RMS02067
RM502068
RM502069
RMS02070
RMS02071
RMS02072
RMS020'13
RMS0207'i
RM502075
RMS02076
RMS02077
RMS02078
RMS02079
-------
FORTRAN
0001
0002
0003
000£i
0005
0006
0007
0008
0009
0010
no 1 1
0012
0013
001ii
V G _EVEL
20
VA_UE2
DA'E = 72006
5UBROU" NE VA_UE2
C
COMMON 15JMI
1 ,MI\RK,LVAR
COMMON IFATI
C 0 to1 M 0 N I F3 U G
Ml ,M2,MCA5E5,Y( 1000) ,XI 1000)
F( 1000) ,LSTORE,KPI5T
I ISaUG,LREAD,LWRIT
C
J = L5TORE
DO 1 I=Ml,M2,LVAR
F(l) =-y(J) + X(J+q).YII+£i) + X(J+l)
FIl+l) = -Y(I+l)-13.95R75 + 1.0392£iS.Y(I)+ 36t.37£i.Y(I+3)
1 +XIJ+7I
F(I+2) =-Y(I+2) + 1.8£i02 + ( .8377.XIJ+3» +(.1116.Y(I+3»
1 + X(J+6)
F(l+) =-Y(I+3) + «YII+l)-Y(I)IIYlI+lI)
F(l+q) =-YII+£i) + (X(J+2)/Y(I+2»
1 CONTINUE
RETURN
END
0/' '/08
~ ~ G E lJtiO 1
RM50208:
RM5020B
RM5020B2
RM5020B3
RMS020SI.f
RM5020S5
RM5020S6
RMS020S7
RM502088
RMS020S9
RMS02090
RMS02091
RM502092
RMS02093
RMS0209£i
RM502095
RMS02096
RMS02097
RM502098
-------
FOR "~AN
00 Jl
0002
0003
OOOLi
0005
0006
0007
0008
0009
0010
001 1
0012
V G LEVEL
20
5UBWU"IIJE VA_UE3
c
V A _UE:
DATE = 72006
Ml ,M2,MCA5ES,Y( 1000) ,XI lOaD)
COMMON 15IMI
1 ,MARK,LVAR
C 0 r1 MaN I FAT I
COMMON I BUG
F( 1000) .L5TORE
I I5BUG.LREAD.LWRIT
c
C
J = L5TORE
DO 1 I=Ml.t12.LVAR
F(I) = -Y(I) + .~2836~X(J+3)+ 2:8303S.Y(I+1) +
F(I+l) = -Y(J+1) + 23.~227+ .9~21S.X( 1+3).X(J+1)
1 + X (J+ ~ )
F(I+2) =-Y(I+2) +
1 + X(J+s)
1 CONTINUE
RETURN
END
137.9922 + .6132S6.Y(I) .X(J+2)
10/~£i/08
.9li81Z.Y(I+2)
;)AG~ 0)01
R~S02099
RMS02100
RMS02101
RMS02102
RMS02103
RMS0210~
RMS0210S
RMS02106
RMS02107
RMS0210B
RMS02109
RMS02110
RMS02111
RM502112
RMS02113
RMS0211£i
RMS0211S
RMS02116
-------
FORTRAN
0001
0002
0003
ODOlf
0005
0006
0007
0008
0009
0010
001 1
0012
0013
0011f
0015
0016
0017
0018
0019
0020
0021
0022
0023
002'+
0025
0026
0027
0028
0029
0030
0031
0032
0033
0031f
V G LEVEL
2(
VA_UEI
DATE = 72(06
Sl ~ROU' NE VA_UEIf
DIMENSION AlIf)
c
COMMON /SIM/ M1 ,M2,MCA5ES,Y( 1000) ,X(1000) ,MARK,LVAR,MKEEP
COMMON /FAT/ FI 1000) ,LSTORE
MM = LSTORE + 1
MT = 5 - r1KEEP
MTI = 2-MT
KD = MKEEP - 1
KDl = KO-7. + 1
DO 20 I=Ml,M2,LVAR
11 = I-MT
MR = MM - MTI
SUM = O~O
6 = ALOG(YIII+8» - ALOG(X(MR+llf»
DO 10 KK=I,KO
K = KK - 1
SUM = SUM + X(MM+K).Y(J+K)
A(KK) = 1.0
IF(X(MM+9).XIMM+K).EQ.O.O) GO TO 9
A(KK) = (XIMM+KD)/X(MM+K».IXIMM+K+MKEEP)/XIMM+KDI»
9 CONTINUE
B = B - X(MM+K+MKEEP) - ALOG(A(KK»
10 FII+K) = -Y(I+K) + AIKK)-y(II+If)
fIII+6) = -y(II+6) + SUM + XIMM+KD)-YII+KD)
M~1 = MM - MT I
FIII+If) =-YIII+If) + EXPIB)
FIJI+5) = -YIII+S) + IXIMM+I0)+Y(II+I0»/II.0+X(MM+II»
fIII+7) = -YIII+7) + Y(II+S)-X(MM+16)
F(tI+8) =-Y(II+8)+XIMM+lS).Y(II+7)
flll+9) = -YIII+9) + X(M~+12).Y(II+7)
fIII+I0) = -Y(II+I0) + XIMM+17)+ X(MM+13).V(It+9)
2U CONTINUE
RETURN
END
10/11/08
)AGE 0001
R~S02117
R ~ .5 0 2 1 I 8
RMS02119
RMS02120
RMS02121
RMS02122
RMS02123
RMS0212'i
RMS0212S
RMS02126
RMSOZ127
Rf1S02128
RMS02129
RMSOZ130
RMS02131
RMSOZ132
RMS02133
RMS02131f
RMS0213S
RMS02136
RMS02137
RMS02138
RMS02139
RMS0211f0
RMS0211f1
RMS0211f2
RMS021lf3
RMS021lfLf
RMS021lfS
RMS021Lf6
RMS021Lf7
RM50211f8
RMS021lf9
RMS021S0
RMS021S1
-------
FOR'RAN
000
r 002
0003
000'+
0005
0006
V G Lt.Vt.L
£0
VALU::5
SUBROUTINE VALUES
COMMON 15IMI 11,M2,MCASES,Y(
1 ,MA~K,LVAR
COMMON IF AT I FCIOOO),LSTORE
IT = 0
RETURN
END
DA.~=72006
10/'+'+/08
:tAGE )JOI
000),X(10CO
RMS021S2
RMS021S3
RMS021S'+
RMS0215S
RMS02156
RMSOZlS7
RMS021S8
-------
rOW "RAN
000
0002
0003
ODOLf
0005
0006
0007
0008
0009
0010
001 1
0012
.V G -::VEL
~o
D::8 JG
)AT:: = 72)06
lC/44/08
:tAG:: : 00
SUBROUTINE )EBUG(Y,"OT._VA~)
DIMENSION Y( 000)
WRITEI6,ZI)
MCASES = LVAR.MCNT
K = U
DO 10 I=I.MCASES,LVAR
K = K + LVAR
10 WRITEI6,20) IYIJ) I J=IIK)
20 FORMATIIH 110Fll.l)
21 FORMATIIHI)
RETURN
END
RMS02159
RMS02160
RMS02161
RMS02162
RMS02163
RMS0216'f
RMS02165
RMS02166
RMS02167
RMS02168
RMS02169
RMS02170
-------
I
I FOn~AN
V G _::VEL
::lAG:: 000
000
0002
0003
0001i
0005
0006
0007
0008
0009
0010
001 1
0012
0013
0011i
0015
0016
0017
0018
0019
C
0020
0021
0022
0023
0021i
0025
0026
0027
0028
0029
20
)A"E = 72(06
S H U-
10/4 4/0~
SUAROU' NE S MUL
c
C
C
~"S0217
R~S02172
RMS02173
RMS0217&f
RMS02175
RMS02176
RMS02177
RMS02178
RMS02179
RMS02180
RMS02181
RMS02182
RMS02183
RMS0218'i
RMS02185
RMS02186
RMS02187
RMS02188
RMS02189
RMS02190
RMS02191
RMS02192
RMS02193
RMS0219'i
IX = }CTAB + LCHECK + NBEFOR RMS0219S
IFIIX-31 69,60,69 RMS02196
60 IFINEWCITI 71,70,71 RMS02197
71 I = IJCIT RMS02198
WRITEILPR,1201 I IAQCRI I ,JI, J=I ,20) ,LPYEAR RMS02199
1 2 0 FORMAT I 1 HI, 111125 X , ' A Q C R ' tI 3 d X tl 9 A 'i ,5 X , , FOR' .I 'i'112 5 X , R M S 0 2 200
1 'CONVERGENCE PROCESS FOR ESTIMATED VALUES WITHOIJT POLLUTION CONTRRMS02201
10LS'1 RMS02202
70 CONTINUE RMS02203
WRITEILPR,50) MARK RMS0220'i
IFIMARK.EQ.2) GO TO 900 RMS02205
50 FORMATI 1H ,11'iX,'.................................................RMS02206
1[[[~..RMS02207
2'/TS,'.',TI20,'.'/T5,'.',2SX,'SIMULTANEOUS SYSTEM NUMBER - ',Il~ RMS02208
MACRO MODEL SIMULATION PROGRAM I - MAINLINE
DIMENSION DELCI000)
DIMENSION IAQCRI92,20)
COMMON ISIMI Ml ,M2,MCASES,yl 1000) ,XI 1000)
1 ,MARK,LVAR
COMMON IFATI FI1000),LSTORE,KPIST,ARAB(10),LCHECK
COMMON ICAT/ AI20,ZO)
CUMMON ISENS/ 1CTAB,NBEFOR,NEWCIT
COMMON I IPRIN IINTT,LPNTT ,JYEAR,IAQCR,IJCIT
COMMON I BUG I IS8UG,KRD,LPR
LSrORE = 0
KRFP = 0
106 FORMAT(ILf)
MCNT = MCASES/LVAR
LMN = LVAR
-------
20
)AT~ = 72006
10/11/(8
FonRAN IV G _~VEL
0030
0031
0032
0033
003'+
003S
0036
0037
n038
0039
00,+0
00,+1
00't2
00,+3
00'+'+
00,+5
00,+6
00,+7
00't8
00't9
0050
OOSI
00S2
00S3
005'+
00S5
00S6
0057
00S8
0059
0060
S P l -
3 " 20,'.'/'.S,'.',"
20,'.')
C
2 WRITECLPR,9) (ARABCJ),J=I,LVAR)
~RITECLPR,10) (Y(J), J=Ml,M2)
9 FORMATIIH ,T5,'.',SX,'ENDOGENOUS VAR.
1 TI20,'-')
10 FORMATIIH ,TS,'.',8X,'ACTUAL VALUE ',3IF9.2,sX),T120,'.')
11 FORMATCIH ,TS,'.',8X,'ITERATION-',I2,2X,3CF9.2,FB.S),T120,'.')
GO TO 69
',3CSX,AIf,'
CONV'),
C
90U CONTINUE
WRITECLPR,99)(ARABIJ),J=I,LVAR)
WRITECLPR,910) CYCJ), J=r11,M2)
99 FORMATCIH ,TS,'.',5X,'ENDOGENOUS VAR.
1 T 1 20, '.' )
91U FORMATC\H ,T5,'.',BX,'ACTUAL VALUE '.5CF9.2,8X),TIZ0,'.')
911 FORMATCIH ,T5,'.',BXt'tTERATION-"I2,2X,S(F9.2~F8.S),TI20,'.')
',SCSX.AIf.'
CONV'),
c
C
69 IF IIT-IMAX) 5,690,690
5 11 = IT + 1
3 CALL GRADN
CALL VERT
CALL VALUE
COMPUTE NEW VALUES AND TEST FOR CONVERGENCE
7 1AS = IT
IT = HIAX
NN = 0
00 If I = Ml,M2
NN = NN + 1
Y Y = Y C I )
YII) = 0.0
KK = 0
00 '+0 J= Ml,M2
KK = KK + 1
'+0 YCI) = YCI) - ACNN,KK).FIJ)
YII) = yY+YCI)
YOIS = AASCYll)/YY)
YOIS = ABSCYOIS-l.)
:lAG~ :00
RMS022J9
RMS02210
RMS02211
RMS02212
RMS02213
RMS0221'+
RMS02215
RMS02216
RMS02217
RMS0221B
RMS02219
RMS02220
RMSfJ2221
RMS02222
RMS02223
RMS0222'i
RMS0222S
RMS02226
RMS02227
RMS02228
RMS02229
RMS02230
RMS02231
RMS02232
RMS02233
RMS0223'+
RMS02235
RMS02236
RMS02237
RMS02238
RMS02239
RMS022'+0
RMS022'+1
RMS022'+l
RMSOl2't3
RMS022'+'+
RMSOl2'tS
RMS022't6
-------
FORTRAN
0061
0062
0063
006'1
0065
0066
0067
0068
0069
0070
0071
0072
0073
007'1
0075
0076
0077
0078
0079
0080
0081
0082
V G _::VEL
20
SIMUL
OAT:: : 72006
10/'+'1/08
PAG:: 000
DEL (I) : Y)IS RMSJ22'f7
IF(YDIS.LT.O.009}GO TO '12 RMS022'+8
IT : IAS RMS022q9
'12 CONTINUE RMS022S0
'I CONTINUE RMS022S1
IF(IX.NE.3) GO TO 69 RMS022S2
688 CONTINUE RMS022S3
IFCMARK.EQ.2) GO TO 673 RMS0225q
~'IRITE(LPRtll) IAS,cY(J),DEL(J), J=Ml,M2) RMS022S5
uO TO 69 RMS02256
673 WRITE(LPR,911) IAS,cY(J) ,DEL(J), J=MI,M2) RMS022S7
689 GO TO 69 RMS022S8
69U IF(IX.NE.3) GO TO 692 RMS022S9
691 WRITE (LPR,105) IAS RMS02260
105 FORMATCIH ,T5,'.',TI20".',/T5,'.',qX,'ITERATIONS TO CONVERGE:', RMS02261
1 12,TI20,'.',/'IX,'................................................RMS02262
2.....................$...........$................................RMS02263
'I..') RMS0226'1
692 CONTINUE RMS02265
If (M2 .EQ. MCASES) RETURN RMS02266
LSTORE = LSTOkE + KPIST RMS02267
Ml = Ml + LMN RMS02268
MZ = MINO(MZ+LMN,MCASES) RMS02269
GO TO 1 RMS02270
END RMS02271
-------
FORTRAN IV G LEVEL
20
GRAON
OAT:: = 72)06
IO//i/i/08
:»AG:: : 00
0001
SUBROU . I liE G ~ADN
RMS02272
RMS02273
RMS0227/i
RMS02275
RMS02276
RMS02277
RMS02278
RMS02279
RMS02280
RMS02281
RMS02282
RMS02283
RMS0228'1
RMS02285
RMS02286
RMS02287
RMS02288
RMS02289
RMS02290
RMS02291
RMS02292
RMS02293
RMS0229Lf
RMS02295
RMS02296
RMS02297
RMS02298
RMS02299
RMS0230U
RMS02301
c
c
c
MACRO SIMULATION I - GRAON
0002
0003
OIMENSION OCI000)
COMMON 15IMI t'1l ,M2,MCASES,yC 1000) ,XC 1000)
1 ,MARK,LVAR
COMMON IFATI FCI000),LSTORE
COMMON ICAT/ AC20,20)
COMMON I BUG I ISBUG,KRO,LPR
106 FORMAT I1H ,1'1)
NN = 0
00 1 I = M 1 , M2
NN = NN + 1
XO = .DOOOI-YII)
IF CXD-.OOOOl) 'f. 'f. 5
Lf XD = .001
5 YII) =YIII + XD
CALL VALUE
D 0 2 J = t~ 1 , M 2
2 DIJ) = FCJ)
YII) = YCI) - 2.0-XO
CALL VALUE
KK = 0
DO 3 J = MI,M2
KK = KK + 1
3 ACKK,NN) = CDCJ)-FCJ))/CXO-2.)
1 YCI) = YCI) + XD
RETURN
END
ODOLf
0005
0006
U007
0008
0009
[JOI0
001 1
0012
0013
001'1
0015
0016
0017
0018
0019
0020
0021
0022
0023
002'+
0025
0026
-------
FORTRAN
0001
0002
0003
OOOLf
0005
0006
0007
0008
0009
DOlO
001 1
0012
onl3
001Lf
0015
0016
0017
0018
V G _~VEL
20
V~RT
DAT~ = 72006
10/1.f1.f/(8
:tAGE 000
SUBROU" NE V ~ n
~"S02302
RMS02303
RMS0230Lf
RMS02305
RMS02306
RM502307
RMS02308
RMS02309
RMS02310
RMS02311
RMS02312
RMS02313
RMS02311.f
RMS02315
RMS02316
RMS02317
Rr'lS02318
RMS02319
RMS02320
RMS02321
RM502322
RMS02323
C
C
C
MACRO SIMULATION I - VERT
COMMON 151M/ Ml ,M2,MCASES,Y( 1000) ,XC 1000)
1 ,NARK,LVAR
COMMON /CAT/ AC20,ZD)
COMMON / BUG / ISBUG,KRD,LPR
7 DO 1 K=I,LVAR
R = I.O/AC!<,K)
ACK,KI = 1.0
DO 2 J=I,LVAR
2 ACK,JI = R.ACK,J)
DO 1 I=I,LVAR
IF IK-I) 3tld
3 AIK = ACI,K)
AIY,K) = 0.0
DO I.f J=l,LVAR
I.f AII,JI = ACY,J) - AIK.ACK,J)
1 CONTINUE
RETURN
END
-------
3.0
PROGRAM lOA (INPUT-OUTPUT ANALYSIS)
3. 1
Introduction
Program lOA contains a 100 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 num.be 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).
If 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.
107
-------
3.2 Job Deck Setup
The Program lOA (Input-Output Analysis) is a 535 statement
FORT RA N 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
/ ):c
Includes the /; Program lOA
mainline and (Input-Output Analysis)
t h r e ~ s ub - ~
rouhnes ({ .
(535 card;;;) Program Deck
Data Input Deck
/ / GO SYSIN DD>:c
Job Control Language
(JCL) Cards
108
-------
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
table
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
109
-------
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,
itis43.
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'l 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 colurnn 8 and so on.
4.
Sector Ordering Cards: 6 cards (for 42 sector 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 I, 4, 7 of the original table
will be aggregated into the first sector of the new table, punch
loon 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 consume by sector.
110
-------
7.
Change in Final Demand by Sector: 42 cards, one card for
each sector in the new I -0 table
Cols. 1- 10
FlO.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 (lOO 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.
111
-------
'0
The program next tabulates and prints the new I-a flow table.
Next, the new I-a coefficients table is generated and printed.
The pro-
gram now calls subroutine MINVR which calculates the Leontief inverse
matrix of the I-a table.
When control is returned from MINVR, the
Leontief inverse matrix and its determinant are printed.
The program
now reads the estimated average consumption by sector, calculates and
prints each sector I s share of the increoased consumption.
Next, a
series of cards, which contain the change in final demand by sector,
are read.
Changes in domestic 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 econOlny is printed, and the program term-
inate s.
The flow chart of Program lOA is given in Figure 3.1 and the pro-
gram listi~g is provided in Section 3. 5.
11~
-------
Sta rt
Read labels
for new sec-
tors matrix
Print the
new 1-0
flow table
Generate the
new 1-0
coefficient
table
Print the
new 1-0
coefficient
table
FIGURE 3. 1
- - ---
113
Aggregate
original 1-0
. sectors into
I new sectors
-------
FIGURE 3. 1 (continued)
Calculate the
Leontieff inverse
of the aggregate
1-0 table
-.-- -
Print the
Leontieff
inverse
matrix
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-
nomy~
~
114
-------
3. 5
Program lOA (Input-Output Analysis)
115
-------
*******~~*******************************************~********
*
""
...
PROGRAM lOA
PREPARED BY CON~AD RESEARCH CORPORATION
APRIL 1971
...
...
1,:
....
,-
*
TtHS PROGRAM READS IN THE lOOX100 INDIIT/OUTPUT TABLE,
RFI111CES IT Ti] A S"IALLER MATRI X,
CALCUL~TFS TH~ NEW ~GGREGATE I/O CQtFFICIENTS
AND TH~ CORRESPONDING LEn~TIEF INVERSE MATRIX,
GENFRt\TES A f\;()DIFIED COi\JSUMp'TIO'~ VfCTOR,
Q~D ESTIvATES THE ~~T ~FFECT OF AIR POLLUTION CnNTROL
O~ TYE NATION~L FCQN1MY
...
...
...
...
1,:
*
*
1,:
*
************************************************************
100 SECTORS I"JCLUOED IN THE INPUT-OUTDUT TABLE
1 LIVESTOCK AND LIVESTOCK PR~nUCTS
2 OTHr:R I\GRICULTURAL PRODUCTS
3 FOqESTRY AND FISHERY PRODUCTS
4 AG~ICULTURAL, FO~EST~Y AND FISH~RY SERVICES
5 I RON t\ 'J f) F F: R P n !\ LL 0 Y 0 f{ E S M IN! 1\1 G
6 NnNFERf{OUS MFTAL ORES MINI"JG
7 COAL MINING
R CRUDE PETROLHP" I\ND f\JATIJRAL GAS
9 STONE AND CLAY MINING AND QUARRYING
10 CH[~ICAL A~D FERTILIZER Ml~ER~L MINING
11 PE~IDE".JTIAL 8iJILDI"IGS
12 PRIV~TE INnUSTRIAL BUllnINGS
13 DTHEQ rRIV~TF ~O~-RESInENTIAL BUILDINGS
14 lTHER PRIVATE CO~ST~UCTION
15 PU~LIC NnN-R[SlnE~TIAL AUILDINGS
16 H I ~ H ~.., t\ Y S
17 I\LL OTHER PUBLIC CO'JSTRUCTION
18 MArNT~~ANCE AND RfDAIR CO~STRUCTION
1<) URf)I\JA!\ICE AND ACCESSORIES
20 ~EI\T PI..f\JIN,;, DRESERVI"JG FRUITS, VfGETABLES, A!\'D SEA FOODS
23 GQ~IN ~ILL PRODUCTS
24 RA'Wl\:'iJ,
1,:
""
,~
*
....
...
~
::-=
""
..'
..,
..-
*
*
*
~,
..-
1,:
OA 00,
IIJA 002
lOA 003
101\ 004
10.1\ 005
IDA 006
lOti 007
IflA 008
IDA 009
lOA 010
InA 011
lOA 012
Fl.A 013
IDA 014
lOA 015
11"1 A 021
lOA 028
IDA 029
101\ 030
IDA 031
lnt.. 032
10/\ 031
IDA 034
Ill.A 035
InA 036
lOA 037
IDA 038
lOA 03Q
lOA 040
101\ 041
IDA 042
ICJA 041
IDA 044
IDA 045
InA 046
lOA 047
lOA 048
InA 049
lOA 050
IDA 051
IDA 052
lOll 051
InA 054
IDA 055
InA 056
!rJA 057
1-;;/ ,e C'
, . I ~ ,:. !
-------
?9 T08ACC'J ,"1 ANt :~CTUR:=S
30 R~lAD AND NARROW FaB~ CS, YARN ANn THREAD ML_S
31 ~ISr. TEXTILE GOJDS A~D FLOO~ CnVERT~GS
32 APPAREL
33 MI)C. FA~RICATEn TEXTILE PROQUCTS
34 llJ'.a3ER AND If/Don PRODUCTS, EXC. COf'nAINERS
3~ ~O~DE~ CO~TAINERS
~6 HnUSEH:-JLO FUR"iITURE
, 7 0 T H r: R F U [{ NIT UR E A. 1\1 D F I X T lJR E S
38 PADFQ ~~O ALLIED PRnOUCTS, EXC. CONTAINERS
3q PAPERR8ARO CO~TAINERS AND ~OXES
4 a DR I NT I i\; G 1\ \j f) rUn L r S I-q N G
It! CHE,vlICl\LS AND SEL ECTED CHE,"1ICAL PRODUCTS
42 PLt:lSTICS A.ND SYNrHFTIC t1tiTERII\LS
41 DRlJGS, CLEI\NIN(;, I\ND TnTLET PREPARATIONS
44 PAINTS ANQ ALLIE~ PRODUCTS
45 PFTROLEU~ RFFT\JI\JG A~D RFLATED INDUSTRIES
4(-, RUf\RFR A"!O r-nC:;Ct:LLt.N!=OUS PLI\STICS PP,Of)UCTS
47 LEi\THER TI\\ji\JP~G A\!n !f\1f1USTRII\L LE"PiER PRODUCTS
lt8 FO:]TWFAR !\ND OTHER LEA THER PRODUCTS
'tC) GLJ\SS -\I\/D GLASS DRODI/CTS
50 STrJNE AND cuw p'
-------
76 A I P. C R AF T AN 0 P I\R T S
77 On~ER TRMIJSD(lQ'ATION EQlIP\1Err"
7~ SCIE~TIFIC ANQ CONTR~llING INSTRUMENTS
7q OPTICAL, CPHTYALMIC AND PHOTOGRAPHIC ECUIPMFNT
1)1) '/01 SCELLMJEnus ~(\\JIJFACTUQ I\lG
81 TOl\j'JC:;Pf1P,fL\TfGN I\~'D \"'I\R~HOU5ING
'p Corv1f"iI/i'JIC.'\T!Or--jS, ~XC. ~ADtJ AND TV BROADCASTING
83 RAnIO A"JO TV BRql\OCASTING
84 E IJ C T R f C lJT I L IT I f S
81) r,!\S UT!L I TIFS
A6 WATER. AND SI\NITAR.Y SERVICES
87 WHfJLESl\LE AND RETAIL TRI\DE
88 FPJANCE .l\ND P.!5U:: A~JCE
89 REA L F Sf" T t: lI'\! I) Q F t\j T 1\ L
90 HOT~LS, PFRSO~AL AND REPAIR SE~VICES EXC. AUTO
Yl BUSI"H:SS SfPVICES
q2 RE~EI\~CH 1\"J0 DEVELOPMENT
93 I\UTO~O~ILE REPAIR I\NO SERVICES
94 AMUSEI'v1::NTS
95 \1Ff")IClIL, FDUCl\TlnNAL SFRVICES, AND NONPROFIT ORr..
96 FEneRAL GOVERNMENT ENTERPRISFS
97 S T !\ T E '\ N I) L 0 (: 1\ L GO V E Q W1 F. NT EN T E R P R I S E S
98 G~~SS fMPORTS OF GOODS I\ND SERVICES
q9 ~USINESS TRAVEL, ENTERTI\IN'v1E~T, AND GIFTS
loa GFFICE SUPPLI=S
DIMENSInN A~(150)
D I 'v1E N S 10\1 C f) L S U ( 10)
DIMENSIO~ RISO,50J,C(Z2,31J ,SUM(31)
DrMENSIG~ L~oF.L(10,5J)
nIMF~SllN STORF(50,3J,P(50,50),SI50,501,VECTOR(SO}
C 0 t.J \1 0 '\J 1\ K ( 1. 0 8, 1 5 0 )
C () 1'-'''' 0 N ;~ ARK, I T EST, I R. 0 ( SO) , I R O\-!, S T E P
C 0 (I; >, O:-J 1\ I 1\1 ( '3 0 0 0) , A ( 1 20 , 1 50 ) , N C , NR , I C , I Q.
CO~~ON ~RD,LP?,KP
Fourvt\LF.NCE (i'.H~, lR), (5 (1, I) ,~,(1, 1), ( 1\K( 1) ,VECTOR( 1»,
1 ( H K ( 1 , 1 ) , p ( 1 , 1 J ) , (II ( 51 , S 1 ) , S T OR E ( 1 , 1 ) )
K~f) : 5
LPR == 6
KP = 7
STF.P = 0.1)'
on 10 1:1,50
10 IROeI) = I
I R ([.:: 0
'v1IJl T== 10
T '" W)= Q
,/If 1;= 1 r)f-\:)
OA
OA
IDA
rot\
101'.
rnA
InA
lOA
101\
I'lA
1011
r r1 A
10/\
FJA
rOo.
I no.
10/\
T1A
10"
101\
rnA
IDA
IDA
IDA
lOll
lOA
rr)~
lOA
InA
IfJ/\
IDA
IIlA
lOll
Ir)A
IDA
10/\
TOti
IQ"
rot.
In"
UJ A
I0A
InA
l'lfl
lOA
1'1/\
J ,'-:
06
107
108
109
liD
111
112
113
114
115
116
1.17
118
119
120
121
122
123
124
12 S
126
127
123
129
130
131
016
017
018
019
020
021
022
023
024
025
026
132
133
134
135
136
137
138
13Q
140
-------
Of) 70 M=., .2
I( K =1
IF (M-12) 21,22,22
22 ~ULT='j
1.'\1]1)=4
J~-fI=540
21 D1 25 J=1,108
I l 0= ( J - 1 ) * ~UL T + 1
IHI=TLO+lflf)f)
25 REA) P
-------
PR ~-S OU' T~F NEW GENERAL FLOW TABLE
"1 = ')
501 00 601 I=l,IC
01 602 J=ldR
'-1 ="1 + 1
602 :3 ( I , J) = A IN ( ~.1 )
601 CllNTP~IJr=
TABULATES THE COLUMN SUMS
O'J 603 J = 1.IC
B(J,I:=(2} = 0.0
DO ,'J03 I = l.I~l
603 U(J,IR2) = 8(J,IRZ)
q( IC.IR) = 0.0
+ B(J,I)
I R 0'11/ = I R 2
ITEST = IC
STEP = 0
CALL PRINT(!)
GENfRATES TYE AGGREGATE I/O TABLE
WRITE(lPR,655) IRl
655 F1R~AT(lHl,~5X,'I~PUT/OUTDUT
I) 0 606 r! = 1 , I C
CO:--J = B(N,!R) + B(N,IR2)
BPJ,IR3) = 0.1)
IF(CClN)6,606,6.
6 00 607 I=l,I~
B(N,I) = ['\('.j,I)/CO\l
627 BCN,IR.3) = rHN,Ip,) + I3(N,I)
607 CONTrNUr=
P.(r-.,;.I,'{2) = ;'HN,IR?)/C1N
606 CO /lJT I ~ U [
TABLE (',12,' - SECTORS)'!n
I R n',4 = I R j
ITEST = IC
STfo = 1
C 1\ L L PR I "J T (R)
I:JRITE(6,6666) IR1,!in,IRl,IRl
6666 FOP'I,I\T(/!1-t5X, "IjIJT[..HH.J-', 12,'+1 = VALUE ADDED (PERC["JTI\CE)'
1 151X,';'W.-J-'.I?,'+2 = PHEQr<1r;OII\TF. P~OIJUCT (PERCF\JT!\GF, SlIM
? T.I <', 1 - '; T) ')' i'r, 1 >: ,;: 'I ;".~ ..17: c; rei 'II' :!
Ot. .89
TOA .90
lOA 191
lOA lq2
InA 193
IDA 194
lOA 195
1'1 A 196
lOA 197
10 A 198
IDA 199
lOA 200
lOA 201
I'H 202
lOA 203
IDA 204
lOA 205
InA 206
TnA 207
IDA 208
lOA 209
InA 210
101\ 211
lOA 212
lOA 213
lOA 214
IDA 215
IDA 216
101\ 217
InA 218
IDA 219
IDA 220
InA 221
IO!\ 222
Inl\ 223
InA 224
Tn A. 225
IDA 226
lOA 227
TOA 228
1QI\ 229
101\ 230
lOA 231
IlJtI 232
Inti 233
[IF P. I (' f\ ?3't
c'
-------
01\
'lA
lOA
l'1 A
lOA
lEONTIEF INVERSF ~~ATRIX OF THE INPUT/OUTP\JT T.AIOA
)'//J lOA
101\
lOA
lnl\
lOA
lOA
InA
lOA
lOr..
trJA
[[)l\
IDA
lOA
IrJA
IrJA
lOA
[ OA
I r):\
.rOA
lOA
I Ill\
In!\
lOA
lOA
lOA
rOA
DOlLARS',5XIOl\
IDA
FJA
lOA
lOA
{nA
Ir)l\
In"
IDA
I il A
lOA
Ifl t\
101\
In"
l
'.
,
THE LEO\!TIFF INVERSE
WPIT[CLP~,742) IRI
742 F (l R tJ.!\ T ( 1 HI, ;? 1 X , ' THE
IBLE C''!3,' SECTOQ,S
on 7701 I=ldCI
rF) 7702 J = 1,ICI
7702 ~(I,J) = -RCI,J)
7701 RC I,I) = 1.0 & !3( I,I)
\1A= ICt
DET" = 0.0
CALL ~I~VR CH,MA,DETAI
WRITE CLPP,773) DETA
7 1'3 FOR 'J A T ( 1 H ,I D!:: T E R "1 I N A \j T = , , F 1 2 . 1 0 )
I R n',oJ = r C 1
ITEST = leI
STFP = 1.0
CALL PRINTCK)
\1nnIFIEO CO\!SUMPTION VECTOR
CON = ESTI~ATED MA~GINAL INC~EASE IN cnNSUMPTION
R{ IC, 2) = Pt~CENTAGE CO~SU~1PTI'JN OF EACH SECTOR
AK( [) = fACH SECTORS SHARE OF THE INCREASED CONSUMPTION
III FORtvlA T C F6. 11
READ(KRD,1111 r.ON
',.JR ITE ( 6, blt 7) C nN
647 ~OR','ATC IHI,18X,' CnNSljr~PTION ~Y SFCTOR'/25X,
l'BENEFIT - CHANGE IN CONSUMPTION'
I ' = , . F 6. 1,' ~1 ILL ION DOL L A R S ' II 3 5 X , ' PER C E NT' , 7 X , 1 I'" ILL ION
2'<)ECTOR n.F.FI~JITION')
DO 641 T=l,pn
A K C I) = CO'\J 1" G ( I C , 1 )
643 !:JRIT~(6,644) I,B(IC,I),AK(II,(LABEl(J,II, J=1,7)
6 44 F n R Jj l\ T ( 1 R X , ' S E C Tn p -, , r 2 , F 1 (t . 6 , 8 X , F 1 0 . 2 , 8 X , 7 A 4 )
I:
I~
~ULTIPLTES THE LEONTTEF INVERSE RY
1. \.1f1I:IFJEfJ cn~!:;I.J'-1OTIO\J VECTf1R
2. CIJNTfUJL C'JST VESTIJR
!'-"
~
'~
..,N = I R 1
STED = ?
I ~ nl., = {R
1- ;- 1<:"
::--i' ,;,(
T I~ ,\\
236
237
238
239
240
241
242
743
2lt4
245
2lt6
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
27f5
276
277
278
279
280
781
...,-, ,-, "'r',
/' 1'" /~
-------
IOA
OA
10./\
lOt!
106.
POLUTION ClNTROL COST - CHANGE IN FINAL OfMANDIOh
IDA
lOA
IDA
lOA
InA
lOti
lOA
IGA
InA
lOti
10"
IDA
I nA.
If}A
IDA
I8A
CONTROL. EFFECT nF COST ON STRUCTUR~InA
IDA
.If")A
lOA
lOA
AIR PC1lLJTION CONTROL. BENEFICIAL EFFECT ON STRLJCTIOfl
ECON:)V,Y'//) lOA
IDA
101\
IDA
IflA
UJfI
lOA
I::JA
rrjA
rOA
lOA
TOI\
1H1,30X,' AIR PIJLUTJON CONTROL. NET FFFECIOA
~ATIONAL ECO~OMY'II} 10/\
TIJA
JrlA
STRLJC IrJA
If) /\
I fl ",
:J,ARI( = 1
r;o TO 444
1000 ~ARK = ~ARK + 1
~rJRITEC6,777}
7 77 F fJ P ..., 1\ T ( 1 H 1 ,2 5 X , ' A I q
J'/I)
REtif)(5,1l> (VECTOP(J), J=1d4N)
71 FOQ'v\/!..T(1CX,'SFCTOR-' ,I2,1=11.3,4X,7A4)
I1 FOR~I\TC~la.2,70X)
f);J 7777 J= 1, fl.1N
'..JPITE(6,7U J,VfCTOR(J),(LABELCI,J), 1=1,7)
i 7777 CONTII>.JUC
I 444 COf\iTIf'!U'
I) 0 2 '3 I K = 1, f" i'J
S T OR f ( I I( , .'vI AR K) = 0.0
1) 'J ? 3 J = 1 , 'v11\!
S C J, I K} ':: !3 ( J, I K) f.' VF C TO R ( J )
STOPE(IK,f\1ARK}= STORI=(IK,MARK) + SCJ,IK)
23 cn"JTI~I)~
IF(""t\P.K-1) 32,31,32
32 vi RITE ( 6, 3:\ )
3~ FO~~:1t\T( IH1,20X,' AIR fJnLUTION
I or NI\TTrJNAL ECONOi.1Y'//)
GO Tn 34
31 WQITEC6,35)
35 F rJ R 1.1 A T ( 1 H 1 , ? 0 X , '
lUQE nF NATrONAl
00 ~6 I C T = 1, MN
I) 0 36 leD = 1, r-11\1
36 PC ICTdCr} = SCICT,ICD}
34 CALL DRINT (5)
I F ( 14 A K K -?) 1 no 0, '3 7, ? 00 0
37 r.1 .1. R K = '3
00381=1,!\.1N
ST[1?I;:(I,3) = STOP[(I,l) + STC'RECI,2)
o (] 1 R J = 1 , '.1 '\J
38 S ( I, J) = P ( I ,J) + S ( I, J)
~v!\ IT E ( 6 , 33 n )
330 F'JQ,\1l\T(
IT O~ STQUCTURE OF
GO Tn 34
2000 CONTI~!LJF.
81 FOR:-1t\T( Hn,?ox,'SUM~~1\RY OF FFFt:CT OF AIR POLUTION CONTROL ON
1 T U I{ F f) F "tA T IrJ N "L E C 0 WH.' Y , I '5 0 X " ( PJ rv1 ILL r 0 N S ) I / /
7 ? l} '( " ,? ( ':I, .x I ( r f F C T r. F ' ) , nr , ' ~J r T ' / 7 C) ~( \
2R3
2R4
28';
2R6
287
288
289
290
291
292
293
?q4
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
31q
320
321
322
321
324
325
326
327
12R
,
':?';
-------
2 3 X " ~:: 1\' EFT' , 3 X , '
ces" ,7X, 'EFFEC.'/)
101\
f)l\
IO~
lOA
IO~
lOA
In/\
InA
IDA
rn/\
IDA
IDA
Ie) 1\
lOA
lOA
IDA
IOl1
IOt\
IDA
TOll.
InA
10/\
BEEN qEAD IN AND STOPED lOA
FI LE C1\RDS IDA
10/\
10/\
10:\
I'JA
lOA
10/\
IDA
COLUMN lOA
IDA
IDA
10/\
r()A
lOA
11;1\
10/\
lOA
10/\
I f) A
FJ /\
lOA
lOA
Inl\
Of] 2222 K=1,3
? 222 S T O~ E ( I q 2 , K) = O. 0
fi~ITE(6,fH)
00222 I=l,\'N
W P, IT E (6, 44 ) 1,( S TOR F. ( I , MAR K ) ,M ARK = 1 , 3 ) , ( lA BEL( J, J) , J = 1 , 7 )
D'l 222 ,<=1,3
ST[lRF=(r~2,K) = STCRE(IR2,KI + STOq,E( I,K)
222 CONTI\jI.Jr:
',.J R IT f ( 6, 4444) ( S T Oq E ( I R? ,NJ ARK) , tv} A Q, K = 1 ,,3 )
4ft r:rJC~'VJAT( 20X, I S':CT(JR-', r 2 ,3F12.2, 5X, 7A4)
4444- FOR ~,11\ T ( '24 X , ' TOT ~ L ' , 3 F 1 2 .2)
STOP
E"JD
SUWHJUT l"JE AGGRf:
[) I ~1 [= "J S I IJ \1 8 ( 1 5 0) , J COM ( 1 0 0) , K C 0:-1 ( 1 5 0 )
C1M,'-IlPJ 9K ( lOB, 150)
c: n ~11" 0 r J v ARK, I T 1: S T , I R 0 ( 50) , r R 0 h', S T F. P
COPt."'!]\! AIN(~OOO) ,A{ 120, 150) ,NC,:\!R,IC,IR
COM~O\l KRD,lDR,KP
PR()GRA~ ASSUVJES MATRIX TJ BE AGGREGATED HAS
BY COLU:-1NS 0\1 r)l SK AW) AL S'1 THAT THE DEF INE
IN THE MAIN PRnGRA~ AQE AS FOLLOWS
NC = rrn I\l = OF COLlJ,''':-JS I N INPUT :v1A TI:~ I X
I\JR = TtlT/\L = m= ROltiS ["'I II\JPUT ~-1ATRIX
IC = TOT1l = OF (OlU~NS IN OUTPUT MATRIX
Iq = TnTAl = OF ROWS IN J~TPUT M~TRIX
** W~RI\JI~G ** THF LAST ROW IN THE nUTPUT ~ATRrx MUST RE
lOR (V!\LU~ ADDED) IN THE INPUT TBl. THE LI\ST
IN THE OJTPUT MATRIX MUST RE COLUMN 108
( C () N S U :-\ P T ION) I NTH E I 1\J rUT Nj AT R I X .
READ (KRD,700) Nr:,N~, !C,!R
\~ R I T E (L P R , 70 a) N C , N R, I C, I R
Df) 5 [ = 1,IC
Dn 5 J = 1,~R
~K(I,J} :: 0.0
5 CONTI"J'JF
700 FOR\1A T (2014)
IF (IC-IR) 11,18,10
10 f"iAX= IC
"JIv1 AX=NC
Gn Tt] I?
it !Iv! !\ X = I R
'J :,1 t\ X '" .'J r~
<~ \'-','
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
V56
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
-.(
-------
'..
.1. R.EAO (KC\9,7~O (JCOt~(I),r=1tIMAX)
R. E AI) K Q I), tno ) K crr,1( I ) tI = . ,NMA )(
1<.=0
") 0 'Y) I = 1 tl r:
J J = .J r. C1:-1 ( I )
DO 35 J=l,JJ
I(=K+1
I I). FC=KC nrv1( K \
rJ0, 16 L=l,N~
36 I')K(I,L) = r3K( 1,L)&A{PEC,L)
35 r:n~!T I'~LJ~
30 CfJ~,TI \!1.1~
1TOP = 1R*1C
f)'l 41 1=1, ITOP
41 AIN(I) = o.n
M=O
on 40 I=ltIC
K=IJ
I) 0 'to J = 1 ,r R
J J =Jcnr.1( J)
~'1= 1.1+ 1
f)fJ 43 L=l,JJ
K=K+l
IF(K-IO'1) 48,47,47
't 7 '..J R I T E ( 6 , 4 g) K
49 FOR:" t\ T( 1 H 1 , , E ~ !:n R . . . T ~ I E f) TO A'J 0 R m-l- , , I 3 " AND THE R E
l'/lx,'r~'=ANrr.lGFUL RQ'.-.JS IN THE INITIAL I/O TABLE'//)
STOP
48 CO NT! ~IJ E
KI(=i,CTvJI K)
41 AINU") = AIN(M)&GK(I,KK)
40 CONTII\!lJF
OUTPUT 1vjllTRIX ,\R~t\"JGED I\J A!t'J RY COLUMNS
QETURN T:) ~1AIN P~O~,Kl\~ TD WRITF. AGGREGf\TED t.,J\TRIX
RE TUR \J
F "J r)
SURROUTiNF. PRI~T(B)
J I MF \J <; I r \J g ( 5!) , 50)
C n W" or....: 3 K ( 1 () R, 1 ') 0 )
c 0 '.~'! '1 \)'.. ARK, I T f S T , I R n ( 50) , I RO IN, S T E P
C D ~ '.1:)'~ .\ I I\J ( 3 0 0 (» , A ( 1 ? 0 , 1 5 0) , N C , ~J R , Ie, I R
cn~~ol\J KRO,LPR,KP
700 F n R \~ 1\ T (1 H . 7 x ,g ( ') X, 'C OL 1 ,I 3)
7 tJ? F fJ R '4 A T (1 X , , p n 'fl' , I ~ , ? X , 1 o!= 1 1 . 0 )
70:; I=OP:1I\T (JX,'Rn'!' ,T3,?X,lOr11.4J
7 (), F i) r~." fI. T (1 X , ' '<. : J '.}' , I 3 , ? X , 10 r 1 1 . -, )
"7 I.~':
F ~~~ F: \i T!: 1 i~(,.~ (' .~,.:
I[)~ 377
nA 378
lOA 379
IQI\ 380
Jr) II 3Al
lOA 3R2
IDA 383
101\ 184
IDA 385
104 386
TrIA 387
lOA 3BB
lOA 389
lOA 390
I[lA 391
Tn III 392
InA 393
lOA 394
lOA 395
lOA 396
101\ 397
lOA 3ga
IO/\ 399
r 0" 400
IOA 401
ARE ONLY 10aWA 402
lOA 403
l'JA 404
10fl 405
lOA 406
IDA 407
IOJ\ 408
IDA 409
InA 410
lOll, 411
IDA 412
lOA 413
IOJ\ 414
IDA 415
IJA 416
IDA 417
InJ\ 418
10/\ 419
lOA 42')
1n1\ 421
f'll\ 4??
f '] i\ ~
-------
i.4
K = 1
~ARK = ( TE~T&B)/9
DO .'q02 J=1 ,"1/\~K
I.J = J*9
F(IJ - JTESTPW5,g04.804
R 04 I J = IT CST
8 () 5 [.J Q I T F (L n R , 7 0 0) (T ~ 0 ( I ) , 1 = I K, 1 J )
IF(ST~P-1.0) 3,3, 332,314
333 Dn 801 v=1,IROW
801 ..IRITF (LPR,702) IPO(M),(R(I,M),I=IK,IJ)
GO TO 9g
334 DO ~22 ~=1,IROW
3 2 ? \.1 R I T E ( L P ~, 70 5) I R 0 ( M) , { P ( I , r-1 ), I = I K, I J )
GO TO 9q
,32 00 811 M=1,IROW
ItJ R I T E (L P R , 703) I RO ( "'1) , ( B ( I ,M ) , I = I K, I J )
g 11 C n f\! T p-j lJ !==
99 CO NT I "ILJ c:
IK = IJ & 1
',~ RITE (l P R , 7 0 1 )
R02 CJ~JTINUF
7 0 1 F n R .~ i\ T ( 1 H 1 )
R[TIJP.N
FNO .
5UR~OUTINE MTNV~ IA,MA,DETA)
OIMFNSlfJN A( 50,50), IR.t 50), ICI 50)
CO "'1:-1IJi'J RK{ 10.'3, 1')0)
C 0 ~:A'" 0 '\J ~"Q. ,<, I T f S T , I R n ( 5 a ) , I R [) w, S T E P
r.O/''''~O:\j ;'IIN{100Q) ,O( 120,15() ,"JC,"JR,ID,'IK
COM~O~ K~O,LPq,KP
DO 11 = 1, '-1A
IR(II=O
1 Ie ( [ ) =0
I)ETA = 1.0
TOL = .nOOODl
5=0.0
R=M/\
*'3EGIN*
2 1=0
J=O
TFST=O.0
IYl ~ 5 K = 1 , MI'.
IF(IR(K»25,39,?5
39 CONTINUE
f)'J ~/-t L= 1, "1.~
[I=IIC(l) 1?4,,?R,?4
, 0, r. 'I '''; r I 'I 111:
f )-
OA
OA
lOA
lOti
1()A
ICJA
InA
IDA
TrJA
IDA
rOA
ICJA
10/\
lOA
IDA
10/\
I'JA
lOA
I () A
10/\
lOA
10/\
I [J 1\
lOA.
FJA
101\
Ifll\
lOA
lOA
IDA
lOA
[()A
If1/\
10/\
IO.A
I nA
IOI\
10/1
lOA
101\
IDA
10/\
lOA
InA
lOA
I rJ 1\
424
425
't26
427
428
429
430
431
432
433
434
435
436
437
438
439
440
4't!
442
443
4',4
445
446
447
l-t/+ 8
449
450
451
452
453
454
455
456
457
lt58
,+ 59
460
461
462
463
464
465
466
467
468
409
',.' 'I.
-------
X=AP'S(fdK,- )
IF (X-TEST ~4,37,37
37 CONTI\JlJ!::
I=K
J=L
P:ST=X
24 CQ~!TI NU[
25 cnl\JTI f\iUF
l~
I-
I:
[
I~
THIS IS THE TEST FOR A OEPENDENT TOLERANCE
(!\ 471
InA 472
lOA 473
10 A 474
I f) A 475
FJl\ 476
lOA 477
lOA 478
10~ 479
IfJA 4RO
lOA 481
IIJA 482
F)/\ 483
lOA 484
lOA 485
lOA 486
101\ 487
lOA 488
lOA 489
lOA 4QO
lOt; 4Q1
J 0,\ 492
lOA It 9 3
InA 494
.101\ 495
IOA 496
InA 497
lOA It 9R
lOA 499
I'lA 500
IIJI\ 501
TOA 502
I CJ 1\ 503
Ifl A 504
lOA 505
lOA 506
InA 507
10/\ 508
lOA 509
IflA 510
lOA 511
lOA 512
10/\ 513
lOA 514
IIJA 515
Inl\ ')16
~ i-I ::, 'ti'
1!=(TOL)73,74,74
73 TOL=A3S(TOl~TEST)
74 CO~T I NUE
c
:: 1;
*Ef\![)~' 5U8/v1XS :t$
I.
PIV=A(I,J)
OETi\=DIV':
-------
I =(K- )4.,16,4.
41 CO~\j- "HI;:
OETA=-[)~TA
811'-tL=1,MA
TF;'-1o=~("",L)
!I( K.L )=fH r ,L)
14 A{I.L)=TH,\P
D(]15L::l,~:A
TEf.t1D=A(L,\i1)
I\{L,/..q=!\(L,I )
15 At L.I )::T[MD
IC('-1)::K
IR(I<)="'i
16 CONTI \lIJF
LEAVE WITH A WARNING OF ERROR
17 PJVFR :: S
RETURt\J
END
OA
IDA
~(1A
InA
I [) fI
rn1\
lOA
lOA
lOll
10"
IDA
I[]~
r 0/\
rnl\
IDA
101\
IrJ/\
lOA
58
5 .9
520
521
522
523
524
525
516
527
528
529
530
'531
532
533
534
535
-------
,.
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 purpose
of the program is to take the output from Program lOA and then distribute
the interregional feedback to each AQCR.
Thes e in te rreg ional feed-
backs are measured in terms of changes in value-added by two-digit
SIC detail for manufacturing industries.
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 5 ections.
128
-------
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)
(/'"
','
~-
l~
109 cards
Data Input Deck
/,',
'.'
Program F BE
/ / Job
/
Job Control Language (JCL)
Cards
129
-------
~-
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.
Cols. 1- 10
Cols. 11-20
F 1 0.4, benefit effects
FI0.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
SICs by AQCR (SIC 20-30) (one card for each AQCR).
Cols.
Cols.
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
51-55
56-60
15, AQCR 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
F5.4, share £01' industry 26
F5.4, share for industry 27
F5.4, share for industry 28
F5.4, share for industry 29
F5.4, share for industry 30
The second 100 cards contain market share of second nine
two-digit SICs by AQCR (SIC 31-39) (one card £01' each AQCR)
130
-------
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
4.4
15, AQCR code
F5.4, share for industry 31
F5.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
Description of the Program Deck
Program F BE 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 multiplie s the cost and benefit for each sector by its re-
spective 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 classification.
The program now punches and writes output
cards containing (1) the benefit and cost to "all other" industries (1 card),
and (2) the cost and benefit to the 20 two-digit SICs (20 cards).
131
-------
.
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.
Finally, the
deck containing the net figures is punched and written, and the pro-
gram terminates.
132
-------
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 "aU other"
1
133
FLOW CHART OF
PROGRAM FBE
-------
1
Write resul-
tant cost and
benefit fig-
Read distribu-
tion matrix
Calculate cost &
benefit share for
each AQCR/SIC
combination
Punch resul-
tant shares
Write
results
End
134
FLOW CHART OF
PROGRAM FBE
(continued)
-------
4.5
Program FBE (Feedback Effect)
135
-------
C
r
'.
::: ~~ ::. ::= ~. =:= -::.. ::::. ~~ ::: =:: ::: ::: ::: :'~ ::: ~::: ~~: :;: :::~ 1,: :::... :::.. ;-;~"/,: :;: ~~':" ~:- ':-: :::: :::: :::: ::': :~: ::': :::: ::: ::: :::: :::: :0:: :::: ::~: ,.;: ::~ .:: :::: :::: :.: :~~ ::::;' ':: ::;: 1,: 1,: ::::: ::: ::=:-;:
:}:
r
C
C
C
C
L ~ r ,'"", :; ;j ~: . J. ), F-
:::
,,-
C r. \. S .. : ',> f S f ,\ Q C H
". r:> ILl t~ -, J
c C :.' r ., (. 'I Tin '!
,,-
"
p" t=:J .\ ~ f- :; ,:. v
~::
,.
P~T:; D~"rl(;''''M,~ C'P,T/lT;!S /' ,)['GU:\:'\1. r,l\r)iq:r Slti'\~:t: I\.I.6T(;>I>:'
i',',I'''; (1\l(:<'L \T(:5 p,r l' T:-~):':-;,I:"\:';!',L FF~r.rl\CV, r:r-:r-!::,~T
T~l 1 "I') 1\('( r
r
'.
,'-
','
C
r
'"
=::
,'.
','
c
**~:~~.*~:***~:~::::*~::*::~~:::::~**~:*:::~~::*~~:*~::~:*~:~:~:**~:*~.~:~:**~:<::~*:::~:***~:*~::**
') r 1 q S I~ r.; X ( It 7), Y ( L, ?), 7 ( 4;' ), IJ ( .2 ), r (1 S T (2 <:) , rq: \J r: (? ') )
!) r ',' F" S J n \1 ~1 1ST" ~ ( 1 (}C , ? i) ), r: r <:, T ~'. /I ( 1 ~ C , ? ('\)
) 1\ T t, C r: S T, ~ r \;" I ') r- ::' r. . f"\, ? (' ::' 0 ." I
r'T :.\!C,!. '\1 -JI <; T'-'C ( 1 j".,? ( )
~ATA V I ? * 0.0 I
1 (' F ,l-"'~ '\ T ( P r 1 (', . I" )
1):-1\'1 (?,H) UIJ), I = 1,L.o?)
):1 ? ;, ! = 1, It I
;') C 1\ r: ( I) , 1 'J) Y ( I ), Z ( r )
'{ ( I) = Y ( I) ::' X ( I )
! ( r) = 7 ( T) ~. ~ ( I )
7 C C(')\JT J ;,' ~ If:
nr] ":\0 J = J.,9
VI]) = \,(1) + vITI
~~ V(~) = VI?) + 7( I)
CroSTfl) = 1(1:1)
"~\Jf(1) = Y(I'))
C0ST(;>' = 7(11)
~ C'.I 1: ( ?) = Y ( 1 1 )
'::'1<;T('3) = .'J ':' l(}n
'3 E ~~ E ( 1) = . ') :~ Y ( ] 2 )
C~ST(4) = .5 * 7(1;)
.) r.: \' elL.,) = . 5 1,. v ( 1;> )
')'1 ltC' T = 13,21
J = 1-8
Cn'=;T(J) = I(I)
') C '! ~ ( J) = Y ( ~ )
4'J ::: n' : T I "; Ie
f)r ::,,' 1= ??,7'i
-:: 'l S T ( 1/+) = C ') S T ( J l+ )
~ r \'::- r 1 ~) = '\ t '.1 ,: { ] I. 1
+ 7 ( T )
+ Y ( I )
~;:. :'-''',!TT'I::
') ., f) r' T = 7 {, , ;>' ~
J = r - 11
C1'~,T(J) = 1(1)
iJ '" I ; 'j
~ ,. r. 1 )
:::
-'-
,,'
;;:-:
::::
'::
':=
'f
.",
','
oJ,
','
: p r- 001.
Cf) r 1')1)2
~IJ ~ ';01
FPC () r It
frr:: A('!C)
FI\F r,of.,
FI\C 0C7
~rH: OOg
FBE 00Q
PH: 010
F 1\ r:: 01 1
Fhf ~12
FP.f: 013
PIE ') l't
F~\ ~ ,}15
FRf .r; 16
rF-i[ (17
FRf- ')1:3
F f1 F: :!]q
t:f1E nZr1
FP. E O?l
Ft\ E n2~
F[i E 073
F'~ r:: 0?4
FR [; 02'1
F::\ E r;7S
r:: [\ E 027
FRE 02q
FBE 029
Fflr: 030
FBi=' 011
q~F :) 37
Flit: (n3
FB F n3/+
FRE 1)3')
F~J: ()l-S
FfI'- 0-31
~~F 818
fl3 E ')3°
F r\ ~ 1)4')
F~r- Old
F~E () 1+ 2
Fr, r '')I.d
P), F I) 4/t
FRE 04'1
FP. C 'lL.of.,
r 11 r r) I, 7
-------
6 C C 'Y" T I 'JlJ':
r, f1 S r ( 1 q) = l ( 2 C)) + 1 ( 3 C )
"3 [ ~J [( 1 '..) = Y ( 2 'i) + y ( ~ C )
:),17; r = "31,32
.J = 1 - 12
CCST(.)) = 2(1)
, R [ \J F ( J) = Y ( I )
7" crl~JT;\I:Jr:
')(' fU' ; = 3,,1,2
I/(}) = V(lI + Y(I)
If ( ?) = \' ( '2) + 7 ( I )
E() r:f~~rri\!ijr:
'.J q IT;:: (f." a ()) If
q (1 F:1 R"1 AT ( 1 H ,? f-- J.'.:' . ~ )
;-Jf.l 11(' r = 1,2:)
J = T + 19
',/ R I T F- ( 7 , } ~ (1) (. : 5 T ( I ) , ~ F '~F ( I ), J
\1 Q I T F ( ~, , C) 0) r "1 S T ( T ), ~j F \ J C ( T )
1(,'! '=)~"1!\T(?H0.1 ,::;r)x, r71
1 1 IJ crJ\J T I I,! II f
1)0 120 r = 1,}1)f\
:>F:'\n(r;,lIS) (i)ISTD1(I,J), J = 1,11)
115 F~QMAT(5X,11FS.4)
1 2 (' en to.J T I ,\'11 F:
f):l 1.?C j = 1, 10,J
':>f:=i\,I)('),llS) (0ISr)l(I,J). J = 1?,7())
'l!" l? ') J = I,? (j
l I <; T r; iJ. ( J , J) = I) r <; T "') P ( r , J 1
nrST::'r\(I,J) = 'lIST!~G(J"J)
1ISTP[IJ,J) = 0rST~~(I,J)
125 cr.'"n r !'JUF
I 1 0 C C \J T ! \J! J F.
i)fl151; 1= 1,110
:It.!. I T F ( ,L), ] L. 0) :J T ~ T I ~ . ~ ( J , 1 ), (n T S T I: n ( I , .J), J =
~:~ !T f' ( 7,14')) rl1 S T ~,\( 1 , 1 ), (0 J S T'~ r. ( I , J ) ,J =
.J q I T r ( ~J, ] tt (";) () J S "< r ( I , J ), J = q, 1 5 ), I
:' ') r T !: ( '7 , 1 4 5 ) ( 1 T ~ 'T " Q, ( ! , J 1, ,J = 9, 1 '5 1, I
-".-: r T t, ( fJ , 1lt I) () J ') T ') 'I ( T , . J ), J = 1 f, , ? 0), I
:,,'P I T F ( 7, 1I~ 7) (,) 1ST '{ r) ( I , J ), .J = 1 f~ ,2 (1 ), I
1 c:; () crJ JT r 'J: I ;-
i) r: 1 f, '! 1 = l, 1 OJ f\
J:< IT;: ( :... li. r;) ,) i S T P :\ ( I , 1 ), (!') I S 1 F f\ ( r , .J), J =
~ r,' I T r: ( 7 , 1 '+ ~:) I)! <; T [:\ '\ ( J , 1 ), (rj 1ST J' fI ( r , .J ), J =
\-}" I T r ( h, 1 t. ,...) (') j S 1 '.' /I ( I , ,J ), .J = (1, I '-' ), I
,J J.~ r T I ( 1, 1/, S) ('J 1ST" ,\ ( J , ,J ), .J = 9, 1 'J ), I
,/ ~). ! T r' ( 1" I I, (» (') 1 ::; 1 'J 1\ ( J , .J ), . j = 1 (, , ? () ), I
! " r f r. ( Ie 1 I, '( I (, ) I c: 1 . .' ( I , ,I ) ,J - 1 I, : ? r \, T
,;-: ;\ E ~J E ( ,J )
':' r.:JSTtJ)
+ i/ r S TP f\ ( r , J )
'~ , £1), I
~ , ,'I), I
FR E 04~
F=R~ Ol+g
FGE 05,)
FBE J51
FR F OS?
F8E 051
FGE 054
F P. F 055
FKF 0'56
r8E 057
FP [- 0118
Fe E OS>:}
FBE ~6,)
FR E 061
FB F:: 062
Fe:: () 6"',
rp F C;,~4
FRE '):) S
FJ3 E 061,
rr~F Of) 7
Ff1E 06>3
FBE 069
FB E 070
FRE 071
F f) E 07?
Ff3E 071
fl~ f 074
F!3r- 075
F £3 r= 076
F!:.E 077
FBf.: 078
FH E 07'-1
FGf 0'-30
FRF ORl
FP, F., OR?
FH[ f)rn
PIt IJP,4
Fr<~ O~')
FP [: ORf,
FPE 087
FRE ()F3H
f: r.\ !: f) ,0. ,~
r: ~\ E ,"")qn
FBf O
-------
F:\!f')
F'3 F ("IQS
FEE f) 96
F')t 097
Fi~ F OYR
F~~ n{'Q
'. 'I .
nu: 1')()
P.i~ 1')1
F r~ c 102
F f~ t:: ]01
F f\ F: 1.04
Hi!= 1')5
FP': IOn
F')f ] 07
FP,f= 10Q
F3 E ] ('9
f:J: ([HC I':' Ir:
1("117:: 1= 1,1:0
.J f'" I T t- ( f), 1 4- (' I D I <) T K :. ( I , 1 ), (l 1ST J r. ( I , ,J ), J = ''''' ri } ,
"~ I? T T [ ( 7 , 1 4 ::;) D I ,) T '< C ( ! , 1 I, ( r) J <) r:;J r ( ! , J ), J = ~,~},
U::: I T ~ ( h, 1 I. J) ('"1 r ST.:' ( ( T . J ). J = CJ, 1 ~J }, T
',I ") ! -; ~ ( 7 , 1 4 ') } ( ) i S T ':> r ( T . I), .J -:: (I, 1 . i ,.. T ( 1 ~ j , 5 F 1 C' . (t , ? ~i X , T -~ )
1 1+ 7 r: (n '.'1', T ( ') r- u: . 1+ , ;( ~ '< , r:l, )
STflf1
-------
APPENDIX
EPA CONTROL COST ESTIMATES
-------
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 source s covered in these estimates are: >:'
Industrial Process
Grain Milling and Handling
Kraft Pulp
Nitric Acid
Sulfuric Acid
Phosphate
Petroleum Refining & Storage
Asphalt Batching
Cement
Gray Iron Foundries
Iron and Steel
Primary Copper
Primary Lead
Primary Zinc
Primary Aluminum
Secondary Non-ferrous Metals
Stationary Combustion
Electric Power
Other Stational Combustion
Solid Waste
SIC2042
SIC2611
SIC2819
SIC2871
SIC2911
SIC2951
SIC3241
SIC3321
SIC3323
SIC3331
SIC3332
SIC3333
SIC3334
SIC3341
>:' See Table 3.3, Volume I of this report.
140
-------
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 OBE.
Thus:
N
C = C
ijt jt
V"t
-2L
yN
jt
where C.. is control cost of industry j in year t in ith AQCR
1Jt
N
C. is control cost of industry j in year t in the nation
Jt
V"t is the value-added of industry j in year t in ith AQCR
1J
VN.
1S 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 ViJ't/y~) has been used,
Jt
because, in some AQCRs, 4- or 3-digit information is not available
to the public. >.'<
):' There is a disclosure problem 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.
141
-------
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.
142
------- |