NBSIR 75-911

A National Network Data

Base System

Richard H. F. Jackson

Applied Mathematics Division

Institute for Basic Standards

Washington, D. C. 20234

September 1975

Technical Report to:

Office of the Secretary

U.S. Department of Transportation

Bureau of Economic Analysis

U.S. Department of Commerce

Office of the Chief of Engineers

U.S. Department of the Army

U S. DEPARTMENT OF COMMERCE

NATIONAL BUREAU OF STANDARDS

NBSIR 75-911

A NATIONAL NETWORK DATA

BASE SYSTEM

Richard H. F. Jackson

Applied Mathematics Division

Institute for Basic Standards

Washington, D. C. 20234

September 1975

Technical Report to:

Office of the Secretary

U.S. Department of Transportation

Bureau of Economic Analysis

U.S. Department of Commerce

Office of the Chief of Engineers

U.S. Department of the Army

U.S. DEPARTMENT OF COMMERCE, Rogers C. B. Morton, Secretary

James A. Baker, III, Under Secretary

Dr. Betsy Ancker-Johnson, Assistant Secretary for Science and Technology

NATIONAL BUREAU OF STANDARDS, Ernest Ambler, Acting Director

ABSTRACT

This report documents a data base "system" created at the National

Bureau of Standards. The National Network Data Base System (NNDBS) provides information on intercity flows of freight and passengers throughout the United States. It consists of a set of FORTRAN programs written

for the NB'S computer (but transportable) and some basic data tapes. In

addition to providing basic data on the modal transportation netv7orks,

and interregional flews of passengers and freight, the NNDBS can make

estimates of market splits among modes, can perform aggregations over certain "zones" in the U.S., and is capable of easy extension to other

uses. This report is intended as a user's guide and includes discussions of the data tapes and each of the programs. Complete listings

and tape formats are also included.

Keywords: Commodity flows, data base, freight mode choice,

market split, modal choice, modal freight flows,

national transport analyses, regional transport

analysis, transportation network

TABLE OF CONTENTS

1. INTRODUCTION 1

2. THE BASIC DATA TAPES--NET70, NET80, AND NET90 4

3. THE MODAL SPLIT PROGRAMS— SPLIT 6

4. THE STATE AGGREGATION PROGRAMS— STAG .10

5. THE BEA ZONE AGGREGATION PROGRAMS--BAG 13

6. CONCLUSION 15

7. BIBLIOGRAPHY 16

8. APPENDIX A: COMPUTER FORMAT OF NET70, NET80, AND NET90 17

9. APPENDIX B: COMPUTER FOR^T OF THE PASSENGER TAPE— NETPS 18

10. APPENDIX C: COMPUTER FORMAT OF THE ALTERNATE IMPEDANCE TAPE— NETTS . . .19

11. APPENDIX D: COMPUTER FORMAT OF THE SPLIT TAPES— SPL70,SPL80, AMD SPL90 .20

•12. APPENDIX E: COMPUTER FORMAT OF THE STAG TAPES— ST670,STG80, AND STG90. .21

13. APPENDIX F: COMPUTER FORiW OF THE BAG TAPES— BAG70,BAG80, AND BAG90 . .22

14. APPENDIX G: LISTING OF THE SPLIT PROGRAMS 23

15. APPENDIX H: LISTING OF THE STAG PROGRAMS 28

16. APPENDIX I: LISTING OF THE BAG PROGRAMS 35

LIST OF FIGURES

FIGURE 1 : THE TAPES AND PROGRAMS OF THE NATIONAL NETWORK DATA BASE SYSTEM . . 2

FIGURE 2: THE MATRIX OF POSSIBLE CONNECTIONS 7

FIGURE 3: NNS ZONE TO STATE MAPPING 11

1. INTROD'JCTIO:

This report docurrents the National Network Data Base System (mNDBS)

frorri a prc;»'5'—er''user po^nt of view. The m!;j3S was created at the

National B-^rea^ c~ Stanci'^ds in order to provide information on flows of

freight and passengc-;-s cn ine rajor intercity transportation network of

the United States, "-"e •c-sscs for developing such a data base, the

rationale used in zr::.:'' : Z''e basic data, and the checks used to

assure validity, are all given in [4J. The present report aims to

provide documentation sufficient to permit ready use of the Data Base

system by any competent crogramimer. It could, therefore, be viewed

as a user's guide to the !,',j3S.

The !;,','D3S is properly called a ''system" because some of the data

available from it are not stored explicitly, but can be obtained by

applying one or more of its programis to one of its basic data tapes.

The alternative to sucn a "systematized data base'' would be actual

storage of all "he various manifestations of the basic data tapes. This

would result in at least 12 separate and nearly-full reels of magnetic

tape. It was felt that the creation and use of so large a conventional

data base would involve too great a risk of crioDiing physical tape problems such as unrecoverable parity errors wnich could completely destroy

the usefulness of a given data base caoe. Although reliability for such

a conventional data base can be obtained through maintenance of comolete

listings and backup tapes, the required costs, not only in terms of datastorage, but also in personnel time to maintain the massive structure,

appeared too high to justify such an approach. It was therefore decided

to create a data base 'system" in the sense described above. This ao- proach would keep at a minimum tne number o^ taoes stored, thus keeping

at a minimum the number of backup tapes and complete listings. The programs are kept on mass storage files.

A pictorial representation of the N'N'DBS is given in Figure 1. From

that representation, it can quickly be seen tnat all tapes are derivable

from the single tape called :;ET70, so that just NET70 and a number of

program.s could be keot as "the system". However, it was decided to keep

three tapes, :;ET70, NETSO, and NET9j as the basic data tapes in the

system 'or a number of reasons. These include the fact_that the program

FPvATAR (see d. 223 of [5]), which produces NETSO ami Nli90 from NET70 by .

a technique known as Fratar orojection (see [2]), stas written for the

IBM 350, and that it need only be run once for eac# tape (NETSO and

NET90) in Question. The cost of converting FRATAR to be acceptable to

the NBS UN I VAC 11 05, and of updating it for use in the NNDBS was_deemied

higher than the extra two-tape-maintenance-costs ircurred by having the

three basic data tapes.

Before mroving on to a discussion of the data m the tapes and to a presentation' of the programs, we note that there ar? two other tape

files which do not appear in Figure 1: NETiS and ^ETPS. NEiTS is a test tape which contains an "alternative network" t3- the one represented

in the other data base tapes. More on NEiTS can be found in the

1

2

discussion of the program SPLIT, and its format is given in Appendix C.

NETPS contains information on intercity passenger movements throughout

the national transportation system. Nothing more than storage is cur- rently planned for NETPS, so that nothing more about it will be presented in the present document except for its computer format, given in

Appendix B.

The programs in the NNDBS were written in FORTRAN for the NBS

UNIVAC 1108,"^ but none are machine dependent. Some will, however, re- quire modification before being transported to another computer. The

modifications required are all related to block transfer of data and occur only in the driver programs discussed in Sections 3, 4, and 5.

Some form of binary block transfer of data should be substituted for

each reference to the utility package NTRAN found in those programs.

Section 2 of this report presents a detailed discussion of the

basic data tapes in the NNDBS, while Sections 3, 4, and 5 take up the

programs used to move from level to level in Figure 1. There are six

Appendices giving tape formats, and three more containing listings of

the programs.

To be most effective, reading of Sections 2, 3, 4, and 5 should be

accompanied by regular reference to the formats and listings given in

the Appendices. A word here on the organization of the listings in

Appendices G, H, and I, is therefore appropriate. Each of those three

appendices contains three separate listings: the first is of the control

stream required for a run; the second, of the driver program in the run;

and the third, of the respective "workhorse" subroutine.

*A1 though program development took place on the NBS computer (a UNIVAC

1108), no endorsement of UNIVAC or UNIVAC equipment is intended. .

2. THE BASIC DATA TAPES--';ET70 , :,'ET80, AND :;ET90

This section discusses the basic data tanes ','ET70, NET80, and

NET90. Since NET80 and :JET90 are directly derivable fro~ :,'ET70 and

their conouter fotTiat is exactly the same, the discussion will center

around NET70. Differences between the 1970 tape (,NET70) and the 1980

and 1990 tapes (>'iET80 and NET90, respectively) are noted at the end of

the section.

NET70 contains information on freight movements throughout the

U. S. and its possessions for the year 1970. The information is stored

as flows anong certain "zones'' in the U. S. These zones, called National fietwork Sinulation (NNS) zones, were mapped out previous to this

project. (For more on N',S zones, see [3, 4].) There 'are 512 such NNS

zones, and each record on NET70 contains data on one-way freight movements between a particular pair of zones. All the data stored on ,NET70

falls into two categories: amount of goods shipped, and the impedances

(distance and time) involved in shipping it.

The total flow of freight fro*^ one zone to another has been broken

down into eight commodity classes, measured in tens:

1) large, high value shicments,

2) small, high value shipments,

3) large, medium value shipments,

4) small, medium value shipments,

5) large, low value shipments,

6) small, low value shipments,

7) large, crude oil shipments,

8) large, refined oil shipments.

Therefore, a tape record contains an entry for each of the above com- modity classes, representing the average daily flow of goods'^ from the

origin NNS zone to the destination NNS zone, of that commodity class an

size. The identification numbers of the origin and destination NNS zones also appear in the record.

The interzonal imoedances are associated with the various modes of

transportation. The NND3S deals with six possible modes of shipping a commodity from one zone to another:

1) truck,

2) rail,

3 ) air,

4) water,

5) crude oil pipeline,

6) refined oil pipeline.

The flovr of goods is neasured in tons.

Each tape record contains an entry for each of the above six modes,

which entry represents the distance in statute miles by each mode

between the zone pair for that record. The record also contains time

entries, namely the times (in minutes) required to traverse each of the

above distances via the associated mode. Each record on NET70, therefore, contains 22 entries: the origin-destination NNS zone numbers, 6

distance entries, 5 time entries, and 8 commodity-class quantity en- tries .

The tape contains one record for every zonal pair that is connected

and has a nonzero flow. When the data tape was first created, the NNS zones were considered to be completely connected--a link was presumed to

exist between each possible pairing. However, it was soon discovered

that many such links had zero flow for all commodity classes. A culling

process was then performed to remove these "zero" records from the tape.*

This cut the number of records on the tape file from 512*512=262,144 to

only 68,227, permitting the complete file to reside on one tape, and con- siderably reducing the time and effort involved for data manipulations.

Some time and distance entries on the tape have been set to a "flag

value" of 99999. This means that the connection via the associated transportation mode for that particular zone pair either does not exist or is so expensive to use as not to be a serious consideration. It is important

to bear in mind that one flagged value does not completely remove a link

from consideration in the "zero flow" sense described earlier. There are

six modes of transportation available between pairs. So long as even one

of these is not flagged, the link is still a valid connection.

NET80 and NET90 are copies of NET70, but with all flows expanded by

Fratar projection so that they are predictions of 1980 and 1990 move- ment levels. The time and distance entries are the same on all three

tapes, the only differences being in the commodity class flow entries.

Appendix A contains a succinct description of the computer format

of the three tapes discussed in this section.

*This step is not irreversible. Zones could be reconnected, with a positi

flow, by use of a simple edit program.

3. THE MODAL-SPLIT PROGRAMS— SPLIT

SPLIT is a subroutine whose inputs are the records from NET70 and

whose outputs are records in which the commodity flows have been allocated

among the available modes of transportation. In other v^ords, whereas

NET70 contains what must be shipped from zone to zone, SPL70 (the output

tape from a modal -split run) shows how much of each commodity will be

shipped by each mode.

The basic idea in a modal split is to compute a set of "weights"

u.., i=l,2,...,8, and j=l,2,...,6, satisfying

6

E u. . = 1 , for all i

j=l

These weights can then be used to distribute commodity i among the possible modes j. The model used to compute the weights in subroutine

SPLIT is given by

\'n 'V] '>

where d. and t. are the distance and time associated with transportation

mode j,'^c. .( .)'^is a cost function dependent on distance, and a, 3, and p are parameters, and e is the base of the natural logarithms. The weights

are obtained by normalizing the co.

for all i

More on the general concept of modal splitting, this particular choice

of a model and the cal ibrations^ performed to choose appropriate values

of a, B, and p can be found in [4]. This brief account is included

here only to provide the programmer/user with background information

on the theoretical foundation for subroutine SPLIT.

The key to understanding the inner workings of SPLIT, is to recall

that the NNDBS involves 8 commodity classes and 6 modes of transportation, and to note that not all modes are available to all commodities.

For example, the only commodity . that can be transported via the crude

oil pipeline is, of course, crude oil. There are other combinations

that are theoretically possible, but in practice are not used because

cost is prohibitive; an example is sending crude oil by air. The com- plete picture of possible connections is given in Figure 2, where a "0"

indicates a usage which is impossible (in the sense outlined above) and

a "1" indicates a valid connection.

6

Transportation Modes

Commodity

Classes Truck Rail Air Water

Crude

Pi pel ine

Refined

Pi pel ine

High large (*) 1 0 1 0 0

High small 1 1 0 0 0

Medium large 1 0 1 0 0

MpH"iiim c:mri 1

1 icu 1 uiii oiiia 1 1 1 0 n n\j

Low large 1 0 1 0 0

Low small 1 0 0 0 0

Crude large 0 0 0 1 1 0

Refined large 1 1 0 1 0 1

h "High" (or " jOmmodity class;

medium" or "large" or

"low") refers to dollar-value

"small" refers to shipment

per unit

size.

for the

FIGURE 2: THE MATRIX OF POSSIBLE CONNECTIONS

li

To facilitate program organization, it was decided to store lj^. • as

a one-dimensional "packed" list (OMEGA) corresponding to the matrix^

given in figure 2, but with all elements whose value is zero removed

from consideration. Since the same storage scheme was applied to a. ., 3.. and p.., there was no problem in matching the correct values of

triDse arr^s v/ith u.. for a modal -sol it computation. However, with this

storage scheme, the Simple row-colunn indexing that would yield com- modity class and transportation mode is not available. That information

is needed to sum across modes for a given commodity, and to retrieve

correct distances for flag-value checking.

The row indexing required in SPLIT is provided by the array COMImD

(see the listing in Appendix G). For any 1 = 1,22, COMI'.'D (I) is the

commodity class reference (row number in Figure 2) for the corresponding

element I of the arrays OMEGA, ALPHA, BETA, and RHO.

It turns out that transportation-mode indexing (column number in

Figure 2) is not required in SPLIT. What js_ required is retrieval, for

each element of OMEGA, of the distance to be traveled by the commodity

quantity in question. This retrieval information is provided by the

array DSTIND (again, see the listing in Appendix G). For 1 = 1,22, DSTTID

(I) is an index into the input record (array IR) where the desired

distance can be found.

The last indexing array to be explained is the one represented by

CSTIND, which is an indexing scheme for the cost functions. For each of

the pipeline refined, pipeline crude, water and air modes of transportation, there is a cost function to be applied to all comimodities

shipped by that mode. Rail and truck, however, have three cost functions, the applicable one depending on whether the shipment is high

value, medium value, or low value, according to a prescribed classification scheme. This yields a total of 10 cost functions, and these

are stored in the array COST. Therefore CSTL'^D (I), for each 1 = 1,22, is an index into COST for retrieving the correct cost function to be used

in the modal -split computation.

It must be expected that cost functions will change over time.

Moreover, their number can also change (e.g., a disaggregation of the

"water" mode might be used for some studies). Given the indexing scheme

for the cost functions outlined above, implementation of any such

changes would be simple and straightforward.

The actual work of "splitting" the flows is now straightforward.

SPLIT works on one record at a time from NET70, NET80 or NET90. First

it computes the costs to be stored in COST, which are based on the

interzonal travel distances by each mode. It then launches into the

computation of the weights >'elying on the indexing presented above.

Having done so, it splits eaihl commodity-class flow among the available

modes by multiplying by the weights. An output record is then created

following the computer format given in Aopendix D, and that record is

passed back to the main program, along with control.

8

The main driving program v.'hich controls the modal-split run is

concerned with reading in fron the input network tape (NET70, iJETSO,

NET90) , mounted on logical unit 7, a block of 100 records, and passing

each record off to subroutine SPLIT. After processing is complete for

the block of 100, the driver writes the ''solitted'' version of the

records onto the outcut SPLIT tape {SPL70, SPLSO, or SPL90) , mounted on logical unit 8, and goes back to read in the next block of 100 records.

The data blocks are transferred by a UNIVAC 1108 software oackage (MTPsAN)

which performs binary transfers of blocks of any size.

The printer output from a modal split run consists of a table

giving accumulations by ton, ton-mile, and revenues, plus a second table

giving sums over all fields of the output tape. Both tables are useful

for checking purposes.

There is another capability that has been provided in the driver

program for SPLIT. It is an option to read an alternate set of impedances from a tape mounted on logical unit 9. This alternate tape

could be uETTS (see Appendix C), or another tape with similar format.

This capability is useful for testing and experimentati on purposes . It

could, for example, be used in discovering the effect of adding and

deleting links from the network. To make use of this capability one simply creates the "alternate impedance tape'', mounts it on unit 9 along

with the normal inpui and output network tapes on units 7 and 8, and

changes the value of the logical variable ALT in the driver program from

.FALSE, to .TRUE.

Note that it is not necessary to carry the comolete set of impedances on the alternate impeaance lape. One may desire to change only

