AED-A235 877

US Army Corps
of Engineers
Hydrologic Engineering Center

GENERALIZED COMPUTER PROGRAM

HEC-4

Monthly Streamfiow Simulation

User's Manual

0 I-1

February 1971

Approved for Public Release. Distribution Unlimited. CPD-4

 HEC-4

Monthly Streamflow Simulation

User's Manual

February 1971

Hydrologic Engineering Center

US Army Corps of Engineers
609 Second Street
Davis, CA 95616-4687

(916) 756-1104 CPD-4

Mr...
.. .~.......
. . ......
... .....
....
... . ... .
 HEC-4

MONTHLY STREAMFLOW SIMULATION

HYDROLOGIC ENGINEERING CENTER

COMPUTER PROGRAM 723-X6-L2340

CONTENTS

Paragraph Paqe

1 ORIGIN OF PROGRAM 1
2 PURPOSE OF PROGRAM 1
3 DESCRIPTION OF EQUIPMENT 1
4 METHODS OF COMPUTATION 1
5 INPUT 7
6 OUTPUT 8
7 OPERATING INSTRUCTIONS 8
8 DEFINITIONS OF TERMS 8
9 PROPCSED FUTURE DEVELOPMENT 8

EXHIBITS

1 DETAILED EXPLANATION OF COMPUTER PROGRAM

2 DESCRIPTION OF CROUT'S METHOD
3 INPUT EXAMPLE
4 OUTPUT EXAMPLE
5 DEFINITIONS
6 SOURCE PROGRAM
7 INPUT DATA
8 SUMMARY OF REQUIRED CARDS

1/85
 HEc -4

MONTHLY STR AMFLOW SIMULATION

HYDROLOGIC ENGINEERING CBTER

COMPUTER PROGRAM 723-X6-L340

1. ORIGIN OF PROGRAM

This program was prepared in The Hydrologic Engineering Center, Corps

of Engineers. Up-to-date information and copies of source statement cards
for various types of computers can be obtained from the Center upon request
by Government and cooperating agencies. Programs are furnished by the
3cv.rament and are ac-ept!d and uzcd 'U,the recipient upon the express unaer-
standing that the United States Government makes no warranties, express or
implied, concerning the accuracy, completeness, reliability, usability, or
suitability for any particular purpose of the information and data contained
in the provrams or furnished in connection therewith, and the United States
shall be under no liability whatsoever to any person by reason of any use
made thereof.

The programs belong to the Government. Therefore, the recipient further

agrees not to assert any proprietary rights therein or to represent the
programs to anyone as other than a Government program.

2. PURPOSE OF PROGRAM

: -ThiR program will analyze monthly streamflows at a number of inter-

related stations to determine their statistical characteristics and will
g7enerate a sequence of hypothetical streamflows of any desired length havinr.
those characteristics. It will reconstitute missing streamflows on the
basis of concurrent flows observed at other locations and will obtain maxi-
mum and minimum quantities for each month and for specified durations in
the recorded, reconstituted and generated flows. It will also use the
-eneralized simulation model for generating monthly streamflows at ungaged
locations based on regional studies. There are many options of using the
program for various related purposes, and it can be used for other variables
such as rainfall, evaporation, and water requirements, alone or in combination.

3. DESCRIPTION OF EQUIPMENT
This program requires a FORTRAN IV compiler, a random number generator
(function RUGEN included, see exhibit 2), and a fairly large memory (64K on
the CDC 6600). Provision is made for use of three scratch tapes, 7 (for
punched output), 8 and 9.

4. METHODS OF COMPUTATION

a. In the statistical analysis portion of this program, the flows

for each calendar month at each station are first incremented by 1 percent
of their calendar-month average in order to prevent infinite negative

'! J
logarithms. This increment is later subtracted. The mean, standard deviation
and skew coefficients for each station and calendar month are then computed.
This involves the following equations:

Xim = log (qim + qi ) (1)

N
Xi Xi,m/N (2)
m--l

S j (Xi-m )2 /(-)(3)
N

gi N I (xi~ x )3 / ((N-1)(N-2)S3) (4)

1m=l i'm i

in which:

X = Logarithm of incremented monthly flow

Q = Monthly recorded streamflow
q = Small increment of flow used to prevent infinite logarithmns
for months of zero flow
= Mean logarithm of incremented monthly flows
N = Total years of record
S = Unbiased estimate of population standard deviation
g = Unbiased estimate of population skew coefficient
i = Month number
m = Year number

b. For each station and month with incomplete record, a search is

made for longer records among the stations used, to find that which will
contribute most toward increasing the reliability of the statistics com-
puted from the incomplete record. The mean and standard deviation are
then adjusted. Equation 5 is used to compute the equivalent record required
to obtain statistics equally reliable to these adjusted statistics and is
the basis for selecting the best record to be used in the adjustment.
Equations 6 and 7 are the adjustment equations.

N N1
1 2
1N 22,N
-,-N
1
I R2

N2

2
 , = , 2 ) RS/S 2
(6)