the highway impedances or only those for rail, or some combination. The

point is that the alternate rape need only carry information that is

changed. Program changes, to accommodaie such partial modifications of

the impedances, are required only in lines 19 and 20 of the driver . program for SPLIT. For examole, suopose it is desired to change the

rail impedances, but keep all ethers fixed. One then creates an alternate tape with 68,227 logical records, each one containing four

words: the origin 'iNS zone, tne destination NNS zcne, the (new) rail

distance, and the (new) rail time. Lines 19 and 2C of the driver program would then be changed from the existing DO loop to the following:

IR(5,J) = AR(3,J)

IR(6,J) = AR(4,J)

9

4. THE STATE AGGREGATION PROGRAMS—STAG

STAG is a subroutine whose inputs are records from one of the

"split tapes" (SPL70, SPL80, or SPL90), and whose outputs are a set of

matrices showing: a) the average daily flow of goods by commodity-class

among the 50 states, the District of Columbia, and the U. S. possessions; and b) the costs of transporting that amount of goods.

Just as with SPLIT, STAG is the major routine which does the aggregation on individual logical records that have been passed to it by

the driver routine. The driver does all the tape handling and the

breakdown of the blocks to logical records. Since STAG'S driver is

similar to SPLIT'S, discussed in Section 3, the rest of the discussion

in this section will be confined to subroutine STAG.

The key to understanding STAG lies in Figure 3, which shows which

NNS zones are mapped into which states. For example, the top line

indicates that all flows originating from NNS zones 1 through 5 should

be aggregated to represent a total flow originating from Maine. The

meaning of the "partial zones" column is straightforward: the specified

percentages of the flows for the given zones are to be aggregated to

make up the flow to or from the associated state. The need for this

column arises because some NNS zones cross state boundaries. The percentages are based on total population figures and were obtained from

Allan Olson of the Bureau of Economic Analysis, U.S. Dept. of Commerce.

The manner in which the mapping scheme is implemented in STAG,

however, is the inverse of the mapping described in Figure 3. That is,

for a given I ranging from 1 to 512 (the NNS zone numbers that appear on the NNDBS tapes), NDX(I) is the state number to which the flows from NNS zone I should be added. If NDX(I) is negative, ABS(NDX{I)) is a pointer

to a row in the arrays PCT and NPCT which respectively contain the

percentages and the states to which the percentage flows are to be

added.

There are 10 matrices that are created by STAG. They represent,

respectively, the volumes of flow and the cost of transporting that

flow, among the states for the following situations:

1. ) high value shipments, both large and small,

2. ) medium value shipments, both large and small,

3. ) low value shipments, both large and small,

4. ) all crude oil shipments, and'

5. ) all refined oil shipments.

As each record is passed into STAG by the driver, the amount to be

added to each one of the 10 matrices is computed in turn and control is

10

STATE

WUMDER STATE WHOLE NNS ZONES PARTIAL NHS ZONES

Maine 1-5 Incl 2 Kpw K^mD 6-9 Incl

3 VpfTTinn t 1 C 1 lllvl 1 W 10-11 Incl Mass. 12-26 incl 5 R. Ii 27 25X of 28

Conn. 29-40 incl 75X of 28 .

7 N. Y. 41-50 Incl; 52-58 incl 8 K.J. 59-65 incl 20X of 66; 212 of 69; 12% of 84

9 Pcnn

51, 67. 68, 70-83 incl B0% of 66; 79% of 69 10 Del. 85 77J of 84

^1 Hd. 86-89 incl 11X of 84; 43% of 90 12 D.C. 28X of 90

13 Hich. 149-167 incl 16% of 113

14 Ohio 91-97 incl; 100-105 incl; 107-112 incl;114 58% of 98; 81% of 106; 84% of 113; 44% of 99; 22%

15 1 ndi ana 115-123 incl; 125-130 incl 3S of 106; 82% of 124; 16% of 205

16 n 1 inoi s 131-134 incl; 136-146 incl; 148 54% of 135; 23% of 323

17 168-181 incl 78% of 296

18 Minn

295, 297-301 incl • 22% of 296; 39% of 323

19 Iowa 302. 304-314 incl 46% of 135; 89% of 303; 16% of 343 20 Ml ^^ouri 315-317 incl; 319-322 incl 652 of 318; 77% of 323 . 21 N. Dak. 324-327 incl en of 328

22 S, Dak. 329-333 incl 23 Ncbras ka 334-342 incl 11% of 303; 84% of 343

24 ICdnsa s 344-352 incl 35% of 318

25 Vi rginia 182-194 incl 292 of 90

26 W. Va. 195-199 incl 42% of 98; 56% of 99; 57% of 200

27 201-204 incl; 205 16% of 106; 18% of 124; 21% of 200; 84% of 205 28 Tcnnes sse 207; 209-211 incl; 213-215 incl 94% of 208; 83% of 212

29 N, Car. 215-229 incl 30 S. Car. 230-235 incl 36% of 239

31 GsoTQ 1 a 241 ;23e-238 incl;240;242-244 incl;246-i:47 incl

17% of 212; 64% of 239; 31% of 245

32 r 1 ui 1 u 248-260 incl 33 Alp hpm;* 261-271 incl 19Z of 245

34 Hiss. 272-275 incl 35 276-285 incl 36 W\t Arlcan KB 1 1 J<i etQ ^O 286-291 incl; 293-294 incl 6X of 208; 65% of 292; 40% of 362 37 vh 1 O llv'llQ 353-351 incl 35% of 292 38 Tsxas 363-399 incl; 412 60% of 362 39 N. Mex. 406-411 incl 40 Ai* i Trt n;* 400-405 incl

41 Montana 413-420 incl 42 Ida ho 421-425 incl 43 WvruT) wjf v*ii inn Illy 426-429 incl 44 - Crt 1 nrA Hn 437-445 incl 45 Utah 430-436 incl 46 Wash. 446-452 incl 14% of 455

47 Wl Cy <J( 1 453-454 incl; 455-462 incl 86% of 455 48 485-490 incl 49 Cal If

463-479 incl; 480-484 incl; en MiasKa 497-505 incl

51 Hawaii 506-509 incl S2 Possessions 510-512 Incl

FIGURE 3: MAPPING FROM NNS ZONE TO STATE.

n

passed to the internal subroutine ADD. ADD does all the work of stepping through the indexing scheme so that the additions are made to the

proper element of each matrix. When each of the 10 matrices has been

properly adjusted for the record at hand, control is returned to the

driver routine which gets the next record to be processed.

When the last record from the input tape file is reached, STAG

closes out the run by computing and printing out row and column sums for

each matrix. The driver then writes the matrices onto logical unit 8

according to the format given in Appendix E.

12

5. THE BEA ZONE AGGREGATION PROGRAMS—BAG

The BAG programs are very similar in form and content to the programs discussed in Section 4. In these programs, however, the aggregation is by BEA zone (see [3, 4] for a discussion of BEA zones), rather than

state, and there are 16 matrices (described belov/) computed rather than

just 10. These two changes required several alterations to the basic

programs written to handle the state aggregations. One such change had

to do with storage: there are 174 BEA zones, so that the 16 matrices

(with row and column sums) require 16*175*175 = 490,000 words just for

data storage. Since that is more than the NBS UNIVAC 1108 can store

"in core", the programs were modified to compute one matrix at a time

with one pass through the data tape, then to rewind and compute another,

and so on.

Another necessary modification concerned the structure of what were

called in STAG the POT and NPCT tables. Recall that these tables provided the indexing and percentages for those zones that crossed state

boundaries. The same situation obtains here, but in BAG, the two tables

have been combined and more information has been stored in the resultant

table.

The following discussion refers to two arrays: PCT and IPCT.

These two arrays have been "equivalenced: so that both integer and real

values can be stored in adjacent words of the same array. Letting I be an NNS zone number, then as before NDX(I) is the BEA zone number to

which the flows to and from NNS zone I should be added. But if NDX(I) < 0,

then letting IAD=-NDX(I), we now have that IPCT(IAD) contains two numbers "packed" into the same word. The leftmost number, retrieved by

letting NL=IPCT( IAD)/1 000, is a pointer back into NDX. The rightmost

number, retrieved by letting NR=IPCT( IAD)-NL*1 000, is the number of BEA zones that are intersected by NNS zone I. There will be NR groups of 3

words each in the PCT table following word IAD. In each group, the

third word is the zone number of one of the BEA zones intersected by NNS zone I. The second word is the proportion of total earnings of NNS zone I that accrued from that BEA zone, and the first word is the proportion of total income of NNS zone I that accrued from the BEA zone. The income proportion is used in the aggregation when the NNS zone in

question appears as a destination zone in the input record from the

split tape. The earnings proportion is used when the zone appears as ari origin in that record.

The DATA statements in BAG that initialize the PCT table were

created by a set of programs not included in the NNDBS. Documentation

of those programs, and the source of the earnings and income proportions

used in the DATA statements, appears in reference [3].

13

The subroutine BAG works similarly to STAG. A logical record is

passed into it, volume or costs are computed and aggregated, and the

internal subroutine ADD is called to step through the indexing scheme

and add the quantities to the proper element of the matrix currently

being computed. The number (1-16) of the matrix being computed is

passed in the calling sequence in the variable M, and the matrix itself

is passed in BMTX.

The matrices computed correspond respectively to the volume and

cost matrices for the same eight commodity classes listed in Section 2.

The only major difference between the driver program for BAG and

that for STAG is that the computation of each matrix in BAG requires one complete pass through the tape. After each matrix is computed, its row

and column sums are computed, it is dumped out in binary form on unit 8,

the input tape (unit 7) is rewound, and the next phase of matrix com- putation is begun.

The computer format of the output tape from a run of BAG is given

in Appendix F.

14

6. CONCLUSION

The National Network Data Base System was designed with one main

goal in mind: to provide a maximum of useful data on the inter-city

freight and passenger flows along the major transportation network of

the United States, while keeping the number of data words actually

stored to a minimum. This approach also permits a great deal of flexibility in data-base usage: a capability to change easily the modal

split equation, to revise cost functions, to aggregate over a slightly

modified index set, to modify the network, and so on.

In fact, this ability to exercise the data-base system and to

derive a set of data different from that built into it has already been

used successfully. A request was made to compute state and BEA zone aggregations for 1972, a year for which no flow data existed in the data

base. What was provided, however, was a set of state and BEA zone growth rates from 1970 to 1972, based on total earnings and total income

for the states and BEA zones. The request was satisfied with ease by

• applying the given growth rates to the row and column sums of the matrices

already produced in the system, normalizing the new sums with respect

to the "expanded grand total" of the matrices, and then using a standard

matrix scaling algorithm (*) to adjust the elements to the new sums.

This is one of the many possible examples of the flexibility of

such a data base "system". Other uses are anticipated.

X*) The program (SCALE) is a p&rt of NBS's Mathematical Programming

Algorithms collection and is based on the procedure outlined on

p. 115 of [1].

15

7. BIBLIOGRAPHY

1. Deming, V/. E., Statistical Adjustment of Data , John Wiley and

Sons, New York, 1938.

2. Fratar, T. J., "Vehicular Trip Distribution by Successive Approximations", Traffic Quarterly , Vol. VIII, No. 1, January 1954,

pp. 53-65.

3. Lozier, D. W., and Stiehler, J. R. , "Machine-Readable Mapping from

National Netv/ork Simulation (NNS) Zones to Bureau of Economic Analysis

(BEA) Areas", NBSIR 75-918, September, 1975.

4. Schofer, R. , "Final Report on the National Network Simulation",

NBSIR 75-912. (in preparation)

5. U. S. Department of Transportation, FHWA Computer Programs for

Urban Transportation Planning , July, 1974.

16

APPENDIX A: COMPUTER FORMAT OF NET70, NET80, AND NET90

DENSITY: 555 f.p.i.

PARITY: even TRACKS: 7 FORMAT: binary

LOGICAL RECORD LENGTH: 22 words

BLOCKING FACTOR: 100

NUMBER OF LOGICAL RECORDS: 68,227

RECORD INTERPRETATION:

WORD MEANING SYMBOL

1

1 Origin NNS zone 0

<>c Destination NNS zone D

5 Distance by truck (mi) 1 u 4 • Time by truck (min.

)

TT

5 Distance by rail RD

6 Time by rail RT

.7 Distance by air AD

8 Time by air AT

9 Distance by v;ater WD

10 Time by water WT

11 Distance by crude pipeline PCD

12 Time by crude pipeline PCT

13 Distance by refined pipeline PRD

14 Time by refined pipeline PRT

15 Large, high value shipment volume (tons) HL

16 Small, high value shipment volume (tons) HS

17 Large, medium value shipment volume (tons) ML

18 Small, medium value shipment volume (tons) MS

19 Large, low value shipment volume (tons) LL

20 Small, low value shipment volume (tons) LS

21 Large, crude oil shipment volume (tons) CL

22 Large, refined oil shipment volume (tons) RL

17

APPENDIX B: COMPUTER FORMAT OF THE PASSENGER TAPE--NETPS

DENSITY: 556 f.p.i.

PARITY: even TRACKS: 7 FORMAT: binary

LOGICAL RECORD LENGTH: 14 words

BLOCKING FACTOR: 100

NUMBER OF LOGICAL RECORDS: 54,128

RECORD INTERPRETATION:

WORD MEANING

1 Origin NNS zone number

2 Destination NNS zone number

3 Distance by auto (mi)

4 Time by auto (min.

)

5 Number of daily auto passengers, non-business purposes

6 Number of daily auto passengers, business purposes

7 Distance by air (mi

)

8 Time by air (min.

)

9 Number of daily air passengers, non-business purposes

10 Number of daily air passengers, business purposes

11 Distance by rail (mi)

12 Time by rai 1 (min.

)

13 Number of daily rail passengers, non-business purposes

14 Number of daily rail passengers, business purposes.

18

APPENDIX C: COMPUTER FORMAT OF THE ALTERNATE IMPEDANCE TAPE— NETTS

DENSITY: 556 f.p.i.

PARITY: even TRACKS: 7 FORMAT: binary

LOGICAL RECORD LENGTH: 14

BLOCKING FACTOR: 100

NUMBER OF LOGICAL RECORDS: 68,227

RECORD INTERPRETATION: Exactly the same as the first 14 words listed

in Appendix A.

19

APPENDIX D: COMPUTER FORI'iAT OF THE SPLIT

TAPES--SPL70, SPL80, AND SPL90

DENSITY: 800 f.p.i.

PARITY: even TRACKS: 7 FORMAT: binary

LOGICAL RECORD LENGTH: 37 words

BLOCKING FACTOR: 100

NUMBER OF LOGICAL RECORDS: 68,227

RECORD INTERPRETATION:

WORD MEANING SYMBOL

1 Origin NiNS zone 0

2 Destination i,NS zone . D

3 Distance by truck (r.i) TD

4 Time by truck (rr,in.) TT

5 Distance by rail ' RD

6 Time by rail RT

7 Distance oy air _ AD

8 Time by air AT

9 Distance by v.'ater WD

10 Time by water VJT 11 Distance by cruae pipeline PCD

12 Time by cruae pipeline PCT

13 Distance by refined pipeline PRD

14 Time by refinec pipeline PRT

15 Tons of nigh value, large shipments sent by truck THL

16 Tons of high value, small shipments sent by truck THS

17 Tons of m.edium value, large snipments sent by truck TML

18 Tons of medium value, small snipments sent by truck TMS

19 Tons of low value, large shipments sent by truck TLL

20 Tons of low value, small shipments sent oy truck TLS

21 Tons of large, refined oil snipments sent by truck TRL

22 Tons of nigh value, large shipments sent by rail RHL

23 Tons of hign value, small shipments sent by rail RHS

24 Tons of medium value, large shipments sent by rail RML

25 Tons of medium value, small shipments s-ent by rail RMS

25 Tons of low value, large si"^,ipments sent by rail _ RLL

27 Tons of low value, small shipmierrts sent by rail ' RLS

28 Tons of large, refined oil shipments sent by rail RRL

29 Tons of nicn value, small snipments sent by air AHS

30 Tons of high value, large snipments sent by waier VJHL 31 Tons of medium value, large shipments sent by water Wi'IL 32 Tons of low value, large shipments sent by water . WMS

33 Tons of large, crude oil shipments sent by water WCL

34 Tons of large, refinec oil snipments sent by water WRL

35 Tons of large, crude oil shipments sent by pipeline PCL

36 Tons of large, refined oil snipments sent by pipeline PRL

37 Sum over fields 15 through 37

20

APPENDIX E: COMPUTER FORMAT OF THE

STAG TAPES--STG70, STG80, AND STG90

DENSITY: 556 f.p.i.

PARITY: even TRACKS: 7

FORMAT: binary , LOGICAL RECORD LENGTH: 2809

BLOCKING FACTOR: 0

NUMBER OF LOGICAL RECORDS: 10

RECORD INTERPRETATION: Each record contains one of the 53 x 53 matrices

created by STAG, written out column by column. They are, in order of

appearance on the tape, state-to-state aggregations of the following:

1) high value, large and small flows (tons).

2) the transport costs of number 1 (dollars),

3) medium value, large and small flows,

4) the transport costs of number 3,

5) low value, large and small flows,

6) the transport costs of number 4,

7) crude oil flows,

8) the transport costs of number 7,

9) refined oil flows,

10) the transport costs of number 8.

Reference Figure 3 for which state is associated with each row and column

number. The row and column sums appear in the last row and column.

21

APPENDIX F: COMPUTER FORMAT OF THE

BAG TAPES--BAG70, BAG80, AND BAG90

DENSITY: 556 f.p.i.

PARITY: even TRACKS: 7 FORMAT: binary

LOGICAL RECORD LENGTH: 30625

BLOCKING FACTOR: 0

NUMBER OF LOGICAL RECORDS: 16

RECORD INTERPRETATION: Each record contains one of the 175 x 175

matrices created by BAG, written out in column order. The matrices are,

in order of appearance on the tape, BEA zone to BEA zone aggregations of

the following:

1.) high value, large flows (tons)

2.) the transport costs of number 1 (dollars),

3.) high value, small flows.

4.) the transport costs of number 3,

5.) medium value, large flows,

6.) the transport costs of number 5,

7.) medium value, small flows.

8.

)

the transport costs of number 7,

9.

)

low value, large flows,

10.) the transport costs of number 7,

11.) low value, small flows.

12,) the transport costs of number 9,

13.) crude oil flows,

14.) the transport costs of number 13,

15.) refined oil flows.

16.) the transport costs of number 15.

22

APPENDIX G: LISTING OF THE

SPLIT PROGRAMS

23

LlSTiMG OF" CONTROL STREAM FOR A RUM OF 5PL

(3RUN,N/R R I CJAK , 3^ 1 ^ I - JACKSn , BE A , 1 , I 00 , n2nnr,

0MSG,W 3613 1-JACK50. I'-'RITE FNAPLE REEL SPL70,

K)H56THANK5....RTC

0A5G , T JM 7 . , (I

, NET70

f3ASG,TJH a , , II , '^PL70'V

|?ASG,AXSUP.

SASG.AX DRV.

Td'FOR.S ORV.'^^LIT,TPF<f;.A

0FOP,S SIIR.SPLIT,TPES.R

QFREE SUR,

RlFRFE DRV.

QXQT

BtHARGE

END ON SITE PRINTOUT ON OCTOPFR 6i 1975 AT I3:M?:'^9

BEA'RUN ( 1 ) . SPL I T ( S

)

i

LISTING OF DRIVER FOR SUbROUTK^E SPLii

PARAMETER MIR=22 r M0R=37 » MAR=1

IMPLICIT INTEGER (A-Z)

LOGICAL ALT

DIMENSION IR(MlR,iOO) rOR(MOR»100) »S(MOR) »AR(MARf 100)

DATA ALT/. FALSE./

WRlTE(6r6)

C**** READ A BLOCK FROM THE INPUT TAPE.

iO CALL NTRAN(7>2»M!R*100# IR(1» 1) f Lp22:)

XF (LeLT.O) GO TO 60

C**** IF r:£CESSARY» READ BLOCK FROM ALTERNATE TAPE.

IF (ALT) CALL NTRAN(9»2»MAR*100f ARdr 1 ) f N»22)

IF (N.LT.O) GO TO 100

-C**** BEGIN LOOP WHICH PASSES OFF EACH LOGICAL RECORD* IN TURNr

C»*** TO THE SUBROUTINE FOR PROCESSING.

DO 50 JsirlOO

IF {.NOT. ALT) GO TO 30

IF (IR(i» J) cNE.ARdp J) .OR. IR(2» J) .NE.AR(2» J) ) GO TO 110

C**** IF ALTERNATE TAPE UP AND RECORDS MATCHr TRANSFER THE IMPEDANCES,

DO 20 islfMAR

20 IRCIf J)=AR(If J)

30 Csc*l

• CALL SPLIT(IR(1» J) f MiRf 0R(1» J) f MOR)

C**** IF FIRST WORD IS 2ER0» LAST BLOCK ON TAPE HAS BEEN REACHED.

IF (IR(l»J).EQ.O) GO TO 900

C**«* SUM ALL FIELDS OF OUTPUT TAPE.

DO t;0 I = lfMOR

tfO S(I)=S(I)+OR(I» J)

50 CONTINUE

C**** WRITE A BLOCK TO THE OUTPUT TAPE.

60 CALL NTRAN{8»l»MOR*100f OR(lfl) »M»22)

IF (M.LT.O) GO TO 90

DO 70 I=1»M0R

DO 70 J=lfl00

70 OR(I»J)=0

GO TO 10

80 WRITE(6»2) L»C

60 TO 900

90 WRITE(6»1> MfC

GO TO 900

100 WRITE(6»7) N»C

GO TO 900

110 WRlTE(6r8)

C**** WRITE OUT LAST BLOCK AND ENDFILE.

900 CALL NTRAN(8»l»MOR*100#OR(l»l)fMf22)

CALL NTRAN(8r9)

WRITE(6»'+) C

WRITE(6r5)

WRITE(6#3) (I»S(I) »I=lfMOR)

1 FORMATC *« ERROR Mr'ISr* Cr»I5)

2 FORHAT(» ** ERROR L=*I5»» C=»I5)

3 F0RMAT(2XI5» 115)

FCRMAT(/1XI10» LOGICAL RECORDS WERE TRANSFERRED.')

5 F0RMAT(//» SUM OVER ALL OUTPUT TAPE FIELDS ///2X • FIELD 12X SUM

)

6 F0RMAT(»1«)

7 FORMATC ** ERROR N='I5»« C=»I5)

8 FORMATC ** ERROR C=»I5»«TAPE RECORDS DO NOT MATCH. )

END

LISTING OF SUbROUTINE SPLIT

SUBROUTINE SPLlT( XRf MIR»OR»MOR)

LOGICAL NFIRST

INTEGER OR » COM I ND > CST IND » DST IND f TYPE t SUMNO r SUMTN » SUMTM

INTEGER SUMRV

DIMENSION IR(MlR) fOR(MOR) »ALPHA(22) r BETA ( 22 ) » PHO ( 2? ) »0MEGA(22)

DIMENSION COST(IO) ^CSTIND (22) »DSTIND(22) f C0MlND(22) rSUM(8)

DIMENSION TYPE(22) »SUMNO ( 22) » SUMTN (22) »SUMTM(22) »SUMRV(22)

DIP4ENSI0N R(3)

DATA NFIRST/«FALSE./

DATA RH0/7e812'*f I5,i4762f8« 1321^10 •1584 f6o988^f» 8,08fi8p 19.0570 »1 .Of

* I,0pie09l,0»l»0»l«0»i.0t-7,tl0fl7> .80l9»l*t^2t^2r".0638»

* l,Or ,768»8e623rl.8/

DATA BETA/- • 9299 » -3. 258«* r 9299f 869f«l,<^726»". 869^-2. 785 r-. 2192

» -aalOt*^ -.2192 1 -9 2288> -•7732 2288 »-1,1186p-1.5626»

* ^~^-.*fa88r-.t*888f -•7t|98>-t 24^66 r -.9052 » -.7928 8^*66/

DATA ALPHA/-e29t4M.#-i.*l29f^?-.29«+4»-.8*f58f-e53i^5p-.8t*58»-1.3691»

* OcO»-.7922»0.0»-.5756»-o2251»-.5756»-.6010p-.1332»

* 0892 0892 1184 P-e 7255 »«. 3275 r -I. 0826 9 -.5963/

DATA TYPE/»THL» f 'THS' » «TML» » 'TMS* t ^TLL' > 'TLS' f 'TRL» r 'RHL*

* 'RHS* r »RML» f 'RMS* # 'RLL* » 'RLS* p »RRL» f » AHS» » «WHL» r * 'WML* » »WLL» f 'WCL* t n/RL' f 'PCL' r 'PRL'/

DATA C0MIND/lr2r3?t*f 5p6p8»l»2f 3>'+p5p6»8r2pip3»5»7»8»7,8/

DATA CSTIND/lf Ip2»2p3f 3p3f t^»4p5»5»6»6p6r7p8p8»8»8»8r9rlO/

DATA DSTIND/3p3f3p3f3»3p3F5»5p5p5r5»5>5f7»9»9»9»9»9»ll»13/

CHECK IF THIS IS LAST LOGICAL RECORD.

IF (IR(l).EQ.O) GO TO 800

IF (NFIRST) GO TO 100

C***« IF FIRST TIME THROUGH* RESET RHO.

DO 10 I=lf22

10 RHO(I)-EXP(RHO(I))

NFlRSTs.TRUE.

C**** COMPUTE COSTS (PENNIES/TON). C(l)-C(10) ARE* IN ORDER:

C**** CTHfCTM»CTLfCRHpCRM»CRL»CA»CW»CCfCR.

100 C0ST(1) = .8817*(-1216. 0+921^. 0*IR(3)**.«I21)

COST(2) = .**961*(-1216.0->92t|.0*IR(3>**.t*21)

C0ST(3) = .U168*(-1216.0+92'+.0*IR(3)**.'»21)

C0ST(t*)=718.+2.96^IR(5)

COST(5)=3'+1. + 1.70*IR(5)

C0ST(6)=21.*t+l#'+7*IR(5)

COST(7)=105'*0. + 18.92*IR(7)

C0ST(8)=296.5+.6t*t^*IR(9)

C0ST(9)=25.8'l+.232tf«IR(ll)

COST(10) = .t|li^8*lR(13)

R(l>=13**9. 5251+3. 5187860*IR( 3)

R(2)=992.«^3256+3.1580756*IR(3)

R(3)=768.85476+2.6590308*IR(3)

DO 150 1=1*3

COST(I)=MlN(MAX(COST(I)fR(I) )»1.17*R(I))

150 CONTINUE

C**** COMPUTE THE UNNORMALIZED WEIGHTS.

200 DO 210 I=:l»22

OMEGA(I)=0.0

JlrDSTINDd)

IF (IR(Jl).GE. 99999. OR. IR(J1+1).GE. 99999) GO TO 210

J2=CSTIND(I)

IF (C0ST(J2) .EQ.0.b.OR.IR(Jl+l) .EQ.O) GO TO 210

OMEGA ( I ) =RHO ( I ) *COST ( J2 ) **BETA ( I ) *FLOAT ( IR ( Jl+1 ) ) **ALPHA ( I

)

210 CONTINUE

LISTING OF SUBROUTINE SPLIT

C  SUM THE WEIGHTS ACROSS MODE TO GET NORMALIZING FACTOR.

DO 310 I=lf8

310 SUM(I)=0,0

DO 520 1=1 »22

JsCOMINDd)

320 SUM(J)=SUM{ J)'f'OMEGA(I)

C** *« NORMALIZE THE WEIGHTS.

DO tJ20 I"lf22

JsCOMXNDd)

IF {SUM(J) .GT.O) GO TO f^iO

OMEGA(I)=0»0

GO TO «*20

<fld 0ME6A(I)=0M£GA(I)/SUM(J)

t^ZO CONTINUE

C***» SPLIT THE SHIPMENTS BY MULTIPLYING THEM BY THE WEIGHTS.

DO 510 I=l»22

jsC0MIND<I)"8-i<*

510 OMEGA(I)=OMEGA(I)*FLOAT(IR{J) ) \

C**** BUILD THE OUTPUT RECORD.

DO 610 lsl,l<*

610 OR(l)=IR(I>

OR(37>SO

DO 620 1=1 »22

0R{I+l«*)=0MEGA(I)'i-.5

620 0R(57)=0R(37)+0R(I-l-ltv)

C**** SUM THE TONSp T0N-MILES» AND REVENUES.

DO 720 I=i»22

JisDSTINDd)

IF <IR(J1) .GE. 99999) GO TO 720

SUMTM < I ) =SUMTM ( I ) +0R ( I + 1 tfr ) * I R ( Jl

)

J2=CSTIND{I)

SUMRV ( I ) =SUMRV ( I ) -e-OR ( I + ii^) *COST ( J2)

SUMTN ( I ) sSUMTN ( I ) +0R ( I + lU

)

SUMN0(I)=SUMN0(I)+1

720 CONTINUE

GO TO 999

C**** OUTPUT THE TABLE OF SUMS.

800 WRITE(6»1)

WRITE(6»2> (TYPE(I) »SUMNO(I) >SUMTN<I) rSUMTMd) tSUMRVd) »I = 1»22)

1 FORMAT (//12Xf » TOTAL 7X « TOTAL 11X» TOTAL 10X» TOTAL' /IX 'TYPE* 7X

* • COUNT « ex » TONS 7X« TON-MILES »7X» REVENUES'//)

2 F0RMAT(2XA3»I12,I12#I16#I15)

999 RETURN

END

27

APPENDIX H: LISTING OF THE

STAG PROGRAMS

28

LISTING OF CONTROL STREAM FOR A RUN OF STAG

QRUNfM/R RICJAK»3613i-JACKSO#BEAp5plOO

QASG^aX DRV.

oasgpax sub,

qasg»a stout,

QASGfTJH 7*fUrSPL70

QUSE SorSTOUT*

QFORpS DRV*STAG»TPFS,A

OFORf»S SUB,STA6»TPFSeB

©FREE DRV,

QFREE SUB,

OXQT

PCOST

29

LISTING OF DRIVER FOR SUBROUTINE STAG

PARAMETER MIR=37»MVD=53

INTEGER C

DIMENSION IR(MIR»100)

DIMENSION V0LS(MVD»MVD»5) rD0LS(MVD»MVD»5)

C**** READ A BLOCK FROM INPUT TAPE.

10 CALL NTRAN(7r2fMIR*100pIR(lpl> rL»22)

IF (L.LToO) GO TO 60

C**** BEGIN LOOP TO PASS OFF EACH LOGICAL RECORD TO SUBROUTINE,

DO 50 J=l»100

CALL STAG{IR(l?J)pVOLSrDOLS)

IF (IRd^J) cEQ»0) GO TO 100

C=C+1

50 CONTINUE

60 TO 10

80 WRITE(6»2) L»C

C**** WRITE OUT AGGREGATED MATRICES ONTO UNIT 8. . 100 DO 110 K=l»5

CALL NTRAN(8rl»MVD*MVDfV0LS{lrlrK) »M»22)

110 CALL NTRAN(8rl»MVD*MVD»D0LS(lrl»K) fM»22)

CALL NTRAN(8»9)

C**** PRINT OUT ROW AND COLUMN SUMS/

DO 150 K=l»5

WRITE(6r3) K

WRITE (6 ft^) (VOLS(I»MVD»K) rI=l»MVD)

WRITE(6»5) K

WRITE (6» 4) (VOLS(MVDr I»K) »I=1»MVD)

WRlTE(6r6) K

WRITE (6»^f) (DOLS(IfMVD»K) f I = 1»MVD)

WRlTE<6r7) K

150 WRITE{6#^> (DOLS(MVD»I»K) irlfMVD)

2 FORMAT(» ERROR L=»I5»» C=»I5)

3 FORMAT(/t ROW SUMS FOR VOLUME MATRIX NUMBER' 13/)

tv FORMAT(10(E13,7) ) 5 FORMAT(/« COLUMN SUMS FOR VOLUME MATRIX NUMBER* 13/)

6 FORMATl/' ROW SUMS FOR COST MATRIX NUMBER* 13/)

7 FORMAT {/ COLUMN SUMS FOR COST MATRIX NUMBER* 13/)

END

30

LISTING Pr Sl.'D^nUTp.'E STAG

SURPOUTTiF 5 T A G ( I ^ , V 0 L S , 0 ^ L 5 ) D I M E N S I 0 '

;

V 0 L S f 5 3 , 5 3 , - ' » D ^ L S ( 3 , 5 3 ,

/ D I '-'E "5 I C

1 C / 7 7 1 D r T / 7 A 3 ) , '\ T ( 26 , 7 1 • 1 P

[) A T A 1 ( C T ( I , ,

K 1 r> ^ T / ] 1 ) , J = 1,3) ) 1 = 1 1 o 1

11 • 2b, 7 r

• / b , 6 « » C ,

e. • ? 0 1 ^ % oT , P P » c

,

7J • 2 1 1 3 , • / V , V , » 0 ,

li 1 "5 3 , t / / , 1 U » • 1 I , 1 1 , C-3 • ^ J j

' ' ) 7 Q \ 7 7 9 . 25 , LO

C Q

• by? 1 ZJ 1 » U 7 • ' - » 7 A

^ 6 , P 0 . 0 . 7 / 1 /I • " ^ J 1 ^ » • Z3 o f 7 A ^ 6 > P P n

O 1 1 u r » - o ,

1 • ' ~- , 27 , Q • 1 6 J 1 -5 1 Q U • 3 J n Pi

1 3 » % 1 c 1 ^ , • 1 , 7 7 A / 1 . 0

1J » b T ?

1 L 1 6 , t ^ c , p P

? -

"7 • *- - > . c; 7 AO, 7 A . 7 1 . 27 , » 1 A 1 O J

1 r

' - J , R ii

, A7 7/ ? •

!" r 1 0 » 7 a •

. 1JO,A A ,00 , n , cJ> D 7 ^ - J

1 7 • 1 / , 7 1 ,

fi n » n

,

Lo • 36, A i-i 7 1 , n , 7 7 » 7 7

-- J > P p , n

,

ao • 6b, 7 / - 6 , 7 CI \ 7 P P > 0

,

Q "7 Q • 7 y , 1 / , 7 7 1-1 1 R P P p /

J A 1 A ( t o r T >

1 p r T /

) , J = 1,3) » I = 7 P 7 A Q O • J J , 7 , P P