S)IS (7)
S, -S = (S' -SI )

The primes indicate the long-period values and those without primes
are based on the sane short period for both stations 1 and 2, and:
N = Length of record
R = Linear correlation coefficient

c. Each individual flow is then converted to a normalized standard

variate, using the following approximation of the Pearson Type III distribution:

tim = (Xim - g4) / Si

Ki' m =6/g [ ((giti,m/2) + 1)1/3 - 1] + /6(9)

t = Pearson Type III standard deviate

K = Normal standard deviate

d. After transforming the flows for all months and stations to normal,
the gross (simple) correlation coefficients R between all pairs of stations
for each current and preceding calendar month are computed by use of the
following formula:

Rii- {l
Rr -- --[ [1 - ( N xXimil 2 Nxx2N
z)/Xi 2_, ,2]n
]
i'm=l i.mm=l i
m--l
(N-l)/(N-2)} (10)

in which:

X = X-X
 e. If there are insufficient simultaneous observations of any pair of
variables to compute a required correlation coefficient, that value must be
estimated. Each missing value is estimated by examining its relationship
to related pairs of values in the current and preceding month by use
of the following formula using i, ,J,and k subscripts to indicate variables
used in the gross correlation.

(1 - Rk2) (1 - R )
Rij R4-

Since, in order to be consistent with the two related correlation coefficients,

the correlation coefficient must lie between the limits given by equation Li,
the lowest upper limit and highest lower limit are established for all related
pairs, and 'We average of these two limits is taken as the estimated correla-
tion coefficient.

f. Monthly streamflows missing from the records of the various stations

are estimated for all stations for each month in turn. Accordingly, when-
ever a missing flow is being reconstituted, there always exists a valid
value for all stations already examined that month and for all remaining
stations in either the current or preceding month. For these remaining
stations, the current value is selected where available; otherwise the
pr-e~iin- value is used. In order to reconstitute the missing value, a
re-ression equation in terms of normal standard variates is computed by
selecting required coefficients from the complete correlation matrix for
that month and solving by the Crout method (See exhibit 1). The missin7
value is computed from this regression equation, introducing a random
(-omponentequal to the nondetermination of the equation, as discussed in
the streamflow generation procedure.

7. It has been found that valid use of the regression technique requires
that all correlation coefficients agree with the data that will be substituted
into the equations and that the correlation coefficients be mutually con-
sistent. Inconsistency in the correlation coefficients causes the dependent
variable to be over-defined and is evidenced by a determination coefficient
greater than 1.0. If this occurs (because of incomplete data), the inde-
pendent variable contributing least to the correlation is dropped, and a
new rerression equation is computed. This process is repeated as necessary
until consistency is reached (which must occur by the time that only one
independent variable remains). In order to make the correlation matrix
consistent with the data matrix, all affected correlation coefficients
are re:omputed after each estimate oF missing data.

h. Normal standard deviates are then converted to flows by use of the

following equations;

I,!~
 tim = (g,/6)(Km - gi/ 6 ) 1 3 -1 2/g (12)

i X + ti,m i (l-
Xim Xt~~

O = AntilogX - qi (14)

imposing the constraint:

Q (15)

i. When the set of flows is complete, all correlation matrices should

be consistent except for truncation errors in the computer, since the data
arrays are complete. Any consistency of matrices obtained in this manner
or of matrices read into the computer will result in determination coefficients
greater than 1.0. If this occurs, consistency of each correlation matrix is
assured by first testing all combinations of triads of correlation coef-
ficients in the current and preceding month for all calendar months using
eqmation 11 and raising the lowest of the three coefficients to obtain
a consistent triad. The test of consistency of each complete matrix is
i-.,e by recomputing the multiple correlation coefficient. If this value
ic greater than 1.0, further adjustment is required. Such further adjust-
.ent is obtained b:i introducing a coefficient, successively smaller by
0.', on the radical in equation 11 and repeating all triad consistency
tests until all matrices are consistent. If consistency is not reached,
coefficients in each inconsistent matrix are moved toward the average
value of all coefficients in that matrix until consistency is reached.

j. Generation of hypothetical streamflows is accomplished by computing

a regression equation, by the Crout method (described in exhibit 1) for
each station and month and then computing streamflows for each station in
turn for one month at a time using the following equation. This process
is started with average values (zero deviation) for all stations in the
'irst month and discarding the first 2 years of generated flows.

I I

K, ,J = IK ,i + K ,l+2
2 ." " + -1 ' , -i + -1Ji,'
K +

R iKl l . . . +Bn~i
K l + - j zi j (16)

5
 in which:

K = Monthly flow logarithm, expressed as a normal standard

deviate
= Beta coefficient computed from correlation matrix
i = Month number
j = Station number
n = Number of' interrelated stations
R = Multiple correlation coefficient
Z = Random number from normal standard population

k. Maximum, minimum and average flows are obtained for the entire
period of flows as recorded and for specified periods of reconstituted
and generated flows by routine search technique.

1. Provision is also included in this program for use of the generalized

model requiring-only 4 generalized coefficients for each station (in place
of 1-) and one generalized correlation coefficient (in place of 12) for each
pair of stations, in addition to identification of wet and dry seasons for
each station. These are defined as follows:

(1) The average value of mean logarithms of flows for the wet
season (0 months). This value plus 0.2 is applied to the middle month
and the avera-e minus 0.1 is applied to the other 2 months.

(2) The average value of mean logarithms of flows for the dry
'
7cazon , mcthz), 'Ai iz api,1:c t', all 3 1,ry .o-ths. !lean Iogarithrs
for months between dry and wet seasons are interpolated linearl,.

() The average standard deviation for all 12 months. This is

anrlied to each of the 12 months.

() The average serial correlation zoefficient for n1l ! mon+'s.

This value minus .15 (but not less than zero) is applied to each wet-season
month, and the value plus .15 (but not more than .98) is applied to each
dry-season month. The average value is applied to all intermediate months.
(5) The average interstation correlation coefficient for all 12
months is applied to each month for that pair of stations.

m. Because of limitations in computer memory size and because of

increasing change of computational instability with larger matrices, the
number of stations usable simultaneously in this program has been limited
to 10. However, the program can reconstitute and generate streamflows for

6
any number of stations in groups of 10 or less. It will ordinarily be
desirable to include one or more stations from earlier groups in each
successive group in order to preserve important correlations. In
.additionto providing flow data for all stations, it is necessary to
lesignate NPASS and to follow each group of flow data with a standard-
,orrmat card with NSTX (number of stations in next pass that were also
uned in preceding passes) and station identification numbers for those
stations. These numbers must be listed in the same sequence as their
lata were arranged in earlier passes. Data for the new stations for the
new pass should then be read. None of these flows can occur in a year
later than the latest year for which flow data occurred in the first pass.

I. As soon as flows are reronstituted for any pass, they are read
onto the flow tape. After statistics are computed from transformed
reconstituted flows, they are read onto the statistics tape (after
ilentification of stations in the nass for future reference). Final
regression equation data for each pass are read onto the same tape at the
same time (for use in -eneration later). For each new pass, the flow
and statistics tapes are searched separately for data for those stations
already used that also occur in the new pass. In order to read and write
intermittentlv and alternatively on the same tapes, it is necessary to
'een traf-! of tape records so as to assure that any read statement does
not read beyond the record mark and so that new write statements occur
at the end of all previous write statements that are to be saved.

o. Once that statistics are put on tape, they are retained throughout
the reconstitution and generation processes. Flows, however, are saved
only for the set of data in which they were reconstituted or generated,
:mntil the last pass for that set is completed. In the generation process,
it is necessary to save the last flow generated for each station in one
set for ,,se as the antecedent flow in starting generation in the next
set. These are saved in the QSTAP array with subscript ISTAP.

INPUT

Intput is summarized in exhibits 7 and 8. Data are entered consecu-

tively- on each card using a simple variety of formats to simplify punching
and handling cards. Computed and generated flows cannot be 1,000,000 units
or larger, and consequently must be expressed in units that cannot exceed
this magnitude. Units should be indicated on one of the 3 header cards.
Column 1 of each card is reserved for card identification. These are icnored
by the "omputer except for the A in column 1 of the first header card, which
is used to identify the first data card. An example of input is given in
exhibit 3. Certain inadequacies of data will abort the job and waste input
cards until the next card with A in column 1 is reached. A card with A
in column 1 followed by 4 blank cards causes the computer to stop.

7
6. O(JTPUT

Printed output includes key input information for job identification

and all results of computations. Generated flows are put on magnetic
tape, and computed statistics are punched on cards in the format usable
later b the program. An example of printed output is given in exhibit h.

7. OPERATING INSTRUCTIONS

Standard FORTRAN IV instructions and random number generator are

required. No sense switches are used.

8. DEFINITIONS CF TERM;
Terms used ti the program are defined in exhibit 5

a* ?ROPOSE.D 'VITLURE DIXELOTPMENT

There are cases where the model used herein does not reproduce histori-aI]
drou,-:hts with reasonable frequency. Consequently, the model is under
Sontinuous study and development. It is requested that any user who finds
en inadequacy or desirable addition or modification notify The Hydrologic
En oieering Center.
 EXHIBIT 1

DETAILED EXPLANATION
OF
COMPUTER PROGRAM

GENERAL

Much of the program is explained by comment cards and definitions

of variables. Supplementary explanation follows, referring to sections
identified with the indicated letter in column 2 of a comment card.

SECTION A

Correlation coefficients, R, and beta coefficients, B, are in double

nrecision for matrix inversion computation, in order to minimize computa-
tional instability. Correlation coefficient, RA, as originally computed
and stored, may be defined in single precision. For computers with wor
len:-,th smaller than ? bits, many other variables in this program should
be in double precision.

When dimensions are chaned, the corresponding variable (starting

with K) should be changed accordingly, as these are used to prevent
exneedence of dimensions. If an excessive subscript is used, the job
will be dumped until a card with A in column 1 is encountered, at which
time a new job is automatically started. If 5 blank cards (with an A
in column 1 of the first) are encountered, the run will be terminated.
Job specification cards are read in this section.

SECTION B

NSTAX is number of columns in correlation matrix. These consist of

,13TA columns for the current-month values and a similar number for ante-
cedent-month values. NSTAA is initial column number for antecedent-month
coefficients. These are computed from NSTA, which is read in if statistics
are to be provided, rather than computed from raw data. If raw data are
to be used, NSTA is defined in the program later and NSTAA and NSTA must be
also. Data for each new pass are processed after transferring back to
statement 42. In the multipass operation, NSTX is the number of stations
used from previous passes and NSTXX is the subscript of the first new
station for the current pass. Station identification for the NSTX stations
must be in the order in which data for those stations were originally used,
because search of data and statistics on tape is made in this order. Flows
for these stations are read from tape IQTAP, and corresponding statistics
from tape ISTAT. Variables LQTAP and LSTAT are used to keep track of tape
position for subsequent writing.
*Provided through the cooperation of the Texas Water Development Board.

IXHhBIT 1
 Months are identified consecutively by the variable M starting with
the month preceeding the first year of data. Some quantities to be
accumulated are initialized. Station combination data are stored for
the purpose of obtaining maximums and minimum (section D) of weighted
flow values later. Tandem stations are identified for cases where a
check on consistency of generated quantities is deemed appropriate.
Station identification numbers are set to a large number so they will
not be undefined. The flow array is filled with -1 values to indicate
missing values. For each station and calendar month, the total flow and
number of recorded values are computed for computing a flow increment
and other statistics later. The minimum flow for each station month
is also computed in order to avoid negative logarithms later.

SECTION C

Station data can be read in random order. Stations are identified

by subscript in the order in which data for each station are first read.
The year subscript is computed. Negative subscripts will occur if data
are for years earlier than the starting year indicated on B card, and
data for these are rejected, with diagnostic printout. The stations
are counted and the flows for each month at each station are counted for
the purpose of computing frequency statistics later. If the number of
stations or years exceeds its dimension limit, the job is aborted. The
number of stations is permanently stored in the NSTNP array for later
identification in multipass operations. The remainder of this section
is self explanatory, except to state that permanent identification
station numbers are given for stations in combination, for tandem stations,
and for consistency-test stations, and subscripts are identified for rapid
computation later.

SECTION D

In this section, maximum and minimum recorded flows for each calendar
month , the water year and for durations of 1, 6, and 54 months, and average
flows are computed for each station and combination. Durations do not
span a break in any record. Quantities are rounded off and printed in
fixed-point format.

SECTION E

The logarithm transform of flows is accomplished here. Missing

values are indicated by an impossibly large numb,.r (the -1 used for
missing flows is a reasonable logarithm and therefore cannot be used
for missing logarithms). Before the log transform, the average flow
for each calendar month at each station is computed and one (costazied
to a minimum of 0.1 flow unit) is added to each flow. If the minimum
observed flow for that station month is negative, that absolute value

2 IXEISIT 1
is also added before the transform. After the logarithm transform,
freq:.ency statistics for each calendar month and station are computed.
An increment needed to convert the logarithws to an approximately normal
distribution is also computed as an alternative future transform. Log-
arithms to the base 10 are used so that statistics are comparable to
other commonly used statistics. A variable IRCON is set to 1 if any
missing values are encountered, so that the flow reconstitution routine
will be called later. A variable INDC is set to 1 if the first approxi-
mation of increments causes any one of the skew coefficients to be
smaller than 0.1 or larger than 0.1. In an optional routine that follows,
the increment for each station and calendar month is adjusted individually
and iteratively (up to 14 trials) until skew is within 0.1 of zero.

Stations with less than three years of data for any calendar month
are deleted, since skew and correlation computations require at least
three items of data.

SECTION f

Correlation matrices are computed here for the purpose of adjusting

frequency statistics for short-record stations. All correlation coef-
ficients are first set to -4.0 in order to identify those not computed
later for lack of sufficient observed data. Then accumulations of the
various quantities required are computed for all items above the main
diagonal in the correlation matrix for each month, using all data common
to the two stations involved. If mor-e than two items of data are avail-
able, the correlation coefficients are computed. Coefficients for the
main diagonal are set to 1.0, and those below the main diagonal are set
equal to their symmetrical element. Coefficients between the current
and preceding month's values are similarly computed. These items con-
stitute an extension of the matrix to the right, which doubles its size,
and the new portion is not necessarily symmetrical. Similar complete
arrays of average values and root-mean-square values for only those log-
arithms common to each pair of stations are found for later use in
adjusting statistics.

A search is then made to determine the station that would be most

useful in adjusting statistics for station months with incomplete record,
and the means and standard deviations are adjusted in accordance with the
following equations:

s = SI (S'2 -s ) R2 S /S
1 1 2 1 2

-X' X. + ( - ) Rl/S2

where primes indicate long-period values, subscripts are 1 for the short-
record station and 2 for the long-record station and,

3 EXHIBIT 1
 X = mean logarithm

S = standard (leviation or' the logarithms

R = correlation coefficient.

An optional check of consistency of standard deviations between

adjacent stations for the same month is next made. This is to assure
that frequency curves do not cross within three standard deviations
from the mean. If there is a conflict, the standard deviation of that
station designated in the input data as the dependent variable is modi-
fied accordingly. All frequency statistics are then printed out.

SECTION G

All flows are next standardized by subtracting the mean and dividing
by the standard deviation for the month and station. An approximate
Pearson Type III transform is then applied as follows:

K = qt +tl) 1/3-
where:
K = normal standard deviate
t = Pearson Type III standard deviate
g = skew coefficient

New correlation matrices are then computed, based on the normalized

variates and using the same standard procedures previously employed for
correlating logarithms. The sign of the correlation coefficient is
preserved, since the coefficient will be used to establish regression
equations. Correlation coefficients are set to zero if the variance
of either variable approximates zero, since the computation of the coef-
ficient is highly unstable and since its use would be of little value.

SECTION H

For Jobs where -orrelation data are given, the portion of the
correlation matrix above the main diagonal for all months and the
entire correlation matrix relating current and preceding month's
values are read, with a different card for each pair of stations.
Values for all 12 months are contained on one card, and the two stations
involved are identified on the same card. An automatic check is made
to assure that cards are in the required order of columns and rows
in the correlation matrix. When generalized statistics are used, only
one correlation coefficient for the entire year is read, but card order
is the same. Symmetrical elements below the main diagonal are then
filled in and values af 1.0 are placed in the main diagonal.

4
EXHIBIT I
 Frequency statistics are then read, 4 cards per station, with 12
monthly values and station identifications on each card. A check is
made of the station order, to assure proper subscripting. When gen-
eralized statistics are used, only one card per station is read, and
this contains the maximum and minimum mean logarithms and the average
standard deviation for the year. The months of maximum and minimum
mean logarithms are also read and converted to corresponding subscripts.
These subscripts will differ from the calendar month number if the year
used in the study does not begin with January.

SECTION I

This section searches for each calendar month the entire correlation
matrix to be the right of the main diagonal for misiing correlation coef-
ficients due to the nonexistence of at least three years of simultaneous
data for the month. As soon as a coefficient between two variables is
identified as missing, a search of the correlation matrix is made to find
established correlation coefficients between each of these variables (i
and j) and any other variable (k). The range within which correlation
between the two variables must lie in order to be mathematically consis-
tent with the correlation with the third variable is established by use
of the following equati ,n:

Rij = Rki RkJ J(l'Rki2) (l-Rk)2

As each successive third variable with established correlation coefficients

is found, the upper limit of Ri1 is constrained to the lowest of all upper
limits computed, and the lower limit is constrained to the highest of all
such lover limits computed. When the entire matrix has been searched the
correlation coefficient is estimated as the average of these two constrained
limits. If this element is above the main diagonal, the value is also
entered for the element symmetrically across the main diagonal. The
search for further missing correlation coefficients is then continued.

SECTION J

Where a correlation matrix is not to be used for reconstituting data

but might be inconsistent, a triad consistency test can be made in this
section. This is done by examining all groups of three related correlation
coefficients, and testing the lowest one to determine whether it is above
minimum constraint established by the equation in the preceding station.
If not, it is raised to that minimum. When this is done, it is possible
that the adjusted coefficient had already been used in another triad test,
and consequently that previous test would need to be repeated. In order
to do this properly, the entire matrix is searched up to 12 NSTA times,
where NSTA is the number of stations, until a complete search reveals no
inconsistent triad (INDC = 0).

EXHIBIT 1
 A coefficient 'A'.' of the radical in the equation is used in order
to obtain complete matrix consisten-y in difficult cases, whenever
possible by this means. A test for overall consistency is made in
section K, and if this fails. FAC is saccessively reduced by 0.2
until overall consistency is reached.

SFTION K

The test for overall consistency oi the correlation matrix for

each month is made by constructing foi each station the correlation
matrix that would be used in flow generation for that station and com-
puting the multiple determination coefficient. If the determination
coefficient of the matrix for any station and any month exceeds 1.0. all
correlation matrices must be reexamined, since some coefficients are
common to two or more matrices. This is done by reducing PAC in the
triad test (section J) by 0.2 and repeating all triad tests. If FAC is
reduced to zero and consistency is not obtained, an index of NCB is set
to 1 and an averaiin routine is .sed For each inconsistent matrix. A
quantity SUM is computed as the average of all correlation coefficients
in that matrix, and each element is modified by multiplying SUM by the
excess of determination coefficient and adding this product to the
uroduct of the complement of this multiplier and the value of the ele-
ment in the inconsistent matrix. The averaged or smoothed values are
rerlaced in the complete matrix for the month, and this requires some
careful manipulation of subscripts. A new computation of determination
coefficient is made and the smoothinc, process is repeated up to nine
times until consistency prevails. If this does not occur, the job is
terminated. When consistency is established all complete matrices
are .irinted out and essential elements are punched if desired.

SECTION L

In reconstititing missing data, a search is made for each month of

record starting with the first for stations that have no record during
that month (Q=T). When one is found, a search of all other stations is
made to determine whether recorded or previously reconstituted flows exist
for the current month or, if not, for the preceding month. If one is
found, it will constitute an independent variable for estimating the
missing value, and its value and pertinent correlation coefficients are
stored in new arrays for computation purposes. The correlation coeffi-
cients with the dependent variable is temporarily stored in the NVAR
(NSTA*I) column to assure that coefficients relating independent vari-
ables which have sufficient array space (they cannot exceed NSTA in
number). A variable ITEMP counts the number of independent variables
(stations for which recorded or reconstituted data are available). It
is incremented after its set of correlation coefficients are stored in
the R array, and is finally used to relocate the correlation coefficients
involving the dependent variable. If no independent variables with data

6 EXHIBIT 1
 are found, as can happen in the first month of record, a correlation is
made with the preceding vaiue for the same station and that preceding
value is arbitrarily set at the average for the month. The regressi6n
equation and deLermination coefficient are then computed using subroutine
CROUT. The variable having the lowest absolute value of correlation with
the dependent variable is identified, and beta coefficients are searched
In order to eliminate all unreasonable coefficients. In the usual case
where the simple correlation coefficient between any variable and the
dependent variable is positive, unreasonable coefficients are assumed to
be those larger than 1.5 or smaller than -.5. In the case where the
variable correlates negatively with the dependent variable, the reason-
able range is -1.5 to 0.'. If an unacceptable coefficient is found, INDC
is set to 1. If this happens or if the determination coefficient does
not lie between 0 and 1.0, the variable with the smallest correlation
coefficient is eliminated, the correlation array reconstructed accordingly,
and the regression equation recomputed. This process is repeated until
all required conditions exist. The missing value is then computed by
use of the regression equation and adding a random component normally
distributed with zero mean and with variance equal to the error variance
of the regression equation.

As soon as the missing value is estimated a search is made for all

established values in the current and preceding month with which it is
to be correlated, ad sums of logarithms, squares, and cross products
are incremented in preparation for recomputing all affected correlation
coefficients. After checking for sufficient (three years) record and
nonzero variance, the correlation coefficient is recomputed. If the
standard deviation of either variable is very small, the correlation
coefficient is set to zero. If the coefficient is above the main dia-
gonal of the correlation matrix, its value is also assigned to sym-
metrical element. Since estimation of a missing value affects correla-
tion coefficients between variables in the current and following month,
which coefficients are stored in a different matrix, this process of
adjusting the correlation coefficient is applied to those values next.

SECTION M

After all flows are reconstituted, the flow tape is read until the
proper position for writing the newly computed flow data on that tape is
reached, and headings are printed for writing flows on the printer later.
Then the standard deviates are converted to flows by reversing the Pearson
type III transform, multiplying by the standard deviation, adding to the
mean and taking the antilogarithm. The increment is then subtracted and
if the resulting value is negative for a variable with zero lower limit,
it is set to zero. In the case of reconstituted flows, the Pearson Type
III transform is constrained so that the excess of the standard deviate
over and above 2.0 is multiplied by a maximum of 0.2 (if the standard

7 EXHIBIT 1

|* _
diviation exceeds 0.). This simply prevents obtaining unreasonably
extreme values due to sampling errors. It is a moderat'on of the
extrapolation rather than an abrupt truncation.

The test for tandem station consistency is next made. and incon-
sistent flows are identified for printout and changed to the limit of
consistency. The downstream flow is aacie consistent witb the swm of
upstream flows. Flows are punched on cards, if desired, printed out,
and written on the flow tape for use in future passes. NQTAP is incre-
mented and represents the total number of" records on the tape.

SECTION N

After convertinri deviates to flows, the frequency statistics are

recomputed in order to aree a: eurately with observed and reconstituted
data. If a consistency, test is rcalled tor, the variable ITRN3 is set
to I and computation is transferred to near the end of section F, where
the test is made ancl the transfer index -auses a return to this portion
of the protrarm. Aj!stel statistics are printed, and the consistent
correlation matrix is rrinted (and, if desired. also punched) by transfer
to section K, usi-in- ITPN., as a return indicator again. The statistics
are then punched. i" desired. P'lows for the snecified station combina-
tions are then comouted.

SECfTION 0

Maximum and minimum recorded flows are computed by transfer to

section S, usin,- ITRN.-=l as a return indicator. The variable ITMP
keeps a record of the remaining years whose maximum and minimum flows
have not been searched -et.

Next, generalized statistics are -omputed, if desired, (if IGNRL

equals two). As indicated, straight averag-es of all 12 monthly correla-
tion coefficients in every category are taken. Means are averaged for
the three wettest consecutive months and the three driest consecutive
months and the seasonal timinc, noted. Standard deviations for all 12
months are averaged. Generalized statistics ac'e then printed out.

Next, generalized statistics read in section H are used to compute

reouired arrays oi statistics. Skew and increments are set to zero.
The mean for the middle month of the wet season is .2 higher than the
wet season averare and means for the other two months are .1 lower.
Means for the dry seasons are uniform, and means for the transition
seasons are interpolated linearly. Correlation coefficients for the
dry season are .1' higher (constrained below .98) than the annual
average, and those for the wet season are .15 low (constrained above
zero). All of these operations are in accord with the generalized
model developed in HEC.

8
EXHIBIT 1
 SECTION P

After obtainih- monthly statistics and correlation matrices, re-

i'ression equations for ehch station and calendar month are computed.
Flows are generated in the station order in which data or statistics
are read and are generated for each month at all stations before
proceeding to the next month. Flcws at each station are correlated
with flows of the antecedent month at that station and at all stations
for which the current month's flows have not yet been generated. For
other stations, fl1ows for the current month are used.

Regression equations are comouted in subroutine CROUT. If any

correlation matrix 'orned is inconsistent (which should not occur at
this stage, except for truncation of computated intermediate variables),
a transfer to section J is effected, and consistency operations performed
on all correlation matri.ces. After such a transfer, all regression
equations must be recomputed, since any correlation coefficient might
have changed. After this, only the beta and alienation coefficients
need be retained, in addition to the freqoency statistics. Tn +.he -,,ti-
pass operation, these are all written on tape I3TST at this point.
3iXTION{.

A routine for nro'ect-i_: histori,.al seruences into the future is

employed here. Values o- ,J RE"V (previous month's deviate) "or each
station is determined as the trans'orm of the flow for the month pre-
ceding the first month specified (by input data) to be generated. The
var-'ele M.1Ais comnutel f'or the subscript o" ; that conforms to the
'irst month of uroj ected "'lows. If the Droiected flow routine is not
to '-e usedl the computer is next set up to yenerate two years of flows,
at the end of whirh synthetic sequences will have a virtually random
start.

In the multipass operation, stations are identified and all neces-

sar : statistics are rontained in the order needed on. tape ISTAT. In an.-
-ass after the first. "lows -,enerated in earlier passes for the same
nernV (the same sec-en-e o" J]ata) must be read from tape I-TAP, and this
tsne nust he rewornd before each pass in order to permit a complete
sear .i.. In any seqtience after the first, the preceding flow for the
'irst month to !e generated is the last flow in the preceding sequence
ur that station, and these are saved in the QTAP array for multipass
operation. If' the multipass feature is not used, all necessary statistics
and 'Jows for -eneratin are in memoryr.
C]rTIO 1

In startin' to renerate flows, a variable JXTMW is used to identify

he ,rear n,,unber o' the "irst -ear of each sequence in the multipass

0 EXHIBIT I

m m
6o .mmmm-m mmmmm im ii
operation. Variab1les AVG and SDV are used to ,-omoute the mean and
standard deviation oF' the deviates "or ear-h flow sequence. These are
later used to 8djiust all deviates so that the means and standard
deviations in ever, -generated sequenc'e will be the same as those of
the historical sequence.

Variables JA and NJ are set up to correspond to the first and

last year of !eneration in each suctessive sequence, depending on
the type of operation. MA has alrea,1. been set up as the subscript
of' ' correspondin- to the first month of flows to be Tenerated (for
use in projecting7 historical flmws re.orded to the current time).
<PR3EN jor each station has 7teen identified as the previous month's
flow for that station. Plows are then enerated for each station,
using stored re'ression equations and a rnndom component. Each generated
flow is immrediatel' entered into the QPRV array, because its preceding
fl.-ow will never a-ain ")c used in that pass.

Tn t'-e multipass routine, flows (as deviates) are written on tape

at the end of ear.h nass, and the last f'low for each station is stored
in the :"STAP arra, 'or use in the next sentuence.
. more than 1 _Gears (an arbitrarili selected length) of flow

are beinf generatel in any! sequen!-e, deviates are adjusted so that their
nean is zero and variance 1.0. Their unadusted mean and standard
deviation a-e printed. Then the,, are transformed to flows, and, if
c.alled. for, consisten"- tests between stations are made. For variables
with zero natural limit, a rheck lor ne-ative values is then made. Flows
are then printed and. if desired, punche(]. Plow combinations are then
cornputed.

SBCTIO 3

Before comptinj meximum anr riinimum values of 7enerated flows,

a nositive value of JX is looked for to assure that flows generated
are not to be discarded (the first two years generated for a random
start). Also, at least NYMG years must have been generated before
maximum and miniium values are -om-pnted (this applies only when the
number of years remaining for 7eneration in the last seauence does
not equal I7%Y(G). Ma9ximum sums are initiated at an extremely large
negative number and minimum sums as an extremely large positive number
(T). Then a routine search of flow stuns Por the snecified durations
at each station is made for the sequence, and results are printed out.
Since this routine is used for reconstituted flows as well as for
generated flows, a transter indicator is used to determine whether the
next step is bac'- to the reconstitution routine or the generation routine.
If the latter, a check is made for the multipass routine. If all passes
are not completed, a transfer to section t, is made. If all passes are
completed for this sequence or if the multipass routine is not being

10 EXHIBIT I
used, a ,he: is -nade o" remai . ii, errs to e 'eerated. I; *rester
than zero, a transfer to section 4r mnade Pf'ter ad.lustin, :-ears yjet
to > e generated . Othei.ie the 'oh is e-ided and a new lob. if an,,,. is
started.

RANDOM NUMBER FUNCTION RNGLN

This random number function is for a binary m;chinv ;and the constants must be
computed according to the number of biLts in an inteoer ,7ord. The numbers
generated are uniformly distributed in the interval 0 to 1.

The function is called from the main program by a statement similar to the
following:

A R Ci., (I.:)

Wherc A is some floati ;,,inL variable name and TX is som, integer variable
name. The argumenL name IX need not be te .T%,, in thie main program and the
function. The argument must be initiol Hoed to zero in the main progra . The
location of the initializing statement is important and depends on tho results
desired. If it is desired to have differejt sets of random numbers for each
of several different sets of computations (jobs) thait are run sequentially on
the same progran, then the argument must be initial ized at the very beginning
of the program and never reinitialized. IT it is permissible to use the same
sequence of random num'ers for each job, tie argument must be initialized at
the beginning of each job. The advanta'o of this latter option occurs when one
of the jobs must be re-run for some minor reason as the same random numbers
will be used and the results will be comparable.

Three constants must be computed by the following equatior.s:

Constant one (Cl) 2( + l ) /2

3

Constant two (C2) = 2" -1

Constant three (C3) = l./2.

Where: B - number of bits in an integer word

The constants for some of the common computers are listed In the following table:

INT J:GER CONSNATS

COPU'fER WORI) Cl C2 C3
GE 200 Series 19 1027 52/4287 0.1907348631-05
GE 400 Series 23 4099 8389607 0.1192092901.-06
IBM 360 Series 31 65539 2147433647 0.465661287E-09
IBK-I 7040 and 35 262147 34359738367 0.2910383046E-10
7090 Series
UNIVAC 1108 if
CDC 6000 Series 48 16777219 281474976710655 0.355271367811-14
11 L(HIBIT l
 April 1960

EXHIBIT 2
Crout's Method

One of the best methods for solving systems of linear equations

on desk calculating maohines was developed by P. D. Crout in 1941.
This method is based on the elimination method, with the oalculations
arrasged in systematio order so as to facilitate their scoomplisbaent
on a desk oaloulator. In this method the ooeffiolents and oonstant
toes of the equations are written in the form of a "matrix," vhioh is
a recatogular arna of quantities arreae in rows and coluns,

The method is best explained by an example. Suppose that in a

multiple oorrelation analysis it is required to solve tLe following
system of linear equations to obtain the unknown values of b 2 , b 3 , b 4
snd b5 .
2 b2 + E,2x b .,- b4 + 5 -
-2x5 x ix
22 b + xxb b5 12
2x 2 +b E 2 b + Ex b + Exx
X2 2 3 3 314 4 3t5 5 1l3
ZY b2 + LExx b3 + 14 b + lxx b 5- EX1'
Lix5b a3x 3 xx b4 + Ex2 b Lx
2 + b 5i ' E4.

For si plioity lot us replace the ooeffioie-'q of the b'e by the letters
p, q, r and a, and the constant terms by t etter t, using subscripts
1, 2, 3 and 4 to denote the respective equa. onss.

p, b2 + q b3 + rl b4 + 1 b5 ti

P2 b 2 + q2 b3 + r2 b4 + s 2 b55 t2

p3 b2 + q3 b3 + r3 b + 83 b5 - t3

p 4 b2 + q4 b3 + r 4 b4 + % b5 - t4

A continuous cheok on the ooputations as they progress may be obtained

by adding to the matrix of the above system a oolmn of u's, such that
u - p + q + r + a + t. The matrix and oheok oolumn are written as
follows#

EXEMBlT 2
I
p1 q, 1 ~ 3

P2 q2 r2 2 t2 U2

P3 q3 3 sr
3 t3

p4 q4 4 %4 t4 %2
The elements pl, q2, r3 and e4 form the "principal diagonal" of
the matrix. Ekamination of the original equatlons shows that the oo-
effioients are symetrioal about the prinaipal diagonl, i.e., q, - P2
rl -p, 'r2 - q, 1,- p4l 2 " q% and 63 - r4 e.
This in oheracteristi of the system of equations to be solved in any
multiple oorrelation analysis. Because of this symmetry, the oompute-
tions are considerably simplified. While the Crout method my be used
to solve any system of liner equations, the oomputational steps given
here are applioable only to those with symetrical coeffioents.

The solution oonsists of two parts, viz., the *amputation of a

"derived matrix" and the "baok solution." Let the derived matrix be
denoted as follows;

INN
3 K 3 5 3 3

P4 Q4 R4 S4 T4 U4

E)XIIT' 2
The elements of the derived matrix are computed an follows

P 1 a Pi P2 = P2 P3, P3 P 4 -P 4

QR _r,_S T U
= u
i 11 p1 1 Pi 1 "

Q - q2 - P2Ql Q 3 "q3- P3Q1 R2

St 2-T12 u2-lP2
" - l S4
S2 - 2 2 2 " "

R r3 -QR -PR R4 r R2-P, S Lt4

33 32 31 r4 3
3 'RRP

T P
' t 3 -T 2 Q 3 -T1IU 3-U21-U1.
33 3R
t-TR-T2Q.-TIP u3 -U2R- U F
5 S R U " 4 R
-s8 -4 4 .U 1
4 4 4 3 -Q 4 S2 - PSl
44
ss-4 - 4 - s -p s

The general pattern of the above computations, which my be applied

to a system containing any number of equations, is as follows&

(1) The first column of the derived matrix is copied from the first
column of the given matrix.

(2) The remaining elements in the first row of the derived matrix
are computed by dividing the corresponding elements in the first row of
the given matrix by the first element in that row.
th
(3) After completing the W row, the remaining elements in the (n+l)
th

oolumn ere oomputed. Such an element (X) equale the corresponding element
of the given matrix minus the product of the element imediately to the left
of (X) by the element immediately above the rinzipal diagonal in the same
colu as (X), minus the product of the second element to the left of (X)
by the second element above the principal diagonal in the sam column a (X),
eto. After each element below the principal diagonal is recorded, and while
that element is still in the oaloulator, it is divided by the element of
the principal diagonal which is in the earn colum. The quotient in the
element whose location is symetrioal to () with respect to the principal
diagonal.
3
tEXIBIT 2
 (4) When the elements in the (n+l>- column and their symmetrical
counterparts have been recorded, the (n+l) t h row will be complete except
for the last two elements, which are next computed. Such an element (x)
equals the corresponding element of the given matrix mimis the produaot
of the element immediately above (X) by the element immediately to the
left of the principal diagonal in the sam row as (X), minus the product
of the second eiement above (X) by the second element to the left of the
principal diagonal in the same row as (X), etc., all divide, by the
elmnt of the principal diagmal in the 1 row as (X).

The check column (U) of the derived matrix serves as a continuous

oheck on the computations in that each element in the column equals
one plus the sum of the elements in the same row to the right of the
principal diagonal. That is,
+
U1 1 + R1 S + T

U2 I + R2 + S2 + T2
1

U - 1+ S + T5

U4" 1+ T4

This cheek should be inAe after completing each row.

The elements of the derived matrix to the right of the principal

diagonal form a system of equations which may now be used to compute
the unknown values of b2 , b V b4 and b5 by successive substitution.
This is known as the "back solution." The computatios are as follows:

b5 = T4

b4 - T3 -S3b
5

b T Sb - R b.
3 -2 2b5 - It

b2 T S1b5 - Rlb 4 - Qlb

3

It is very important that the computations be carried to a suffi-

cient number of digits, both in computing the coefficients and onstant
terms of the original equations, and in computing the elements of the
derived matrix. It is possible for relatively smal errors in the
coefficients and constant texias of the original equations to result in
relatively large errors in the computed solutions of the umkznna. The

4
MXIBIT 2
greatest souroe of error in computing the elements of the derived matrix
arises from the loss of leading signifioant digits by subtraction. This
must be gmrded against and can be done by carrying the computations to
more figures than the data. As a general rule, it is reoommendod Vkat
the coeffioients and constant terms of the original equations be carried
to a sufficient number of decimals to produce at least five significant
digits in the smallest quantity, and that the elements of the derived
matrix be carried to one more decimal than this, but to not less than
six significant digits.

5
EXHIBIT 2
 aooi--4~ormC * 0 0 0 Qfl( O 00 0 0
- aP%* a0 *Nt.
Y J O* .Y ** * .Y tf-N
- N . I-L 00m -'* UC n

tinh'-P 0 'O 't 0O4e

*.k *l in In6ka0 0 -t
- .n a 0 * * COO.- in 0 ow *0 0 in a;-.
- ON 4 -ow% f-Ni M. * 4-
CU LN)A -

0 0 *'0O'4 4 i . -Of 0~4

~N
M o440 -p
"~' 0 ~ 0- - 4 I
Nu - N
1 NO ON tJP- Al.. 0 04(MP.- pi os .
soat 0 M '* &
0O~sf~.4 -0~
%~0N 40
(

N -z
M. .4 -6 N -I
0n -n r- i I

zI) a; 4 0 N0 M0 OM sg 0 E) 0 ... * NI
Z "-4Cinb~i to M CN42 zn in13.Oi
M Ajm( 0 (Y .1 N In .4 't NM 4.. 14 m
Anjo. M zz o
() . t I -F'*-.amN. x .9
t) * 04N4- M'4t -

CO6 * 0 0 N
M.4 -0 OD*.
)v- Nr-U .n0 0 * k C)-0
COflU) - NP-inV 0 0 6 0 6 : C;00.-
*9 9)**M4)'E")0I tIflOO 1 - * . *t'm4t')mc L~n 0 C
00%) (v-44
N N)- 39 -- 0'D ." .4. '- -V
I-..JJ -P UI.- (7k t- _j Z n in (1 Lfl I" O
4c L. Z .4.j %0 D O000e.
CCe q U N '0 C O0 000
0 3E4 - 1-- -0*
'O 0 0 0 * c W'D~ oN'O
OZU -( 0 * o a 4b a 00
a . .*f'C'OLIn.%-N4C 4W *&
9 *M.OOMN -N 't
LL))
I.-u *N. - I- w a 0'( m'N N.4
t ..
#-d 3 04* "'oM Q I--U) Z 04
2 0M
(Ai C3 Mo Q 4m * 0 0 0- Wf#A Mi 0 QU' -I M
Z C0 f-NP4 1 0 0 *P)N'OE'M.- rNt 0 o * a *0D
rtn - * M tn t
M-u .4
04-PP- P.' 'r'*P-n "-1s12a

o z *4MGOt-e # 9 o0 1 40w
-M
0- 4 .o
-4 NI- N ly 0P. N 4 NP- c Ok I*- .44 N 4
4PI .- Onnte go*-100 4'OP- tnfltIn. C 0 *0
*CMI
* a a a * y ' E ')MOD0 0 o *
0 f3
a *'''f,%C
f, r-_ . - t) 0 0 .i r-
,a"N
m
(04 - 4
P 4 in
M %ain' NoP-0*i4nr- 4 Lfl%'OP- 4 nNo f. - f 0--
'
~
o ooooooa0 cocoeco coa
CpOpC
-4 o 4.4444"- 4 4--- -0 -4 -4- 4---4"-4-4P-4.4
0C)0-0P 0@ .- " r- p- f-coco0 a --
-444---4 -- 4P-4-4--4 -40 -D 4-"--4- -4--4
444 "O "0P4
I*Xw X"0-0 04.4 II e- Z I~-ZI
 *~w 0 NO . .ft is (P
intU

P. kA %D 4* p~~f~4
i % 0 .O* 0 00om w~

,; 0 N0P .sgNOa-4
o 0 . . S 1l0 0 0 E0~ 0

'1- NoL i f
4 04 *04 00 f M, C
4tf0
* ofq-N-0&f 4 0 -oNm ll 00 *0
;-4--o a; ;---
- - rt -Q a) Nn
cocoy
-GNU-~~~4 NoDs @GNa0I
4..D
P_4P~'. N
ON4'EY0. CY .Y ._ D *l 4%C
Ny
to4 0 4O-0C N DM0 YA
- cc0
i ODI- cot~A) fn.
Go.. 42
. D . cc
- 0 fG.00
.4 . .O~l7 0000 * C30
*0oa
tn flIn DU t- to tn t 40 C
.o Z0C L 4% D% 9-4 00 %0ODD

z UfLn C.
QzI.- Lflukn 4. z - .04aN C
o) 4. z z z *N m r-bm 0 4 0 .9 . .Dz *Nt1
Nl~ 0 m 17a.
41 Ck .$ %DOD0
C0 4 1.- - ODO
-%D0. 40 - * 0 0
044 44NmN'4m U0 .0 4 - 4N@N 4 tc) 0
_j 0 o .j0 0 -10x Y IO~-coC)
I o"Z 14 *" 0 1- 0 s~n A c Cpo4Dtaf.4 0-....
I 0 a * * .40Q00 0 1-NM * a . 0
AD OD0 1C-- a0
O cOr -
x4 " '4N4 - M - M4 X 4c"" '0 4N I-)(bM(4-q M X * *
hfMN 4 M 00 M- -M I 0
-4
,1i I- N 4 M 0 I _jJ-cp
Wx w Itr JDM mx M 4 com N XIA. 00 a3
l0
- .0 ) NOY 0 OD o o 0 P-m-" LJ Nm0-44 00 o0 0 o P.-e "00 C aa0
4A 0 0 0 0. D - 0 0 flC
)z * &* * .1 - 00-4 (I0UW0 401
I*I g'4
m~ 6nr) o mLf l I* M fn (IlfL
(0J-)4 0C I U 6 0
it.CL 4.4N m 4 - Is0 .4." (m - 3040
40C 40 oc00.
I- J 3Ut in 'A'lG
a- t % i-_ w fln cmm % in- 01 #- 4iU) cm I-Nomr-
4LL.0 0 C0N0 0 e o oee 4b.. 0 0 CD 000O5O.& . 4 L. - 0 0 OF-I 4
D x .1 * e0 0 0 f.NO'%at- 0 x_3 f Z C; L; Z*N0CIDOi-*-t1X-I--
AtLL moc
'O Lfnq-4 4
co00Lq - N 4 4 .9 S *
w-a 4N- -y 1- w 4N'I-4
Ui
- O004 -4'o0t Q I. -a c404 -N'o M co Lai - NC 00
fr-Ut~ 0 OMt'%OM bU 0
0W... ~ ~* lE' 0 * an
%D NIt"okin inLf
"-4 .a . I~N.I' i.- 4 ** .4-
4
~ N0tI~~~

1-- 00f "r--f.- -4 I-i C04 - r- P- 0%- 4 #.- -3-C

Z C04in 0 r) 0S0 ZOC .4 M C M0 0 zZA 4II ln-N
0 OM -Y OU 0 M -Y 0 k y74.ek
o) 00
N4 m_ -?k C 1 _ -4r o0(04
N.Ln M 4 0 N N n 0. P4 N D00*.Pk
~
M_.~~~ ~ ~ Td74 1-4
r . .f'o .

m 1i40atQ P tf't PkPkk4 el %04 k04fl 0

-41.
----- 4-4-4 -44
a-4 44-4- - -4 *fF_4Q
cacao---1 00 0404- -4 Cb 0 0-r-o
.0. 4 " W --- 4 -4444444 -4 a4a4a4.."a.
-444----444 444----44p. 44---

ML ~ IuII40UIIII..40
 * 0000000.. 0 0 0 )

0 P-4
Q%1,k OD Ln 1- M 0 o 'N
(y,4Lf
0%
(7 M 0%0 (.4"
, a-lJQ'OM .VODI-n S I-

0%0%t~3%0%COm..4t.. E 4(~
0-4
0 0 0 0 0 0 0 * 0 * *
0 **

CO
7113 * C C - r I --D 0
jmN
a, 010 il ~ (3 Y p-.NE~00 ~ N% 4 L
* 0 0 00 00 0 0 0. 0 0
N

~~00
* 0 0 0
% m fn n 0- so M 0> co 0 *0Lo

I N I N 1
> wL
0%M -- M W (%cy ej0 f. n 0 0 T
10M-4N0'0 o%00U,Mtn emao 00 40a, 0 z L
*~~~ 0 *** * a & 4 1 Z-
I N
r.4 mo
N z

* * 0 0 0 ~00 0 a . * x 4U

0~~ 0% t- -

-TU n aoa
~- 0l~J(A0
L .0o
I~ r-' Np-%D 00*OD I...- . 0-
- a. m m rn (n m 0'J- 4'- 0---4

0 -~ - 0 0
N% or't~.-
0 ~ *m~LA 4LAO4'
1-
- 03 0 * Co *n *DF- *, _* (P a
-A - S-

LA0- 00%0'

LAco tao 40 0 z N
%0t

*' 0 0 D CC OD N 0 0K 0 j 0nLI M m0 0 0 S

'--' ( N 0D(K -D0N.. o'

* S0 00 0 0 0 0 0 00 0 0 0

C*
CD CDI -D 0 -P 0r 4
-4'-N40 % U~o..r~o,.
%~r~) -40Lf -qL r) L

0% O'O''O~~No *-44N *NLA ~X L ) 4 O ' '

 ina
a,~* 0P o 4C.4
opN a,-N I C -a
4n a,04o o O nN m 0 .

1. S.', . w 0. -4i C P4N 0'O

0fta..O 00'Ca 0

N0 ON mM4 004 U"

0~- fnfl a~

0
f-s 47
044t f*Om n ' ' -

of a-- C. In

o-,a f" 0. I

13 P.4t. IA 0 '4

w . .n . 4.
%

a47

-40~~o 40M-
F-'N t-i-
-
lo 4 0 t 0. *lMN4
ki4 . NNM" O t n c

IA 5 4 4 4 4'.% 4

M Qa 4)Ia on -. I.. 4t4 ".j

.0 '
N' InaN r-7P-' r- in
cm C.4 O0 .00.44 f"5P. 04N I.P. Na
04 NO'(P' .4
-4 4 N Nt

cm A# (.

%0 -4' 0 aN O PN4 N

w ww .4 o ui It N
-4 N a
C. 9 C A4 - u
8--In 0z . w
4A 1- 3.t O4S ' 3- a4 at4' V. M Z I.-

7,-A

411K~PA 44P4 II ITNIl

 0 N
r-4.N So0 C"2

.. 4 C.
'o D 4.

o ,m , coy

v I. *0 .a
NN0. m

14fj.
'a 0m 0'~ 0z

-0 41 N

04 ,A 41

-~ 'JN ~ tn '~r.f r II ..

-. . . 4 1N4 N
-Aif U.4 K v a NJ r0 P

b..3 ~~ ~ C' 0'

Ge

.4 ;9 .4c

FACI 2EMUIT
 4N~~ iner-

NA-

w *j
4 .0 ~ t-
4 1-**
* -4 W
* r- a.

10O' IVNI e0w

in- N-n4)0
-

M.4t4 "rWQUU%. .04.-'

4
Its~4
w* 9%w
P44- 5 * -t'.,

* 4 4 0 4 . " -4 Nu 0-N

Wl go I" ff VlOi
.4 dJ'

ONMO 00...0
N40

OQ4
- N r0

piN %F N0~ NI -

r-P A .AVA
a- It~ t$.. -
-. .4 0i N-IiI C4 . C-i4 N fm N

04 -4 - -

tv- r4 0 10iN c no%

4 . -y0
c 4 -4 I -4 0 -4:"

,n4 0 1-1o
...
-
~~Ii~~~~~
, .4.4)
ft W... 4
*
b, a-%4-f"
w intIa
N
N
Op C #A~.N
0 mC

"d -v1nW -a 4. .in4.4- N -

At- 44 - -- m C

44

InC. #>-4 -4

its
1.1 1 Az 4 u k
W% 1--W1$
ft,~~~' c 4 U W . 4 UZ1

:u 44 4A 4

'aBZ PAC
 04W.04% 0 - 4t
't rA4 4'tdNom
in 4.4 IA

0 .Ea F * .M
0 %

~~4 ~ - C.

N I9I

IA0I

31..
NK I
%-
J a N

0 -0$.

v 6CCL 00
in 4 do fA.4

np 0
$A 4 K Cie

-44NJ ge 0inrn 0 -

PAO XtZ 44

LaN
 to olV~ -4004
6-1. OP c c~f"~AA'. .4 ni 4 0V PI
q. on A 24 *0
N0~
e4

tu .4 .4 CPf-

QN oft 04 .d99

C4~ CR1
t4' ar
4I 4v1.4

tv~ 4"- .r

t4. A0 of 0.4-l A4
0 -4 m ta -# 0 (P n s4 co .
.OrOC4 eq N~I m
OP@4" .efdO- m

ra W% i I a- C44

em
eq

& .0 4o10 11%

fa a 0 U

to)

lu
J.
ryON
N41P--
f -

4N4 '4-Ca1)q
441d'@C4- .4 E% Cm .0 . -

-'t InI4gg- lN.

41~D4@-4.
O4'I?.4 E.41444A
.4 vi =.

PAN1-4
 t.. ~~~
i i .
c0r4 ~ -4 0
C-44
x nQ n W%4I

l"
AdQ~ 4) m

44 In4. Ct

-4 F- 0 a (p -

cv0
30 ey4V0

'o P-4inOO
t m4. v 'D -

~e-00At in %A
.4 ~ ~ P C.44 gl

N V

.401 tnc m 0 mtof--#1- 1 i f

-21A-41A~ 0 f~ ~ ~In.F.~. InN4.4

10 I

eq VN M -IN

IN InNf. mmOIU1410 P

m* WN0
m m14l
jg 04C

-4 V4
43 U'

c~c n "1o -0.F

4 4

-4 -1I G

)hI'd A 'N0 N-
)I. bU

-
10 a - . 4 I -I
AL -W.. *1 -
0 W

ME 6 IAl
 4

10 140fINI
-- 0 0 N4.4 A0 0 fn n f" 0 u%4~d
am 0 .4 Ia MM.4.4 NOV

md% f" -* f- m m
4 in in 0SO f"* nA
0 0 0 .:

@4NO NoSIA on coo N QOA

N a 0
0.4fn 4 f" a0.0#
0. A*o
2i .4 . . .i . .4
0.4 f"l 0 in
9-0 #n " o i
In4 F. I nI~oI in IAvr 040 ' m%

ISi N a P.

z CN 0 Ia W 4 I enJI N N0J
w6

4 , -6 W *.4l. NIA* .0dI~ ONfA 1"4

in* aI-C In -*%f 4.in Ar .4 -t..

o" -4 N4 .4 ini p4 44 N 4 q% a 04"

44 0 I'd0 4 N, .4 V4 .1 0 .N 4
0 Iu

of % f f A 1lII em 004p In en m lf--

2a N
0n- U)0
in 0 0I
a4 NAN4 us O -o 40

o..-a l P-.4 .4 .4

.A z4 ot. .A
0CI 0j 0Ad 0 Pc
04ON.4. 444. -44401N4f

.tI11 -4 * .. -4 fA C
04 on. t14 p 4.

zd. o. wl z w is I.-W .5- 2-.49b- zo- a

-10 on0A 4.14 u0.04 4
0 .3u . 3o a z4

144
.4 5-W.

40MR 7 1.4.4T 441

 a-~d rI IQA
b-a-
# f n

E 0 fo

a...~o-.a N u 09.0
- 4 In

,u u-

.40
r- In 70

4
w
Lu Nt c: .tN
P w-4n i .

O(A~~~Nc n(i%4V

4aD. .. a l P. 0I'd'
N It P Ir
ccNlc 49%u 4o f

1-V~~. r4 -
O-4-40

In' U 04jf 4~ * 0
oAON 4 *rsn m r .4 -0
r4 I" e4

I-e

XO-0 go q 4 j.
EU~u At-O-N f- m vA
IV 4- c4 4- cm " in on r4
44 In

3 je4~b-tJ I -

M j A .. 4- ' -
1- :-, -.. v 'y jI .4 ma

.,j C) c . %" -jI- -D

13 a-
wC. Li 7- M s
-- ro v -- I. -u 4 4~ xuuI
-1J IL 4 0 uj '3 !: La 0i at " t
- e . 0IA4 *0 % .~ 43 v.. mVI' tk'I~ c
Li .4 -kAL

4 0P.U #T .4 cmc

It IL

I P"~ NO 1 4
 zooa..
KOIA I^

Co.s x*

oz:

4
1 1

z .4

6-.

us 0

IIt

af 0 C9

In 4n3 .

-0II Cl.m

ac
IIN m N P t

o4 O 0 1- 0-I

. PAO 9 SMU 4
 'Om 00C

,o 40

m m0 41 41

I,-- I n I

'10 10-

4D Q -. - N

w0 0 cwCCDt
OF-N N~J 'WN

lA ~ 1~NO A .r fa co 0AN

* o .1 3 in
a. 14 -f' t I
NI
" 10 1 t In. ~4

z *.7W 'A no Ao.

t- NaU - p 4 Ns
U " Nm

.0~I .S -s

104N -'f 4.4, N40,'.

 -0.0 0'N

04'

M-la

'.a In Mo

In N

'ju
'AI 0 cc N N

IN4

'0 0

C.) 0 -

4.- I.-

cc -.- I

jj 4- j b- r-.0-w%
IC~ tA 0
Q , In
4 I- -o
'U OtWk

ILAo 00 a. z

IA- '7 - Zl
-~~~4 -L - 4- A .

- -44-4

U.Laz ~
 K. a

f- O~mftdO

M.4

00 C
Z e

*~*-4-4-do *N 4N

40 -00104 Wj

4A
P- i1i - % i %r

40~~ on N e 0y n
'0 -4-

.41 4@4u~ 0.t-'.

wC00
13C tn14 4n

In

z~ at

ki I-

NM 0

z Ix

on.,*
f- -I -t4 O
.7 % C4 C .4 -4 -4 -

0.0

44 1

N9 -0140 C

II-oc I.N
* N 3 c woN

0A 12EI K
 4m 1 fft . 4 4q0

~~~4-~~ ae M. i. r
.i .4

. j On -
.4 on
%A~ .

0" dm4mm
#"I

.4~V n04e0 -4

I~W%

IV
- ra A 0 P P

'.-tN OD". -%W

f- w%0DA

.~~ o3 mf.4

r C4 0

~ ~
0~ ,4 cO-
IIn

In 4

ct 4A

11m

43 .O. IL -f- -0
- ~ N

14 Q.-a4

&. w. w
-4 P.n I.r
1
SL. - - 0 0 N 0.
'A 4

11~.0-A %
a~M 3 a

1doo

Kw z
.441.~ 416

PACA Ii Xftryl 4
 0" in

4%.
N 2b
PtNp~owNN0m
a0 A N
N 0 0 .i

tea K MI 44 00

.4N 4

r-N4,r4..4..o I0 I-4'

C4

I" M

M-. cc (1-
-0.4IOIMM
0.P0ca mCjr

- O II %14 M Nn

'o0-' obe 0
A0IU
MI
lt. IM I IA.
f4 (n f" 0

IVI *.'tl

m'0o4aF-.t0 f"Z

IL n p0 UN %A MIN

m.4 M4 4 tn U. Ln '
co 4 a 0

ot..4 I NN . 01
0

144V4 NU Na

C"4J 44 t t 'o0tn

WLC .4 q 0' N I40 CO N Nt

a"
CYf

kS ILh IA
ac

usC 0i K00
4L 4 I- N 30 .e ~ J ~ . 0
'J 4t t" " t
i Q6e -4
to'U ..

PAN 1A 121111 is
 0 o

1 9 0.

oe~

-'"
4em 10 'A:91.4
m.' OQA 000
AD mo d on .4
S - - - I 44
co N 10 0-4** -*
.4

ac. Nl M4 A 4$N'
30 1"'Am Iti'0 N A14A 4
t t- '"
M"A 4
on an N0 4 -f N

C .4b
c to a
imU C, F-N M k"V M 0

wmu
4 Nb O e Ob N4 a,
It N cmN
*I-

In
' ;AO.. 'A.4N
m% r %
'Ae. % 4C,

W% ID illt -t 4n w 4. 0,
.3 0%
N 0 P-4 . C
ON I n V%71 .'t
nS. AN0 0 C~.4 N
CP 11 #APN 0% 4AN13Anc o A -rN
O2t- t..
NO a r 4 ccNy in47N, N0 N %
040 4 43 f r- d ( 4 9 C,**V

CMM . PNNN
Fa. 0,
-j NC %J4 It*7.N.
k4 NI~ N '- N 43

P
o - in . *

4OS X -4 .4 N IV N4
'A t. E10
'A &L 'A
oZ .3 0 1 w 3.4 m A
at
4% 0' Wk4N -AID -%.0M -% O%
-C4Con 7
'J. N .4
I- . . 14)
. C **

14 In 1- N.
jA li j
N Or V. 4 ab
W%
*
N NI ~'A OP M
4IN
4 wN J
-t 4;PNC
IN

WI 11O' it i 'A w 1o*a$2u"9;l 09 1 u

N -
r
VI ~ 'AO..'A I
.4
0 ' 444 *

*~A
13 ITN 4.*~
 4I- 4 0 % 4. a (P

o 0- ta 0t

04 - o9.4 o

f-O cm

so ~ in~ .4f-m 04
IM.4iOC
10.CU.

10
'C'J i~~

0-00N00 0Q-.
04 in in
f4.4 4

,4 w~t.
m WWAJ
CinC
V% .- 40-NY
$OO ~ 0' ~ L *.4 N .4' If%n.N

co- 400'

-o o 0 04 M00t m0
r4-4400 .4 N

40m~r 04-

M to -4Cfnt4 m
Q I De- .
4 40
0-~~- 0-N O"N * O

'i4A I-"0.a 0lO 'i

m -4 .
U m . 4P.- OD .
'At 0 0 %n.

-o. . CP.
r" enp -0f

t z:4 t Q. C, or
13 u OD ~N4 -o 4 C r

eU.i 0 in tnW

- .,
..i
EX
AI
1.-
'r'
:3
JI
cc.
.4 .40L0
fn
0. 00
~O
- r - r T -

U. Ap g. mt-'0
tn -t.

p- 0
-, P tAZ '.3 , 4 0% I-1'A' 0' ~ 0 I. 0
Et

.5 7 z s- ;p>x$
0- f. .4 o 1- 0 .jt U w4
W .9
4WW .4 aUWX
ev deI.
4 - W0 a

41 0.-

~ ~
VI ~
I4 ~ .4 ~o . ~ .4 ~ 4 0.4.' -0
v "a
4 49

IPAW 26 2=131T 4
 A on A

C~.0

~ U)

0
1-J
o N

ev -

Lie

0 0u
414

-9-

.. 4 I li i.

PAC 17 MI4I 4
 N It #AJ*
4 ~i O4

* 0 .. 4

in ev .0N40taf

to 0P ~ ~ 41 .q ~ .0k ~ ~ - 4D

N 4 - n

'4) - 4. a(i
-.

- ~- . % 0 4( 4-( % . .

cnq 0 w Q- r ,% O f.'U 'n.

ly - T f
--

Nvq in P 4i

of. wD
10 Q4 ( "-

~~L" r- - - -n . 0..a'f. 4j

1' 1 7
-4

-J 3 - , t '
U.v

44 -*~Im ~
C-- 0% * Li44--

421 '7 44
W o u P( a ~
(85 ~ ~ ( a-~c~ ;
0%-c.i.
A- -- -4 p-4 * -44444 (-4.--4

1M3
16 inm 404
 CA v

or . ow 4

wle 0 0o- an wm.w -Wc C t- k'l9

a.. 0.e a' Oe * fu c cee *- *PV
0- a o; 6 a

oil'
~oa..~
*-Vk .0 oo0 In
fE'Oi 40
9 n
le--*
0a- V.O
O
l
V
00
@' 4
x 040. 04- *IE
0-0
I- a - a 4 c-oo- - 5-4

4 . 40.1 m NIPN w .0 o, 9 r.Ew

I.- P n IV in a c P-a on. " pqIo n p
449..0

C)e
Oecf. # Of P
ag -4 -
IL
9wwl4,N 4t wl fu-0 Vl toN qf
zIDl1
Jo inuG .0 0050
44 .. *c c. e . .

1~- 0.0 0-
in0 em fu U
41 0 0- v AS
zC -0~ :IVl' 0 0

tv 0'5E 0.E P--4.i44

E 0000 - a 0 t1
ID v
i- 6 DOf In so . *e c
-0aGG ml

z to.*a, -W n *0 4 G 00. 4 co40

WN or In0-5 Vo05i
WA . 0.e5 4q on -4n-p
S 040 In4 414o14 PIGU VI
It41i Ac;
%
G- NO'a inU 45 4 1
-v Iw fw

410 Im 0q
43I1051
xEI.
VIO
on u8-1 P
00 P 4
0111. ~ p

0 c nE aAI No N- N

.0 6

f- k ckO 44 &,N on a -w5- t-fl Al 0 0 3,0Ia fl

z I--
1-5 XI 0.50 04

us- 4- Q .

* ~ 4 ~.1- 14. I5W

9004~~' c0 0. 41

N
4I ~ 41.
~ N
Go
1
G -
C3X4
- u r Iv

4190 PACE 19-DI

M1L O e OIVi0040 09 0
 40 ~ ~P D0 ~ s.ino ~1.n ~ C
l-EW *U l

ui S-a- VE.W

on a I-u r~

1 ~ inA
~ ~
uin. ~ doQrm ~ e v 4
fu 49e ' *mmw4D
.
.0t.Il.4.0iI~ e a-.0 0 0

40. L" in con in a 17 C.

h f owg l0 u Dr i
In-V~,l- ID0~ w M l 04D 0 0v *" 0D
40 a
4w in. .0 MC
O 'D q
t.j 40n0.u

10~ 49
zoloOs 4w a

In0. in tnain1
M* Iii IItfi fno
rn-o3
Rw in ~ a~.. W

V) M 0,4nPP D
0in v P

*03

*0v -604 NJl a-A . It-PIt

$- a
-a ni f S Inna
bo0
fu I44
Vi M 0l .NEo i
3L 10 xW Zis4ae 1MI UP do 04 f on M .
Mon m

or a- 04 Pt9I- in O

T bo .4
"4-~ on
V*0O . . ..
.4 ~ ~x ~ ~ ~ 0 aS .t W q 40

a so a
a 2 0 09 w so
NI4. Nn ow cJ P.0hub 0 .
41 z ~ n 00,

*44 -* -C 4C CC 0

OZ a- 2= U - a-
.4 . .4 .4 .4 . a V

N.~~~PG N*. 3MST i . 20 ~I 4.3

 aft- - -k

C4 43

C00 Upsfl

1111
c4 I l - - '

.C Al

- Wt in

zhi

Ce i

ol V

14 4L in IV P p

z - 6 B- o -on
0~ x X 0
z V in

zz
C, or z h f I l 4 @1 B
11 z IL *- e
O Z 3
om w l1 -a
J- a U aX.. 0 a
OK 4A K tb

o~ ~~~~
.- . * op-

4r4
9% a an inIP ci N
a1 a 0 : orU

U- 4.0 00
9 . .91 4. -
0 0 MW 3.4a
we
300. ve:II vI OP.
.. a~~~ .49 Wl4 n %.a, Inh.

.4 04 ft
00 4 AL.
X 0 : 2

2~~PG 44 480.1 44BT 21 4

 omnow 0
MOM m

XIlMP m 0W

04P
Wm D4 4F -v
a.,. a3< .

,Ble
U' f U)
@fWk
, w'

sW 0in8

ma0
0@O aOM .

.4 ,4 .4

a noo m
.4-d .4

4)0

o . 8 F -.40
0 w

.44 4; m
C,
W
II- w2
do c
do
fu0f.
O
fn on
Z. .

L t
Q 0Q

S, - 0 I
40 41 in m on.

6- ac

0 0-
I- 0 0 ..
el - w -

I
CD S_ IL a- v
OL a z C2W

2 lD Pa
01D us
:1 01 3 5
IuI 0. W III
A Li cue
5-
1- 6- - -5

1 nC
40 a or wk .4 in2-0
0-c V0rb
a-0 .0 M

0 0 0 g X K;
U~I .4 V #A 0.4 -a* on

lot . .I
I .
IIII. . . . . . . . . . . . . ....
pPAC 22K
 a a oe a v4
4rK

4Si*G v-Uop e. 0 3c &ni#

I ou C3"i

a one-

~
0~~~~ ~~
04w1mi u-aI
m ~ .0

mI- mfaw- .4l molm -,4fmoi i

.. 4 -44- , MA 4p -W v a- Mio

4ma0"a
00momr-W~a~ie *aat 0i
.- D,-mne- a mitw -0&- rA
Dct-M o
wo 0wn.-
am0 ue-1da.-
m in

w v n.Al .01Im0 a. i a ~ 1 ~ l oam pNd..Inw00

. ~ ~ f-- MMW D~~~iV

I-DQ 0 24D

0 I'CO

ovuOmewMfVI Lommulimmo0m Z
2-

Si .0 IV

NWM- ~ ~ -lm- AS.4DNO

bl. .40 4

IL .q 0 f

aM - f - C-4 Inv
CY~~~~p
V. M- -NM0f 4r0 f mI
OA
i
of or
aO~w.w~ie~IL
AS~p00I~Q~iaJ
4 an 0- 00 m c ma o ASiA-
aaij a 1 - - m

IL j aS
F i i .N

w a.

- - t~-VA~--
.. - - b.-4
$ 2--

PAC 23 MMDI
 I.-

al
cz.0 O4O 0 4

'm

za s l 1.1 0 0 -Nr - i

4'-

ixN ~
in$AA NomP-( Nr
00

I-

4O i-'m3Pe Uo~temotm

-4 4i44.4 in f

I.-A

1
3w P.1 -I .r F . . A00 yr

PAG 24 WMsI
 's InC U

M.
. IV5
44 4~
ri0N~. 0-
.. d
i-~~~~~~ roe v- m..

NO.4 0-i(A

ab n r4 0'M

70 P -t W 4.
(7 N-o -oA
a ~~10M4n^ CY 47'n0M.
4

41 ~~ ~' ~ ~M InMr
~ 4 u .t p m0' -"
10UNI '0 ,j

05.0x
0- F inj In' 44 In

-4 - - Oj

-0~~L tN f N .0Niw m

In IN I 'A F

fn .,ql
N a CV

..L tJ p-N

14 C

n an -r t
u
'A

L- V%

- --- -- 4 F

PACE 25 ~. EXHBI 41i.

 -70a

CL -

I
..
ft
9 5 4 - .

in-I
X~
=

A.a
fwt

4p

I.. j
0 0 f

j 0 A
'a t- i a
..

.j -
aa
0 a0Q 0*09

3.491 Ub 40 64 cc~
xw0
a - 0
2 a N20 a U ~IlIO@I0
S- 3- 0. i
w -W
0 0 64" 4wf o
b4 at .r I.M I6-
ft. w "m..-z e Al.4 IOO IEIw-0
fto 0- 0-
uz a1 - 4t
x isa 39 a- IL 2 -0
z w4
-800 " a C W 4 -
wI u. Do-9 ina a o VI 400
a 0 6
aa I- 0 0.. -
0 9- P o- 0
*9 C 0 0 0 0 D404 m. 4
a-~ 0 a A9 F.
in WI.- av
qh. 21 x z 0

lo.u LO W
z z u~o~e, "a h -
zuuo~m

ftI~f .~ aMI 26 4
 a~UP 0 f M

0. 4W =0 n

P- ia- a4 MUM44 mmo MO

w-~ ~~ ~ A~ n So4
4OI

*,wi~p1-5140 U m-0M
IV
10 ~
nMVlt0wso&o ~ ~ fto ~ ~ CO 0
M,*~~
", V -4 A "

04o
0 t4 o"J n .f -- lu- .40 4 -A 009E 4

4r v~
o 4S d
I" go CPi o

I..MI0~J~~l~W~.. 0 514D
I VE
46 N fuAAmI OI
w F o QonWP. @S*PaP
51.4M.. 5on O

*pS -0Ih-514 404S

In inIm I

QV ft 0 IVQ ft-Ai
.4O 0I

*-c.9-vI..o~~mfl4qlq400 15 VU

-m C4 0.

MEM51EV1d4
44*G 2I7.
 EIIT 5
DaFIITI(S- 7235-42340

AC1 - Alienation coefficient for station 1

AC2 - Alienation coefficient for station 2
AC3 - Alienation coefficient for station 3
ADJ - Plus sign indicates value smaller than upstream sum by
tandem test
ADJI - Equal sign indicates value adjusted by tandem test
ALCFT(IK) - Alienation coeffiient array
ALOG - Computer library function of natural logarithm
ANLOG - Number of logarltm
ANYRS - Number of years of record
AV(I,K) - Muan logarithm
AVG(I,K) - Average of the generated deviates
AWG(I) - Average monthly flow for a station
AVMN(I) - Average logarithm of flow for minimum 3 consecutive months
AVMX(I) - Average logarithm of flow for maximum 5 consecutive months
B(L) - Beta coefficient
BETA(I,KL) - Beta coefficient for generation equation
BLANK - Blank space
CROUr - Program subroutine to solve simultaneous equations
CSTAC(KXK) - Coefficient by which flows are multiplied before adding in
a combination
DABS - Computer library function of absolute value of double precision
number
DQ (IK) - Increment of flow
DTRMC - Determination coefficient
E - Letter E indicates estlmated value
FAC - Temporary factor
I - Index for calendar month
IA - Indicator in column 1 of first card for each Job
IANAL - Indicator, positive value calls for analysis
IENDF - End of file indicator
IGNRL - Indicator, + 2 calls for computing generalized statistics
and + 1 or + 2 calls for using generalized statistics for
generating flows
IMN(I) - Month sequence number of last month of 3 driest consecutive
months
INH - Calendar month number for first month of water year
IX(I) - Month sequence number of last month of 5 wettest consecutive
months
INDC - Transfer indicator
IP - Month number for precedin month
IPASS - Sequence number of pass (subset of stations)
IPC!Q - Indicator, positive value calls for writing discharges on tape
IPCHS - Indicator, positive value calls for punching statistics
IQ(I) - Fixed-point conversion of flow values

IXIImT 5

UL -
 IQYAP - Tape number for storing flown
IRCON - Indicator, positive value calls for flow reconstitution
ISKz - Positive. value calls for varying flow Increment (DQ) to
mw skew zero.
ISr(x,L) - Sequence moer of upstream station for tandem test
ISTA(IC) - Station nmbter
ISAC(OC,K) - Station nmbter in a combination,
I3I'AN - Temporary station nmbter
ISTAP - Statics sequence niumber for all passes
ISTAT - Yap. number for stor ig statistics
ISTII(L) - Station nmbter of domsre tandem statics.
IITri(J L) -Station number of upstream tandem station
ISI'X(LS Station nmter of independent statics for consistenices test
iarr(L) -Station number of dependent statics for commistences test
ITEMP -Temporary variable
IYMP -Temporary variable
ITNPP -Temporary variable
ITP - Temporary variable
ITMN - Transfer indicator
IX - Temporary variation of I
IlK - Argument for random number function
IYR - Number of current year
IYRA - First year of data
IYRPJ - Year of start of flow projection
J - Index for year
JA - Sequence number of projection year
JTT4P(L) - Matrix column number
JTP - Matrix column number
.11 - Temporary variation of J
J)XI'P - Temporary variation of J
K - Index for station
KM - Dimension limit for number of consecutive months
KPASS - Dimension limit for numiber of passes
KSA - Dimension limit for total number of stations
KSTAC(KX,K) - Index number of station In a combination
KDI'AP - Dimension limit for total number of sations
1(1 - Temporary variation of K or combination sequence
KYR - Dimension limit for number of consecutive years
L - Index for related station
LA - Temporary variation of L
W1'AP - 1.imber of records up to present position on tape IQiTAP
LSlA! - lumber of records up to present position on tap. IODAT
LTJ4P(L) - Matrix row number
LYP - Matrix row number
LTRA - Letter A
LX - Temporary variation of L
LYRPJ - last year of each projection
M - Serial number of month
MA - Sequence number of month of projected flow

2
uMWIiT5
 MO(I) - Calendar month number
MPASS - Temporary counter for number of passes
MTHPJ - Calendar month of start of each projection
MXRCS - Number of years in each period for which maximum and minimum
recorded and reconstituted flows are desired
N - Serial number of period of flows
NC - Counter to prevent continuous looping
NCA - Counter to prevent continuous looping
NCAB(I,K,L) - Number of values and cross products used to compute correlation
coefficients
NCB - Transfer indicator
NCOMB - Number of combinations of stations max. and min. quantities
are to be computed
NCSTY - Number of consistency tests
NINDP - Number of independent variables in regression study
NJ - Number of years in computation sequence
NLOG(I,K) - Number of logarithms used to compute frequency statistics
NMNMX - Number of months following dry season and preceding wet season
hnMfXN - Number of months following wet season and preceding dry season
NPASS - Total number of passes in job
NPROJ - Number of projections of future flows from present conditions
NQ - Counter for number of flows
NQTAP - Total number of records saved on tape IQTAP
NSMX(L) - Number of upstream stations in tandem test
NSTA - Number of stations in analysis
NSTAA - NSTA + 1
NSrAC(KX) - Number of stations in a combination
NSTAT - Total number of records saved on tape ISrAT
NSTAX - NSTA + NSTA
NSTNP(I) - Number of stations in a particular pass
NSTX - Number of stations in current pass that occurred in preceding
passes
NSTXX - NSTX + 1
NSUM(K) - Number of stations upstream from a station for tandem test
NTNDM - Number of tandem tests
NVAR - Total number of variable in regression study
NYMXG - Number of years of generated flows in each period for which
maximum and minimum flows are desired
NYRG - Total number of years of generated flows
NYRS - Number of years of recorded flows
Q(M,K) - Monthly flow
QM(I) - Monthly flow
QYIN(I,K) - Minimum flow
QPREV(I) - Flow for previous month
QR(M,K) - Identification symbol
QSTAP(I) - Temporary storage of QPREV
R(K,L) - Correlation coefficient in a given matrix
RA(IK,L) - Correlation coefficient
RAV(K,L) - Average correlation coefficient for 12 calendar months

3
EXHIBIT 5

..
....
..
RMAX - Maximum consistent correlation coefficient
RMIN - Minimum consistent correlation coefficient
RNGEN(IXX) - Program random number function
Ri - Correlation coefficient being tested
R2 - Correlation coefficient being tested
R3 - Correlation coefficient being tested
SD(I,K) - Standard devaition of logarithms for calendar month
SDAV(K) - Average standard deviation for 12 consecutive months
SDV(I,K) - Standard deviation of the generated deviates
SKEW(I,K) - Skew coefficient of logarithms for calendar month
SMQ(J,K) - Maximum or minimum flow for month or duration
SQA(I,K,L) - Sum of squares of first variable
SQB(I,K,L) - Sum of squares of second variable
SUm - Average correlation coefficient of matrix
SUMA(I,K,L) - Sum of first variable
SUMB(I,K,L) - Sum of second variable
T - Large positive constant
TEMP - Temporary variable
TMP - Temporary variable
TMPA - Temporary variable
TMPB - Temporary variable
TMPP - Temporary variable
TP - Temporary variable
X(I) - Value of independent variable in regression equation
XINCR(I) - Iteration value for flow increment
XPAB(I,K,L) - Sum of cross products of first and second variables

EXIBIT 5
C 7F3*Vb-L?340 O WTvWLv 'yqEAPvL', SIIAT0 J4C- of ff (11# Rev 1970
CA a* **.* * L~THuAN, FIJ4CTICKS At OG# OA4!1 * , * ft * * * t * * ft * 1002
C PpOGRA4StrMIJi C4CUT,RGFk we &ff CtIMENTS IN PINGIN 1003
C INDEIXES IuChLENDAQ 'IONH JNVAR Xa3TA LoNELAYED $?A ,,AV"C0331*1 1O0"T" 1004
C 0pm3n 1006
*Otll),i(10.17.1007000

pZS7A(ok DT~mCN"P9.L8 0)I (?)0AB1l1?a

*NLFORM(.1l) 1(120. ) ,c1020000I. 0 O~tV~O ~N1O~l~

3C OPHTM100).PA12.a foa
no rn;"T2.1O)3,sNQC30.OA) 102212,spce1.o)SML2.~8
90 IPRNALIX.1O1?0)?E0)hC1),at(2 10113)
60NORMT1([email protected],'S ) 1024C)II(.OI
ToIIL PRECIGION
(tir,263 1026
10 Fnai'ATCIH1) 2F,3 10219
40 FnpMAT (11.71.1916o) f0200l
IOC FPoPNA(IX#,Ib.I2!p 0 1028
110 FORMAT (x.A3qA4pl0A4) 1021
120 FMiRMAT (1X.1J.4,p26) I0 300
130 FORMT034GNAIE TATISx.FT.O.918.)$TO RAV 1031
160 FflUMAT736 .1 S(1WVP Ath MANO0.IH 102
q0FIAT (N.7128614 1034
I(C
OWAP29 .0.61! 1 t)20

KYRA8*11 1039
KIn2YAPU9* 1040

KSTAZIO 1016

C WASTE CARDS UNTIL AN' A IN COLUPN It FIN$T TITLE CARO 1043

C to CARO A to 1044
ISO REAO(5,110) IA#, (M#0#41)MUI.20) 1
If (!A.NE.LTRA) GO' TO 110t94
wRtTE(6qo 10)4
~EA01.0) (S~4tM.KP3.20.Ku,3)1650
t4UIE(6i0) (Sea(N,).MI.?O.KU.3)1051
Ct CARo a CAAO C to 1052000
IiEAnc(.3O.) YVtRAIMNTWJANAL, AQCS,1KYRG.NVfNG#KPASSPCNQ. IPCNSNSTA t103
I NCfMI;,NT Nil M,"C.1T, t~l~t,#1.OPRflJ, TI rPJpm~*PjtLyRPJ 1091
C VF.RmZKATE WITH S PLANk CAR03, AN. A IN CfL I CF FIRST 1056
ITONIAMAL,'.7q0 1057
IF(?TM,'.GT.O)GG O ISO
160
STOP 1059
140 WOTE (6,*1701 hYNSN3TAptKCflV4N.PAS3 1060
ITO FOMAT (119i4 rIIMEWNO ENCECEt #5%g4WNYNSI4*1X$4"mNTA9I3pSX. 1061000
iSMKICftMA. I3.5xr.MtP44S.113) 1042000
Sof TO ISO 1063
sac WRITE(b.19l) 1064
19C FOIWATC/IO8N typA M"KTN.JAAL 'INaCi 1IRG hyl$XS NPASS IPCI-O IPCHS 1065
I KSTA mCapfia NVNO04 NCS7V YOMSL. PROJ ZYNI'.10iPJ LYRPJ 1 1066

IF (LYIPJ.IYRA.GIE.KYMI 60 TO' 160 1069

cq * * * SET CONSTANTS 1070
tws**1071
t4SIAA6'.ITA01 107?
NSTAI"T'1A4005A 1073
Tmeqqqqqq 0 1075
!Y3*sIt'Y 6
 p457 ~a01077
IPASSI 10ll
REWIND IOTA? 1060
9ST4T&O 1061

R4liIIO TOTAP 10,83

NOTARUO 1084
LOTAPRO loss
Do 195 J8l*EPASS
NTNON(J)80
an To 170 106
C SAVE 3TATIONS FROM PREVIOUS PASSES IF MEC93S*O 1087
100 IPA3IPA&3+I, l086
~4ITEC8.10) tool
If (IPA$S.61.KPASS) s0 To 160 1090
c . CARD J *4 10l1
WRITCE.i!2e) IASSIST(XhWIINTX)1093
ate PFORnAT5SHOPASS #13/2H8M $?AM FROM PREVIOUS PAWSS v1016) 1096
NSTYXXNSTN.I 1095
REWIND IqYAP 1098
LOTAPSO 1097
REWIND 18747 1098
READ CISTAT) 1100
LRTAT21 1101
I7EmpxsjSy"P(HPh8S) 1103
ITMPP20 1104
onO 5 2~1,RSyx 1105
220 RFnIOVAP)TTNP,(0CNK).Pu1,ITP) 1108
L4YAPLO!AP+1 1107
IF .fT4TA(K).NE.IT"~P1 go TO 220 W
230 IWPPNTPPP.1 1109
IF0"Pt-T1E:P 7"n
To240 1110

LSA3LST?*11113
IF (!11P.EQ.!STA(X)) 60 to 250 111*
Gn Tol 230 1115
240 REA"CTSTAT) 1116
LATATULSTAT#1 1117
HPASS84PA38+1 1118
tTEMP32J8TNPp(PA33) M19
JV.PRO 1120
Go Tr 23fl 111
250 CflPKTtk'iE 1122
On ?60 XmI.NsTN 1113
On P60 ?vist2 11es
260 fLOq(I.E)uNVRS 1128
270 IF(!AN'AL.6T90) WSTAWNSTI 1127
0n 24f0 TwIP12 11e6
4o0(1)3J1PNTH* 1119
IF(PO(I.LT.t3)G#0 TO 100 1130
HM ane)~?(!)*12 1131
260 CONTINUE 1138
I7(NCPMR.Lf.0) go1TO 320 1133
C IOE4.TIFY STATION COPRZNATIO4S 1134
On 300 NXlOCOMA 1135
14 CARD 0 ** 1138
R!*D5.3)ITP CItAC(.L)LSI.TPI1137
WmITE (60290) .TIAfNLLUZP)13
290 FORMAT USM CoNBalass" ITA.ISZS) 1139000

mEo9.@ TEP(SA(NLIASeLIt) CARO 1 * 1141

WNIUIT wits
 V w:JTrE(1-, lt',I STAC (oL. P43)L1ITOy)
S1j' FflQM4T (7V1 eAT!CV.X,11JFM.3) 1144000
JPC TFCP4YTND'.LF.0) GO TO 350 1145
MTN04 CrP1AS)?'TNDN
00 540 Ly1-1.NTNDM 1146
gE00 3) YSTN (LW) .11~'P.(ISTTCLXL1,L.1,1TMPI it&&
WQITEC6.33()) iXaISTN(LIO, (JSTTLLX,L),LSI, yyhpJ 1149
*IIC FnR'AAT (/13H YANO~FM ~L, 3 &,o~w8qAiSAi.x 1150400
315NLJPSTRFA4 STACS),IO15) fissmnn
.340 KS~IiCLXI=TTKP 1152
350 I(TPh8X.F..1)Y40 1153
0n 3S0 K=",i1vWpK3TA 1 154
4sU"(K, TPAsi)MO
ZI1TA (K)1fln0O0I 5156
C 7Ij1TJATE -1, NO RECIRO FOR ALL FLOW6 117
')a 360 Mhl.KM 1156
AfI(MOK)U.-10 1159
am 370 1.1,12 1160

00 IN 11. 1163
37C CONTINLIE 1164
34C Cr'NTIt~IE 1165
YF(NCSTyLE.0) 60 TO 420 1166
~4R!TE(6p. 1 1167
34C FflQ4AT(/30yv8HSTAT~flNs/I7w CnNSJS1'ENCY TE8Tg5X#23M'INnEPENOENT OF 1166
1'Ex 0ENT) 1169
01 400) Lz1,NCSTV 1170
C .. CARO a *0 1171
.9EF4 (5,30) YarXCL)T3TY(L) 1172
4"0 IqlTT(b.JII01 ~ L#TT(L)oI.1TYCL) 1173

VIC IF(IANAL.Lh..")G0 TO 57 1175

CC 000 q AD AND PROCH1i3 I STATION-YEAR OF OATA I*4A a * 1076
C .,CARO N * 1177
qleC QFA0(5.O TTNY.U9,1$) 11783
C 0*CAhO 1 1179
C RLA K CiAD IkOTCA7ES 9wD OF FLOh DATA 118C
1F~tSTAN.L?.1151 TO SrO 1101

C A' qTqN SOUSCRIPT TO STATION fls

On 40 IK:NS1'1,NSTA 1184
IT(f3TAN.Fj.13TACW))GC TO 460 11o5
oiAc CflNTPIJF 116
a!O NSAAN3TA*5 ts
IF(N3TA..K3TAJ GO TO 160 118720
KINST! lice
ZSlTAf)*19VAHI t189
C A;INSUf3SC41PT TO YEAR I140
44"C JaXYqvYP4 I191

ZFJ,~0*JLEPii, GO 71) 480 1193

W101TE t6.'%7n)1yq 1194
(670 FMPMAT (/I'iH UNACCEPTIHLF VEAR,IS) 1195000
'rflTO 150 1196
C S~fl.E FLOWS Ji~ 3TAT1C'N AND MO.;ITM ARRAY 1197
"Jo,11.1 1
Ma 1194
rfI 490) ?21.12 t199
~1 1200
SGO ril. 490 1201

1203
4qC CIJ.TINUF 1206
'i' To) 430 1207
5nC NSTAAuN.A7hOt 1208
IF (NVPI.O7,AY4.DQ. T7A.NCO?,B.GT.fI'A) 90 TO 160 1209
F(PJSTA.LE.fl) GD Ta 160 1210

.3. C~tI18TT 6
 rip(NCOMS.LE.O)GO TO Sgo120
C IDENiTIFY 3TA SLBSCA1PTS FOR STAB IN COMBINATIONS 1a14
nf 530 9XVUI.cNCe lass
TYPEpSTACtKX, IPASS)
LV80 1217
on 520 Lwl,!TP lia
IVfMPS!STAC(K~oL) 1219
no SID KvI.NSTA 1220
IFC!STA(K).NE.ITE4P)GO To 510 1221
LYNLX*1 1228
AITAC1KX.LX. IPASBSK
GO To 520 1aae
SIC CO'KTINUf 1225
S2C CMNTINUE JIM2
C REDUC~E STATIONS To THOSE IDENTIFIABLE 1227
N4YhCCKV, IPA3S)2LX
530 CON~TINUE 1229
C IDENTIFY STA?IOCS IN TANDEM 1230
54c rtNTNOM.LE.0) GO TO 600 1231
Pr)1540 LX~1.*iTNOM- 1232
On 550 XuINSTA 1233
tFtISTACK).EIJ.I3TN(LX)) Go To 560 1234
9;!C Cnh?rNIJE 1235
SAC TSINCLX)xK 1236
NAUM(K. yPASs)=N3"XCLX)
ITvaNSMxtL.%) 1230
On 960 t*u14rMP 1239
Do 97e) I(W:jq3TA 1240
IF(15TACkW).Fa.ISTTfL~,L)) GO TO 38o 1241
570 CO~NTINUE 1242
SAC IATfK.LsIPASs)zKX
59C CnN!TINUE 1244
C 10E~irIFY PAIRS CF STATIONS FOR CON31STENCY TESTS 1245
60C IP(NCSTY.LE.0) GO TO 630 1246
.jr 040 tL,Nc3TY 1247
OM 630 Kxt,43TA I248
1VC7STACK).Eu.ISTXCL)) c-c TO &10 tagq
tFCI8TA(K).En*ISTY(L)) GO TO 620 1250
iml TO 630 12s1
610 !ITX(L)uRc t252
GM O7 30 1253
62C ZATY(L)XK 12S4
630 CiINT!NUE 1255
bac cf"xTINuE 1256
690 ITMp3NSTA~jCtPM 1257
CD A A 0 14AXCAND MIN RECOiNDEn VOLUMES 0 tOt*0f t0f t001258
C TNITIATc SUMS 1259
LOfl 79n KXNSTXXITmp 1200
AvfOf)KXmO* 1261

)m, A60 1010ts 1263

60)C 1 M2(I . X)*-T 1264
.)I 676 T ul6 p3 O 1265
6CSmntIK)T 1266
0
7T4Pui~ t267
42t 1269
piz0 3270
On ?fn JU.p~jv7S 1271
pm 770 721012 127
I1273
12PI +1
K:fV 1 1274
IF(K.LC.NSTA3GG TO 700 3275
C CompuitE COm8INED I;LO'd8 1276
-(YSXNSTA I27Y
ITPU=4STAC(KYD IPASS)
l (P, K)6. 1279
P't 490 Lnt#ITP 1200
ITf4Pv9STACfKWe, IPA3I)
C CU14111PEU FLOb PIASING 1ama
~aIEP~O*1~6.f"K)E~u. GO TO 600 12M

EXH4IIT 6 4
 Gn TO 690o 1285
0C 12.%u..
as
bqc CANYThue 124?
?0c I~lMK.E- GO TO 710 1248
C STAwiT NEW ACCUMULATIONS WHEN F40W MISSING 1369
oJso 1290
TmpmO. 1211
YMPASO, 129Z
GO To 770 1293
?10 TEMPX0(fl.K) 1294
C I-MONTH FLOWS Iq

tVCSMOL28,K) .rT.TEMP)SMC(28,K)uTCMP 1299

C b.MflNTM FLOWS 1300
FMP&UTMPA+TFMP 1302
1 Cw.6)7ehO.730,72O 1303
72c TPvTvP.)M-6,K) 1304
?3c 1PCTmP.LT,.4n(29,K))8Sflt29,K)QTMP 1305
1I lIMP .GT.1.,( IU K) 3SM)C14pK)aTWP 130b
C 5Lj-inrNTH FLOkS 1307
IF (W.44) 760.750,740 1305
74C T !PkSTMPAfl(M-54,K) 1309
15C IF CTMPA.LT.Sr'inC3O,PK) ).1t0C31,X)zTPA 1310
IFCMPA.OT.0 ,0(15,K))8MQ((15,K)uTMPA 1311
c AIEFAGE FLOW 1312
76C AVGQC)(K)2VGQCX).TEMP t313
Nnskl~l1311
770 CONTINUE ISIS
nrINUE
r~0 1316
TFlMhUNO 1311
AV~tl K2A VG (k) /TED4P 1318
7 CnNYz'.uE 1319
XTE(6 800)1320
00OC vriguAT(,/34if MAXI?41'I VrLLP6 (IFRECORDED FLOW) 03?1
*'PITE (6,510) m(Mfl Cl), 12) 1322
Ble FOnmAT ($H AT&#12I7,3SH S1MO b.MO 54*NO AV 140) 1323000
ITPPz'N3yA*NCflmo 1324
.)r 431) Kz'iITwX,rTHP 13a5
lTPAvrfl(X)*.5 I26
j~p lftt1.,15 132?
82C 'f'P0C ~ . 1328
At~4c (.f0 ~A~) (aI.:,5ae" 1329
S'~C
A HPACW I',11,28I, 1330
i~iT C6,~0)1331
1)SC FmcoATIC/1A P1NIMUM VOLLJE3f) 1332
.4U7TFC 6r8lC) (C CI), I1, 12) 1333
rrfl 470 Xkl~3yX'A,ITMP 1334
Dm Rho 1:1,13 1335
67C .1E6O")5AWTf)I11 133?
CE * * * * * * C014P1TE F'NECJLEKCY OTATXSTICS 0 *0 a 1338
icI1E C,pit) 1339
APC Fll&;4&T (1214W FA0UjF'JCV STATISTIC3) 1360
prrE6,R0)c~(I)T~l,2) 3aj
hqC Fnr-mAT (/t4; 37* ITEM',!7pZ1A) 1342000
c "rS5qfG FL.OW POEM!NG FIRST RECORO MO?4T14 1343
,In 100 MaPiNTVV,N3T4 1304
9(c :J(1.KX)UT 1343
!UCOWUO 1346
17fhIP 2 NITA 1306100
17 (jEHP.GT.N$TA) GO TC i1.0 1347100
IF(w.Lf.N3Tv) GOl TO 942 13'7ma
91C nn q2 0 Ta1,12 1348
TFMPww4 OG ~1,A() 1349

!0M1W&D(,) *.,) G(,KuD(IIc.QN?. 1312

92C COP.TINUE 1313

 1354
glC An 44M yal12 1359
LV (1 * )sav)* 1356
smc1.K)s30. 1357
SKEW IN)UO. 1358
TmPUN 1359

400 Cn~fINUE 1361

942 al~~ 1362
On 970 j.1,WY4S 1363
On 960 1.1,12 1364
Hmmm I 1365
tF(0CM#K).E0."le) GO TO 950 1366
C REPLACE FLOW ARRAY WITH LOG ARRAY 1361
TFMPAOG(gC(i.)*fl(I))/2.3025I51 1368
r(0,KuEMP 1369
IF(K.LE.NSTX) GO To 960 1369-2
C SUM, SnUARES, Ahn CUBES 1370
Lv (1*K 3uAV CT K) .TEMP 1371
$fA(r,K)2sfl(I,K)+TEiP*TEPP 1372
SxEA(rN)3SXEW(1,K)TEMP*TE4P*FMP 1373
GmTO 960 1374
C MITSSING FLOWS EQ2UATED TO T 1375
95C ji(m.K)=T 1376

9Lc CnklTINUE 1378

970 C"n.YINUE 1379
!~W.E. ~Gfl TO 1180 1379o2
!wTW~nO1380
n~llono 121,12 1361
T7PL0G(DK) 1302
IFCEi4P.LT.3.)Gl TO 1120 13R3
T,,P:Av CT.'C) 1384
LV CTK)mTfl~TEmP 1385
lF(5PlCIpK).LE.O.)G0 TO 980 1386
T~w~SO C. ()1387

TF(,AflCI.,).LtE.0.) GO TO '3P0 1388.2*

S')
51341 I, j:5) C ~) *
tFQr(,K.T.GOl) 2 Tr 990 1390
S'F .C! ,K):(TF *TIP*1(EW (y,K)-3.TEPT PTMPA*2.,T0*TPTMP) 1391
1/(TrmPCTE1P*P1.)e(TEIF-t2.)a9n(TDK)**3) 1392
Irf.;iFw I, r) .(.T. f-.1) . [IQ SXEW( I,K) .GT..1) INDC21 1393
IF ClX1E.4T,PK) *I7*3.) SxEh (TK)s3*

G1Tr) 1000) 1394

91C inC?,K)S0. 1393
qqc 'qEw(7.w)UO. 13%6
iomc CflPjo 1397
.41 1398
1:(N.rT.llGO TC 1060 tl99
..
rTFC6,11)ISAtKdCAVC?.K),I31,123 1400
toic Fn;."AT CI'.t5,8i4 mcL&F.12FO.3) t40100O
toC FlkmAT (?V.?,iqTn rEv,12Fui.31 1403000

to C rjr-4%T (0,uSC 1Pf.)1405000

1
iA'C "IU-47 ($l%,h~IoCuT.F7.2,ltF8.2) 1407000
1f,' p61C~C~,Q!11)10
10AC lr(N.AT (q.514 1, TO 18 14100
JF(?n0C.LE.nj f;- 7 1180 1411
C T'E FLLl;1dNq ACUjTPkE WILL ADJL3T THE INCREM4ENT TC 1412
C TRY Tn CdTATIN ZERC SKEW 1413
C CHAI%C~E THE FCLLC'wTmW. 3747 T9 I8'CZs To ACTIVATE 1414

IF(TAKZ.LE.0) GO TO 1180 1416

lyps-11 1411
,n ti10 t31.1 1411
1421 ?* I 1419
on 1080 JaloonYs 1420

EVH6
T
6*6
 YV~f(f),I().1O#T) Go To 1070 1422
T",Fxrl (pip ) 1423
J(M.K)UlO."*TMP nfo(I*K) t424
Goi Tn 10ao 1425
107c AR(,M.K . %426
JOAO CM'yr~je 14P7
T~Mf2SKw (IiK)1426
XVTMPfT.~1)AJ~1MPL.,)00 TO 1110 1429
1VC71NP) iaq0#11011100 130
100t I(,~)0GIK)2 1431
GM TA 1110 1432
I IOC 0A(T,X)XUOC,K)*XTNCRt(I) 1433
111e Cro1'tNut 1434
Gn TAi 930 1435
C * * * 4 * tnELlFTE 3T&T1GNS WtT14 LESS THANi I YEARS OF DATA'*a * 143b
I1120 .,LJT6(6.ll3O)TSTA(K) t'137
113C Fnr'~A7 C/44 374,1b.281M DVLETED, INSUFFICIENT DATA) 14138000
~. '~T~t .~ A.!t439
N9T4ASSAI 1440
NR~TA ~ ,N
S7+N 3 TA 1441
IF(K.GT.NSTA)GC TA 1180 1442
C QEfluCE 311~8SCRIPT9 OF SU03SEOUENT 3?ATTON8 1443
On 1170 KXSK.N3yA 1444
1.sTA(Kx)mISTACXWa1) 1445
imj 144b
an 1130 JV1,NYHS 1414
An 1100 1M1 1 12 1446

r~i ~6A TjV214159

ll CPTN)UE
p~ NCC~1 1453
'.1 G (tK'A)NLAG(IK*1)1452

WIC c'".rpUE 1456

nnfl 910 1457
CntNiIF
110Ccn 1436
TFCTsjCP?.LE.d)) GO TO 1370 14160
CF* * * * * A~jtjqlmENT Of FRECUENCY STATISTICS TO LIP4G TE41N 11 * *a 0

Or 1 191) KaltiSTA 1463

oA 1190 Lml,NSTAX t6
NC44N(1 * K,L) cA 1465
3tk&!1.KL)MO. 1416b
S1.I"IC!,I(,)0O* 1467
SrA tT.K.L0. 1460
S1,;rT#9,L)V0. (469
wP An~(I,, L.) =A 1470
af!(*,K *1)Z . 171
t
11qc rmkTItjE 1472
On 12?0 Ku1,NSTA 1473
,(Vcx~l1474
'11 1 1473
nn 12PQ jaiol'iS 1476
D00 1210 141012 14177
TI MP05 (P* ) 1479
1FrFrP.Er).T) Gil to 1210 14810
(3n 1200) LUM)(,NTAW 1461
LVGfL-%ST4 1402
IFRXL.T.1) TMPUIJ~tfL) 14aj
JF (1.%pr.T.O)
TMI'2A1(fit eLy) t464
IFfT$40.E0.T) GM TA 1200 1485

St)4CI.4L)~ll.JqCTNeL)Tpp1486

IF(L.GT.NiTA) GO TO taoo 14qa

ova EVMTIYST 6
 NCA4(1,LK)I2C&A3(!,K,L) 1493
Silt# ,L)1494
fi(I L tK I 30MR (t
.Sr'AI(I# ,1*K 313 AI #)(L) 1495
SI AF(I#* .,K) 23) K i1 K L) 149
l2~C Cfl~I~uE1498
1210 CnKTINUE 1499
I220 CnNTpIJUE
1500
t2CcIN'C 1501
[In 2K4 K*IvN3TA 1502
Dm 1260 Ixt#12 1504
4A(I,K.K)zj. 1905
,n 1250 LzKK.NSTAN I506
tFCNCA03UK#L).LE.2) GO TO 1230 1507
TF:PZCAB (I .KPL.. 09
T74ns)q I q '1) 1510
T'4PACqUHA(r.K.L) 1512
I (T'PP-LT.O.) TlIPP0O. 1514
ma (TK.L)T~lP...51515
T-Pf'z TP.T0P4i.*/TEmP) /TEMP 1516
kF (TMPPp.LT.O ) TmPP3Q, 1517
51C!,KL)Tm4PP**.5 Is19
(FeT'4P.LE.(i.) G;1 TO 123M
1511
!0 %32WP1nf(I.L)TPAT)IPf/TFMP 1522
r 1523
tF('4ALT3)Tmpau.Tmps
1520
T'fP.':y 'PA AT 'PA / 1 -P
Isas
IfP(TmPA .1T.0.) T44PA:0* 1527
.A(f,K,L)TH4PI*T~lPA**.5
12
L A =I.

tf (L.ILE.N9TA) GO TOJ1233

fr(' P.Lr.1) -ATP212

t. -A-
1235 rF(Io(I.K).LT..0001.OR.SOCZTP,L&).LT..0001) GO TO 1230
".T 1240
1529

1530

1534
1?M. I'.TTNUF
1535
1~AK=1,0TA
ITh 153?
1539
L 1540
(L.K.ACKL~.~J)GO TO 1270 154a1
1502
'k~'GfT*L).L..PLUG(IX)) Gr Tn 1270 1543
T'IrI2t.CAPC?,K,LJ 1544
"'r'=LPGTL)1545
fV(TP.Lt.TE'4P) GO TO 127?0 1547
Lva t. 1344
re~varp159
127'C CM' T!4UP 1551
I )I.LE.0)) rO To 1280 19
f'oiCri
T" 0!AA(KL)/SnR(3I,K,LX) 15,13
1550
r"'AgSII'4A C? K#LW) TMPF
1355

Evp~T8Trr 6 .8.
 T"FuSup4Iq CIKpLX) /T'4P0 15
As ,K.):TMPA4+(AV(TLX).TINPI.1aQA(1,KLx)aymp M
3r CTo K)8OA (I vK LX) +(SO(1tLXf) 3(l1S (I#K.L0 )*AI#CK,L U *a2*1iP 1598
C LOJUST SYAND&II OEVTATON3 FOR COJNSISTENCI 1560
C (NC.E0 TO; 340 5

KCISTY CLX) t564

LUXISTY cLir) 1565
00 1321 t1112
1566
TFHP=PUAVCT,9).AYCJ,L))'3. 1567
IFCAVCIX).rTAYCX,L)) GC TO 1300
1568
TFPzTEM0+SflCIN)
1F(1,0C1.L).LTflmP) Gt, TO 1310 1569
1570
IFCInCXDL) - TEmP) t120,1320.1310 157e
1300 TFmpzyEmp~gnCI,1C)
CF(5n(IL).GT.TEmP) GO TO 1310 1573
1Fv~a:5nC,K)a2. - TEMiP 1574
1575
IF(qfftL.rE*TEPP) GO 701320 1576
1310 3'C,L)RTEMP
t3PC C,.,TINUE 1577
133c C0m.INuF
157a
JF1FTR4~S.f5y.) 1579
GM TO 28e0 1sma
134C 1F(TkOC.LF.n.AKfl.NC3T'Y.LE.O) GO TO 137 1501
-1)1TF (b. 1350) 1542
t3SC rnrATC/39H FRE9iiENCY STATITIC3 AFTER AOJU3TMENTS
I1583
(,n 1360 Xm1IoNSTA 1585
dIQIl'EU,.1OIO)ISTACK),CAVCTK),Zu1l,12) t586

1589
136C COPKTINuE
jgqO
CG * * * * * TOsAIS~FlP To. STANOI)DIZE0 WARIATES we* a aa aa
13 7 C '1' 1420 Kvl1iSTA
'Is1 1592
1593
On 11410 Jxu14YR3 59
.0~m+It~ 3,12
OR(Hx~zSANK1596
1F(n(NK).EQ).T)GO TO 1400 19
IFUqnTK).Erj.O.)Gtj TO 1390 1600
CEAANTYPE III TRANSPORM lO
IFCSkEW(.K7.FG.O.)G0 TO tiOO 60
TF P.5*KEl(I#K)JaOCMsK)+1
*yvx
1604
160S
IF(TENP.E.n.2GO TO 1380 1606
TV 'PM-TEP
1607
* 1604
1 3PC
1609
G1i TO 14lO0
1610
I.ICr~.Tu
11cC CfII.TTNF 1611
IGICCnKYNUE1613 1IE

C It a ** * Ci'~tIT SUMS3 OF srtJANES AND CR03S PRI1OICIS a ai a a11

00' l4Ob kU1.-STA
1616
rn 1030 L3t.NSTAX 1618
.IHAC!K,L):o
1620
3"'A low* )ZO,1620
in~~f(TKPL)M0*
1622
10 '( T * K.1)0.
1624
IU4C (.Ko)mO,
162!
4&(Toxx~sl,1626
145C CO'.?IvuE
1628

-9- EvHRIT 6
 )Iri 1540 KZI,NSY4 1629
KVZK~t,630
,-*1 1631
nfl '460 J:1.NYHS 163a
n 1470 IC10i2 1633
.
va 1 1634
yrmpacnM#K) 1635
IF(TENPsE0.T)GC TO 1470 163b
DO 14160 LzWXNSTAX 1637
c SUBSC$41PT3 EXCEED1ING p.574 RELATE TO PRECEDINGl MONTH 1638
LY*L.NSTA 1639
IF(LXLT.1) TmPxQC'4pLl 1640

IFCTHP.En.T)Gl) To 146i! 1612

c Cr~i'iT A"O JiI5 ChLY REC04floD PAIRS I64
(I XKt.) auvAI (1,Kpt) 41
l~A Ib44
IIA( I PKA )300 A(IsK 1)4TD'P teas
~ ' p1646
31 ,I ,I)4 fi(7,K, L 'i I~4 (T
7 1647
3rJAu*i: A (T,K,L)*YEIPPTENP

IF(L.rIT.N5TA) GO TO Id4bC 1610

JJ (TK,L,K)=J.14Aj (1,KaL) 1654

SmA (K.L,K):Sjn (r,kPL) 1654
wq(I .L.K) cxPAlCIe(Kl
,LJ165
1460 CnTINUF Wy~
14~7C rmNTTYNE tb58
141AC CflA.TTNIUE 1659
C . * * a CuMPUTE CUH4EL4TIAlN COEFFICIENT3 a * * aa1660~
On 1530 T= tI I 16bl
'10 1520 LZKX,NSTAX 16012
LV=LNsTA 1663
C EL14INATE PAIRS OITP' LE3S TpiAN 3 yR3 UATA 1664
IvCvCAB(!,A,L).LE.2) GO TO 1510 1665
TFNPZNCAPT.I(.L) 16405
TMP~(.ljA (,K,).5JMACK.~a5UACIKLrEM)*C(IKL)SJ t 67
pK ~L)
*SU43 I PK L) /TE~lP) 1668i
C ELTrITNATE PA1;45 WITH ZEkO VARIAN~CE PRODUCT 1669)
!F(TF'P.LE.O.) CS1 TO 1500 1670
TmP01231 * 167
T"IPAXXPAf3C!,K ,L) .5IJ~'A (1,K.LIaSUI4ICIKL)/TEMP 1672
C '4FTA1N ALG~t1qA1C SIGN 1673
7F(T 'PA.LT.0.)Tmr~z.TmPP 1674
=TMP A* T45, A/ Tf!~P
T'-P A: 1675
:;'(VI4,LET!)*TMc'Aaa.5 1676
yp~i! 1677
LAL 1676
IL C *.i.TA) rfl TO 149 1680
1 1-1
IFCtTP'.LT.1) ITPc1a 1641
LAI X 16ma

( ?fl 1910 1684

i S 0C4 A( T KL) a(. 1665
SilC:rfI..0T.K'TA) GO TO le 16fl
I~
*1.,K) UR4A(Tv ,KL) 1697
15?C ClIPYINIIE 16p8
I 53C C? TIN iJZ 189
154C nKTINIJJ Ith90
inr TMl 2170 1691

156C FfRMATC/ISN DATA OUT CF CgROEq) 1693

r,0, TO 190 1694
C;4 * a * - o at PE40 C0U'Rt'LATir% COEFICIENTS 1695
1SCOm 163o KUIMSTA 1696
19 (i(.EG,1)Gf TO 1800 1697
ITONK1t 1690
On 1590 L21-1TP 1699
C CuRRI4CT IMfNTh CCOPELAYXON 1700

EV~rUI?*t.1
C A* CARD L *01701
;AO(5.7o)iT'4PITE4P. (R&(TKL).zul12) 1701
*" Vf(iL) 378(1 *K. LI 1703
IFCTCNWL.FQ.I11TAY(K)ITVP 1704
IF(TTMP.NEI~sTA(K)IO TV! 1550 1705
Ic(1TEMP.?JE.ISTA(L))GC TO 1350 170*
C-i 15AO 71112 1707
1380 flA (T.L K)NRA (!,K,L) 170a
1590 CflhINUE 1709
C PRECEDOING MONTH CORRELATION 1710
1600 LV:'k-S7AA 1711
IF fIGNPL.EO.1) LXQNSTA*K 1712
LAMNSTX1713
IF UIGNRL.EJ.I) LAgLX 1714
rr1 1610 LGLX,LA 1715
I!P21.-ISTA 171b
c a CARD K( OR M A* 1717
M79
IF(k.EQ*1fl8TA(ITP)=ITE0P 1720
IF (IGPJRL.Efg.l) PAYCKpX)9QA Cl ,'L) IM2
XF(ITMP.NE.IA(K))f',O TI! 1550 1?72
JFI9PN.-9AIP); TO' 1550 1723
1610 CrNITTNUE 1724
011 16?0 71.12 1725

103C CrM\TINUE 1727

C a * * * READ FREQUENCY aTATTSTICq A Aaaaaa aaaaaaaa 1728
01h 1640 XK:1,NiTA 1729
C ARO84 N C4 0 A* 1730

rP(TTP.NF.TSqTA(K))Cfl TO 1550 1732

C GENF'WALIZEO STATISTICS ON ON~E CARO PER 3TATInN 1733
AV44V K)ZAV (1.K) 1734
A Vtlk~(K)Z3A V(2, K) 1735
S~~Av(K:AVCK)1737

1739.2
IF (IXK) L.)i1X(X)m:PYC),12 1740
IF (1M')tc) .LT.X) I~iw1(K):IpN(K)+%2 1741

C aaCARD P 1a
743

?F(TTP.PJE.XSTA(K))t;(J TO 1550 1745

C aaCARD 0 aa1746
qEAf(,Sf)Tp(.qF~K(.).1l.12) 1747
rF(TTP.NE.I'iTA(K)).o Yf 1550 1748
C aaCAkO 4 1749

1F(TTP.NE.15TA(K~))flA TO 1550 1751

thao C?, TINUE 1752
Cy EFITIMATFJ M1351KG V.'WIELATION CtI!FFICIENTS 1733
1650 IFfI'1NPL.F,, 1lG(I Tf 3'2e 17,54
tF(NT4.LF.IGO. TO 2311) t755
00),1720 rm1,12 1756
zI:-'I. 173?
TF(TM.LT.I)101I2 Il58
01 1710 K1.&N9TA 1799
ITP~l(+l1760
It' It)O LUTS 1761
C L A~n K ClOR4ELATII!N POSSIBLY m17ds7NG 1768
!F4A!K,).F.(G) O TO 1700 1763
~i'*Au1 *1764
4M1N.1 *1765
C iw SEARCHES ALL QELATEO CORRELATIONS EXCEPT FOLLChING mTm 176*
Or 1640 LXN1.NSTAV 1767
Irt ~ 1b90
oTO.) 1768
IF(L.FO.LV)Gn 70I 1690 1769
TFNPPACIXPLX) 1770
IF(L.LE.NATA)GC TO IbbO 1771

.11, EV),rsyy *
 IF(LK.LE.NSTA)Gn1 TO 1670 1772
c SOT10 L AND LX P~rSethr PIECINGI NOW 11
T11E4PNLXftN8TA 1775
Tt'P*RA TP#EIV1PIT(Npq) 177
Gn TO 1680 1778
c L REPRESEPITS CLRMFNT 14O*4YT I778

Gil T0 1680 t170

c LX AND NOT L NEPOFIENT8 COJRRFNt MONTH 1782
1610 T&4P~q4(1.LX#0 7)
t6AC IF (Tm0.TEMP.L7..2.0)Gc rr b 1783
T4pAU(..TE.P*TEmp)*(..TP~Pe7,P)) F8
IF (TM0A-LT.0.) T4PAm* %IRS
Tb.PAvy"PA*.5 174b
TbPP TN P.TEN P.7 $PA 1787
1'rCTfiPR.Lr.RiiAx) mAxmThPB 188
Tep4%TmPRTM*A.TfMPA. 1791
IF T490TRisRm~~~ 1190
169C Cfl\IINUE 1791
c AVER~AGE 11HALLEST P~AY AND LARGE31 MEN C0NS1STENT VALUE 1792

IF CL.LE.NSTA)PACJLJ()MOAC!,KL) 1794
tVIC CribTINuE t793
tIC CflNTtN,1E 1796
172C C'ip.rTNUE 1757
G1 Tn 2311~ 1798
Cj
* a TEST Fr~f1 TRIAb CC?4JS7TENCY a * * S* 1719
173C NC300 1800
174C FAC=I. 1801
~ ~ ~411802
I-;UNCA.LTNATA*1Z) Go Te 1790 1803

In TO LS0 1805
I75C NCIX 180b
t 7 4 C 1"no1807
1b0
07E' 111 18019
tPIT-1 1810

C Ko Li AND 0X SEARChb ALL 9ELATfn TRIUS OF COMPEL CCEFS l812

1F~Q
!A2o mioNTA, Ms1
IT "PD M.1 1814
m ISIO LZDTPONSTAX 1813
V (L.Et2.-4STA)G TO 1611D 1816
LA41.LNSIA tat?
,4I ib.A(IoK.) 1il18
I~wcL~l A19
nr 1400 LY.2!TP,?3TAX 1820
I 1~ 2LXN9 TA1821

(L.L.H5T)RS~(ILLW11823
C 5lT~4 L AND LV RFPPESEN? PRECEOINto MOMN 1424
t(L.0,T.NSTA)q3:RaCIP.LAeIT~t4P) 1825
C RAI~k Lni'VST CCEF21C!ET47 IF tNC%)N313TENT 1826

Le t(1 ..- aR* 1829

IF(41.GT.02) GC TO t770 1830

I-.1 PZ*p3.ACP.ACJ*FAC 1632

IVfb~4b.L.1*I NZ-1 1633
I~W.GMP1Gn Tn 1800 1834
I' nrt 11133
P~ t~EL)QNIN1836
IF (I.LE.NSTA) RA(1,L.K)IURPIN 1831
Gn TO 1800 tills
1770 IFI1?,.GT.03) GIn TO 1780 1839
kmpIbmaj*3wAC1*AC3OFAC 1640

ININS01- 1653

6P
EW,4?PIPIN TO*010as
 ;A CTt,LX)uRti!N 1044
IF' (LI.1-E.NSTA) 181.41O5MP
(;n TO 1600 1646
t7A~ IJ~c~1.2.AC*AC~FAC1847

1FC03.G9.R11J) GO TO 1660 149

tNOr*1 log0
IF (L.GT.NSTA) GO TO 1790 18s1

IF (LY.LE.NSYA) RAClLX#L)DRMINi 1853

(fl To 1p00 1814
1790 PA(PLA.IT6MP)21Rm1 1855
RA (!PITE0PPLA)24mIN 186
I1IPC CnTINUE 18537
1510 CrnhfImue 1451
tA~r COnTINIJE 1659
la3c COnTINUE 1860

IF(NC.LE.N3TA*I?) GO TO 1850 1862

,ItTE(6.18a01 1863
1~an FnRiI&T324 CtRRELATTON FATRIX ?NCOP.SISrENT) 1864

On
18 r(rNMC.Ef).1) GO TO 1760 18046
r * * * * * TEIT FORJ OVEP-ALL CMNSISTVICY a * a 1867
IfEMPSO 186
C"1 TM 1870 1809
i36C tTEMPSI 1870
C WHFP1 ITEMPU1D CUPFNT MIONTH uSFD FOR~ ALL %'DFPE4CFNT 3TAS 1571
C nTHE,7'I5E. PRCC MTH USED FOR CUeAENT AND SUBSEQUENT ATAS 187
1J70 J4TOfPZNSA 1573
.vAQZNINDP,1 1874
?m;160o
0.1.12 1875
Tv=Tr.1 1876

C CONSFR'JCY CUP'PLETE COPREL MATRIX FOR4 EACH MONTH AND 8Th 1878
4
5M 2150 Kml.p!ST 18?9
C L IS Rflti NU,4eEq, J T3 COLUMN NUMOIR t8'0
,)n P020 L=1,.JSTA 1841
LYZL+NSTA 1ose
Inr 1900 JX10)STA 1883
Jv=J+N3TA 1864
TFrig *) 18l80,1920.1960 tops
IAl"cT;V(I~ 090. 0100.1900 M8b
t8qf q(LJ) 2 JLECA(1,L,J)) 1847006
LTHP(L)CL load
3T&PCJ)2J 1889
Gr 170 t 1800
1QnC IFC(fl.MP) 1910,1910,189' 1841
191C O(L.J) 2 OEFiLECNA(I,L#JX)) 1892000
0T10MXL 11193
JTP' (J)Zjw 1h894
n 1970
fr 19

!1 7(LJ) 2 ORLE(QA(I#J#LX)) 1047000

LymP(L)NJ 1898
J1TP~P (J) Z L 1849q
~r' yO 1970 1900
t94s .tL.J) a OLE(RACPPLJ)) f96100t)
LT 'P (1 )KL Me0
J7IPCJ)UjE 1403
Gm' TM tqln t904
1n!C ' T fFPP
8ir) P,4r0,
I 1930 1905
194C fFCTTF4P) J9e'0#1920e1A48q 1906
197C I(JL)GP(LPJJ 1407
1450 CMTTNUE 1913e
L11'PL~xX1909
C iPECIAL SIJSSCRIPI FOR DEPENDENT VARZADLE 1910
IF (L-0) 144M,23iC.2000 1911
194C R(L.,NqTSA) a nPLECA(#x#L)) 1912000
Gm Tn Rojo 1914
,fIF (1TEMPGT*0) no TO 1990 1915

-13o Ey"IftI? 6
 ?olCRfINSAA)* O1LECA(IM.Lt))1916000
JTAO(NJTAA):LX 1917
20iC CAANTNUP 1918
c tNATRIX COJSI3TINT IF CORREL DOES NOT EXCEen t~o 1919
NoO 1920
N. C2 ft1921
P03c CALL. OVYCR) 1923
c *==uCZuuau 1924
1FfiATR11C.LE.l.) GO TO 2130 iqq.,
WRI (6,0f405 K, I#KOtYRMC 1926
P04C Fn9MAT (/36H4 IN~CONSISTENT CnRREL MATRIX AOJUSTEI),314,p12.3) 19V700
c NITH~DRAW 1928v1931
FACuFAC..2 193a
IFCFAC.Gr...12GO TO 1750 1933
,0'
(%a 1934
tnEN* 1935
IFCN.nT.10) GO TO 150 1930
3tImmo. t1137
On ?nAO' LCIPtl!NOP
1938
On ?o70 LXxt,NVAR 1939
IF(L.E').LX) GO TA 2070 1940
TM~a4 CL.LX)1941
~5U*TM
# Pp 1942
207C CnKTIUE 1943
2PC COSTN~pUf 1944
TFMP=N1NDP*NtNOP 3965
611' S U /TEM P 1946
TE 4 PSI)TQP- 1947
Ir(T~!mP.nT..1) TENP2.1 1948
TmiPm1.-TEMP 19'49
On P120 LM1,NlNnp t950
JT,-I b1
1951
On P110 LX41TPNVAR 95
qCL-LY) a 0LF.CTMPP*TP StJm*TEMP) 193300
IP(LY.L-E.NXNOP) P('.v#L)-:4 CL#LX) 1954
LTPZLTPI'fL) ls
JTP=JTtlP (LX) 15
!Vr(.TI.LE.NSTA) GM TO 2100 1957
IF(lTP.Lt*NSTA) GO TO 2090 1958
LTFtLTP-N3TA 1939
J'I'=JTP.NSTA
1960
WA(T0,t.TP, JTP) = L ,LX) 1961
R4&(lPJTPpLTP!=P(LLX) 1962
?O9C ITM~5LTP 1964
L7PZ,51P 1965
JT~'aTTmP 1966
TF C.TP.LF.NarA) kA(1,.JT0L7P)=R(L.,Lx) 19#18
Plic CnN T I :tiE 141,9
21pc CrNTINUS 1970
;nrTn 2fl30 1971
;)13C I;CnTL;,C.GE.O.) GO TO 2140 1972
-Q ITYF (,S 71) 1 , K, 0T~t!C tq73
r)T'- "C20, 1974
Pl)jC 1F(oNCA.GT.fl) GO TO 1740 197!
;).I Im7NTtiJe 5970
?16c Cfl%TINiJE
MY7
Ir fYTF'4P.El.O) GO TO 1860 3978
TFITTRNS.EQ.2) Go TA 31CO 1979
217C .ep1TyEt, 1 03 1900
C * * a PRINT CORRELATION MATRIX aa ****aa**a*aa 19o~i
IF (TYTA,.L! *n) wGITE (6,24O)PO~) t9p]
?1AC FIRNAT (11394 RAN. C"lNRELATflnw COFFFICIENS F0R MAN?N,13) 1984000
IF(CTTROJS.GT.0) 14t.ITE(6#?1903 Mori5) 19#k!
P19C Fni "AT (//404 Cfl'jSZSTFNT CVQAELATTON MAThIX FOR ?4ONTH.13) 598,000
"Q1TkU6,2?f0O) CI9V(K),1c11.M8TA) 19P7
pIoo ~Fnl~hAT (/3X,3H3TA,18I7) 1948000o
w~1T8C6,a21) 5949
2210 !AWkATC201,19H '4174 CURiOEWI "OI)1910

EvWMTIT 6 14a
 flo n.ao xu,,4 1 qq j
401'"9(60223n) SAtKtIPSA ,2
M~&~
P230 Fnpm47 tt,5~73
q93 000
1990
,;0T~F6.2I~) PkECEDTING g4ON7 AT A8OVE STATION) 1999000
224C iFl-,OA7 (20%,3iITH411
IP3'157A,: 1996
On ;250 K23I143TA 1997
P250 '-:DTClS,1A(K)(RA1,(.L),LUITPDNaTAX) 199S
A2AC crWTI1UE l9q9
ZF(TANAL.LF.O) r' To 3tnC 2000
1F(TIRNS*AF*Q) or TO lbsfl Pool
IF(IPCHNI.LE.nl Gos To 2870 2002
C PU'ICH f5E~.rNTAL ELEmENTS OF MATRIX 2003
t'f P30o xwj,)* 7& 2004
IF o8C
rU 70EQt 2005

P2C "I1 A 2 90 LU*NSTAA,0.TAX 2009

4 M010
1fF PzL.NSTA
P2C Qlikfor O 2011
rn rn 2850 2013
CL * 0 * * * * ECIk~TjTUTE MTS31N. QATA II It II II IN **UU 2(114
P310 F(TANAL.LE.0) GO TO 31PC 2015
19 (TQCSK.LE.6) rO To 2610 2016
,:VAP2N3TA*1 P017
"421 2016
D', pb00 32184YRS P019

2021
2023
01 PSAO 1K31,WqTA 2024
IFCGC(m.N).IF.T) Un To 2580 ?025
c FL3wo r(.;4F.LA?irh oi&TRiv FOR EACH MISSING FLO% 202b
Ntkpso Lt4T 202V
LYZI*NSTA 2029
TIvCOML).Nk.?) 6n TO 2320 2030
117(n M-IL).F'I.T) GO To 2 39f) 2031
Nn7NflNnOp* 2032
lyc"P~yhop2033
x( rnv~q~m-vL)2034
4(PsTNOP.k.VAhi) I O1LE(RAAI,#,LV)) 2035000
T 2330
70 203b
232C Nr~.OPaNTP.Op*l 2037
ZpIT pikP 2038
xr~ ~r')EOC.L)2039
2 0~ECUAI,14L))2040000
.131C (',P0TtDlP .0 P041M00
IrrL.E0.NSTA) GO 10 2JqC A042

0r 23JA( LANITP,NrTA a044

JIII1114TA &1045
1F(i(m.L).fQ.r) r.O To 2.150 2006
Cn.tA.~~n)ro ?C ?340 20'17
~Y9.~Pu~ff'4'*12048
~e~Ifl~T~uP
IIL(RI(,LA))2009000 2
nn TO 2370 20l0
113JC 1(c1LaE.)GO TC 238n 201

,u~bplr.jI4Pt) a rI1LE(RA(IvL.JX)) 2053000

Gn TO 2370 2054
P35C IF2M.&.C3T)G TC 1360 2055
2056
~4tN~Jfl.
OL~fO CIL' .x))2057000
!E~P)*
Gn TO~ 2370 25
a36C TC 1300
400MlL).(.~Q 1019
IT"~PIITE"P*1 2060
a rB.ECA(Iv#LtLA))
2(TnI1p .7001000
c 00 SYH'PI!TRICAL ELEPEIKTS 1081

also y"IbTT b
 a37C Rf1EMPmNItJDP)2R(NINDPo, hP) 26
P3AC cnh?pkuE Z064
p3qC CnN?!NUE 2065
IrCN1INOP.GY.e0) GO TO 2400 2066
?N0~pot 2061
Y(I)S0. 2066
Pria1) 8 1.000 2069000
LX3S.+4STA 2070
R(I-14VAfl) R D0LECQAC1 1KLX)) 2071000
140C TTEHOUNNOP,1 2071'
(inl?d LNl.MINDP 2073
2411C !7fL.ZTfP1=(LNvAR) 2070
C 222:2X8389 2071
24a2C CALL CpOUT (q) 207b
C m=vu2u2077
ITE04PUNINOP41 2076
TF'HPu1 2079
Ip'nezo 2080
nln P0d40 Lx1,JfNnP 2081
TmPUDAAS(R(L, !TEMP)) 2082
IP(T"P.GY.TEpiP) GO TO 2430 2083
?FPPxTmtP 2084
I?P*L 20.S5
P43 GC To 2440 2067

CMNO 0
Ifa 2086
TI:CJ1C.GT.0) GO TO 2450 2090
f~rTtCL..A~..TICG.. 6o TO 2510 2091
C IF MAt~jX INCCPNS13TEN7, OiMIT VARIABLE WITH LEAST 2092
c CoqIIELATION 21093
2'490 IT!'PUNINOP-1 P094
IF(TTP.GT.IT4P) GOr To 240 zq
D)n 2'io LQITI'. ITMP 2096
Dn P'IAO LAMI.TTEMP 2097
if1-,1A)zR(L+1 ,LA) 2098
2-17C 1YtIAUY(L+1) 2099
T 2100
C'~ 11 ?500 L1,lT 1'P

pac(L-L.A)2R(LvLA$1) 2102
P501 (.nNTINUE 2103
-1TNfl I 11p 2104
f'n TOi 2420 2105
C ADDt RANIDCM CrMPCN.VNT TO PRESERVE VARIANCE 210b
p5Ir frmpml20 2107
Q, ?20 Lxtoh 2108
T0~?.~IN;EN(XX)?109
21?C :Ep~NENIX 2110
C C0I4;)UTE FLOW 2111
AL~f..OO~C2*.52112
~ ~ :TEM, AL2113

Z!33C TPP3FMPOl(L)*X(L) 2115

"I".0TEMP2116
1::14 t ) 2 2118
C ADD NEW VALUE TC Sumis nF SOUARES AND CROSS PRODUCTS 2119
Or 4560 L81PNSTAX a22
I( L.El.K) Gn TO 2560 2121
C SuflSCRIPTS EwCEEDING NSTA RELATE TO PRECEOING IIUSTH 212,2
IVml,.NSTA 2123
IF C'.LT.1) ?T'OCM.L) 2124
IFCLX.r?.0) YmPzrt".1,Lv) 2125
!F CTMP.EO.T) GO TO 2560 212b
C C0jII1T AND Usk CpLY QECOROEO PAIRS 21a7

11~A (,ItX#L)uSUMA CIlcL)4*TP 1319

2130
I" ACf.K,L)zt,IIpC,K,L)+4PTMP

IF(LogToNSTA) GO TO 2560 2134

owI4THtT 0116
 .11,4A (I.L.x)xs0PF1 (I#K#L) P136
'%j' (I L. K)03tGA (TI KL) 2139
XFAtiCT.L.K)=XPA1(I,K.L) 2140
c RECri1I2LITE C0ARdELATyflN COEFFICIENTS TO INCLUOE NE% DATA 2141
P340 IF(Nr.A5(I.,L).LE.2) GO TO 2560 2102
TFMP8NCA9CI#KtL) 5
2143
TMPZtSOA(IK.L)SUlA (y9 KL)*SIJMACEL)iTEP)*C QSC!,x.L)SUmi 2144
I(IX.L).SUM8C!,K,L)1TEMP) 2145
C EL1'41hATE PAIRS WITH ZERO VARIANCE PROUUCT 2t46
IF(TMP.LE*O.) GO TO MO6 2147

TMPAUNPA4 C!,KeL)-3UMA(ZigpoL2.8Ufl8(I.KL)/7E"P P149

C RETArp, ALGEBRAIC 1Gb' 2100
IF (T'fFA.LT.O.) TmPvUTWmP 2151
T4'PA=TMPA*TmPA/TPP 2.152
PAcT,K.L)=THf,Tf1PA**.5 P153
ITPUI 2154
LAML 2155
!V(L-.LE.NATA) Gfl To 2550 2156
TF CTTP.LT.1) ITPxl8 2158
LAZI V 0
2159
255C IFUPr(IK) .LT..OO11.OU.SO(TTP,LA).LT..00 13 RA(I#K#L)uO. P160
?tCL .GT.?dqTA, Gn TU 2560 2161
.4A(Tp1i.X)2RACTKoL) 2162
?s~c cr~kTINULE 2163
IT~p=NvQ3*12+t 2164
TFCm.(,E.!TMP) GO TO *580 P165
TC P:O (1,K) 2166b
On' r7O Ls.,PjITA 2167
T'41 C(MIp ) pi160
TF(T?tP.EnT) GO TO 2570 2169
L~2~.8?AP170
IF (TrP,flT.12) ITP:1 2172
tA"(ITP,L.LX) :bCA4 CITP,L.Ly) *1 2173

SIA (ITP.L.LX)m30A (ITP,LKX.1TP*TM4P 2170

SOP (TT0.L,LXx):SQJ CZ(17,t V)#7P*P 217
XPAR(TTP,L.Lx):xPA8(1TP.L,LX)*TP*TMP 2170
I1:('A!)CITPL,LvJLE.2) GO TO 257n~ 2179
TF,'0:P C.AFiCTpDL .LX) 2100
T"P.F:(Or)ACITPL,LX)-SLJ*A(ITP,t ,LX)hsumA(TTP,L,LX)/TEIP)* 2101
1Cr ITT,,.-lm)IT t ),8*UIF4ITPLLX)/TEM40) 2182
JF(TM4P.LF.*. GOf TO 2570 2183
rpi ;1 * 214
ThPPAXPA0i~tlL ,Lx)-3M* CTP,L,LX),MUM8(ZTPLeLX)TE4P 2145
0
TF (T',PA .L. .) T ,Pl-TpMPB 218b
Tm~PA=T A* TP4PA /14P 2147
6:Af!TP,LLx)=T#41p8&TA41PA,..5 218
IF(Ir)CIo) .LT..O'OCI.OD.80C(IVPuL).LT..O0I) RA(ITPLLX)8O. 2189
'T70 CM'TfktlE ?190
p~p cnrT 2191
plqc C"IvTwuIE 2192
POVMCCnkTTNUE 2193
2051C rF(TAb4AL.LE*f) fqO TO 3100O 2190
Cm # * * * CaNVEQT 3TA14 fCARO l.'yTATEs TO FLUO ' *a. 2195
IF C,".'8.LE.I) GU TO 263C 2196
7v~pu~jQS.1* 12197

?FCIOqTAP. 0ITAP) GO TC 2.30 2199

NV& CTUTAP) 2200
P62C LATAP2LOTAP*1 2201
P63C 4PI'TV(hot0) 2202
WRITE L6,2640) 2203
P640 FO-'4AT(3w RFCCQOEO AND gFCOWITITIJTEO FLOW$) 2204

265C Ffrh-Artu P48,1t3i 2206 S

EUNIOIT.T*I

*17
?P S
 1Frx GT.NKTX1 *IT~ v6b: 2209
P66C FIQ4AT C/114 S1'A YfAiflImfjA6x*5TOTAt.I 2210600
fl 1211
'In ;1760 J11IDIYRA 2212
TTpuo 2213
On ?720 Tutle 2214
!4--+ 221!5
TF4Pm(MpX)2236
TwfPc3KEW0~i,I() 2217
3
I" (TMP.NE.O.) TEMPE(TYIP.(TrMP.TrP/p,.)/6..1.)"* ". )*Rsiy"P 221
15(flP(i,K).NF*C) Go Tr 2b~n 2219
IFCFMP.GY.2..AND.SDC!,A) .f?..3I y~w~mZ.,(TEiiP.)*.3/SOCIrX) 2220
IF(CuP.LT.-.OlO1.OR.TI'P.GT..0001) "PXp(-2.)1TMP
2
1F (SX"W9Tp) ififD bq0,2ht0Q 2222
267C 1F(rEmP.GT.T-1P) TEPIN371P 2223
1;n TO 2490 R224
264C TFCTFPP.LT.TiP) TE"PuTiP0 2225
2690 T"-PmTfMP*'%D(Te()+AVCIpK) 2226
-10' KIN tlO.*T-P.0a (I K) 2227
2229

IF(FTMP.LE.01 GO TO 2710
TF40xn. 2232
10 7?0 L31IT P2233

27flC TFPTm~ntX 2235

tF('fMvXI.G?.TEmP) GO TC 2710 2236
,A er *ADJ 2237
?vfn0(a,N).NE.E) GO TC 2710 2238
Z"(T)XA0J1 2239
K I SE14P 240
2YtC ri2241
;?72C ITF'=TTP+10(1) 220a

-v .LE.hATX) Gfl 73 27r6O 22 at

!F('PCHfO.IF..O16n TO 274" 2245

273C F r'v~fAT(214, iPTb) 2247

;!? dC ' .TTE(6plKIi ( I, I I;(I Q IUP
75!)8,1A I Ip 1)O) 2249
276C C~f"TINUE 2250
1w(P!0AS.LE.1) rf TO 2765 2250-1

P tr-ffQCOIf.LE,.o) r TO 2l8I0 2?53

0 0 ? 254
e 0 * * *RECfl4 i1TC ?4F.AN Air 3TANDAFIO DEVIATION
or R770 tal.12 2255

I , xI 0. 2257
???C -j~(t,'()3O 2258
."1 2259
~7q~
7O *0r92lN S 2260
, = -- 0 12262
?PPALOGCQel,X).D)Qe1.K))*.4342q45 2261
5fT()UKj2A I,*TEmP*0 2264

P79C r..",r INuE 22147

I- "MOO) Tuttle 2268
IF-M1AVCZ.X) P269
?F-A=SCI.X) ??70
Tp*) 2211
IP (Tm0.LT*Os) ,"i'OO, 2272
Av (?.W)NTfMP/kNTIR 2273
Sn~j,)8TM**.52271
T-p=SKEw(If%) 227S
F ~ (,W) O *1276
I~~,(,)LE.90 60 TO 2800 2971

CVNI4?tT 6 *left
 5 IW,X,(UANYRS,2T.P3.*AYP0TEMW*TMPA2,OTE4Pi*JI 2?74
I r~ygCAVl31.)(py5P.*SO(I.K)-.3) 0279
PaflC c"'JNtIe 2240
281C Cnk~TI'NUE 2241
LOTAPSNOTAP 26

IFIRC .0) GO Tn(!I(

2930 . 12)
,OTNCTE C.0f) ro Tt,),!3I,12 2281
C PQI")T ADUN8STEDT CAEVIENCA TATlN C 22A61~
a82C
4pfTE(610 228
flfl 217f
TFl62W
C01 PUNCT/3H
ADJ CUCTEY STATINCSTAICSTI 2?31
rn t*60 XUNSTYX,N15TA 232
w~rTE (,101)
!YAO )(A T,(IX*
tulu2) 2303

APIT (6102) (
n~tx),al,2) 237

C84 cO'4PIirE as
2296!iTINFLW
C T~~ 2309
289C IF (ICHS*LE.) GO TO 240 2300
I HP INriS1PGEC 3AIT 231
DP! ?86 K=?45TlPv3T 2302
I ~ ~9~ ~~1,~~tW~2303
~ VA 2314
I (' K) 0 *2306

?Ahr -,nTIP4[E 2320

PC co4m.TINCOBIATONFLW 2301
C 2322

267I("KCS.LE.0) GG TO 2910 2310

j~TP LE.0 CU T5243 2311

PA90x XxI , ic rmpt2313

.(=KlNsTA23314

. K
Ar P0 C. 313
2T'l

On -JA
0:0 TT 2335

c o*A aa CO~xml~i'()+(i#ITGEEI)*CSTATXLISTCRA33) I400 a23

C10 CV.IAG CfQPLTIN OFFCI 2320

PqrC
CM*q5N 23021
Lc .T 2303
Cn MAX ANDL MIN PE*SI DF1W 2323
9ir';=O 0 11 2324
TYPNSU7.,L 2300

P92I(L.E.V uI,( TO 2q) 23.2

~90~ AV tI.L)*RV(KL*TP29
l4V(.L~~AV(,L~12.230

.14. EwiPPXPC5
233
 w'nTTE(6.70)!STA(K),!STACL),RAV(KL) P350

C AVE91AGE LOG!; FCR VbF7 ANDO OQY 3FASONS 2

AVPXK(X)6AV(11,IitAVC12,K)$AV(,() 2353
AVI'N(K)M&VNXCK) 235!
TMP8AV(12.x)+AV(1,K)+AV(2#N) 2357
IF(AVHX(K).GV.TiP)un yu ?960 A3158
AVPV(K)flTMP 2359
IMI(XI022360
Gnf r0 2970 2361
796c AVYNCK)=IMP 2362
Z.N (K)tA. 25313
c IND AVENA(UE STANDOARD DEVIATION 2364
2y7c snav(x):snci1,x)#s$ic2,x) 2365
Dm 299'j 133DI ?366
V'~A(K)U0AV K).SUK)2367
~MpZVCT.~K)AV(I1,UA~tK)2368
lrAm()rFTPG O 29P0 P369
AvMV~(I()"m 2370
tIXI(KINT 2371
2qJC IF1AVmN()AF.TMP)r,0 TO 2990 P372
AVMNJ(X) CrMP P373
IIS(K3a 1 2370
;-iQc clK rrNIJF P375
Avymx (w ) mAX ( )~/3* 2376
AvMI()vAv'!(K)/.5. 2377
ilfV (K) uvAV (K) /ti2 P378
30('c CrnKT!N(IF 2379
I-P!TEfb,100) 23P0
fl3010 X2101TA 2381
JYPZIX(K)232
IT~'!M11(K)P363
3('I ~TF.(,t20IST(K),A(K,&VCK)S0AVK),0CIP),M(ITP) 384
C *~ * *AIPitYrJENERALIZED STATISTICS.* Al 0Aa** 23F5
il?C 1F(tfNkLj-.i-)Gfl To 31rjl 238b
e~ 3041) KCjIWSTA ;S387
<Vrwe N3TA 2348
C rIEfRPE0IATE eCkTH~S 2389
2390
~ 2392
ly"l 3040 TclIi2 2393
C STI~jOAR0 DEVIATION LJNIFCRM, SKEW ZERO 2390
ix I * K) O 2395

,f(T,K)ft5AV(K) P397
10 1030 L21.'JSTA 2398
C ZE~n CnIREL&TIC% WITH OTSE4 STLTIONS ANO PPECEnIkG MONTH 23q%
LY:LNSTA2400
~& (~kL)~O.2401
IF(L.rE.X)r,n TO 3030 P402
C I!NIFn,1'4 SERIAL CCRREL TITERMEDIATE MCNT143 AND INTER-STA 2403
;A(T#*L) 2!4AV v(eL) 2004
;A(TL.K~qA~jjL) 405
3mic CnlTTNLiF 2006
26(7 ,K,K'i):'qV(K,K) A4407
N A(7, K* K)~ *2408
AndC Ld'NTINUF P409
4
c ' fA~i AND SERIAL C1,RREL, WE? &NO PRV SEA504S a410
tMPNRAVCK.IX) .15 pa11
frP,2T142. *3 241a
fF (?mpGT..qA;Ympa.98 241J
IF (TE"P.LT.O) TEMPxOS 24t4
tpaIKX(K) 241!
tv(7?P.X)2AV!NX(K)..1 2416
3(TTP~k.Kx)%TEHP
1 2417

IFr(ITP.LT.1) 172212 2419

AV TTP.K) 2AVlX (K)..2 2420
2A~l~~x~xwzTEAPaila

EX'47E!' &20
 2
IYP:,4x().. 2422
UF(TTP.Lr.1)!TP*1YP+12 2423
AV (T'rP.K) :AVIIX(K) ..1 2424
t.tA (ITP, KkX) --T"P ;0425
ITPOIMN (X 24a26
AVCT1'P.X)XAV4N() 2427
RA(TTPK,KX~mn7P 2428
lTP21MN(IV).t 2429
!FCYTP.Lr.1) iTPUIP 2430
AV (T4,I()AV-lN(IC) 2431
RA CTTP.X.Kk)zTMP 2432
ITPUINdN(g)-2 P433
IF(TTP.T.IjTIPt*1? 8434
Av(T?0.K~AVmAJ(X) ?433
4 (YTP,.Klx)m?tp 2436
C HESJ3S FCR MO'kTIN F(ILLOWING WET 3EA8ON 2437
!F(NHivM?.LT.1)Gfl TO 30bfl A438
ITPZIXCK) 439

T~lm0u(&AVMV(K)-..1AYMN(K))/1EKP 2441
011 1050 !X2I,NMKMHN P442
T4P=T'i 2443

101C AV (Z,X)=Av(ITP,K) .TEMPTYqP 2446

c "r4AS Ffl4 tiOIhTI5 FCLLCwING DvY 6EA30.4 244?
106C IF(N~lNMV.LT.j)Ca T(O 309r, 2448
1TpriNa(K) P449
T F0t00NMN X*1 2050
TFmP. fAVNm(K) -.1 -AVmN R) )/hEfP 2'451
Dn 3070 IX,NNp~X 2452
T4PCIV 2433
I=PN (Id.rX 2454

1,17C AV (I.X2AV(lTP~,K),TEmP*TMP 2456

10AC Cr'kTlNUE 2.457
309C IrKPL90 t1458
19CONSO 2499
rf' TO 1730) 2460
l31)C IF(I YrG'.LF.fl.AND).NP CJ.LE.0,AN0.NPA3S.LF,1) GO Tn 15o 24161
CP * * * * * FL(Iii GENEqATTOk EGLIATIONS 0 *0**0o 2462
NT '0 P2 NSTA 2463
,VAR1k,4TA*1 2464
nr 3200 Y21PI2 2465
1P--T.1 2466
IF (IP.LT.1) IP2 2467
i,, jqqo k1,113TA 2468

c Cf N FA.ATIONS IN CURQEKT PONTY4 ?470

IF IL.GE.P) (r! TO 3120 2471
K:(l,NVAQ) a f)'4L(QAIpKL)~ P072000
-m .4110 LA=L.,'SA 2473
LV'TL-A+IPTA 2474
IF (LA.LT.K) R(LoLA) a ltBL(4A(I,L.LA)) 2475000
IF CLA.GF.A) -I(LLA) ObLtUFA(!,L,LX)) 241600
311C i*aLA,L):mQ(L#LA) 2477
ri Tn 3100 2478
C Cr!P-4ELATInN3 kt7k PPECEDNG 14ONTN 2479
31PC LvzL*N3.4TA 2440
J(LNVA41) 2 tn"LE(AA(1p,~LV) 21181000
0rMItI0 LA=LNSA 2442
41
4fL.L ) u fMLF(QA(IPvI.,LA)) 24A3000
'13C '4(LA.i.)2Q(L#LA) 21084
114C COITINUE 248S
C m2222 2486
CLL CAOLICl1) 2487
c x2musu:333 2488
flo M50 LulpISTA 2489
313C lFTA (l,',L)n4(0) 2490
IF(flTp'iC.Lt.t.) GO TO 3570 2491
w$IE C6,3160)1I,K.OYRIC 2492
316C Fnqt''' CJ0 INCO~iT3TF'.y crRREL MATRIX FOR ta,13g44 KUlp 249J00

-210 EXIBIT b
 t'li4 orR'iSnF603) ?494000