t1 • o b , 7 7 c / 7 , , P P , p

,

7 7

1 6 , 7 7 7 'I — - 1 P 1 c

,

J

7 O • 37, A 1 • o I , 7 1 A 1 , P P

q

1 / • 1 c ,

1 o H U • '!} ', 7 7 A J » ^ p

» 0 , n n • ' - J 7 / On, A P- JO, 7 R P P , n ,

t U ' 6 , , R A U7 7/ J , P . c/

I 7 6 7 8 Q7

11 U A 1 H (

i ! A ( T 1 T — 1 , 1 _ ) /

tI 1 , 1 ,

1 1

1 J

! ,

1 t 7 7 , 2 , 7

3 , n7 , M , y , 7 ^ , M , H , M , M , 5 , - 1 » 6 6 » a ^ > ^ 1 6 , 6 ,

/o , P , 6 , 6 , 6 cD 7 , 7 , / , 7 , 7 , 7 7 , 7 , 7 7

zo 9 , 7 , 7 , 7 , 7 , 7/

, 7 , 7 , R 7 8 , p > Q 0 P , _— 7/ , 9 , 9 , - 3

rS 9 , 9 , oV - oV . 9 , Q Q

, 9 , 9 Q , o 9 , o 9 , - M , 1 0 , 1 1 5 1 1 ,

I 1 ,

1 1 — ,

1 ^4

, 1 ^ , 1 a , 14, 1 y I > - 6 , - 7 1 U

1I 1

' ^ 1 1 5 M »

1 u _ a

! » 1 ,

1 Ll

1 ^ I Ll M » I » -

» ! » 1 b , 1 ' »

1 5 , 1 5 ,

1 5 1 ^ » 15, 15, 15, -in. 1 c

1 S , ! 5 ,

J r

1 b 1 r 16, 16, 16, 1 1 »

t A 1 6 ,

! 6 , ! 6 I A , 16, 16, 1 A 16, 1 A t A 1 6 , ! 6 ,

1 7 1 7 O 13, 1 3 , 1 3 , 13, 1 3 ,

I 7 1 3 , 1 3 ,

\ 7

1 i 1 7 7 13, 13, 13, 13, ! 3 ,

1 7

1 - 1 1 3 ,

' 7 ,

1 7 I '

1 7 17, 1 7 , 17, 17, 1 7 ,

! 7 ,1 7 ,! 7 ,

1 7

1 '

1 7 9 17, 25 , 7=^ , 25 , 25 , ' 5 , 25 , 25 , Z b 7 C / / P / DATA ( MO ¥ { I ) , 1 = 1 9 1 , 3 B n ) / 25, 2'=', 25, 25, 76 , 2 6 , 2 ^ , 2 6 , 26 1 7

1 27 , 27 , 27 , ?7 , -13, 2 7 , 7 R ,

.- I , 2*^ 7 R A ^ , •y

e- 28,-15, 2P, 2^, 7 R , 2^ , 29 , 7 9 , 7 9 7 n 3 29 , 2 9, 79, 79, 2 9 , 2'' , 2 9 , 2 9 , 20 P » M 30, 30, 3n, 30, 3 0 , 3 1 . 3 1 , 3 1 , - ! 6 3 1 , 31
-A
LISTINC, or SUBROUTINE STAG

C • « • *

3 1 » 3 1 3 1 » 3 1 - 1 7 , 3 1 » 3 1, 3 ? , 3 2 3 7 1 7 . 3 2 ) 3 ? i 3 ? 3 7 » 3 ? 1 1 7

.

3 ' , 7 1 1 -5 - »

•5 -1 J J »

-i -i -5 J

7 7 7 7 7 -J 3 3 , aa -7 ^ s 3 S . 3 S s . 9 3 5 I 3 ^ » 3 5 » 3 S 1 3

Cj 3 A » 3 6, 3 6, 3 6, 07 A0 .

J 3 6 1 1 a ,-19 , 18, t ft . 1 ft .

1

1 9 11 97 t Q ,1

1 97 . 1

! p7 ,

1 P

£. 1 9 - i 0 1 7 n 7

'->, 7 0 - 7 1 7 0 /C7 r1 ,

-i•J ^ <i J 7 1

t /- 1 7 1 — 7 7 "? 7 * - ) 7 7 n 7 7. 7 7 9 9. 7 3 c 0 7 3 7 3 7 77.

^ -7 1 7 - - 77 .f 7 M 7 M 7 4 7 4 , 7 4 , 7 4.

o ' 1 / ^ »

-) 7 -5 ' J-1 7 7 7

•1 7 7 7 , 0 '

, 7 7 7 7 7 7 3 ' , /

-7 ~ / T » - " » ^ "

,5 0

7" D 7 Q

, -5 ri

, J n , 3 ft , Qo

"i Q 3 H f ^ Q J Q T n

1 - 7 0 Jo, 7 Q 70/ J H /

[) A T A XNil/ ( I ) ^ ft I . S 1 7 1 9 ft . 3 ft . 3 ft 3 ft f 3 8, » ft . ft . ^ 0 f 0 . •J \

"K ft . ft 7 ft 3 ft ft .

1

1 • J »

ti

I n'J I 'IP, ^ n 7 Q 1 0 7 , -J ^ ; 7 P 7 0 J T , 0 0V J 4 1

T 1 ) ^ 1 ^4 1 J 'I 4 1 , 11. a 1 4 1 ~ < •

f 7 J

* ^ 1 M 7 M 7 ^4 7 q 3 ' -J " 7 J

1 J JI - V J M S q C, q q C n T >

^4 4 ^ 4 'Of 4 A , LO "Of 1 0 J ^ 7 a 7 — 9/COA 44 7 7 U 7 u 7 il 7 7 it 7 U 7 " 7

J " 7 q 9 , i4

' Q7 , u 0' * U P ' ^ f oo M17Q 1 L| 0 M 9 q 9 i4 0 f 17 , ' ' 4 9 ^4 P Q u' 0' f U P' ' > 4 9 " 7 M n 4 ft "Of 4 0 « n J n

0 » 0 » 0 , sn , sn

,

sn , SO ,

1 SO , so

,

so

,

sn , sn , S 1 , SI, 1 » S I , S 7 , Z 5 ? /

I F { 1 R ( 1 ) . EO GO TO 1 on

CTH=.ftft|7-*(-17l6.0 + '774,r«T'=(3)**.47!)

CTM=.4^61«(-1716.^+9?4.0«IR(3)»«.4?1)

CTL=.41ftft»(-l?lA,0+924.O*TP(3)*».47l)

RTH=1349.S7Sl+3.5187ftA''*IR(3)

RTM=99?.4 3 2S6+3.IS807SA*IR(3)

RTL = 76" . «S47A + 2 . 6S9P30«* I P/ 3 ) CTH = MI^'(^1AX(CTH,PTH) , ! . 17*RTH)

C T M = M I ^M f

' A X ( C T M , R T f-l ) , I . 1 7 P T ) CTL = MI'"(MAX(CTL,PTL) , I .17«RTL)

CRH=7 18.+7.9A*IR(S)

CRM = 34 l. + 1.7C*Tf?(S)

CRL=?!.4+).47»IR(S)

CA=10S40,+!8.O7»IP(7)

CW=296.5+.644*IR{9)

CPC = ?S.P4-f.?324*IR( ! 1 ) CPR=:,414ft«lR( 13)

C»»*» DO COMPUTATIONS FOR HI r,M VALU"^, lARGF ftMD «^KA1L.

00L={ IR(?3) + IR(2?) )*CRH+( IP( I^) + IP( 16) ) •CTH+Ic{3P)«C'-'+If(29)*CA

V0L=Tf(23)+IR{27)*IR( JS)+1R( 1A)+IR(30)+IR(?9)

CALL A DD ( 1 ) C**«« DO COMPUTATIONS FOR M«^n VALUE, LARGF AMn SMALL.

DOL=( I''(24) + IP(7^) ),cRM+( TPf |7) + i9( ift) )»CTm+ip{3i )«CW

V0L=IR(2 4) + 1R(2S) + IR( I7) + IR( 1P) + IR(3 1 ) CALL AOD ( 2

)

C»»»» DO COMPUTATION'S FOR LO^ VALUE, LAPGF AND Sf-'ALL.

D0L=(I''(2f) + IR(?7))«CRl+(l''{I<=) + IRf20))«CTL + ]f^(32)«C'/.'

V0L = 1P{2A) + IR(77) + IR(19) + IR(2'^) + IR(32)

32

LlSiniC, or sup r.'OiiT T rjr t A

Ce»«* DO CC'-'^IITATIO''^ FOP C^'Jn?:.

D0L=I^M33)'»C''+I^(3B)«C^C