ITANSUZ 49
i~TO 1730 2496
117C !F(flTQMC.(GE.O.) nC TO 3114,0
.r.JTE (bp70)! .,DTPNC

lIAC '[-CFTC I.Km(I .-O)TIMC)aa.S 2500

l19C CnkTTNUE 20
32nCCnNTNUE2502
(C
a tNTU * *
FLa * aE4,0 250i
IF(NPAS.I.tE.1) GLI TV 3240 ?505
321C IF(LSTAT.Eg.NSTAT3 GO Tr 320 2505
(ISTA( ) 2506
LqATtLATAT*1 2507
i-,, TMl 3e,10 2508
12?C 250
ITATUNSTAT+1 25t0
IqATUN5iTAT P511
NI 1230 Kru.N~ITA 251.2

1 !4ET& (1 ,KCL) .La1,PSTA),ALCFT(I'K) .IXl12) 2514

1.3TATONSTAT 251b
TF(TiAS5.LT..PA3S) GO Tr 200 2517
324C *1"I 2)518
:OAASU1 2519
2520
2521
!F CNR4fJ.LF.0) t;O TO 3310 P522
C~a* * * * P,^nj[.TEfn FLCW 1E'ENCES aaa* aaaaaaa aa 252 3
32tC .A=TYRPJ.TYRA P52 a
,IZLYRJ.Tyqh25253
Ir. j , 526
TPjMTHPJ-IMNTH.j 2527
IF(ITP.NE.0) GO TO 3260 2528
1 T p * 122529
3P6! (IrP.LT IPCITP1-T 2 2531
7' 1290
q (l'TA 2533
!~
P,'),
5~ st To 3280 253a

2536
!F iSKFWTT,).FPn..) GL TO 32 qO 2537
T I.5~~*j aKFW (IT P, X) a P H V(CK)*I1 2538
T .*~c 2539
r-(FOrFM)r TO 3210 2540
rf-'vf-EMP) 2501
'T'" (-TAP) 2542
327( -,; rvf()cb*T Pa*TE'Paa(1*#3.).1)/AKEi(ITPM)*8r(Ew(rTMK).I6. P50i
'11 3290
T 25044
3?pr P545

jiztvppj.1 2547
C N 2.!)'UE'CE kr.o M4 8 MON.T1 NO,# JX d YEAR HC. 2548
311C '130-1 ;)5a9
T) 3330 sso
C ITAv'T 41TP4 lqrP OFVTATION AT ALL STATIONS3 ;.551
3 P512
331C 53P0 KcI,K. TA
1321, I"Ev()n. 2b53
C CE-q4ATE 2 YEARS FOR D!SCARDING ?3

2
JVX* 25196
333C IF '"ASqLF.j) Gl; TO 3uCC P5
lciIPkS5.GT.I) %7JO 3340 2538
T~ ISTAT
7'F.' 21559
.ITAPSO R560
(9IAPSQ 25131
334C Pr Io IOTAP 25362
L.IT'P90 .3563
PF~I(ISTAT)'~qTXX,%5Th, (I3TA*(KI@NSA) 251#4
.'src~l~w~l 565
IC04STI.LE.0) GO TO 3380 ?5&6

EVHI4HT 82
6
 Ic(!PASS.Lf.t) GO Tfl 3360 26
139L G.FAn(onyAp) ryE"P#,CQP',K) p,PoITP) 2570
Lq)TAP8L0TAP.1 2571
1(tTTFMP.NF.I3TA(K)) GU T173350 ps72
33hf :rAD(!8TAT, 0 AC,)SlIKKE(g K,)nI6Lu ?S73
l#K'5YA).ALCFT(IK)pXU1,12) 2574
337C Cn TINUE 2511
13$C "'I 339V KZNTYXN3TA F576
ISTAP81STAP.1 ?I??
tFCN.GT.O) )vfJEV(K)xQ9TAP(j9TAP) 2578
33;C :rA)CtjT AT) i f, p(Av(r 4) pSfln K 3XE w CI K)D QCIK) CE TA I K L I L 1 2579
1,-j3TA),AjLcFr(T#K)#j:1#!?) ?Sao
C* ** * * rE'NFfOTE. CORAFLAME 3TANDARI) OEUIATE '* . a *2541
:jr1isCHM C"(pASSj.Ej.
) NMZX2

-T n~cfluTND 'i( IP AS5)

134i20 KutI'.ITA 2583
:1r) '4410 121.12 2S64
&~G 1.1(30.254S

ilitcC n.TT1NI 2540

j F ti. j_ , 0) rn TO 3440 2589
.Q[TE(b,10) 2590
,I)!TE (6,31130 N 2591
3uJ3C rriAT (21H ';FPERATED FLCwS FOR PFIIOD,13) 2592000
IF ('NPA~q.GT.1 ) wR!TEC6,26560) IPA3S 2593
344" 1:510 JZJANj 2594
!=I?*f-tl~l2595
'r '0012,i2 2596
.3597
IF(NSTY.LE.0) Gil To 3460 P598
1(; 34501 K-0j,'ST.% 2599

TV (H.LPF.Ph) rc To' 3500 6

1490 KZIf'T'AA,v9T4 2602
c 1WAN~fle COMPON~hT 2603
Y~F' :0 * 2604
111-70 L3 s 2605
TP=FF P. '4PN
tF%(.C XX) 26C6
3i ?C TVP=T(>0-i )iAN~(,X) 2607
?.-f:TFmp.*%LCFTCZ,'K) p2618
Jm' f4AO Lc1,'STA Pb09
34 7:T ; E w0 T AC I , K pL)*Q1E V (L) 2610
EMP IX) Avr( ,K)+
',vf; 611

* rO-,KI:TEMP 2613
w;' F (K)= TEM 2614
. I, I'.hiNJE 2615
3~' ' TINLJE 2b6

tF~053L.)Go TO 330f 2618

St(L0lT A 0. ,).'O9T 4P) GO0 TC 3530 2619
* F n (T(nTAP) 2620
I.'!TAF'UL0TAVP* 2621
;n Ff 3520 2622
353C TTP=Nlts.!1l 2b?3
IqT 2(x~qTI'..'RTA*NSTX 2624

V 7APII I(
JI UI( TA(K, CQCP ,K).M22.1P) 2b26
YAPPJ()T AP+ 1 2629
iI fhPajsTAV+1 2620
:F (15TAP.IrT.1(STAP) 0)0 71! 160 2630
3IIf'.ITAP(I8T*Ph:!((IyP,K) 2430
355C I"r MN J J A #1 3
i.'o
M-..TXNSTA
WO 2634
2633 a
W(NJ*JVT47rT0) WQZTE(6,266O) CpOMP)151012)

n e)zAvr;,(j,U() ANLCQ
4'ir 63
5nvr!.X1(3,,V(x,a.AVuZI,K)..l.AL6)/PLCGa)tto 2b6

.p3-?
 j v2Jy T4'2
1660 jajA.NJ
ill) 2641