VOL=IP(33) + If?("^^)

C»««» DO CnnPUTATIO'''? FO^" PE'^lN'^n,

D0L=I^(?8)«CRL + 1^(?l)»rTL+lP('»M)«C"+lP(3,'>)«CPP

V0L=IR(2")+I^(21 )+I=(3")+IP(36)

CALL A ( 5 ) GO TO 2nc

10 0D01inK=i,S

D0I1CJ=1,52.

DO 1 1 D I = 1 , B2

V0L5(S3,53,K)=vnLS(?3,F3,^)+vr^LS( I ,J,K)

V0L5{T ,5 3,'<')=VPLS( I ,53,l^)+vnL^( I ,J,K)

V0LS(S3,J,'<')=vnLS(53,J,K)+V0L^( T ,J,<)

D0LS(F3,^3,K)=0OL5(53,S3,.<)+'^nL?(T,J,i^)

DCLSf ! ,?3,l<)=^^'LS( I ,';3,^')+nOL"^( ! ,J,<)

1 IC D0L5(?3, J,K)='-PLS(53,J,i^)+nOL?( I ,J,K)

200 RFTU?^J

C**** SURROUTI'ir TO n'TER^RET 5UBSCPIPTR A^Jn no THE AG'^-^EGATIOf,',

SURROLITIMF ADr^fAi<)

IMPLICIT INTEGER (A)

A 1 = I P ( I ) A I = ^.' n X f A I ) A J = I R ( ? ) A J = NDX ( A J

)

IF(AI)3 0,9r!j,JD

10 IF (AJ) 'S0,9C'",7?

C»**» AGGREGATE FOP CASC '^F ^JO S^'SA'5.

20 VOLS(AI ,AJ,A<)=V0LS(AI ,AJ,A>^WVOL

D0L5(AI ,AJ,AK)=00L5(AI ,AJ,AK)+DOL

G 0 T 0 R o P 30 IF (AJ) cr,9 0'","r,

C»*«* AGGRFG6TF FOR CASE OF ORG AN S^^'SA,

'J0D050A1-=1,3

A2 = NPCT ( - A I , A 1 ) !F (A2.FQ.n) GO TO oqq

V0L5(A2,AJ,A<)=V0L5(A7,AJ,a<)+PCT(-AI ,M )*VOl

?C D0LS(A7,AJ,AK )=^0LS(A2,AJ,AK)4.PCT(-AI ,Al )*r^OL

60 TO 999

T**»« AGGREGi^TE FO^ CASE OF ^FST AN 5*^5A.

60 DO 70 A 1 = 1 ,

3

A2=nPCT(-4j,A!)

1F(A2.FQ.C)G0T099<'

VOLS{M ,A?,AK)=V0L5(AI ,A2,AK5+PCT(-AJ,Al )*VOL

70 D0LS(AI ,A2,AK)=D0L5(AI ,A7,AK)+PCT(-AJ,Al )«DOL

GO TO 99?

C#»*» AGGREGATE FOR CASE OF RPTH O^G AND DF5T AN 5 MSA.

RC DO 90 A 1 = ! ,

3

A2 = NPCT ( - A I , A 1 ) IF(A?,E0,n)G0T0-9PO

D09lIA3=1,3

AH = rjPCT(-AJ,A3)

IF (AM.FO.O) GO TO "0

V0LS(A2,A'4,AK)=vnL5(A7,A'4,AK:)+PCT(-AI,A!)«f'CTf-AJ,A3)«'.'nL

90 D0L5(A7,AM,AK)=nnLS{A2,AM,AK)+PCT(-Al ,Al )«PCT(-AJ,A3)«00L

GO TO ^99

LTSTiwl OF SUftH-OilTlME 5 'l A G

9C0 WR I TE ( 6 , 1 ) A I , A J , AK , T R

1 FOf?MAT(' fTPROP** W { I '? ] \ 0 ) ) 999PETURN

END

END 0N5ITE PPINTOUT Orj OCTCBER 6, 197S AT 13:M8:S7

BEA*SUB(1).5TAG(5)

34

APPENDIX I: LISTING OF THE

BAG PROGRAMS

35

LISTING OF CONTROL STREAM FOR A RUN OF BAG

CRUNrN RICJAKf 36131-JACKSO»BEAr i+0/D20 00 » 125

BASGfAX DRV.

CASGrAX SUB,

OASGrTJH 7.»U*SPL70

QASGrA BAGOUT,

QU5E 8.»BA60UT»

CFORfS DRV.BAG^TPFS.A

GFORfS SUB«BAG»TPF$,B

0FREE DRV.

QFREE SUB.

QXQT

QCOST

36

LISTING OF DRIVER ROUTINE FOR SUBROUTINE BAG

PARAMETER MIR=37» MB=175 » MB1=MB-1

INTEGER C

DIMENSION IR{MIR»100)

DIMENSION B(MBf MB)

C**** LOOP TO PROCESS ALL 16 MATRICES.

DO 300 MAT=lfl6

0=0

C**** READ A BLOCK FROM THE INPUT TAPE.

10 CALL NTRAN(7f2rMIR*100»IR(l»l) »Lf 22)

IF (L.LT.O) GO TO 80

C**** BEGIN LOOP TO PASS OF EACH LOGICAL RECORD TO SUBROUTINE.

DO 50 J=l»100

IF (IR(IfJ) .EQ»0) GO TO 100

C=C+1

50 CALL BAG(IR(1»J)»MAT»B)

60 TO 10

80 WRlTE(6r2) LrC

2 FORMAT (» ** ERROR L=:»I5»» C=»I5)

C  COMPUTE ROW AN COLUMN SUMS.

100 DO 110 I=1#MB1

DO 110 J=1»MB1

B(I»MB)=B(I>MB)+B(I» J)

110 B(MBr J)=B(MB» J)+B(I» J)

B(MBrMB)=0.0

DO 120 I=1»MB1

120 B(MBrMB)=B(MBrMB)+B(I»MB)

C** ** REWIND THE INPUT TAPE AND WRITE OUT THE MATRIX.

CALL NTRAN(7»10)

CALL NTRAN(8rl»30625rBfM»22)

WRITE<6»3) MAT

WRITE(6»1) (B(I»MB) »I=1>MB)

WRITE(6»4) MAT

WRITE(6»1) (B(MB» J) » J=lfMB)

C**** ZERO OUT THE MATRIX FOR THE NEXT PASS.

DO 200 I=1»MB

DO 200 J=1»MB

200 B(I»J)=0.0

1 FORMAT(10E13.7)

3 FORMAT(/» ROW SUMS FOR MATRIX NUMBER' 13/)

I* FORMAT(/» COL SUMS FOR MATRIX NUMBER* 13/)

300 CONTINUE

END

37

LISTING OF SDRROUTTME RAG

SURROIITIMF BAG(IR,M,PMTX)

DIME M 510 M IR(37) ,PCT( 19! 3) ,IFCT( 1913)

DIMFMSIOM fjnx { S ! 2 ) , PMT X ( 1 7^ , 1 7B ) EOUIVALFMCF (PCT.TPCT)

DATA ( M n X { I ) . T = 1 9 J • , HO)/

I 1 » - 1 2 9 2 , 2 , - P , 3 , -15, q ,

, 3 «9 2 - 2B » ^ q , q , q , q , q , q , q , q' ,9

«-? —' - M 7^ .

9 - M 9 , q , -56 , q , -63 , - 7 0 5 . • ^ 9 u1 c

1 1 t; •9 — 77. / / , ij 9 - P M , t M . 1 q /

A T A T - q 1 t ^ -» / /

t

I —~ 9 1 7 I A A . C> 9 7 . 7 . 7 . 8 9 9 9 9 J ye.. 1 n - 9 R . 17. J ^ , 1 '.0 J A . 1^9 J '4 . 1 1

.

1 q , 1 1 1 q . >J 1 q - t 1 R 1 c

* - 9

1 c. . 1 !j J 1 V) , I J 9 1 0 f1 7 . -12 5.

q t1 o7 t 3 . — 13 7. 1 1 / J -1^7. A A . 6 6 , 6 6 1

1 Aij 1 6 . 5 1 6 , - 1 A 2 , I 6 , 1 5 9 1 7 9 1 7 9 1 7 , - 1 7 2 » 1 8 I 1 R , AU /> 8 1 6 7 9 6 P , 6 8, A R . -17 9. 6 6 ,

> A A ,

, -18 6,

7 6 M , A 7 . 6 2. 6 3 1 A 3 ,

» -217,

6 9 - 2 2 7 » 7 0 , -237 , - 2 M M , 7 ^ 9 - 7 5 1 , 6 1 1 6 1 , -7^7,

9 - 7 7 7 » 6 C , - 2 R 7 , 55 , 5 6 , - 300 , 5 9 , - 3 1 0 , 7 6 » 77,

» 7 7 » " ? » -317, -330 , 7 9 , - 3 '^ 0 , 7 R , 7 8 1 -350 I 5 8,

1 B 7 , 5 7 , 5 7 , -373 , - 3 R 0 ,

1 1 M . - 390 ^ -MHO I - M 1 0 , 2 7 M , 74 , 7 1 , 70 9 7 1 , 7 1 , -M 20 , 7 2 » 7 2 , 7 1,

3 7 , -M27

,

7 3 , 7 3 , - M 3 M , M 1 , -MMR , -M"^ 5 , -462, -469,

q -'18 2 f » 8 M , RM , 8 H , -M 95 , P 3 , -502 , -^12, -522 , -S32 » 1 7 » -539 , 2 2 , 2 2 , - 5 ^4 9 , 2 I , -556 , - 5 A 3 > -573,

2C > 20 , -5Rn , 5 ) , -590 ,

I *' 9 -597 , -6*^7 , 5 2 1 5 2 , 7 — oA ^9 nu — A 3 — A q 7 . C, q — 656 , - A A 9 . q A . — A ft 7 q' 9 9 — A 9 — A O . - 7 n A J 0 .

9 — 7 1 — 7 7 n . 7 7. -7 77, 7 A - 7 q 0 , QT — 7 c, r 7 1 7' — 9 — 7 A r — 7 A 7 '7 J7 9 — ^ / -i 7 q - 7 R 1 /

A T A f M n Y f T ^ T — 7 > 1 , 1 - / / 1 q 7 r> » /

1I — 7 Q 1 — 7 Q P 7 0 ^ " 9 — p 1 1 7 P — R 7 7 q q"9 — R0 ,T7 q~ > 7 » J ^ — P q q o — n c: 7 q 7' ^ 9 — C 7 O " ^ 9

<4 1

' I , — 0 R n 7 R — P R 7 3 q t J ' 9

«. n o u J r> f J ' 9 _ 0 r> <i .> 9 1 '4 > A 7 A -.3 "i A0 — 971 > J7 97 J q• — 9 2 R - 9 3 R . q n . • - 9 — 9 c, 731,1 q C "9 — 9AM. q 7 q 5

.

- 9 7 1 « Q P 1 , _ 9 9 q I ' 9 -10 0 7, 1 R . - 1 p 7 3 . IIP. - 1 P 3 0 . 13 9.

A 1 3 o _ J n q n , 1 -> J 9 —10 4 7. 1 7 . - 1 1 A 1 . 117. 117,

7 1 1 R .

•- 1 0 R R , R 7 . R 7 , - 1 inn. -117 1, -113 1, 9 1

9

g 1 n 3 . — 1 1V A0 1 . 1 , - 1 1 A Pt 1 1 0 n ,

1 O^ . R I .

9 -I 1 7 R . - 1 19 1 . 1 1 7 1 1 8 0 , 9 .-1714. IDA. — 1 7 7 7 . — 1 7 q 111. Ill, 111. Ill, 111. Ill' — 1 7 q q . 1 1 A - 1 ^ 7 - 1 7 A ^ . 1III, 1 q - 9 . - 1 7 7 q — 1 7 R 7 . - 1 7 9 q . c 7 1 n p , - 1 3 0 1 ,

1• - - 1 3 1 1 . 99. 1 0 7 . -13 7 1 . - 1 1 2 R . _ I 3 3 p 10 2, - 1 3 q 5 .

1 (J o $ 1 -) > -> f

1 n 7 1 n 0 — 1 A 7 1 Jo-', 111 — 1

"» 7 7 1

1 !

1 n'J — 171 J 00 7r ' J\ 1 n 9 1 1 9 — 1

1 -JVC", 1 9 A ^ 1 q 0 1 — 1^19 1 7 n — 1 q T 7 1 7 1 — 1 ^4 9

q 1 I .

9 1^19 - 1 q q A 1 .1 * / ,

1 q 1 .

1 ^ • 9 - t q A A . — 1 q R c.

1 I -J J

t q *4

.

9 1 1 ' 9

1 q q 1 q 7 — 1 n c;

1 ~ u r> ,

11 q1 7^ ,1

1 '7 9' 9 -1^17. 1 o 1 ^ I 1 7 R 1 ^7 7' 9

f, 1 27 9 » 9 -15 2 7. _ 1 c, 7 9 . 1 7 1 — 1 c, 1 . 1 7 7 .

1 ' <i 9 17 7, 17 7, - 1 5 A q . 7 1 2 3 . - 1 c, 7 1 . f

I » / i f

1 7 q 1 7 q _ j q Q f)

,

1 7 A . - 1 A 0 , 1 7 7 1 A 2 . Q - 1 A n 7 1 A 3 . 1 A 7 110/, A 7 . 1 A 7 — 1lot', A 1 q . — 1 A 7 1 — 1 A 7 T . - 1 A 7 CI

1 0 J

~ 1 q 6

,

9 I'D 1 ' ' ' • - 1 A q 9 i4 7 " , — 1 A C 9 — 1 A A A Q C, 9 Q — 1 A 7 9 / DATA ( NO y ( I ) . T = q 7 1 . ^ 1 7 1 /

1

1 - tI Ao Ro Ao - 1 A 9 J n ^ 0 $ — 1 7 1 A — 1 7 7 9 — 1 7 •> A — 1 7 q 7

1 ' ' J 9 — 1 7 c n — 1 7 Q 7 1It)!, ^ 1 2 1 J 1 9

1 1 — 1 7 A 7 — 1

1 7' 7' q' ,

1

' q' 9 9 — 1 7 q , 1

1 q1 R- , 1 q R 3 -1797 9 M7 ,

i H7 , - 1 ROq , 1 M9 , - 1 R 1 1 ,1 5M , - I 8 1 R , - 1 R25

,

1^'=>, - 1 B3S 9 1

''3 ^ 9 157, 1 "^7 , 1 57 , - I PM2

,

1 56 ,

1 r p

, 1 5R , 158,

5 1 -^7 , - 1 RM9 , - 1 , - 1 R63

,

1 6« , 17 1, 1*1, 167, - 1 R 7 3 , 6 -1886 9 17 1, » 7 1 , 1 6^ , - 1 R93 ,

I 66 , -1900, 1 A5 , 1 65 , 165,

7 1 65 9 165, 1 6M , 1 A5 , 16 1, 16 1, 160, 1^0, 1 60 , - i.''n7

,

fl n 9 0 , 0 , 0, 0 , 0 9 172, 1 '2 , 172, 1 7 2 , 38
? I-
LTSTIMG OF SURROIJTINE PAG

9 1 12 172, 1 72, 172, 172, 1 • 1 7M , 1 7M/

DATA ( PCT ( I ) T = 1 , 7 ) /002002 ,

1 .0792 , . OA 26 , 0'"' 1 , 2. 9?0M , . 937H , 002/

DATA ( PCT ( I ) »

I = 8 ,

1 4 ) /006no2

,

1.2723 , . 3296 , 00 2 , 2.7277 , . 6 7 0 M , 004/

DATA ( PC T ( I ) » I = 1 5 , 24 ) /008r03

,

1 . 093S ,.0755, 0 3 , 2 . 1 B M ^ » . 1 6 '4 1 , 005 , 3. 7220 , . 76nq , 004 /

DATA [ PC T ( I ) . I = 25 , 34 ) /O 1 1003 ,

1 . 1 900 ,.1700 , 006 , 2 . 6 1 OG , . 5900 , 003 , 3 .2000 , . 2 ') 0 0 , 005/

DATA PCT { I

)

. I = 35 , 4 1 ) / 0 2 1 002

,

1.566^, . 5 M 1 0 , 00w 5,—'7 2 . M 33 1 , . M590 , 0^4 /

DATA PC T ( I ) , I = 42 , 48 ) /023f^02 ,

1 . 807R , . P 2 7 M , O'^ 6 , 2. 1 922

,

.17 26, 005/

DATA PC T ( I ) , I = 4 9 , 55 ) /024n02

,

1 . 007 1 , . r?7r

,

- - » 2.9929, .9930, 005 /

DATA ( P C T { I ) . 1 = 56 , 621/026 0 02,

1 . 0 3 1 r , .0330,

2 . 9A9C , .9670, 004 /

DATA ( PC T ( I ) » I = 63 , 69 1/028002,

1 . fl98M , .8929, 0 0 ^4 , 2.1016,» . 1 J 7 1 , 0'"'5/

DATA 1 PC T ( I ) » I = 70 , 761/0 2900?,

1.2836, . 3609 , C04 , 2. 7 I 6M ,. 639 1 , 0P5 /

DATA ( PCT ( I ) » T = 77 , 83 ) /035'^02 ,

1 . M966 ,, . 505 7 , 2 . 5n3M . M 9 M 3 , 0 14/

DATA PC T ( I ) 8 4 , 90 ) /O38002 ,

I . 1 03M ,.1245, 0 !^ 5 . 2. 896A , . 8755 , 0 14/

DATA PCT ( I

)

* I = 9 1 , 97 1/04 1002,

1.8537",, . 8823 , 0 0 6 , 2. 1 463 ,.1177, 00 7/

DATA PC T ( I ) , I = 98 , 1 07 ) /052003

,

1.3012 , . 3009 , 007 , 2. M25

1

, . 4 6 M 8 , 0 12,

3. 2737 ,.2342, 008 /

DATA PC T ( I ) » I = ! 08 , 117) /054nQ3

,

1 .6756 ,.6715, 0 12,

2.2195 ,.2313, 007 , 3. I 0H9 , .097?, 00 6/

DATA PC T ( I ) t 1 = 118, I 24 ) /06 2P02 , 1.8770 ,.9156, 0 14,

2. 1 230 ,.0044, 0 15/

DATA PCT ( I ) . I = 125, 13 1 ) /070r»02 ,

1 . 72M I , . 7453 , 0 15,

2.2759 , . 2547 , 0 13/

DATA [ PCT f I ) . I = 132, 14 1 1 / O 7 300 3

,

173, 173, 173, 174,

39

LISTING OF SIIBROUTIMF: BAG

I . 3 R90 , . 3 3 ? 0 .012,

2 t S 30 8

,

. 5BP I ,011,

3 . OBO ? , . 0 P 7 9 1 0 1 3 /

DATA ( P C T ( I

1 . 6 6 ? . 6 « 6 6 ,011,

2 . 1 6S 1 , .1579 , 0 A 6 , 3 . 172 1 .1556 ,009/

DATA PCT f I 1,1 =

1 . 3 39" . 3 ?5 1 , 066 , 2 . 297r , . 3n09 (010,

3.363 1 , . 3 7 0 1 0 6 7 /

DATA [ PC T f I 1,1 =

1 . 1 G 6 1 .117 1 ,013,

2. 570!? , . 5903 ,016,

3 . 3239 , . 2926 ,015/

DATA PCT ( I . I =

1 . 3 1 7 fi . 3 5 7 P , 066 , 2. 682? .612! ,017/

DATA PC T ( I I =

1 .9659, . 9 t; q 3 ,068 , 2 . 0 3 M 1 , . i:^M57 , 066 /

DATA PC T f I » T =

! . 3757 , . 3953 , 068 , 2.62^3, .6017 ,070/

DATA ' PCT ( I »

I =

1 . 9623, . 9 ? 7 5 ,061,

2.0377, . 0775 ,066/

DATA ' PC T { I » I =

1 .0 229, . 8521 ,061,

2.1771 , .117 6 , 0 ^ ? /

DATA PCT ( I 1,1 =

1 . 3 9 1 r , . 3 7 1 B ,062 , 2 . 3 5 'I M , . 3676 1 " 6 1 . 3 . 2 5 M 6 , .2606 ,052/

DATA PCT ( I > , I =

I . 265n . 2 B 1 5 ,069 , 2.5555, .5515 , "^6 3 , 3 . 17 9 6 .167*^ ,061/

DATA PC T ( I , I =

I . M9M0 , .15^5 ,068 , 2 . 1 H 7 B » . I 2 M 2 1 6 . 3.3583 ,.1161 ,061/

DATA PC T ( I , I =

1 .5782 , ."^99? ,075,

2 . M 2 1 B , . 1 0 0 P ,070/

DATA PC T ( I , I =

I . 3 ?B 1 , . 2560 1 - -J % 2.6719 , .7110 ,076/

DATA [ PC T ( I 1,1 =

1 .212 2 , . 2 1 6 n 1 0 7 5 , 2. 2397 , . 256 2

3 . 091

B

,.0710 » 059 , M . 389? , .1269 ,060 , 5 . 03H 1 , . 0 3 " 0 ,076/

DATA PC T ( 1 1,1 =

1 . 2026 ,.1^31 ,067 , 2.7171 ,.7416 1 06 1 , 3 .0803 , . 0650 ,060/

112, 15 1) /07H003 ,

152, 1 6 I ) /075003 ,

162, 17 1) /0B2P03

,

172, 178)/0BB0 02,

1 79 , 185) /097002 ,

186 , 1 92 ) / 1 Q0n02

,

193, 199)/) 02002

,

2 00, 206)/101P02,

207 , 7!6)/105003,

217, 2261/110^03,

227 , 236 ) / 1 1 2nC3 ,

237 , 213 ) / 1 1 100?

,

211, 2 5 0/115002,

251, 266)/117"05,

267, 276)/120n03,

40

II-^TI^JG HF S«'PP0UT1ME BAG

DATA ( PC T ( ! . I = 2 7 7 1.718^, . 7 A 3 « Q L n

2 . 1 « 2 6 . 16 0 8 , 0 5

'J

3 . 0 9 9 r . r 7 M 0 6?

DATA ( PCT ( I 1 ^ = 2 8 7

1 . 1 2 R 6 , . c 1 2.638°, . 7 1 2 0 5 5 3 . 0 7 5 P , , P c u 1 0 6

H . 15 6 7, .13 16 0^0

DATA ( PCT ( !

> 1 = 3 0 0

1 . 6 2 M 9 , .6^57 0 5°

2.1313, .^973 05 6 3 . 2 M 3 ^ , . 2 7 0 , 0 6 C DATA PC T ( I » I = 3 1 0

1.15 7 5 . 1 0 5 n c Q

2.8425 .85^5 , 077

DATA PCT ( !

T -

3 1 7

1 . 2763 .20 = 2 2 . 2276 , . 2750 077

3.063^, . 0 « « 5 (08!

'^ . M 3 2 7 , . 'I 3 8 3 0 P 2 DATA ( PC T ( I . I = 330

1.15 0 2, .12 3 6

1-170 2.7717, , 0

I a q 0 7 7

3 . 073! , , P c, A 5 p 7 a DATA PCT t I >

I = 3 u n

1 . 7M5? , . 9027 0 7 8

2 . 17 8 3, .1^3'^ ! 1 3

3 . D7 5 3 , n ^ ^ r ,079

DATA ( PCT ( I > I = 35 0

1.391?, ."18° pen

2 . 6 0 R 9 .5811 ,077

DATA p C T ( I » I = 3 5 7

1 . 1 2 3 .12 7 7 ,056

2 . 2^1 6 3 . 2 5 8 » 1 1 3.3888 , ^ (05^

^.1215 ,.1232 ,055

5 . 1 IT . 0 " 5 P 1 0

' 7 DATA P C T ( I 373

1 . 8^455 .8-455 » 1 1 2 . 1 5M5 . 1 5 5 ,057

DATA P C T ( I , I = 3 8 0

1 . 1 7 6

=5 .17^5 , ! 1 5 2.6528 .6663 . 1 1 3 . 1703 . 1 ^ ^ 3 1 ^ 5 DATA ( PC T ( I » I = 3 9 0

1 . 5059

,

.5415 » M 3 2 . 1 9 , .17^5 > ! 1 3.3092, . 3 3 30 ,^57

DATA PCT ( I 1 I = MOO

1 . 1 6 M f .13 2 0 ,073

2.230°, . 1 p p q , 0 7 M 3 . 6 0 6 , . 6 7 6 0 7 A DATA ( P C T ( I 1 I = y 1 0

1 . 6 1 9 P . 62 3 3 , 0 7

'4 2 . 2 M 3 2 , 2 4 5 M n 0

3 . 137^ .13 13 , 0 7 A / DATA PC T ( I > . I = M 20

2

a 6 )/ 1 2 1 0 0 3 ,

29O)/!2 3"0'4,

309)7126 0 03,

316)/128"C7,

^7C)/13300'^,

33° ) / 1 3'^'^03 ,

3M9)/136-03,

356)/13°^0 2,

3 7 2 ) / I M 1 G 5 ,

379)/l«5^07,

3«^)/lM6 0 03»

3PP)/!'^^'^03,

"19)/15''"'C3,

y7A)/157P07,

41



L 1 SI I fiG OF SUB'^niJT T Ml HAG

DATA ( PCT ( I

. I = 563

1 . 6Q 1 0 , . 7 M 9 n • 019,

2. 2n33

,

. 1 A ,018,

3 . 1 1^7 , . 0 « 1 1 ,071/

DATA ( PCT ( I , T = 573

1 . 3 R M 3 , .3130 . 020 , 2.6157, .6 870 ,019/

DATA ( PCT ( I , 1 = 5 8 0

1 . 0 1 ^4 P , . 0 1 O'J .051,

2.1379, . 1 O't ° ,075 , 3. 8M7 3 , .8078 ,070/

DATA ( ^CT ( I . I = 5 90

1.7292, . 7G92 ,05 2,

2. 270P

,

.2908 ,051/

DATA { PCT ( I 1,1 = 597

! . 8 1 3 A , . 8 3 6 « . '?65 , 2.12AM, .117 1 .052,

3 . O/SO*^ , .046 1 , 066 /

DATA ( PCT ( I , ' = 607

1.3715, .3719 05 2 , 2.0173, ,0170 .065 , 3 . Mn?2 , . 5379

,

06^^,

M . 1 220

,

. ""773 , 0 6 6/

DATA ( °CT ( I , I = A20

1.5070, . 5 4 6 3 0^3,

2.1919, .10 9 9 062 , 3. 1 325 , .16 18 05? , M . 1 » 86 , .10 1^ , 0 5 M /

DATA ( PCT ( I , I = 633

1 . Mn32 , . 39R5 . 050 , 2. 2970 , . ? ^"^6 . 0 - 3 , 3 . 299^^ , .3159 '"^

'-' 2 /

DATA ( PCT ( T , I = AM 3

1 . M 1 2M , . 3P 1 7 0^3,

2 . 1 7M M , .1233 0 4 9,

3 . H 367 , , M 7 « 7 1 5 M , M . 0766 , . C2'^3 ,050/

DATA ( PC^ ( I ,' = 656

1.2890, . 27M^ . M 5 , 2.0606, , 0 3 7 . 05M , 3.2690

,

.3191 . 0 'J 9 , M . 3 8 I 3 , .36^6 ,055/

DATA ( PCT ( I » I = 6 6 9

1 . 077 1 , .0136 -'go

2.8131, .8196 , r q A , 3.06^9, . 05 77 , r 'J 7 , M . 095C , , 1 n 0 r: ,115/

DATA ( PCT ( I ) , I = 682

1 . 890^, , . 9r 2 M , C , 2.1094, .0^76 , 0'J7/

DATA ( PCT ( I » , ! = 6^9

1 . 1 7 M , .1756 » 0 ^-

, 2.6332, . 6 f r, n ,049,

3 . 1957 , .19 3 6 , 050 / DATA ( PCT ( 1 1 . ' = 6 9 9

1 . 9823 , . 9 7 fi

, OM 8 , 2.0177, .0177 ,044/

DATA { PCT ( I > , T = 706

572)/l89 0C3,

579)/l90n0?,

5np ) / 1 93'^03 ,

506 ) / ! 95002 ,

6061/197 103,

619)/l9e00M,

632 1/20 1 OOM ,

642)/? 0 3003,

6'^5)/20Mr'nq,

6 6 8 ) / 70 6 0 0 M ,

6^ 1) /207'^OM »

6« n ) /709ro2

,

6«8 ) /? 1 1003 ,

705 ) /2 1 7n02 ,

717)/213 1". 02,

43

LISTING or SUpROUTpiE

I » f^l \ 7 1 . 6670 » C 4 P , 2.3993 , . 3 330 ,050/

DATA PC T ( T ) , I = 7 13

1 . 2fl 1 'J

, . 250 1 ,050,

2.7196 ,7499 ,051/

DATA P C T ( I ) ,

1 = 7 2n

1 .0687 ,0789 ,049,

2.9313 » • 9 ? I 7 ,027 / DATA P C T ( I ) ,

1 = 7 2 7

! . 7 7 6 M ,7916 ,0 76,

2.141? . 1 " 1 5 , C 75 , 3 .0653 ,0579 ,0 77 , H . 0 1 65 , n n Q r ,0 79 /

DATA PCT ( I ) , I = 7'^0

1 • 5 1 23

,

.5432 1 0 7 6 , 2 • 1 8 M 9 , ,1^61 (0 73,

3 . 302 r>

,

. 2 7 9 ,025/

DATA ( ^CT ( I 1,1= 750

1.672", . 6360 .025 , 2.0430, . 0307 ,070,

3.3342, , 3 3 32 ,07 3 / DATA ( P C T ( I , 1 = 7 6 0

1.8259, . 84 25 0 2 3 , 2.1741, .15 7 5 0 2 2 /

DATA ( PCT ( I , I = 7 6 7

1 . 7 2 5 f. , . 7402 02^ , 2.2745, .2^9^ 0 7 3/

DATA ( P C T ( I , I = 774

1 . 6884, . 6 6 3 2 » 023 , 2.31 16, . 3 3 6 P , 024 /

DATA ( PC T ( I , I = 7 9 1 1.6=52, . 6559 ,030 , 2. 2632

,

. 263 1 . 029 , 3.0817, . 0 « 1 1 ,026/

DATA PCT ( I ) , 1 = 7 9]

1.5707 , 7 5 9 ,031,

2.4293, .57"! ,029 /

DATA PCT ( I 1,1= 799

1.1464, .1526 .032 , 2.5426 ,.5355 ,03 1,

3.2759,, . 2973 .029 , 4.0351 . 0 2 ,033/

DATA PC T ( I , T = 1 1 1.5139 , 5 :> 9 6 • 028 , 2.2032, ,3734 ,03! , 3.1232, . 1 r> P 7 . 076 , 4.0319, . 0726 . 03 7 , 5.nt(7fi , . 0 3 A 2 ,029/

DATA ( PCT ( I , I = 9 2 7 1.7615, . 7 3 65 . 044 , 2. 23S5 . 2635 ,04 8 / DATA P C T ( I , T = 9 3 4 1.2169, .231' .042,

2.6365, , 6 5 ° 7 .044,

3 . 14 6 6 ,0 3 7 / DATA PCT ( I I , I = P 4 4 1.5959 , . 6 r 1 4 ,03 3 , 2.2034 , . 20np • 0 3 2 ,

7i9)/2'ic,rc2,

7 76 ) /7 1 6002 ,

739)/2lPnD4,

749 ) /72C'^03 ,

759 ) /222'^C3 ,

766 ) /725n02 ,

773 ) / 2 2 6^02 ,

790 ) /27 onC2 ,

7 9 ) / 2 3 n 0 3 ,

797)/731'"0 2,

P!0)/232^0 4,

9 7 6)/?34-'0 5,

P 3 3 ) / 2 36O02 ,

P43)/73Pn03,

956)/74r'''^4,

44

LISTUJG OF 5U PRO IJ TIME PAG

^ 1 "7 D 1 1 O 1 9 n ti ? U " i 1 M . 0 ? 2 6 , o 1 lit " 3 M /

DATA ( P C T ( I ) , I = «57 ,

1 .0932, , 0990 , 0 3 3 , 2 . M 5 5 , . M 6 1 1 , 0 1 » 3 . I 70S , . 1 4 M 9 , 0 M 2 , H • 2 8 1 8 , .2950, 0 3'*/

DATA ( P C T ( I ) , ! = 8 70 ,

! • 1 2 , .5120, 0 4 M , 2.1627, • 1 6 5 M , 0 4 2,

3. 2961 , .3726, 0 4 3/

DATA ( P C T ( I \ , T = P 8 G ,

I • 80H 1 , •7^97, 0 4 1,

2 . 1959, .2003, 0 'J 3 /

DATA ( P C T ( I 1 , I = 8 8 7 ,

I • 7 M 2 7 , • ' h , C 3 ^ , 2 . 2 S 7 3 , • 2 H 0 , [ ' 3 t* /

DATA ( P C T ( I , T = 8 9 4 ,

1 • M 1 8 P , . M B 7 6 , 0 3 4 , ? • 1 60 1 , . 1 u 7 7 , 0 3 7 , 3 . M 2 1 1 , • 4 0 M 7 , 0 3 5/

DATA 1 P C T ( I 1 , I = 904 ,

1 .2752 1 .3379, 0 3 5 , 2.5991 , .5170, 0 37 , 3.1357, • 1 7 5 1 , 0 3 6/

DATA ( P C T ( I 1 , I = 9 14,

1 . 7003

,

. 6 M 0 6 , 0 37 , 2.2997

,

. 3 5 ° 4 1 0 3 6/

DATA ( P C T { I , I = "21,

I .6 203, .6060, 0 3 8 , 2 . 3797 , . 3 9 M 0 , 0 3 9/

DATA P C T ( I , I = 9 2 8 ,

1 . 8826, , R o 9 ^4

» 0 0 , 2.031 1 , . 0758 .13 7,

3.0862 , . 0 7

'J 8 (039/

DATA P C T ( I , T = 9 3 8 ,

! • 3909 , .3888 0 " 0 , 2.1987 , . M6« .045,

3 . 1658 , . 1 P I 0 1 1 3 6 , H » 2 '^ 6 , . 2 8 3 ,13 7/

DATA [ P C T ( I 1,1= 95 1,

1 . 5371 , • 5 'J 9 7 ,045 , 2.2 805 , . 2779 , 04 3 , 3.0261 , . 0 ! 7 ,044 , M . 15 6 0 ,.15^5 , 040/

DATA C T ( I ) , I = 9 6 4 ,

1 . 3637 ,.3190 ,045 , 2.6363 , . 6 P 1 0 ,047/

DATA P C T ( I ) , I = 9 7 1,

1.10 3 5 , .0779 ,045 , 2.5786 , .56 25 ,048 , 3 . 3 6 7 9 , . 3 5 ° 6 ,''47/

DATA P C T ( I ) , I = 9 8 1,

1 .0859 , . 9 2 M , 047 , 2 . M 1 9 6 ,.385" , 046 , 3 . M 397 , s 7 I 8 » 045 , M . 05Mn , .04 99 ,134/

DATA ( PCT {

I

) , I = 994 ,

869 ) /24 2004 ,

879 ) /244a03 ,

R86 ) /247"02 ,

8"3 ) /249002 ,

903 ) /25 1 nc3

,

9 13) /254003 ,

970 ) /255'~C2,

977>/259r0 2,

9 3 7)/261003,

950 ) /263''04 ,

963 ) /265004 ,

970)/2 6Rno2,

9P0)/27ir0 3,

993 ) /277n04

,

I 0 0 6 ) / 7 7 3 0 4 , 45

LISTING OF 5 U PRO U TINE RAG

1 . ? A6n , . 2"^0 1 ,

1 35 , 2. ?9M5 . 39(S 3 , 134,

3 . 1 M3 1 ,,.1502, 04 5 , M . 2n6M , . 2033 , 136/

DATA ( PCT ( I ) , I = 1 on7

,

1022) /275005 ,

I . 2Q I 1 , . 2A39

,

138,

2 . 079C

,

. 06^1 , 1 35 , 3. 2MM I . 232 1 , I 36 , M . 38 1 n

,

. M2 ! 7 , 1 37 , 5.01 MP

,

.0131, 134/

DATA ( PCT ( I ) , I = 1 C23

,

1 n?o ) / 277^02 , 1.7767, . 7S96 , 138,

2.2231

,

. 2404 , 1-^9/

DATA ( PCT { T ) ,1=1030, 10391/27900 3,

1 . 6q3?

,

.6411, 1 39 , 2. 307 I , .3100, 133,

3. 0M96

,

.0489, 138/

DATA PCT ( I ) ,1=1040, 1046) /282002 ,

1 . , 133,

2.1358 1.0945, 138/

DATA PCT ( I ) ,1=1047, 1053) /284002 ,

1 . A23A

,

. 648 ! , 132,

2. 376M

,

.3519, 133/

DATA ( PCT( I ) , I = ! 0 5 4 , 1060) /286002

,

1 . 8977

,

,.9165, 134,

2.1023,,.0835, 117/

DATA ( PCT ( I ) ,1=1061, 1073)/287r04,

1.1034,, . 0994 , 04 6,

2 . 3 3 1 P , . 3450 , 134,

3 . 2 0 5 Q . 1877, 117,

M .3606 ,.3678, 13 1/

DATA PCT ( T ) ,1=1074, 1 ':'8o ) /29o^a2

,

1 .572^ , . 5452 , 117,

2 . 2 8 G , . 4548 , 04 6/

DATA [ PCT ( I ) , I = 1 O'^ 1 , 10871/29100 2,

! . 3 1 7M , . 2 8P8 , 117,

2.6826 ,.7102, 046/

DATA ( PCT ( I ) ,1=1088, 1 1GO)/293004,

1 . 2702

,

, . 2?«>5 , 116,

2.6 157 ,.6711, 1 '9,

3. OM 1 8

,

, . 0383 , 118,

M .072? ,.0610, 1 1 7/

DATA PCT ( I ) ,1=110 1, 1107) /294002 , 1.7598 , . 7525 , 117,

2.2402 , . 2475 , 118/

DATA PCT ( I ) ,1=1 108, 11 20) /297"C4 ,

1 . 1 755 ,.1644, 0<'7 , 2. 1 385 ,.1246, "9 1 , 3 . H 3 5 8 , . 45-^6 , 092 ,

*4 .2502 , . 2603 , OP 7/

DATA [ P C T ( -I ) , T = 1 1 2 1 , 1130 ) /298n03 , 1.8238 , . 8^07

,

09 1,

2.0636 , .077? , 087 , 3.1127 » . 1 1 9 1 , 0^7 /

DATA [ PCT ( I ) ,1=1131. 1137) /299r02

,

1 .6732 , . 6 A94 , 09 1 , 2.3268 , » 3 306 , HQ?/

DATA PCT { I ) ,1=1138, 1 1 47 ) /30 1 003 , 46

. y c A 7 ,091,

2 • 3 R 7 3 . M 2 1 , ^ 0 , 3 • 0 Q 9 6 ) , 0 0 A 2 ,0^9/

n A T A P C T f I ) . I = 1

. ^ A 9 D

, J 0 q ? . 3 M 7 7 . , 3 3 n (

!

• 1 ri n R

. n Q A M 0 Q <? /

DATA P f T ( I J 1 — 1 1,8070. . 892 1 » 1 " 5 » 2,1 12?) . 13 7 9 , 1 /

DATA ( P C T ( I

T - 1 >

*

. M 7 3 , n 8 1 » 2 • H B 6 5 , • M 6 7 0 > 1 " 5 » 3, 097"^ . 1 " 9 A ,11-13/

DATA ( ° C T f I I » I = I

I • 1 6 H 3 , .15 9 3 ,

1 7 , 2,2276, , 1 9 o 2 , 1 A , 3 . 7 M 1 . (11'-*,

*j , 0 P 4 n , Q p 3 7

-DATA ( " C T ( I

1 • 1 6 , .13 3 1 0 8 0 ,

^ . n A 9 9 . . r A Q c; • 0 " J . I 'J ^ 1

• 6 5 1 0 , . A A " 7 1 -1 A /

DATA ( 9 C T { I J I ~ 1

1 . A 7 u .5778,. p R n . t! 0 f 2,3776, , 3 a I 9 (079,

3.10^1 . r a ' 1 w - 1 / DATA ( p C T ( I . T = I 1 1 — 1

. ^4 6 7 7

7 - ? <5 A P .• w D J 7' ill-},

J • J _3 3 O 1 - I T t: . 1 -3 J J . P 7 9

q , pe^ 1 1 , , Q q 5 1 (080/

DATA ( P C T { I . T = 1

1 . 7 1 7 , 3 ^ ^ 1 ^ A 2,6788, • A M s 7 I 7 / DATA 1 PCT f I 1 .1 = 1

1 .705^ , . 7 2 ° 6 112,

2 . 1 3 9 A . ! 1 3 .

•\ . icq/. 114/ II'/

DATA 1 p r T f I ,1 — 1

1 . 3 =; 1 , -J n n, Q 1 1 1 7 . 7 A 7 A . 77 7 7 t 1 7 3.1 1 A 1 . , n 9 "J M 111", 1 A . M . 7 A A 1 . . 3 7 2 >^ 1 1 « / DATA ( P C T ( I J I — 1

1 . 9 M 6 7 , , 9 c; 6 q ! ! 6 , 2,0533, . n ti 3 A , 1 1 M /

DATA i P C T ( I »

'

. A 7 3 1 1 1 . f 2 • 3 3 G 7 , . 3 n f. M

3 . n 7 P 7 .0743 11'''/ 1 1 A / DATA 1 P r T f T ,1-1

1 , 7 L| Q 9 J .76^5 (092,

n 0 7 3.0735 , . "A73 .093 , «».n52 ! , . 0 =5 n , 09 A /

l'^S,IlA0)/3n7nOM,

1 6 1 , 1 1 A7 ) /30'4n02 ,

16 8,11771/305^0 3,

178,1100) /308004

,

1 = I , I 2 0 3 ) / 3 0 o r c ,

204, 1213)/3l 1^03,

2M, 1276)/3 l2rnq,

?27,1733)/31M00 2,

7 3'i,12'^3)/3l5003,

2'4M,i256)/3loroH,

757, l2A3)/32ino7,

76M , 1273) / 32?'"03 ,

77M, 12O6)/ 3 25'^0^,

47

LISTIfiG nr SUBfnUTlME BAG

DATA ( PC T ( I 1 I = 1 287 9 1293) / 326002

,

1 .9417 • 9 M 7 3

r) p 6

» 0 7 7 DATA ( P C T ( I -1 = 1 9 1 — I 7 9 4. 1 300) /377007.

I • 6 I 3 3 , . 5 P ? ,097 9 ? . 3 IR A 7 , • 1 7 6 , 096 /

DATA ( F' C T ( I 1r .9 !I = 1A w) ^ 1 9 1 3 1 C ) /33oro3 9

\ * -J \ c. L - ! 7 7 . 0 9 O 1 U 7 9

f. « J o o o . . 3 0 , 0 9 fl 9 3 • 1 M 9 n , . ! R 7 B (10 0/

DATA f p r T f T - T = 1 9 J 1 j 1 1 91 3701 /3370n3.

1 •"Job 4 7 4 n O P P

"> C n R 7 C, 1 D 7 1 n 1 I J J 9

"k. n C 7 Q n r: 7 D 1 n o / DATA ( P C T ( 1 . T = 1 3 7 1. 1 377) /33A007.

I • ! 72

1

, . t AS4 , 1079

2,8779 , 8 3 q A 10 3/

DATA ( P C T ( I 1 = 1 3289 1337)/337O03,

t • *4 4 Q . U £4 p A 1

1 r»1 / -J3 9

c: • U D O L? . n q Q 4 . 1 7 J .•

*4 O 7 ! 1 / ^ 1 10 7/

DATA ( p c T ( I : 9 I = 1 3 3 8 9 13'<M)/33R"029

1 . R A 1 9 , • R 7 7 3 1 1 ? 2 9

*- • loot .17 7 7. 1 4 fl /

DATA ( P C T ( I - T = 1 9 i 1 3 4 "^i .

1 5C4) / 34non3.

1 . <9 7 A . c, Q n /; , 1 0 R 7 . 7 0 c; c. 3 14 0 10 7,

7.1 1 A R . n Q Q 4 1Iff r- 7 /

DATA - T = 1 1 3A1 ) /3470n7.

1 . 19^9 . 1 4 1; c; . 1 n 7 - 7 . R n 1 . Q c q c. f • n n " rj 1 1 0 R /

DATA ( P C T ( T -911T = 1 J7 A0 z7 .

9

1 "^71 1 / 34c, on 3.

\ . R 1 1 • O J 1 D . P c 7 p 1 1 1

I 1 9 7 . n M o » r 4 R R 110. I 1 w 9 (10 9/

DATA ( P C T f T - T = 1 7 7 . 1 3R1 1 /347nn3.

1 . ^ 7 7 1 - ^ 4 7 C, 1It

1 A" 9

, . 4 ! 7 . 1 1 1 - i 1 1 9 3 . n 4 A c, 110/

DATA ( P C T ( I -1 = 1 9 I I ^ R'> f-7 -

9

1 *J7^0 3RR)f 1/349n07. -JI7'1.'Z 9

1 . 7 7 9 - 7 1 P A .110 1 1 9 7 . 7 7 4 i 9 • ' " 1

' 1 P P / ( I 7 / n A T A f p r T f I 91T —— 1

1 7J n p7 9

1 . R 7 I 7 • P 3 7 A 1 9

7 . 1 A R 7 . I A 7 4 9 • 1 O / ^ 110/

DATA ( P C T ( I 9 T = 1 3 9 6 9 140?)/354O07.

1 .7997 . 7 A A 3 9 • ' O O J ( 111

I 97

1 7 . 7nn 3 9 .• 7 JT J 7/ 1 1 A / DATA f p r T f I \ V 111 T = 1 9 1 ~ I 4" .J130 9 l41RW3t;cr>nEi.

9 • '/ 7 1 \. 1 7 0

7 . 1 > 7 1 n p*A 9 • 1

'

. ^ o 17 7

- 3 . 17 0 0 -.1709 .131 0 1-1 9 - 1 iT A . 1

1 1I R" 9 -

c: P P 9 III'/ 119/

n A T A r p r T r t T — 1 » 1 — 1 uHIV, 1 P i4*>i l"Jl \)/Jt)o' / *c:Anri4 Ui9

1 * 7 1 n ^ A D 7 ( 11 7r 0U 9 2. 1 69R 9.1731 » 1 2 1 3 . 020M 9 .01 73 .127,

48

LISTING or SUP OUT IMF RAG

H. 09 9(, ,.1272, 118/

DMA ( PCT ( I ) ,1=1432, !43R)/358 0 02,

I . ,.^311, 1 20 , 2 . SS 1 2 , . 56 " 9 , 12 1/

DATA { PCT ( I ) , I = 1 4 3P , 1 445 ) / 360r>02,

1 . 6 1 6P , . 6C'^5 , 120,

2. 3fl3? ,.3955, 119/

DATA { ^CT ( I ) , I = 1 M M 6 , !4 55)/363n0 3,

I . 272a ,.2645, 13 1,

2. 2S7A , . 2«65 , 1 27 , 3 . M697 , .

'^ f'^O , 130/ - DATA ( PCT ( I) ,1=1456, 1465) /364003

,

1.7369 ,.763 1

,

13C,

2. 1 3^"^ , . 1 1 2 1 «

1 40 , 3.1286 ,.1248, 14 1/

DATA ( PCT ( I ) ,1=1466, I48q)/368n06,

1.0231 , . n?67 , 124,

2 . 086S ,.0687, 14 1,

3. M206 ,.5076, 1 28 , M. 3 1 87 ,. 286 1 ,

1 29,

5. OMMS ,. 0338 , 1 3D , 6 . 1066 , . 0772 , 127/

DATA ( PCT ( I ) ,1=1485, 1497) / 369004 , 1.1231 ,.1320, 143,

2.0888 , . 06^45 , 1 29 , 3 . M8M 3 , . 5065 , 142,

*4 . 3037 , . 2970 , 14 1/

DATA ( PCT ( I ) , T = 1 4 98 , l5nq)/3 7!'^0 2.

1.7959 ,.7925, 143,

2.2Q4 I , . 2075 , 14 4/

DATA ( PCT f I ) » T = 1 5^5 ,

1 5 1 1 )/ 3 7 5 0 0 2

,

1 . 828 A ,.8199, 14?,

2. 1 7 1 H ,.1801, 143/

DATA { PCT ( 1 ) ,1=15 12, 1521 )/37PCD3,

1.5135 ,.4539, 1 , 2 . 3 1 H 8 , . 3383 , 129,

3.1717 , . 2n78 , 126/

DATA ( t^CT ( 1 ) ,1=1522, 1578) /383CC2

,

1.2193 ,.1733, 12 1,

2.7807 , . 8 ?67 ,

1 27/

DATA {°CT ^ I ) ,1=1529, 1 55C ) /3S4007 , 1.115 1 , . C'^57 , 128,

2. 1 Ml M , . 1 30M ,

I 25 , 3.5373 ,.6059, I '^5 , M.0383 , .0235 , 1 ^40 ,

- 5.1163 ,.1087, 1 27 , 6 . 0?M I , .0220 , ! 26 , 7.0275 , .0737, 12 1/

DATA ( PCT ( I ) ,1=1551, 1563)/386n04,

1 . '457H ,.4177, 121,

2.1153 ,.1195, 122,

3. 357n , . 3^24 , 123,

M .0703 , .0704 , 1 25/

DATA ( PC T ( I ) , T = 1 564 , 157 0/ 39010?,

1.7555 , . 7 26n , 123,

2 . 2MM5 , . 2732

,

122/

DATA ( P C T ( I ) ,1=1^71, 1 5P9 ) /39?"06

1.1269 ,.1073, I 3n ,

LlSTIMfi or SUBPQUTI'lf: ^ A G

tl U 0 1 , • M M H I » 1 / I 1 3.0207 , . 0 1 P 2 , I ? 6 , M . 2SB 7 , .2^79 > 1 7 5 , 5 • 0 B » . 0 P 5 P t

1 7 3 , 6 • 0 7 n 1 . 0927 1 1 7 7 /

DATA ( r C T ( I ) , I =^ I 5 9 0 , 15''9)/ 3 96r^03,

1 » 5 3 S '^ ,.5510 , 1 7 M , Z • 3 7 U 7

^ / /I o

, • 3 6 ^4 ^ »

1 ^ 7 , 3 • 0 9 M M , . n 8 M 1 ,17 6/

DATA PCI { I 1,1 = 1 6 0 0, 16n6)/398~02,

1 • 2 1 9 7 , .2306 , 1 M 5 , 2.7803 , . 7694 , 1 7 '4 /

DATA ( PC T ( I ) , T = ! 607 , 1613)/'4D10C7,

I . M967 , . 5 r n 1 ,16 3,

2.503 3 , , y Q 9 9 ,16 7/

DATA [ P C T ( I 1,1 = 1 6 1 M , 1620) / '4C6 0 02,

1.7959 , . P 6 r 6 ,1^5,

2 . 20 H 1 , • 1 3 9 M ,1^6/

0 A T A " C T ( J ) , T = 1 6 2 1, l627)/'4n7nQ2,

1 . 6 7 2 M , « 6 7 " P ,1^6,

2.3776 , .3202 , 1 ^4 9 /

DATA P C T ( I > I = 1 6 2 8 , 163'4)/4r. 8002,

1 . 827 6 , . P I 7 7 » 1 6 , - 2 . 1 7 2 ,1 . 1 P 7 3 ,> 1 7 7 /

DATA 9 C T ( I ,1 = 1 6 3 5 , !6'44)/Mn9nQ3^

1 . 8598 , . 8 .3 3 9 1 1 7 2 , 2 . 0 8 9 M , .073'^ 1

^4 5 , 3 . 050 » . 0 3 1 , 1 " 6 /

DATA ° C T ( 1 ,1 = 1 6 " 5 , 1651 )/'4l30?2,

1 . 6793 « 6 7 3 P 15 3,

2.3 207 .3762 , 0^ ^ /

DATA P C T ( I .1 = 1 657 , !65P)/HlM'-!0 2,

I .7705 , . 7 6 P ,15 3,

2.2295 , . 2352 . n o ^4 / DATA P C T ( I » I = ! 6 5 9, 1665)/416G02,

1 •9^6'^ . 953 2 , C ^ M , 2.0536, . 0 M 6 8 5 /

DATA P C T ( I »

I = 1 666 , 167P)/417'*'04,

1 . M 66 M ,'4519,,093 , 2 . 2 9 M 0 . 3 C 7 0 ,095 , 3.21 19, . 2 1 « 7 ,,094 , M . 02 77

,

,.0779, 0 9 6/

DATA P C T ( I 1 = 1 6 7 9, !685)/'42C00?,

I . 1 2 I C ,

, . 1 7 f> 7 , n 9 4 , 2.8790, . 8 7 1 P , 095 /

DATA P C T ( I .1 = 1 6 « 6 , 1692)/'47IOC2,

1.6337, .581 !

, 1 5 M , 2.3663, , M 1 q 9 , 15 7/

DATA ( P C T ( I 1 T = ! 6 3 , !7P2)/^27CC3,

1 . 1 3 Q ? , . 1 M M 9 , 159,

2.6362, .6^32, 15 2,

3.2255, .2170, t 5 M /

DATA ( PCT ( I ) .1 = 1 703 , 17 15) /423C'0'4 ,

1 . 60*^ 9 , .6192, 15 2,

2 . 0 M 7 P , . 0378 , 1^1,

^ r" 3.2705 1 • 2 6 3 , 15 4,

M . 076P

,

. 0566

,

159/

DATA ( PCT ( I ) .1 = 1 7 16, 1 77P ) /42'4P0M ,

LISTING OF SURROMTINE RAG

1 . 7 A 1 I * O /L O I % • J /I _ n J 15 2,

^ • 6 1 V 6 » • 6 3 1 T » ICO 1 S V , 3*0339 , . 0 ? 9 R , 1 5 M , M • 0 2 0 M ,.0179, 15 1/

DATA ( P C T ( I ) ,1 = 1 7 2 9 , 1 735 ) /425O02 t

1 • 96 I 3 , .9636, 15 9,

2*0307 , . 0 3 6 ^ t 1 5 M /

DATA ( P C T ( I ) , T = 1 7 3 6 , !7M2)/M26n02,

1.9035 , .9130, 0 9 5 , 2 . 09 A , . 0970, 10 0/

DATA ( P C T ( I ) ,1 = 1 7M 3 , I 7M9 ) / M270 02,

1 • 86H ? , . R 5 5 M , 150,

2 • 1 3 5 8 , I M 6 ,

! n u /

r\ A T A U A 1 Ar D ^ T / T \

I r C T ( 1 ) 75C , I7b6)/M28r'u2f

1 • 9 ^4 5 ft , • 9 ^4 7 M , 1 5 C , 2 • 0 B H M ,.0576, 1 1 /

DATA ( P C T ( I ) ,1 = 1 757 , 1766)/M29 0 03,

1 • 3 M 7 P , • 3 5 M 3 , 15 0,

2 • S 5 1 ? , . 5 6 n 1 , 15 1,

3 • 1009 , . 0 P 5 6 , 15 2/

DATA (

- C T ( I ) » I = I 7 6 7 , 1773)/M35'^n2,

1 • M 7 20 , .^692 , 16 1,

2» 57B0 , . 5 3 P , 15 1/

DATA ( PC T ( T ) ,1 = 1 774 , 17P3)/M36003,

1 • 2229 , . 2223 , 15 1,

2 2 2 6 7 , . 1 R 1 3 , 16 1,

3 • B50 , .5064, I M 9 /

DATA ( PC T { I ) ,1 = 1 7PM, 1796)/43prOM,

1 .0959 , .0600 , 1 ^4 7 , 2 • 7 3 9 M , . 7 6 9 ft

, IMP,

3 • 1 5 7 o , . 1 6 4 S ,

t /I n

1 'l V ,

£1 r> ."^ ^ 1 H • U J 7 1 , . U C ^ 6 , 15 0/

DATA { P C T ( I ) ,1 = 1 797 , IP03)/4Mir!02»

1 # 7372 , . 7 1 P ,1 ^4 7 , 2« 2A28 , . 2 3 P 2 , IMP/

DATA ( PC T ( I ) ,1 = 1 ROM, 1P10)/MMM0C2,

1.963? , . 9^,26 ,1 M 7 , 2.0365 , . 0 4 7 M , 1^9/

DATA ( PC T ( I ) ,1 = 1 fill, 1P17)/4M600 2,

1 . 2 M 0 6 ,. 2351 , 1 5 M , 2 . 759M , . 7 6 '4 9 , 155/

DATA' ( PC T ( I ) ,1 = 1 P 1 8 , 1824) /MM8002 ,

1 . 5 3 ft P , .5072, 1 5 4 , 2 • M 6 1 2 , . M P 2 8 , 15 6/

DATA ( P C T ( I ) 1 I = 1 P 25 , 1 83M ) /449r03

,

1 . 6 6 M , .6609, 15 6,

Z . 3 0 0 R , . 3 (3 3 P ,

1 5 M , 3 • 0 3 H 3 , . 0 3 5 , 15 7/

Data ( C T ( I ) , T = 1 P 3 5 , 18 4 1 ) / 4 s I 0 0 2 ,

1 .8391 , • P I M P , 15 5,

2 . 16 0 9 , . f P 5 2 , 15 7/

data ( P C T ( 1 ) , T = 1 P M 2 , IP48)/M56O02,

1 . 9 5 6 M , .9559 , 15 7,

2 • 0 H 3 6 , • 0 M M 1 , I 6 /

data ( PC T ( I ) ,1 = 1 P M 9 ,

1 P55 ) /463002

,

1 .5363 ,.539?, 170,

2.^637 , . M608 , 169/

DATA ( PCT ( I ) .1 = 1 P56 , 1 862 ) /4640n2 , 51

1 . 6039 , . 6 P o 4 , 1 6" , 2.396 1

,

. 3 1 , 1 6° /

DATA ( PCT ( I ) ,

I = 1

1 . 399P , . 4 >=. '4 7 , 1 A P ,

. M c M 9 , ! 7 1 , 3 . 052S

,

.0604,

I

6 9/

DATA ( PTT ( I ) , 1 = 1

1 . 376S , .3^46, 1 6 7 , 2. 3?2« , .26-^7, 1 A P , 3.2727, . 3C56 , 1 A A , '4.02a 1 , .0331,1 65/

DATA ( ^CT ( 1 ) , I = 1

1 . 654 1 , .7100,1 67 , 2. 3H5'' , .2900, 1 7 1 /

DATA ( PCT { T ) , I = !

1 .0933 , . oa Ro , 1 65 , 2 . 9G67 , . 9 1 so , 1 66/

DATA ( °C T ( I ) , ^ = 1 1.3499, . 4079 , 1 A 6 , 2.650 1 , .5971,1 65/

DATA ( T ( I ) , I = I 1.9997, .

'J « 7 4 , 1 60 , 2 . n J 03

,

.0126,1 6 1 /

ccccccc cr c S^FX I ^ r M I c H MA

R63,lP,7?)/465'^03,

P73, 1P«5)/47P"04,

PP6 , 1 P02 ) /47 1 '^0^ ,

R93 , 1 89=5 ) /475002 ,

O00,19r.6)/477'^C2,

of^7,19!3)/49n'''02,

ccccccccccccccrcccccrcccccrcccrcccr

CCCCC S^FXl^y '/MICH MATRIX IS TO BE CO"°WTEn CCCr

CCCCCCCCcCCCrCcCCCCCCCCCcrcCfCCCrCC^CCCCrCCCcCCCr

GO TO f 1 0 , 2 r , 3 5r ,160-^),

CCCCCCC

T fi T I ON P cc r CCCCCCC

(

1

5 ) + T R ( 22

op ( VOL )

CCCCC 00 C!

100 VOL=

CALL

F T U » ^J CCCCCCC CCCCCrccCCCCCC

HQ CO'--^UT A T I OV' FOR

r

HIGH VALUE - LA^GE,

cc^ccrcrcccccrcccccccrccc

cc^cccccccccccccccccccccc

ccccc'-rcccccccccccccccccc

! 0 no, i)nn,i7nn,i3Q,;,

CCCCCrrrCCCCCCCrccCC CCCCC

'mE ccccccccc^^ccccccccccc

CCCCCfCCCCCCCi

CrccCCCCCr

WTG^ yALijr - LAf'GE, COST

ccccccccccccrrcrcccccrcc-.:

crcccccccccrrcccccccccc.

cccccccrcccrcccccccccccc.

CTM= ,

n

1 7 * ( - I 2 ! 6 • + 9 24 . I 9 ( 3 ) . 4 2 1 ) R T H = 1 3 4 9 . 525 1+3 . 5 1 p 7 8 6n* ( 3 ) C T H = M I

f

1 (

' A X ( C TH » T H ) J

1 . 1 7 ^ T H ) C9H=7

1

n • + -> . 9A« I ( 5 ) C V = 296 • 5* . A a q • I p f 9 ) DOL= I P ( 1 5 ) «CTM + I D ( 7 ? ) *CP I R { 3 " ) « C

"' CALL A 0 D ( nOL ) RETURN'

ccccccccccccccccccccccccccrcrcccrcccccccrcccfcccr

CCCCC 00 COm'^UTATIO'J for high VALUE - SmALL, VOL"

ccccccccccccccccccccccccccccccccccccrccccccrcccr

300 VOL= I ( 1 6 ) + I R ( 23 ) + I R ( 2^ ) CALLAnn(VOL)

RETURN

cccccccccccccccccccccrccccccrcccccccccrr^cccccccc

CCCCC 00 COmPMTATIO'! HIGH VALUE - SmALL, C^S^

ccccccccccccrccccccccccccccccccccccrccccccccccccr

400 CTH=,nP17«(-i:'!6.+9''q.«IR{^).«.q7| ) RTH=13'J9.525l + 3.5187«6n«lP(3)

CTH = 'MM(maX(CTH,oth),1.!7*RTH)

CRH = 7 !«.+2.9A«IP(5)

52

CCCCCrcri

mE CCCCrccCCCCCCCCCCCCCC

cc^cccccccccrccccccccccc

CCCi

ccccccccccccccccccccccc

ccrccccccccccccrcccccccc

LISTIMG OF SURROtJTlME RAG

CCCCC

500

CCCCC

CCC

600

C A= I CSMO , + la.''7»IR(7)

D0L=lf?(l^)*CTH+TR(23)*CRH+lR{?9)»CA

CALL Ann ( DOL ) RFTUPf)

no co^'^UT^TIorJ foR "ed value - large,

cccccccccccccc

VOl= I f'(17) + I^(2'<) + IR(3n

call AHDCVOL)

PETU^'fJ

DO CO'-if'UTATIO'J FOR ^FD VALUE - LARGE,

VOLUME rcCCCrcCCCCCCCCCCCCCCCC

rccrccccccccc^ccccccccccccc

rcccccrcrcccrrccrcccccrcccr

rOST cccccrcccccccccc ccccccccc

CTm=.'1''A1«(-1''1a.+92".*IR(3)*«

RTM = 99?,q3256-»-3. 15nn7B<5*IR(3)

CTM = Mn(MAX(CTM,PTM),i,|7»RTM)

CRM=3M1.+1.7«TP(5)

CW=296 . . AMM* I P ( 9 ) DOL=IR( 17)*CT"'+IR{24)*CRM + IR(3 1 )*C'"'

CALL ADDfDOL)

R E T U R rj ccccccccccrcccccccccccccccccrcccccrr

ccccccccc

M2 ! )

CCCCC DO FOR

700

CCCCC

CCCCC

CCCCC

COMMUTATION]

r CCCCC

VOL= I R ( 1 R ) + I R ( 2S

CALL ADD(VOL)

RETUPh'

r DO C 0 M ij T A T I 0 f! FOR

ccccccccccccccccccccccccccccr

800 CTM=.iJ9fel*(-1216.+92'J.*IR{3)*

RTM = 992.M32S6 + 3. lSBn7c,6*IR(3)

C T M = M I N ( K A X ( C T M , R T M ) , 1 . 1 7 * R T M ) CRM = 3M 1 . + 1 . 7* T ( B ) DOL=IR( in)*CTN' + IR(2?)*CRM

CALL ADD(DOL)

RE TURK'

cccccccccccccccccccccccccccccccccc

DO COr-pUTATIO'] FOR 10"" VALUE - LARGE,

cccccccccccrrcccccccccccccccccccccccccccr

900 VOL=,I P ( 1 ^ ) + I R f 26 ) + I R ( 32 ) call" A'^D(VnL)

R E T U R ^.|

CCCCC DO COMPUTATION FOP LOW VALUE -

1000 CTL = .M1P«^»(-1216.+92M.*TR{3)*«.*42) ) RTL = 76R.P5*<76 + 2.659n.-^CP*IR(3)

CTL = ^MN(MAX{CTL,PTL),!.17*f?TL)

CRL=2I .H+l .H7*IR(5)

CW = 296.5+iA4M»IP(9).

D0L=1P(I9)»CTL + IR(26)*CPL+IR(3 2)*C"^

CALL ADD ( DOL

R E T U R M ccccccccccccr

CCCCC DO COMMUTATION FOR LO".' VALUE - SMALL,

1100 VOL= I P ( 20 ) + I P ( 27

53

CCCCCCCCCC

cccccccGccccccccccccrcccccccccccccc

MFO value - SMALL, rOST C C C C C C C C C C C C

CCcCCCrCCCCCr

M 2 n

rccrcccccrcccccccccrccci.

vOLU^'F rcccGcccccrrcc

''CCCCCCCCCi

E, COST

CCCCC'

:ccccc

ccccccccccccccc.

ccccccccccc

CCCCC

VOLUME

c

LISTING OF SUrROUTiNe RAG

CALL AHD ( VOL

)

ccccc

CTL=.M18P*(-I7IA.+92'^.*IR(3)**

RTL = 76a.85M76 + ?.6.59n3rip*IR(3)

CTL=MIM(MAX(CTL,RTL),I.17*RTL)

CRL = 2 I .M+J.q7*IR(^)

D0L=IP(2n)»CTL+IR(27)*CRL

CALL ADDfPOL)

RETURN

. , .

: r DO COMf'UTATIOM FOR CRUHE - LARGE, V 0 L U ^' E

1300 V0L=IR(3 3)+IR{3S)

CALL AnD(VOL)

RETURN

CCCCCCCCCCCCCCC

CCCCC, DO COMPUTATION FOR

CCCCCCCCCCCCCCCCCCCC

IMOn CPC=7S.PM

CW=29f.5+

DOL= I R ( 33

CALL AOD(

RETURN

CCCCC DO COmPUT

CCCCCCCCCCCCCCC

1 BOG VOL= I R ( 2 ! CALL A n D

{

RETURN

CCCCC DO CO PUT

cccccccccccccc

160n CTL=.M!RR

R T L = 7 A R , P, C T L = M I N ( M CRL = ? I . +

Ci'i/ = 2''6.S + CPR= . q 1 q R DOL= I P ( 2

1

CALL ADD

(

RETURN

5URR0UT I

N

IMPLICIT

CCCCC

CCCCC

CCCCC

ASSIGN TO

WHICH COR

THT 7PMFS

C R U n E

C +.237M*IR( 1 1

)

. * \ R { 9

)

) * C 'V + I R ( 3 5 ) * C P C DOL )

ATION FOR f^EFINED - LARGE, VOLUmE CCCCCCCi

) + I R ( 2 R ) + I R ( 3 4 ) + I R ( 3 6 ) VOL )

ATION FOR REFINED - LARGE, cr\<,l CCCCCCCrrr

cccccccccccccccccccccccccccccccccrccccccCfi

*{-1216.+97M.*IR(3)»*.M2I

)

S M 7 6 + 2 . A S 9 9 3 0 R * I R ( 3 ) A X ( CTL , RTL ) , 1 . I 7*RTL

)

1 . q7* I R ( B )

. AMM I R ( 9 ) * I R ( 1 3 ) )«CTL + IR(2fl)»CRL+lR(3M)«C'^'+IR(.3A)»CPR

DOL )

E ADn { X )

I

^.' T E G F P ( A )

AO A^I^ AD, RESPECTIVELY, THE VALUES IN THE NDX

RESPOND TO THE ORIGIN AND r>ESTl NATION NN5 70NFS

LIE BET VEEN ZONES 1 AND S \ INCLUSIVE) CCCCCCC

A0= I R (

M

IF (A0.GT.S12) GO TO 6000

IF (AO.LT.M GO TO 6000

AO=NDX(AO)

AD= I R ( 2

)

IF {AD,GT.S12) GO TO 6I0J

IF (AD.LT.l) GO TO A!1C

54

LISTING or SUBROUTINE RAG

AD = fJnX ( AD ) IF ( AO ) ?nno ,/«,?^n , i 7on

1700 IF (AD) I OQO , A ?nn , inoo

CCCCCCCCCCCCr

CCCCC IN THF FOLLOWING COMMENT STATEMENT

CCCCC 0=0R1GIM NN5 70NE; n=DESTl NATION

T

CCCCC AGGREGATE FOR CASE OF ROTH 0 K P E

18P0 PMTX(AO,An)=RMTy(A0,AD)+X

GO TO 9999

CCCCC AGGREGATE FQR CASE Op n IN SINGLE

S THE rjOTATlON "'ILL ^E! CCCCCCCCCCC

NNS 70 Nf^ CCCCCrcCCCCCCC CCCCCCCCCCC

C ' r I

r-J T I P F L Y V.' I T H I M SI rg GLE RE A ZO rjES CCC

REA 70kjE ^[^r^ D Im MULTIPLE ZONES CC

r

1 900 A 3 = I C T ( - A D ) /

IF (A3.NE.IR(

AM = NOO{ IPCT(-

DO 1 9Sn K= I ,

A

ADK= I PCT ( 3»KPCTDK=PCT ( 3*K

RMTX ( AO , ADK )

GO TO ^^999

IF (AD) 2200

,

AGGREGATr FOR

A1=IPCT(-A0)/

IF ( A I , N E . I R

(

A2 = !'10P( IPCT(-

DO 2 1 5C J= I ,

A

AO J= I rCT ( 1» JPCT0J=PCT(3*J

BMTX(A0J,A0)=

60 TO 9999

CCC

AGGREGATE POR

100 0

2)) GO TO ^HOO

A n ) , 1 C n ) 4AD ) -AD-2 ) RMJX ( AO , ADK ) +PCTDK'*X

A 2 0 0 , 2 1 ? n

r CASE OF

^ ^ ^ r

1 000

n ) GO TO

AO) , 1000)

2

A 0 ) - AO- t ) M T X ( A 0 J , crcccccrc

CASE OF

0 IN MULTIPLE ZONES

C C C C C C r

6300

A rj D 0 I N

' S I N G L E ^ E A ZONE C C rccrcccrcccccrcccccccccc

AD ) +PCTO J»X

rcccrccccccc

nOTH 0 f- D Ir

c c c r c c c c r c c c c c r i N U 1. T I I E 7 O r s cccccccccccci

CCCCCCCCCCC

A1=IPCT(-A0)/1 nnn

A3=IPCT(-An)/^^^^

IF (Al.NE.IR(l)) 60 TO 63 00

IF (A3.NE.1R(2)) GO TO 6M00

A2 = M0n ( I PCT ( -AO ) , 1 COG

)

AM = MOD( IPCT(-AD) , lOOD)

DO 22Sn j= I , A2

A0J=IPCT(3*J-An)

PCTOJ=^CT ( 3* J-AO- I

)

DO 22SO K=I,AM

ADK=IPCT(3*K-An)

PCTDK=PCT ( 3*K-An-2

)

RMTX(AOj,AnK)=PMTX(A0J,AnK)+PCT0J*PCT0<*X

60 TO «7999

CCCcCCCrCcCCCCC

THE RPMAjr^nER OF THIS

cccccrcccccccccccccccccc

PROGRAM C0^JSI5TS OF FRHOR M^^SSAGPS CCCCCCCCCCCC

6000 W'R I TE ( A , 600 I ) 1 ( 1 ) , T ( 2 )

' 6001 FORMAT(» c-RR0P,RAO irjPllT RECORD.

-MO) 55

U N K N 0 N ORIGIN. T P ( 1 ) , I P ( 2 ) I

• 2

LISTING nr SlJpROUTTME RAG

GO TO 999^

6|n0 I T=- ( 6 , 6 1 n 1 ) IR(1),!P(?)

610 1 FORMAT (• EPPOR.RAD INPUT R.ECOf'O, UNKfjni"M D F 5 T I N A T I 0 N , IP(1),IR(

2 )

: ' 7 I 1 0 ) GO TO 9QQ9

6200 '.'.'RITEI^.AZ^l) I9(l),iP(2),A0,AD

6201 F0PMAT(' ERROf'jRAn N'DX MAPPpJr. ,IP(1),IP(2),A0,AD:MI10)

60 TO P99?

6300 ''.'PITE{6,63ni) 19(1), 1^(2), AO, ^C, A 1, A3

630 1 FORMAT{«EPPOP,BAr> PCT ••^APPIflG (Al-T^(]n. IR(!),IP(2),AO,An,Al,A3

+•6110)

GO TO ^999

6HC0 /RITE(6,6'^01) I^(l),!9f2),A0,An,Al,A3

6M01 F0RMAT(«EPP0R,PAn PCT MAP!=lrJG (A3-IR(?)). IP(l),IRf2),AO,An,Al,A3

+ • 6 n n ) GO TO 9979

9999RETUP^l

ENO

END ON SITE P'='P'TOUT ON OCTOPFP 6, 197S A"^ 13:'^9:'^7

3EA*SUB(l).f^AG(!4)

I

I

56

i)

v^A^HlC DATA

SH c

1. rriii.ir ATiON ok Ki-.p()Rr \o

:;3si-i 7 5-911

3. : ; e( 1 }Mi-ni <i Accc^-.

.il LL; and Si:aTITLH

A Ma::icnal ^Tetvork Dana Base System

5. I'uSl i.:.iiu'ii D.itc Seotember* 107 ^5

6. I'l-rlnf min>; Or);.ini ion f . J'-

7. AUTHOR ^ Richard K.F. Jackson2. rf rriin^ C-r ;an. ?.irp'^':: -

9. PERFORMiNu ORG AN iZAl" iON XAS'F .VNO ADI3H!:?5

'^ATiOnAu 3uSHAJ of STAinDaSDS

DEPARTMENT Or COMMERCE

VASHINGTCN. D.C. 2Q234

10. Pr J ici. : i a s'< . ork L'ni: No.

11. ' -oinraL t, Grant No.

1 2. Sponsoring Crsanlzdtion Nane anj Cc-plete AJdresi I'Slreet, City, S;ale. ZIP)

U. S. Departnent of Transportation

U. S. Departnent of Cornnerce

U. S. Department of the Anny

13, I'vpe of Rc-por: i Fer^oc

CovereiiFinal

- "

.

. - 14. ^ponst^rin^ Agency CcJ?-

15. SLPPI.E.MKNTARY NOTES

16. AM-^l.'<\C I . f.-l llyit- Mi^rtl •ir U'ii [.i-niii .nmrr.irv r,[ mnnl 1 1:^) 1 1 i r ,in I in I unn^i 1 1 nn . I ! rin' 'irT-n I inf.in'l'-s n s i tin 1 1 1 c.inl hi bl ittiir-iitliv «•/ lilfrnlitrri -.iirvv, ni'-iili,»\ i( Ifn- ) This report documents a data base "system" created at the National Bureau

of Standards. The National Netvork Data Base System (NND3S) provides information

on ilows of freight and passengers throughout the transportation system of the

United States. It consists of a set of FORTRAN -programs -\^it-ten--for -the-NBS --

U^^'IVAC 1103 (but transportable) and some basic data tapes. In addition to

providing basic data on the transportation network, the NNDBS can produce modal

splits, aggregations over certain "zones" in the U.S., and is capable of easy extension to other uses. This report is intended as a user's guide and includes

discussions of -the data tapes and each of the programs. Complete listings and

tape formats are also included.

17. KEY ^0RD5 (six to iv-e'ive rr.lries: a .'/ ca; order; capitalize only the lirst letter of the fiml key word imlz^ss a proper

ncjfTie; sepdrflfed fay sem;co/ons; Commodity Flows J Data Base; Freight Mode Choice; Market Split* Modal Choice; Modal Freight; National Transport Analysis; Regional Transport

Analysis; Transportation Network.

13. AVAiLABtLITY __X_ L'p.l imued

[ For Official Distribu:ion. Do No; Release to N'TIS

: Or.:er From Sup. of Dnc, U.S. Gnvernir.e-.t H-;n:int; Office

i.\^hingtpn, D.C. ZO-iOZ. SD Cat. No. CH

i Order From National Technical Information Service (NTIS)

Spr-npfield, Virginia 221'^ I

19. SF.CLRITY CLASS

(THiS REPORT;

UNCLASSIFIED

21. NO. OF PaGFS

60

20. SFCiURlTY CLASS

(THIS PAGE)

UNCLASSIFIED

22. Price
.0 7'


NBSIR 75-912

(This number has been canceled in

order to allow the records of

Office of Tech. Pubs, to be brought

up to date. A new number will be

assigned when the report has been

approved by WERE)