JN*J+% ph26'
IF (JX.LE.0) GO TO 3660 2&44

on 3690 Talolp

IF (MjEF.MAl SC Ta 3640 2b6

c rU*ANiSvFOM TO LCG PFA4S3tw TYPE LII VAPIA7E (FLOO) ?b09

f jt,,)V (2CV',K)-AVG(CToK) ) IND V(1#9) 2651

IF (TPE(. GOC TO) 36rC IaTD~ 652a

I~ K~vgt,1Y) 3,380,.3902655

F (TMP.GT.T!)Y?4PzEPP
t51 265b
5n1 TO 3810 265?
TMPUEMP 2650
15q(' tF(1I4P.1.T.TE"P)
qn To 3610 ?659

3~C TU!M,~O2681
?I~FT~i*~0
K)?b62 I~h)4AYCT

ZF(TTMI'.Lf.0) Gn TO 363C
TFV'0S0 * 266b
LvZrST(KsLolT4 33
j3?C rFv~mE1P+(4(.X1 bb9

3
~((
AO
VO110. QWt~CTK)~ Q(MtK) 0. 2071

P614
!6QC Cr"NTINIIF
NITHOREV 2615
c 261b
I(t 31atTP
- ITF (6, 100) 1314 X)*Jyo(117(l)plt# 3,1 2617
I~fti'CW.L.C')G TI) 36 20111
2079
v'iTff,c??f0) LaTA(K), jx*gf(0C1),Tl1i2)
3b"C CONTINIUC WbO
37C C 0 . T INJEF 26A
1Ff ,raj.LE..1) r 0 TO 3720 2b8
5710) JUJA,*4j ?6A3

Cfl~4PwTF Cr4TJN'ATlfN FLaO.3 2881

c

,=x+S TA 26L*9
1TP=NTI%7h(KA, JP&33)
(!4,)CO.2691
T ?b9?
-i !SAO SLZ,T '

N(A ?S
I C;hI TK
1).
I.L
269b
37nc C-i'kul
j71r Crmr.VtIUF 2b51
in lzia ?b,48
37?C f;,:<.LT.N' 1jn1 c.0
sb0
IFi4MxGtr.))GO TOJ Van
C 1. 44X A~Sn mjN GkKERATfO FLnw$ 0B00*0 2700
1r(JK:L./IG- TO 1970 2701
INUFFICtEftT 210a
C 3919PAXMIN IF NE"ANIN0 VEA43
IFVYXE.V?.O.A. .NjoL.k~.Yo XG)CO TO 136l 2103

171c T1U4,I0 ~2105 arGb

On 3800 K3J0TwxolYP 210?
r.Pi ALV.OAi
I4 PC lot?, NAV MO l3s 6-140 14, 14-MO 15
n3,00
t27i 01
 3740 i'41)(1.K)ENT 2709
C M4IN 0ALENDAp4 lit P..27, 11IN mOj20g.m 1g9, 54vmu 30 2710
[in~ 3750 1218.30 Rill1
3750 31"OCTvK).T 2712
C TMP R A.14OP ?CO.P 2 434.Mq VOLUME, 7MPP a I-MO) 2713
TFI'Dafl* 2719
TMiPEO, 2715
AVG0(Kx)Uoo 2718
Ntlan2717
IF (JTRNS.OT.0) mm(N.1) *1EXNC312*1 2719
,rl 37q0 JU1IJ 2720
~fln170 1.1.12 2721
N~2I11 272
HVIA+
2723
14PA: 0CM, K) 2724
yGl(lk)XAVG0CPO*TM'h 2725
N'12 N a#1. 2726
I V (TA4PA. GT . SH~I
(J. w) 1 3m 01 9 ) mTM pa 2727
1F (MP A LT .V9IJ(X CK) ).SMOC!)X) TMP A 2728
FCT14A.GT.3ltfl13,K))St 0C(13.K) TpA 2729
IFCTPA.L.I~if)(S~i))3t0(2#K~nMPA2730
T14P2TNP#YMPA 2731
TF,"P2 EMP+THP~'A 2732
IP~e9.LT.8,GO TO 3760 2733
TriP2ywr'.e)(HM6,) P734
2735
IFU%.LT.56)r0 TO 3710 2737
TF'PTE~p (*v5'J.K)?739

*.n~ TO) 37A0 2741

0 Jl-AC 1t'C4K) XTMP P7 42
7177' ri-il I . K) 13TE,iP e743
37AC Cr''NYIIUE .37440
fi. T INuj
379C!r P745
c AVEQAGE MONT4LY FLOW 2746
TV!'PZN0 P747
'4rV(KEAVGqCK) /TFIiP 27?48
3#Oj I-IPNTIN.JE 2749
4QfTF (4,10) 2750
IF(IR?4.Gf.))AITE~o3PO)PINi 751
3II r's1f4AT (1271H *AXImJ4 VULUPFS FORl pE'4IO,1303H OFpta. 27500
I4-1H YEARS f]F 9CONOEO Ahnf PECflNSjTUTkC FLOWS) 2753fl00
IF(?TN3.LE.li4QTE(6,3PiOC)N.Nj P754
' 14AT (12114 44YIMUM VOLL14ES FCR PEPICJO,13,3M OF,14, el'55006
12%k YFARI OF 5Y~ITHFTIC FLrtS) 275800
.UITF~~plol 757
)n 144 X3NSTXX#ZYP P759
iT ,,J3tAVrfO(K)+,5 2760
lo 31'3o 121,I'5 2761

364C C-f INUE 2760

k f ha fl0) ;!7865
'541.0 KNTV'\,LTP 2767

38AC C!N I I NLIE 2771

c TRApSFFR L&Cx TO RECOWITYTUTED FLOW$ ?772
IF(TTQNS.GT.Ol,3 TO 292(t 2773
3070 . a NyMXG 2774
rTO 3 190 P775
3880 uiYR
aj 2776
loan jF(,lP&SS.LE.j) GO TO 3900 2777
2778
60 TO 3310 2779
1V(TPAss.LE.NIA.%s) GO TC 3J40 21 80

-25. Evs-IdYT 6
 1l'A33I 2781
c Go to ,4k~ oa20
3911 ' IV(NYRG.LI.') G13 TO 150 V7A3
rfNj .GT.W1:?c)NJ=NykG ;P784
,YPAUNY146-s:j a761
nmr Tn 33no 2746
END 2781
S319B00UTTP. CPPLUT(9!() 1001

nnLhIALF PJFCI31ICN R,13,RX 1003

Crw'4r]h 0TRTiC,CNINDP.B 1004
..vAPvNXNnP*1 1005
gm a ul1N1INOP 1006

,'r) to K=1.NVAq 1008

if ri.K)m~y(j~v) 1009
2c crmkTZNIjF 1010
Tv(qINflP.rnV.1)SO Tfl 30 101l
O~m~(1)aiC1)1013

~ :~ o XO~JVAN1016
C 1 FI4IJ MATRIXK/~(
*1 1015
If1 40 K-u2, JV A4 1018
AO )(2KNTO 1019
On 60 JzW.*,T'jDP 1020
,,, rso I~1.Ir 1021
I022

yrCJ.EG.K) !Z'j Tn 60 1024

.(KJ)2C.K)/R1CKtK) 1025
r. ("T'vrywgjf 1026
'' 70 :1TTP 1027

7r wN F A i LV-R R( 1029
~'fclzn8Af ;4(ix K)) 1030
IrrcC'P.rt"..)(00) rGc an 1031
2TL!CU1 1
t03-3
', TUr tot033
PC, r; &VAR )f,? (W.VAic) /P ( ~K) 1030
(a ~ * PACK SOLUTION.aaaAaaaaaaa 1035
(N~h)1y=12"T11P,4AR)1036
'n 100 Ix2,;4P'DP 1031
j.5%VAq~- I 1038
t: I * 11039

'40 2,019 1041

'n)+L1042

1 11fIO%
P.U I 006
tTTD
a 0cr; Nf) 1005
C 10) NJjI I: i145ir F II AR!V 'AH 1002
Ct r'7F~g
Ct 0 t AS
'4-C+ICjI*4X(.NI 104
1009
' j 1)~ J YKQ'. E~OT R~N toot
-,-rTr(IN N(TTA) (

c ~ N 'v 1Ef.cflP lA RFAqYIG

Q? AIN 1002
C a-I-~
,aOAFa&~rP
aIRF 1r(2aaC(PRtYj?)a3 INf a a
a3a 00

c * a a aFc'*i:./2aai)aa 1011
C OnyfF~TkP BL1tr TC. 7twE INEE!H. OION~
, lw"EQE 1012
c NYLFE C' Nhrqto 1013~hFP
c ANTANIA r,11V121/ v- r M 0LClrjE)AT.N 101a

I :!U 1013

e!,4livf 6 *ab"
 Trr.~m 1 S71Pj~lot?
i~ TYy&~c~Njlots

11N-,~F Pis IV 1021

FCou3.3S5271lb1 8E-t4 1022
NG PtNZ$N Sf. *l'FC~C
R NJ 10171

END 1021
t026

ESdT6
 MXUBIT 7

INPUT DATA 723-X-L2340

CARD VARIABLE CO0MMMS

A Three title cards, first must have A in column 1.

B First specification card.

1. IYRA - Earliest year of record at any station.

2. INMMH - Calendar month number of first month of water year.
3. IANAL - Indicator, positive value calls for statistical
analysis routines.
4. IRS - Number of years in each period of recorded and re-
constituted flows for which maximum and minimum values
are to be obtained, dimensioned for 100.
5. NYRG - Total number of years of hypothetical flows to be
generated.
6. NYMXG - Number of years in each period of generated flows
which maximum and minimum values are to be obtained,
dimensinned for 100.
7. NPASS - Number of consecutive passes, each pass consisting of
a new group of stations which can be correlated with
specified stations in previous passes, dimensioned
for 5.
8. IPCIQ - Indicator, positive value calls for writing recorded
and reconstituted flows and generated flows on Tape 7.
9. IPCHS - Indicator, positive value calls for writing stat-
istics on Tape 7.
10. NSTA - Number of stations at which flows are to be generated,
not required if flow data are supplied. NSTA + NCOHB
(C-i) dimensioned for 10.
C Second specification card.
1. NCOMB - Number of combinations of stations, the totals of
which are used to obtain maximum and minimum flows,
dimensioned for 2. If positive, provide D and R cards.
2. NTNDM - Number of tandem situations, compares sum of monthly
values of upstream stations with downstream station
and adjusts if value Is less than sum and that station's
value has been estimated or generated, dimension for 10.
If positive, provide F card.
3. NCSTY Number of consistency tests. Adjusts standard deviation
of a dependent station in tandem with an independent
station to prevent frequency curves from crossing,
dimensioned for 10. If positive, provide C card.

EXMBIT 7
CAo VAMANZ 0aer's
C (Cont'd)
4. IGM - Indicator, + 1 calls for reading generalized stat-
istics and using for generation, + 2 calls for cow-
Vitin generalized statistics from flov data and using
for eneration.
S. NPW)J - Number of projections of future flows from present
conditions, usually 0.
6. IYRPJ - Year of start of each projection.
7. M)HPJ - Calendar month of start of each projection.
8. L RPJ - Lst yea of each jrojection, number of recorded and
reconstituted years plus number of projected years
dimensioned for 100.

D Identification of combination, NCOB (C-i) sets of D

and N cards.

1. INSAC - Number of stations in this combination, dimensioned

for 10.
2. ISTAC - Station number (!SrAC values).

E Combining coefficients, MCOB (C-i) sets of D and

9 cards.

1. NSTAC - Sme as D-1.

2. CSTAC - Coefficient of flow used for adding, corresponds to
respective items in D-2.

F Identification of tandem situation, NDM (C-2) cards.

1. I&M - Station number of downstrea station.
2. NMX - Number of upstream stations, dimensioned for 10.
3. IS"T - Station number of upstream station (NSMX values).

G Identification of consistency test, NCSTY (C-3) cards.

1. ISTX - Independent station number.

2. I lY - Dependent station number.

H Flow data, cards in any order, omit if (B-3) is

LANAI
not positive follow all flow data cards by 1 blank
card (I card).
1. Cole 2-4, Station number
2. Cols 5-8, Year number.

2
EnHBIT 7
CARD VARIABLE COMMENTS

H (Cont'd)

3. Cols 9-14, 15-20, etc., Flow in desired units. Units should be

selected so generated flows will not exceed 999,999. Use -1 for
missing record. If record for entire year is missing, omit card
for that year.

Card blank after Col 1 to indicate end of flow data,

omit if IANAL (B-3) is not positive.

J Identification of stations in previous passes to be

used in current pass, supply only if NPASS (B-7) is
greater than 1. The variables NCOMB, NTNDM, and
NCSTY apply to the current pass only.
1. NCOMB - Number of combinations of stations, the totals of
which are used to obtain maximum and minimum flows,
dimensioned for 2. If positive, provide D and E cards.
2. NTNDM - Number of tandem situations, compares sum of monthly
values of upstream stations with downstream station
and adjusts if value is less than sum and that station's
value has been estimated or generated, dimension for 10.
If positive, provide F card.
3. NCSTY . Number of consistency tests. Adjusts standard deviation
of a dependent station in tandem with an independent
station to prevent frequency curves from crossing,
dimensioned for 10. If positive, provide G card.
4. NSrX - Number of stations from previous passes which are to
be used with the additional data in current pass as a
means of maintaining consistent flows between groups of
stations, number of stations from previous passes
plus number of new stations dimensioned for 10.
5. rSTA - Station number of station in a previous pass which is
to be used in current pass (NSTX values). Must be
in same order as stations first appear.

Note: Flow data for current pass supplied As described for H card and
follow data with a blank card (I card), supply NPASS-l sets
of J, H, and I cards (also D,E,F, and G, if necessary) when
NPASS greater than 1.
K Preceding-month correlation coefficients for first
station, omit if IANAL (B-3) is positive (NSTA cards).

1. ISTA(K) - Cols 2-4, Number of first station.

2. ISTA(L) - Cols 5-8, Number of station from 1 to NSTA (B-1O) on
successive cards. If IGNRL (C-h) = 1, only first
card is used.

3. RA(I,K,LX) - Cols 9-14, 15-20, etc., Correlation coefficients for

successive months between flows at first station and
preceding-month flows at stations from 1 to NSTA (B-10)
on separate cards. If IGNRL (C-4) = 1, only general-
ized coefficient (in cols 9-14) is given.

3 EXHIBIT 7
CARD VARIABLE COMMENTS

L* Current-month correlation coefficients, omit if IANAL

(B-3) is positive, (NSTA-1) pairs of L and M cards.

1. ISTA(K) - Cols 2-4, Number of station, progressing from K = 2

through NSTA (B-l0) stations on different sets of L
and M cards.
2. ISTA(L) - Cols 5-8, Number of station, progressing on different
cards through all stations from L = 1 to K-1.
3. RA(I,K,L) - Cols 9-14, 15-20, etc., Correlation coefficient for
each successive calendar month between flows at station
K and concurrent flows at station L (12 items). If
IGNRL (C-4) = 1, only generalized coefficient in cols
9-14 is given.

M* Preceding-month correlation coefficients for remaining

stations, omit if IANAL (B-3) is positive. Paired with
L card.

1. ISTA(K) - Cols 2-4, Same station number as on corresponding L

card (L-1).
2. ISTA(L) - Cols 5-8, Number of station, progressing in same order
on different cards through all stations from L = 1 to
NSTA (B-lO). If IGNRL (C-4) = 1, only card with L = K
is used.
3. RA(I,K,LX)- Cols 9-14, 15-20, etc., Correlation coefficient for each
successive calendar month between flows at station K
and flows in preceding month at station L (12 items).
If IGNRL (C-4) = 1, only generalized coefficient in
Cols 9-14 is given.

N Generalized frequency statistics, omit if IANAL (B-3) is

positive or IGNRL (C-4) does not equal 1.

1. TSTA(K) - Cols 2-8, Station number for NSTA (B-l) stations on

successive cards in same order as supplied by L cards
(L-l).
2. AVMX(K) - Cols 9-14, Average mean logarithm for wet season (3
months).
3. AVMN(K) - Cols 15-20, Average mean logarithm for dry season (3
months ).
4. SDAV(K) - Cols 21-26, Average standard deviation for the 12 months.

* Sets of L and M cards are required for each station from K - 2 to NSTA.

EXHIBIT 7
CARD VARIABLE COMMENTS

N (Cont'd)

5. MOMX(K) - Calendar number of last month of wet season.

6. MOMN(K) - Calendar number of last month of dry season.

0 Mean logarithms, omit if IANAL (B-3) is positive or IGNRL

(C-4) equals 1.

1. ISrA(K) - Same as (M-1).

2. AV(I,K) - Cols 9-14, 15-20, etc., Mean logarithms for successive
calendar months.

P Standard deviations, omit If IANAL (B-3) is positive or

IGNRL (C-4) equals 1.

1. ISTA(K) - Seine as (M-1).

2. SD(I,K) - Cols 9-14, 15-20, etc., Standard deviations for suc-
cessive calendar months.
Skew coefficients, omit if IANAL (B-3) is positive or
IGNRL (C-4) equals 1.

1. ISTA(K) - Same as (M-l).

2. SKEW(I,K) - Cols 9-14, 15-20, etc., Skew coefficients for suc-
cessive calendar months.

R Flow increments, omit if IANAL (B-3) is positive or

IGNRL (C-4) equals 1.

1. TSrA(K) - Same as (M-l).

2. DQ(I,K) - Cols 9-14, 15-20, etc., Flow increments for successive
calendar months.

Five blank cards with A in Col 1 of first should follow last Job.

Note: Cards K through R are not required if cards H and I are supplied. Cards
K through R are as punched by computer when IPCHS is positive.

EXHIIBI r 7
 EIBIT 8

SUIUARY OF REQUIRED CAMDS

723-x6-L23.O

I See Icard)
NPASS-1 iiZ (Ie H I I)I-
sets (1)(e Rcad
J COMB NTN4DM NCSTY NSTX I STA fSTA T.. NS X Val a
Mi I lakcard, after all ta,

(1) H UT Q Q Q Q Q (12 Values) format 214, 12F6.0

NCMT cards G I
1ST!

IVUM cards F ISSNSMK 1ST? MM?...18 au

NCO (PTA CSTAC ICSTAC I ..... IWAC aluer

st[D" ( ST"Rn
IMAC IISTAC1 .. 0.. 1 NSTAC I Valm?
C NCOMB NMD NCff! IaNRL NpROJ rm]RJ KWJ LYRPJ
B DWH rAKAL mKRw8 " IAS IB'
NYR I11 IP I STA:

A WT ITLE1 CARD

A TTLE CARD

Notes:

(1) SupplY Only if ZANAL 0B3) Is Positive. Papeat U card tar each
station-year of data before supplying I card.

EXHIBIT 8
 SUbMY OF REQUIRED CARDS
Continued
7'23-X6-L231eO

R DI NI
IST le
Q IT SKE SKW...(2Vle)oiti IGIL(4

NSTA sets (4)

(4) N A AVMX AVMN SDAV MOMvc M.gIM

L TA K) IST-A L CRA 1,K,L) HA(2,K,L) ..... RA(12,K,L)

(2)(4) K TA 1 ISTA(L) HA 1,K,LX) RA(2,K,LX) .... (12,KLX)

(2) L designates correlation with current month and LX with preceding month. If
IGNRL.(C4) =1, only one (generalized) coefficient is given following station
numbers on each card and only 1 K and M card is used for each K station, with
L = K. Use sam format as H card.
(3) Repeat set of L and M cards for each K station except first.
(4) Omit if IANaI (B3) is positive.

2
MIB~hr 8