PB161598

ecltmcai vjote

Moulder laboratories

92c. 97

TECHNIQUES FOR COMPUTING

REFRACTION OF RADIO WAVES

IN THE TROPOSPHERE

BY

E.J.DUTTON AND G.D.THAYER

I. S. DEPARTMENT OF COMMERCE

IATIONAL BUREAU OF STANDARDS
THE NATIONAL BUREAU OF STANDARDS

Functions and Activities

The functions of the National Bureau of Standards are set forth in the Act of Congress, March 3, 1901, as amended by Congress in Public Law 619, 1950. These include the development and maintenance of the national standards of measurement and the provision of means and methods for making measurements consistent

with these standards; the determination of physical constants and properties of materials; the development of methods and instruments for testing materials, devices, and structures; advisory services to government agencies on scientific and technical problems; invention and development of devices to serve special needs of the Government; and the development of standard practices, codes, and specifications. The work includes basic

and applied research, development, engineering, instrumentation, testing, evaluation, calibration services,

and various consultation and information services. Research projects are also performed for other government

agencies when the work relates to and supplements the basic program of the Bureau or when the Bureau's

unique competence is required. The scope of activities is suggested by the listing of divisions and sections on the inside of the back cover.

Publications

The results of the Bureau's research are published either in the Bureau's own series of publications or

in the journals of professional and scientific societies. The Bureau itself publishes three periodicals avail- able from the Government Printing Office: The Journal of Research, published in four separate sections,

presents complete scientific and technical papers; the Technical News Bulletin presents summary and pre- liminary reports on work in progress; and Basic Radio Propagation Predictions provides data for del- the best frequencies to use for radio communications throughout the world. There are also five series of nonperiodical puolications: Monographs, Applied Mathematics Series, Handbooks, Miscellaneous Publications,

and Technical Notes. A complete listing of the Bureau's publications can be found in National Bureau of Standards Cir< 460, Publications of the National Bureau of Standards, 1901 to June 1947 (S1.2S), and the Supplement to National Bureau of Standards Circular 460, July 1947 to June 1957 (SI. 50), and Miscellaneous Publication 2

: July 1957 to June 1960 (Includes Titles of Papers Published in Outside Journals 1950 to 1959) (12.25); avail- able from the Superintendent of Documents, Government Printing Office, Washington 25, D. C.
NATIONAL BUREAU OF STANDARDS

technical Ntote

97

October 17, 1961

TECHNIQUES FOR COMPUTING REFRACTION OF RADIO

WAVES IN THE TROPOSPHERE

by

E. J. Dutton and G. D. Thayer

NBS Technical Notes are designed to supplement the Bureau's regular publications program. They provide a means for making available scientific data that are of transient or limited interest. Technical Notes may be

listed or referred to in the open literature. They are for sale by the Office of Technical Services, U. S. Department of Commerce, Washington 25, D. C.

DISTRIBUTED BY

UNITED STATES DEPARTMENT OF COMMERCE

OFFICE OF TECHNICAL SERVICES

WASHINGTON 25, D. C.

Price $ 1.50

TABLE OF CONTENTS

PAGE

ABSTRACT

1. INTRODUCTION 1

2. LIMITATIONS TO RADIO RAY -TRACING 3

3. AN APPROXIMATION FOR HIGH INITIAL ELEVATION

ANGLES 4

4. THE STATISTICAL METHOD 5

5. SCHULKIN'S METHOD 6

6. THE FOUR-THIRDS EARTH MODEL 8

7. THE EXPONENTIAL MODEL 10

8. THE INITIAL GRADIENT CORRECTION METHOD 12

9. THE DEPARTURES-FROM-NORMAL METHOD 13

10. A GRAPHICAL METHOD 16

11. SAMPLE CALCULATIONS 16

REFERENCES 26

APPENDIX I. 27

APPENDIX II. Tables of Refraction Variables for the

Exponential Reference Atmosphere

Tables of Coefficients, Standard Errors of

Estimate and Correlation Coefficients for use in

the Statistical Method 36

Refraction Variables in Various N(h) Atmospheres

41

Initial N Gradients, AN , in the C . R. P. L.

Exponential Atmosphere. 49

Schulkin's Method Sample Computations 50

Departures Method Sample Computations 54

TECHNIQUES FOR COMPUTING REFRACTION OF RADIO

WAVES IN THE TROPOSPHERE

by

E. J. Dutton and G. D. Thayer

ABSTRACT

Eight methods of computing atmospheric refraction of radio rays

are discussed with appropriate theoretical background. These methods

are:

(1

(2

(3

(4

(5

(6

(7

(8

The high-angle, or astronomical, refraction case

The statistical method

The low-angle, or terrestrial, refraction case (Schulkin's method)

The four-thirds earth model

The exponential model

The initial gradient correction method

The departures -from-normal method

A graphical method (Weisbrod's and Anderson's method).

Sample computations are included for each of the above methods

TECHNIQUES FOR COMPUTING REFRACTION OF RADIO

WAVES IN THE TROPOSPHERE

by

E. J. Dutton and G. D. Thayer

1. INTRODUCTION

If a radio ray is propagated in free space, where there is no

atmosphere, the path followed by the ray is a straight line. However,

a ray that is propagated through the earth's atmosphere encounters

variations in atmospheric refractive index along its trajectory that

cause the ray-path to become curved. The geometry of this situation

is shown in figure 1. Figure 1 defines the variables of interest. The

total angular refraction of the ray-path between two points is designated

by the Greek letter t, and is commonly called the "bending" of the ray.

The atmospheric radio refractive index, n, always has values slightly

greater than unity near the earth's surface (e. g. , 1. 0003), and approaches

unity with increasing height. Thus ray paths usually have a curvature

that is concave downward, as shown in figure 1; for this reason downward

bending is usually defined as being positive.

If it is assumed that the refractive index is a function only of

height above the surface of a smooth, spherical earth (i. e. , it is assumed

that the refractive index structure is horizontally homogeneous), then the

path of a radio ray will obey Snell's law for polar co-ordinates:

n r cos = n r cos ; (1)

the geometry and variables used with this equation are shown in figure 2.

With this assumption t may be obtained from the following integral:

n 2

, 2

, = - ( cot e — , (2)

1, 2 J n

v e
FIGURE 1. Geometry of the Refraction of Radio Waves
FIGURE la. Differential Geometry of Radio Ray Refraction.
FIGURE 2. Bending Geometry on a Spherical Earth with

Concentric Layers
-2-

which can be derived as shown in Appendix I, and also by Smart [ 1931] .

The elevation angle error, e , is an important quantity to the radar

engineer since it is a measure of the difference between the apparent

elevation angle, , to a target, as indicated by radar, and the true

elevation angle. Under the same assumption made previously € is given

as a function of t, n, and G by

_ cos t - sin t (tan 6) - — ,

e = Arctan J £5 >- . (3) — tan 6 - sin t - cos t tan 6

n o

s

The apparent range to a target, as indicated by a radar, is defined

as an integrated function of n along the ray path,

R h

R~ = j

ndR = J Hn-e-- < 4 >

However, the maximum range error (R minus the true range) likely to

be encountered is only about ZOO meters, hence the evaluation of (4) is

not of great importance unless one is dealing with an interferometer or

phase-measuring system.

The preceding material should suffice to show the importance of

radio ray bending in radar systems evaluation and allied types of radio

propagation work. Unfortunately, the integral for t (2) cannot be evaluated directly without a knowledge of the behavior of n as a function of

height. Consequently, the approach of the many workers in this field

has been along two distinct lines: the use of numerical integration

techniques and approximation methods to evaluate t without full knowledge of n as a function of height, and the construction of model n- atmospheres in order to evaluate average atmospheric refraction.
-3-

The following sections are devoted to a discussion of these methods.

2. LIMITATIONS TO RADIO RAY -TRACING

The user should keep in mind that the equations given in the

preceding section are subject to the following restrictions of ray-tracing:

(1) The refractive index should not change appreciably in a wavelength.

(2) The fractional change in the spacing between neighboring rays

(initially parallel) must be small in a wavelength.

Condition (1) will be violated if there is a discontinuity in the

refractive index (which will not occur in nature), or if the gradient of

refractive index, dn/dr, is very large, in which case condition (2) will

also be violated. Condition (1) should be satisfied if

(dn/ dh) per km < Q N kc

where refractivity, N, is defined as N = (n-1) x 10 and f is the carrier

frequency in kilocycles [Bean and Thayer, 1959]. Condition (2) is a basic requirement resulting from Fermat's principle for geometrical

optics. An atmospheric condition for which both conditions (1) and (2)

are violated is known as "trapping" of a ray, and it can occur whenever

a layer of refractive index exists with a vertical decrease of N greater

than 157 N-units per kilometer. A layer of this type is called a "duct",

and the mode of propagation through such a layer is similar to that of a waveguide [Booker and Walkinshaw, 1946] . Taking into account refractive index gradients, a cutoff frequency may be derived for waveguidelike propagation through a ducting layer [Kerr, 1951] .

In addition to the above limitations, it should be remembered that

the postulate of horizontal homogeneity, made in order to use equation

(1), is not realized under actual atmospheric conditions; some degree
-4-

of horizontal inhomogeneity is always present.

3. AN APPROXIMATION FOR HIGH INITIAL ELEVATION ANGLES

A method may be derived for determining ray-bending from a knowledge only of n at the end points of the ray path, if it is assumed

that the initial elevation angle is large. Equation (2) in terms of refrac

tivity, N, is equal to

n,, e,

1,2 " j cot dN • 10

-6 (5)

N

l ,e

1 assuming n = 1 in the denominator. Integration by parts yields:

N2' 6

2

1,2 'j cot 6 dN • 10

N

r

e

i

N cot • 10

-6

N2' 9

2

N

i- e

i I

9

2- N

2

N

de- 10

-6

er

nx Sln 9

(6)

Note that the fraction, N/ sin 0, becomes smaller with increasing 6 for

values of 6 close to 90 . If point 1 is taken at the surface, then 6=0 1 o and NIs= N . Then for = 10°, N = and 0^ = ir/ 2, the last term of

o 2 2

(6) amounts to only 3. 5 percent of the entire equation, and for the same

values of N_ and 0_ but with = 87 mr (~ 5 ) the second term of (6) is

L L o

still relatively small (~ 10 percent). Thus it would seem reasonable to

assume that for

> 87 mr (-5°),

the bending, t , between the surface and any point, r, is given sufficiently

well by

1,2

N cot x 10

-6

N , r r

N , s o
-5-

or t , = N cot 9 x 10" -N cot x 10~ . (8) 1,2 s o r r

The term -N cot 8 x 10 is practically constant and small with resr r pect to the first term, for a given value of 9 and r, in the range

> 87 mr. Thus t is seen to be essentially a linear function of o 1, c. N in the range ^ 87 mr. For bending through the entire atmosphere,

s o

(to a point where N = 0), and for > 87 mr, (8) reduces to

t = N cot x 10" . (9)

s o

o For initial elevation angles less than about 5 the errors inherent

in this method exceed 10 percent (except near the surface) and rise quite

rapidly with decreasing .

4. THE STATISTICAL METHOD

Another method for determining high-angle bending is the statistical

linear regression technique developed by Bean, Thayer, and Cahoon

[ 1959}- It has been found that for normal conditions and all heights the

right-hand integral of (6) is approximately a linear function of N (0 ,

r constant) for > 17 mr (~ 1 ) and that the second term of (8) tends to o be constant. Thus (6) reduces to a linear equation,

t = bN + a, (10)

1 , c. s

where b and a are constants (as in tables I - IX) and N is the surface

s refractivity.

The form of (10) is very attractive, since it implies two things:

1) t may be predicted with some accuracy as a function only

of N (surface height and constant), a parameter which may be

s o observed from simple surface measurements of the common
-6-

meteorological elements of temperature, pressure, and humidity.

2) The simple linear form of the equation indicates that, given a large number of observed t versus N values for many values

1 y Ct S of h and 8 , the expected (or best estimate) values of b and a can

o be obtained by the standard method of statistical linear regression.

This is what was done to obtain values listed in tables I - IX.

Tables I - IX also show the values of the standard error of estimate, SE, to be expected in predicting the bending, and the correlation

coefficients, r, for the data used in predicting the lines. Linear interpolation can be used between the heights given to obtain a particular

height that is not listed in the tables. For more accurate results, plot

the values of t from the tables (for desired N ) against height, and then

plot the values of the standard error of estimate on the same graph.

Then connect these points with a smooth curve. This will permit one to

read the t value and the SE value directly for a given height.

5. SCHULKIN'S METHOD

Schulkin has presented a relatively simple, numerical integration

method of calculating bending for N-profiles obtained from ordinary significant-level radiosonde (or "RAOB") data [Schulkin, 1952]. The

N-profile obtained from the RAOB data consists of a series of values of

N for different heights; one then assigns to N(h) a linear variation with

height in between the tabulated profile points, so that the resulting N

versus height profile is that of a series of interconnected linear segments

Under this assumption, (2) is integrable over each separate linear

N-segment of the profile (after dropping the n term in the denominator,

which can result in an error of no more than 0. 04 percent in the result),

yielding the following result:
-7-

2(n

l " n

2

) AT

l,2

(rad)=

"J

cot0dn = Tanei+ tan8? n

1 ,0

i

or

2(N - N ) x 10" 3

At (mr)= —

• (11)

1 , 2

v tan + tan

1 Cm

For the conditions stated above, this result is accurate to within

0. 04 percent or better of the true value of At , an accuracy that is

1 > c. usually better than necessary. Thus it is possible to simplify (11) further

by substituting 6 for tan 0; this introduces an additional error that is less

than 1 per cent if is under 10 (~ 175 mr). Now (1 1) becomes

2(N

X - N

2

) AT

i,2 (mr) =

e +0 ' Vnd0

2

mmr - {1Z)

1 Cm (mr) (mr)

where may be determined from (33), Appendix I.

The bending for the whole profile can now be obtained by summing

up the At for each pair of profile levels:

. .. y

2(N

k

- N

^i»

T

n

(mr)= Z 9

k + 9 k+| • (U)

k=0 (mk

r) (Sfr

1)

This is Schulkin's result. The degree of approximation of (13) is

quite high, and thus most recent "improved" methods of calculating

t will reduce to Schulkin's result for the accuracy obtainable from RAOB

or other similar data. Thus, provided that the N-profile is known, (13)

is the most useful form for computing bending (for all practical purposes)

that should concern the communications or radar engineer. Soine other

methods have been published which are actually the same as Schulkin's,
-8-

but have some additional desirable features; e.g. , the method of Anderson

[ 1958] employs a graphical approach to avoid the extraction of square

roots to obtain 0, . k

6. THE FOUR-THIRDS EARTH MODEL

Perhaps the earliest attempt to utilize a model of atmospheric

refractive index for the solution of problems in microwave radio propagation dates back to 1933, when Schelling, Burrows, and Ferrell [ 1933]

published their discovery that radio propagation through an atmosphere

with a constant refractive index gradient of -/^a, where "a" is the radius

of the earth, was equivalent, for purposes of calculation, to radio propagation over an airless world of radius 4a/ 3. This was a great simplification, since it meant that for the calculation of radio field strengths,

etc. , the atmospheric refractive index could be ignored provided that

4a/ 3 was entered in the calculations instead of "a" wherever it appeared.

This method quickly became known as the "four -thirds -earth, " and has

formed the backbone, until very recently, of radio refraction calculations

since its introduction.

The 4/3 earth method, as originally proposed, suffers from two

serious shortcomings, only one of which may be overcome by use of this

kind of a model. They are as follows:

a) The gradient of refractive index near the earth's surface that

is implied by the ratio 4/3 ( —4 ON/ km) is valid only for certain

areas and at certain times, e.g., temperate areas in winter; the

gradient implied is less than average for temperate climates in

summer, always much below average for tropical climates, and

greater than average for arctic climates.
-9

b) The gradient of refractive index implied by the 4/ 3 earth model

is nearly constant, decreasing with height at uniform rate, and

thus the values of refractive index implied quickly reach unrealistically low values; free space value (N = 0) is attained at about

eight -kilometer -height.

The first of these drawbacks may be avoided by a simple modification of the original 4/ 3 earth theory. All that is required is to pick a value of the "effective earth's radius factor", e.g., 4/3, which is consistent with the meteorological data that are available for the area under

consideration. Hence, a location that has a normal gradient of refractivity

near the surface, of -100 N-units/ kilometer , would have an associated

effective earth's radius factor of 1 1/ 4, and the effective earth's radius

for this location would be 11a/ 4, or about 17,500 km.

The shortcoming of the 4/ 3 earth model listed under "b)" above is

an objection to the effective earth's radius theory in general, and hence

cannot be avoided by a change in the size of the effective earth's radius

factor (except by making the factor a function of height).

With the above considerations the following recommendation is

made: when dealing with problems involving ground-to-ground communications systems or other types of low-altitude radio propagation problems where the ray paths involved do not exceed one, or at most two,

kilometers above the earth's surface, the effective earth's radius method

should be used to solve the associated refraction problems. The user

should refer to the tables in Appendix II, where effective earth's radius

factors are tabulated along with other refractivity variables. Table A-l

may be entered with N and table A-2 may be entered with AN(N ) subtracted

from the N value at one kilometer above the surface. In both these tables
-10-

Probable Errors when using Effective Earth's Radius Model for a Ray

2 , with = 0; using h ~ d / 2 ka.

True Height of

Ray for an

Exponential Profile:

Calculated Height and Percent Errors

For Normal Conditions

dN/dh = -50/ km

For Superrefraction:

dN/dh = -100/ km

1.000 km 0. 987 km, 1. 3% error ~ 0. 95 km, 5% error

2.000 km 1.950 km, 2. 5% error ~1.8 km, 10% error

linear interpolation will suffice for any practical problem. The variables

listed in these tables are for the exponential model of N(h) that is covered

in the following subsection.

When the effective earth's radius treatment is used, height is

calculated as a function of distance, for a ray with 9 =0, with the

equation h = d / 2ka, where d is the distance, k is the effective earth's

radius factor, and a is the true radius of the earth (~6373 km). The table

above will serve as a guide to the errors likely to be incurred when using

this equation, assuming as a true atmosphere an exponential N(h) profile

as given in the following subsection.

7. THE EXPONENTIAL MODEL

An exponentially decreasing refractive index in the troposphere

has been recognized for some time [Bauer, Mason, and Wilson, 1958;

Anderson, 1958] . Recently Bean and Thayer [ 1959] introduced an exponential model for N(h) based on an analysis of observed profiles from many

climatic areas (mostly in the U.S. ). With this model, the value of N as a function of height is given by the equation

N(h) = N exp{-c h} , s e (14)
-11-

where c is a function only of N , and is thus a constant for any given

e s profile, and h is the height above the surface. The quantity c can be

related to N and AN by

s

C

e = 4v^ [ ' U5)

so that a relationship between N and AN would fix each exponential

s profile of form (14) as a function of the single variable, N . A brief

description of the development of such a relationship is given in the

following paragraphs

.

If it is assumed that N(h) is indeed an exponential function of

height, then the gradient of N(h) would also be an exponential function of

height. The most extensive amount of data with which to evaluate the

coefficients in the exponential is that of AN (the value of N at one kilometer

minus the surface value, N ) which has received wide application in radio

s propagation problems. Thus one would expect ~ = kiexp{-k

2 h} (16)

to take the form

AN = k exp{-k },

J. c*

for the special case of Ah = h = 1 kilometer. Examination of available

AN data reveals that k is dependent upon N , i.e., the higher the surface

value of N the greater the expected drop in N over one kilometer. Further

examination indicates that

k = k N , 2 3 s and the resultant equation,

AN = k exp{-k N } , (17)

1 is
-12-

may be solved by least squares. The least squares determination is

facilitated by converting (17) to the form

In | AN | = -k N + Ink , (18)

J 5 X

that is, expressing the natural logarithm of AN as a linear function of

N . The values of k, and k„ are established from 888 sets of 8-year

s 1 3

means of AN and N from 45 U.S. weather stations. The results of this

s study are shown graphically in figure 3, and the least squares exponential fit of AN and N is given by

AN = -7. 32 exp{0.005577N } . (19)

With this equation the C.R.P. L. Exponential Reference Atmosphere

[Bean and Thayer, 1959a] is determined; the profiles are completely

defined by equations (14), (15), and (19).

Ray tracings have been computed for this model covering more

than the normal range of N , and the results found in Tables X through

XVII may be used to predict t for any normal combination of Nso,6 , and height [Bean and Thayer, 1959a] .

The exponential atmosphere is considered to be an adequate solution

to the bending problem for any larger than about 10 milliradians and

all heights above one kilometer.

8. THE INITIAL GRADIENT CORRECTION METHOD

The importance of the initial gradient in radio propagation, where

the initial elevation angle of a ray path is near zero, has long been recognized. For example, if dN/dh = -l/a (the reciprocal of the earth's radius), then

the equation for t is indeterminate, an expression of the fact that the ray path

remains at a constant height above the earth's surface. This is called
""Sin>3

c. • A

'• ii

'.

r

yv*: •

£:&-'. •.

vi.

A- • ; : v*\ • • :\ '-

'. ' \ \ 1 •\ o o o O CD 00 o o

CO o o o

rO o m

cj

—

OO

rOo

CD

rOOOCO

rOOin

rOO

roOrOrOO

CM

rOO

rOO

rOO

CMOCO

CMO

CMO

CO

CMOinCMO

M"CMOrOCM

HSoomcoomouo

\Zco

a,3Co• •-I m

CD

CDuV

at

f<0 wo

l<I
-13-

ducting, or trapping of the radio ray. The effect of anomalous initial

N-gradients on ray propagation at elevation angles near zero, and for

gradients less than ducting, (| dN/ dh | < 157 N units/ km, or dN/ dh > -157 N units/ km), may also be quite large. A method has been developed

for correcting the predicted refraction (from the exponential reference

atmosphere) to account for anomalous initial N-gradients, assuming that

the actual value of the initial gradient is known [Bean and Thayer, 1959b] . The result is

t, = t (n , e ) + [t (n *, e ) -t._.(n , e )], (20) h h s o 100 s o 100 s o where t, (N ) = t at height h, for the exponential reference atmosphere

corresponding to NSo, and N # is the N for the exponential reference

S atmosphere that has the same initial gradient as the observed initial

gradient; t is t at a height of 100 meters.

This procedure has the effect of correcting the predicted bending

by assuming that the observed initial gradient exists throughout a surface layer 100 meters thick, calculating the bending at the top of the

100-meter-thick layer, and then assuming that the atmosphere behaves

according to the exponential reference profile corresponding to the observed value of N for all heights above 100 meters. This approach

has proved quite successful in predicting r for initial elevation angles

under 10 milliradians, and will, of course, predict trapping when it

occurs.

9. THE DEPARTURES-FROM-NORMAJL METHOD

A method of calculating bending by the use of the exponential

model of N(h) together with an observed N(h) profile is given by Bean

and Dutton [ I960] . This method is primarily intended to point out the

difference between actual ray -bending and the average bending that is
-14-

predicted by the exponential N(h) profile and is a powerful method of

identifying air mass refraction effects.

The exponential model described in subsection 7, can be expected

to represent average refractivity profile characteristics at any given

location, but it cannot be expected to depict accurately any single refractivity profile selected at random, even though it may occasionally do so.

In order to study the differences between individual observed N(h) profiles and the mean profiles predicted by the exponential model, a variable

called the A-unit has been developed; it is defined simply as the sum of

the observed N at any height, h, and the drop in N from the surface value,

N , to the height, h, which is predicted by the exponential profile for the

s given value of N .

Thus

A(N , h) = N(h) + N (l-exp{-c h}). (21)

s s e Thus (21) adds to N(h) the average decrease of N with height, so

that if a particular profile should happen, by coincidence, to be the same

as the corresponding exponential profile, the value of A (N , h) for this

s profile would be equal to N for all heights. The above analysis shows

that the difference between A(N , h) from N , 6A(N h), is a measure of the

s s s departure of N(h) from the normal, exponential profile:

6A(N ,h) = A(N , h) - N = N(h) - N exp{-c h} . {ZZ)

s s s s e

It seems logical that the application of the A unit to bending would indicate

the departures of bending from normal, in some way, just as it indicates

departures of refractivity, N, from normal. This is indeed the case as can be seen in figure 4, where for an N = 313. exponential atmosphere.

s A(3l3.0,h) is plotted on one set of graphs for various typical air masses,

and the corresponding bending departures from normal are shown in the
i ro

ro

s

uo

o

o

>Q

do

•i-t +->o

0)

ni

en

cu

05

CO

•r-l

CJ

o

CO

O

i

w

D

Oi—

(w>|) 1H9I3H
-15 -

second set of graphs corresponding to the same air masses. Obviously,

the bending departures between layers are highly analogous to the A unit

variation. It can be seen from figure 4 that the similarity exists, although

it is less, for higher initial elevation angles. The similiarity also decreases with increasing height, owing to the fact that the bending departures from normal are an integrated effect and at low initial elevation

angles, are more sensitive to N-variations at the lower heights. This

causes an apparent damping of the bending departures from normal at

greater heights. However, the A-unit variation is not similarly influenced;

hence a loss of similarity arises at large heights above the earth's

surface.

If (21) is differentiated and substituted into (2) the following equation

results:

k

h N

k+

T

0,h = TN

(h) + I

- 8. +9, [AA(N

s

)] Xl0

" 6

' (23)

(rad) (ra

Sd) k=0 (rid) frtk)

N

k

where

AA(N ) = AN(h) + A[N { 1 - exp( -c h)} ] = AN(h) + N c exp(-ch)Ah

s s e see'

t (h) is the value of T tabulated for various atmospheres in tables X - s XVIII, 9. and 9, are in milliradians and must be from the N exponential k k+1 s - atmosphere used. AA(N ) is obtained from subtraction of the A value at

s layer level, k, from the value of A at layer, k+1 . The A value may be

obtained by adding any given N(h) value, obtained from RAOB or other

similar data, to a value of N [ 1-exp {-ch} ] for the same height which may

be obtained from figure 5. Since t (h) has been calculated only for a few

s of the exponential atmospheres, these being the N =209.9, 252.9, 289.9,

3 13.9, 344.5, 377.2, 494.9, and 459 . 9 atmospheres, one of these
m-pp;

LUX

(/)

</)

LU

>

o

ICL

X

o>

I

c/>

w>i Nl V1H9I3H

in

IT)

* " w

.xz ri o

LO

o

Q.

D

X h

<D

1 o

o 2
-16-

atmospheres must be used in the calculation of bending by the departures

method. The selection of the particular atmosphere to be used is based

on the value of the gradient of N, dN/dh, between the surface of the

earth and the first layer considered. In table XVIII are shown the ranges

of the gradient for the choice of a particular exponential atmosphere.

10. A GRAPHICAL METHOD

Weisbrod and Anderson [ 1959] present a handy graphical method

for computing refraction in the troposphere. Rewriting and enlarging

( 1 1), one obtains

T(mr) = 1 500 (ton 6, + tan 9, „) ' (24)

where r will be the total bending through n layers. Terms for the denominator can be determined from figure 6. Equation (24) is essentially

Schulkin's result with only the approximation, tan 6 = 9 , for small

angles omitted.

The procedure in using figure 6 is as follows: Enter on the left

margin at the appropriate N - N(h). Proceed horizontally to the proper

height^ h, interpolating between curves if necessary. Use the solid

height curves when N - N(h) is positive and the dashed curves when N

- N(h) is negative. Then proceed vertically to the assumed 9 and read

500 tan 9 along the right margin.

11. SAMPLE CALCULATIONS

The following problem will serve to illustrate the application of

the various methods of calculating bending.
1000

1000

GRAPHIC REPRESENTATION OF

SNELL'S LAW FOR FINDING

500 TAN

h in thousands of feet

6 in mr. o

irvrr in i ii -u-j:

J I000

FIGURE 6
-17-

A particular daily set of RAOB readings from Truk in the Caroline

Islands yields the following data:

height above the N value

surface (km) (N units)

0.000 400.0 = N

0.340 365.0

S 0.950 333.5

3.060 237.0

4.340 196.5

5.090 173.0

5.300 172.0

5.940 155.0

6.250 152.0

7. 180 134.

7.617 125.5

9.660 98.0

10.870 85.0

What is the total bending up to the 3. 270 km level at initial elevation

angles of 0, 10 mr, 52.4 mr (3°) and 261.8 mr (15°) by (a) Schulkin's

approach, (b) the exponential model, (c) the initial gradient method, (d)

the departures from normal method, (e) the use of regression lines, and

(f) the graphical method of Weisbrod and Anderson? Since the gradient

between the ground and the first layer is

AN 365. - 400. , Ah"

= 6~M0 = -102.9 N units/ km,

and this is a decrease of N per km that is less than the -157 N units/ km

required for ducting, no surface duct is present. However, should a surface duct have been present, it would have been necessary to calculate

the angle of penetration,

9 = \l2[N

s

- N

h

-156.9 (Ah)(in km)] ,
-18-

to find the smallest initial elevation angle that yields a non-trapped ray.

Any initial elevation angle less than 9 c annot be used in bending calculations.

(a) Schulkin's approach of (13) yields the results shown in table

XIX for mr, table XX for 10 mr, table XXI for 52.4 mr and table XXII

for 261.8 mr, where 0,,-, is determined from (33) in Appendix I using

r, = a + h, , and a is the radius of the earth, k k

a = 6370 km.

It should be remembered that 0=0, 10, 52.4, or 26 1.8 mr only for the

first-level calculation, and that thereafter 0, is equal to the 0, ,, computed

k k+1

for the preceding layer, e.g. , for the second layer of table XIX (9 = Omr),

0. =6.15 mr, which is the 9. , , calculated for the first layer, k k+1 '

(b) The exponential model solution may be found by using tables X

through XVII. Interpolation will usually be necessary for N , 9 , and so

height; this interpolation may be done linearly. In practice, one of these

three variables will often be close enough to a tabulated value that interpolation will not be necessary, thus reducing from 7 to 3 the number of

interpolations necessary. Since in the problem for N =404.9, h=10.0

s km and 9 =10 mr

o

T

n in n (in \ = i5 -°84 mr,

0, 10.0 (10 mr)

and at h = 20. 00 km, 9 =10 mr.

o

T0,20.00(10mr)

= 15 - 946mr'

and thus by linear interpolation for h = 10.870 km, 9 = 10 mr.

o

T

0, 10.870, ,10 mr, = 15 ' 084 + «".9*6 - .5.084, ^g^I^gg

= 15. 159 mr.
-19-

Similarly for N = 377. 2 in the exponential tables,

s x .

, = 13. I20mr.

0, 10.870, (10 mr)

Again using linear interpolation, but now between the N = 377. 2 and

N = 344. 5 atmospheres, the desired value of t at 3. 270 km for N = 360.

s s and 8 =10 mr is obtained,

o Thus

400. - 377.2

T

0,10.870 (1 mr)= 13 - 120 + (15159 - 13120) 404.9-

= 14. 798 mr.

For the 9 =0,52.4, and 26 1 . 8 mr cases, by similar calculations, using

o linear interpolation:

T

> 10.870,(0mr)

= 2l - 386mr >

T

0, 10.870,(52.4 mr)

= 5 ' 8l6mr'

T

0, 10. 870, (261.8 mr)

= lZ70mr -

(c) The initial gradient correction method may be used if one determines the N * which corresponds to the observed initial gradient and

then applies (20). The initial N gradient is -102.9 —:

, which, as km

can be seen from table XVIII, corresponds to the N 450.0 exponential

atmosphere. Therefore, using the exponential tables of Bean and Thayer

[ 1959] and (20) to determine the bending for the 6 =0 mr case, one finds

by linear interpolation

T

10, 000 (0)

= T

10, 000 (4 °°- °" ° mr) + I T

100 < 450 ' °' Omr) - T

!00 H00. 0, Omr)]

= 21. 309 mr + [5.908 - 3.65 7] mr = 23.560 mr.
-20-

The t (400. 0, mr) is determined by linear interpolation between the

404. 9 and 377.2 atmospheres . At h = 20. km as given in the tables

T20,000(Omr)

= T

20,000 (4000 ' 0mr) + [T

100 (450 ' ' 0mr) - T

100

<400 ' - 0mr)

= 22. 191 + [5. 908 - 3.65 7] mr

= 24. 442 mr .

10, 870 - 10, 000

Hence by linear interpolation

T

10.870 (0 mr) = "' 560 + [24.442 - 23.560]. ^ ^ _ ^^

= 23. 637 mr.

The bendings for = 10, 52.4, and 26 1 . 8 mr are as given below:

8o

= 10iT 10.870(10)

=15 - 053mr

8o

= 52 - 4;T

10, 870 (52.4, = 5864mr

eo

= 261 - 8;T

10,870 (261.8)

='-280 mr.

(d) To use the departures -from-normal method of determining

bending it is first necessary to know the atmosphere which must be used

for the calculation. In the problem,

dN i / - TT" , = 102. 9 N units/ km,

dh initial

which is within the range of the N = 450. exponential atmosphere, as

can be seen from table XVIII. Thus one will use table XVII to determine

the 0's and the t's in the N = 450. exponential atmosphere, or one can

s use the exponential atmosphere tables of Bean and Thayer and (33) in

Appendix I.
-21-

For an N = 450. atmosphere:

s TTM in x (10. 870) = 30.776 mr, N (0 mr)

s

TTM Mn ,(10.870) = 19-414 mr, N (10 mr)

s

t

tvt 1^7 A ,(10.870) = 7.024 mr, N (52.4 mr)

s

t

kt ha, a ,(10.870) = 1.506 mr. N (26 1 . 8 mr)

s Equation (33) of Appendix I should be used for the interpolation in preference to linear interpolation, although, if no tables or other facilities

are present at the engineering site for easy acquisition of square roots,

linear interpolation will suffice. Proceeding in table XVII with (33) of

Appendix I for the first layer at h = . 340 km and the 9 =0 mr case:

'Sr7~ 2U. - r ) ~~ 0=9or+ x 10 - 2(N si- N)

o

= 6.388 mr. The remaining 8's for the various la/ers are shown in

table XXIII. To determine the value of A at the bottom and top of the

layer, one makes use of (21) or figure 5. First, however, one must

determine the value of c in (21) to be used. Usually interpolation will

be necessary in table XVIII, but in the N = 450.0 case it is not possible,

s and thus the straight N = 450.0 exponential atmosphere values are used.

From (2 1)

A(N , h) = N(h)+N [ l-exp(-ch)],

o S and for the layer running from h = to h = 0. 340 km, figure 5 yields

377. 2 [ l-exp(-cO)] = 0. 0,
-11-

and 450. 0[ l-exp(-c x 0. 340)] = 32. 8,

and therefore, A(450. 0, 0) = 400. + = 400. 0,

and A(450.0, 0. 340) = 365. + 32. 8 = 397.8,

whence

AA = A(450. 0, 0. 340) - A(450. 0, 0) = 397. 8 - 400. = -2. 2 N units

Therefore, the departure term of (23):

-z r V ^wAA(tvJ '

becomes

365.

r 1 -2. 2 = + 0.689 mr. 0+6. 388

400.

The remaining calculations are tabulated in table XXIII for the 9 =0 mr

case, in table XXIV for the 0=10 mr case, in table XXV for the

o =52.4 mr case, and in table XXVI for the 6 = 26 1 . 8 mr case. The

o o sum of the departures for the mr case is

k

» -z r

. Nk+ i

I ) 6, + 9, ,

i AA (N = - 5 . 335 mr.

k4

k k+1

L

s

\

N

k
-23-

(e) Determination of the bending is required in part (e) of the

problem by using regression lines. By (10), using table VII and VIII,

it is found for the 9=0 mr case that at 10.0 km (from table VII)

o t = (0. 1149)(400.0)- 18.5627 ± 7.5227

= 27. 3973 ± 7.5227 mr,

and at 20.0 km (from table VIII)

T

0, 20.0 =(0. 1165)(400.0) - 17.9573 ±7.5131

= 28.6427 ± 7. 5131 mr.

Thus, by linear interpolation

t

i.-2

=to.io.87= 27 - 3973 ' (28 - 6427 - 27 - 3973 Coo:to:oo * 7 - s227

? (7.5227 -7.5131) ".87-10.00

20. 00 - 10. 00

T = 27. 5056 ± 7. 5218 mr

1 1 &

Similarly for the remaining 0's,

t. „._ . = 13.9548 ± 0.9701 mr,

1 , c { 1 mr

)

t, ,.„ , = 5. 2186 ± 0.0817 mr , 1, 2(53. 4 mr)

T

! ?/?A! « * = 1-2695 ± 0.0158 mr.

1 , Z(Z6 1 . 8 mr)

(f) Determination of the bending by means of the graphical method

of Weisbrod and Anderson yields, from figure 6, for 500 tan for the

first layer:

h(m) At = mr At 9 ^ 10 mr At =52.4mr At 9 = 261.8 mr

o o o 0.000 0.0 5.0 26.2 134.0

0.340 3.0 5.8 27.0 134.0
-24-

which yields for the bending in the first layer.

At =0 mr At =10 mr At = 52 mr At = 26 1 . 8 mr

o o o o 11.67 mr 3.24mr 0.66 mr 0.13mr

Similarly, the bending for the entire profile may be obtained, and shown

to be At =0 mr At = 10 mr At = 52.4 mr At = 26 1 . 8 mr

o o o o 24.42 14.00 mr 5.32 mr 1.18mr

The answers to the several parts of the problem are summarized

in the table which follows on page 25. Bending values for the assumed

profile, from a method which exponentially interpolated layers between

given layers and then integrated between resulting layers, assuming only

a linear decrease of refractivity between interpolated layers, are included

for the sake of comparison. The computations were performed on a digital

computer.

The reason that the answers to part (e) vary so radically from the

remaining answers for the =0 mr case and not so much for the =261.8

o o mr case is the fact that the accuracy of the regression line method increases

with increasing initial elevation angle, . It must be remembered that the

statistical regression technique, like the exponential model, is an adequate

solution to the bending problem for all 's larger than about 10 mr, and

all heights above one kilometer.

The reason that the answers in part (f) and part (a) agree more

closely than with any other of the answers is because (24) is, as mentioned

before, Schulkin's result with only the approximation, tan = for

small angles, omitted. For this individual profile the bending obtained

from an exponential atmosphere does not give particularly accurate bendings; however, for 22 five-year mean refractivity profiles, figure 8 shows

that exponential bending predicts accurately within 1 percent of the average

bending for these five-year means. Figure 7 shows the r.m. s. error in
i

i 1

.. _._

L

— *fl s* SS-i CM — C\J i*

.§:

->l

9 co

CD o \

II£

"_i

ii ^ ~

o >

• Ct>

" -1 r

degression

1

1

— -

Regression

—

^ /

1

- <

L fc_

<b *

«^j ; /

1

• v; 1 *

f S-i

, 5

II

»5 • j

•

s- . // h

5V/ |<

£ /

£ * » /

Co <*; 5 / \3 v;

*J

ll \

v

/**

>

> ki

i . - y / 1

~f^J/

'i r s / // <oki

\ c ..-•' i t

..•*'* y .-••• /

.••'" 'H-

< v- k CD

C\J S

C\J

C\J OC\J 00 CD CJ O 00 CD CJ O

00O

ph- 0) a

01M

J3 o H

CD

do

•r-l O

O in

\i) ~z_ s < Pi Q be

£2 < •IH o a:

sf ^^ Tl _j e _i

Oh ^ o oro ?

M

O

M

o w

Q>

co O a

C\J

o

eo

M

O

Oh

r-W

SIAId %
en

UJ

\-

O

_J

si

i

o

•rHu

<u A

Cu

VIO

O

43

4->

O

•i-t o

U

<D Pi

M-lo

o

Vi

•iH

oU

00W

o

CVJ

LO LO

SNViavamii/\i ni j-
-25-

Summary Table of Refraction Results for the Sample Computation

Problem Method Bending Bending Bending Bending

Part used inmrat inmrat inmrat inmrat

G = Omr 6 =10mr 9 =52.4mr 6 =261. 8mr

o o o o

a. Schulkin's

Method 24.248 14.008 5.341 1.196

b. ExponentialModel 21.386 14.798 5.816 1.270

c Initial

Gradient

Correction

Method 23.637 15.053 5.864 1.280

d. Departures from

Normal

Method 25.441 14.858 5.350 1.143

e. Statistical Regression

Method 27.506 13.955 5.2186 1.2695

±7.522 ±0.9701 ±0.0817 ±0.0158

f. Graphical

Method 24.42 14.00 5.32 1.168

Comparison

(exponential

layer interpolation) bending 24.171 14.104 5.343 1.168

predicting bending at various heights as a per cent of mean bending (not

including super -refraction)

.

In summary, it is recommended that the communications engineer

either use the statistical regression technique or the exponential tables of

Bean and Thayer [ 1959] without interpolation (i. e. , pick the values of height,

Nso, and that are closest to the given parameters) for a quick and facile

bending result, keeping in mind the restrictions on these methods. However,

as mentioned before, use of Schulkin's method is recommended if accuracy

is the primary incentive.
-26-

REFERENCES

Anderson, L. J. , Tropospheric bending of radio waves, Trans. AGU 39»

208-212 (April, 1958).

Bauer, J. R. , W. C. Mason, and F. A. Wilson, Radio refraction in a cool exponential atmosphere, Tech. Report No. 186, Lincoln

Laboratory, MIT (August 27, 1958).

Bean, B. R. , and B. A. Cahoon, On the use of surface weather observations to predict the total atmospheric bending of radio waves at

small elevation angles, Proc . IRE 45_, No. 11, 1545-1546(1957).

Bean, B. R. , G. D. Thayer, and B. A. Cahoon, Methods of predicting

the atmospheric bending of radio rays, NBS Tech. Note No. 44,

(1959)-

Bean, B. R. , and G. D. Thayer, Models of the atmospheric radio refractive

index, Proc. IRE 4 7, No. 5, 740-755 (1959a).

Bean, B. R. , and G. D. Thayer, CRPL exponential reference atmosphere,

NBS Monograph No. 4 (October 29, 1959b).

Bean, B. R. , and E. J. Dutton, On the calculation of non-standard bending of radio waves, J. Research of NBS 63D, (Radio Prop.) 259-263

(May- June, I960).

Booker, H. G. , and W. Walkinshaw, The mode theory of tropospheric

refraction and its relation to wave guides and diffraction, pp. 80-

126 of Meteorological Factors in Radio Wave Propagation. Report

of a conference held on 8 April 1946 at the Royal Institution, London,

by the Physical Society and The Royal Meteorological Society.

Kerr, D. E. , Propagation of Short Radio Waves, pp. 9-22 (McGraw-Hill

Book Company, Inc., New York, 1951).

Schelling, J. D. , C. R. Burrows, and E. B. Ferrell, Ultra-short-wave

propagation, Proc. IRE 2 L 429-463 (March, 1933).

Schulkin, M. , Average radio-ray refraction in the lower atmosphere,

Proc. IRE 40, 554-561 (May, 1952).

Smart, W. M. , Spherical Astronomy, Chapter III (Cambridge University

Press, London, 1931).

Weisbrod, S. , and L. J. Anderson, Simple methods for computing tropospheric and ionospheric refractive effects on radio waves, Proc.

IRE 47, No. 10 (October, 1959).
27

APPENDIX I

The approximate relation between 9 and 9 is derived here. This

relation holds for small increments of height and small 9's. The relationship was used in making all sample computations in preceding sections

Since for small 9's

9 Q

cos 9 = 1 - -j- and cos 9

z = 1 - — , (25)

and knowing6 .„ ,, r = r + Ah (figure 2)

then substituting in (1) yields

^2 2

G

2

G n

2

( ri + Ah) i^l -- )= n^ (\ - -± , (26)

2 2 2 or 9 9 9

Vl+ n.2

Ah-n

2

r

1 --n

2

Ah -\ = v -n r -j . (27)

Dividing by r

n Ah n 9

2

9

2

9

2

, 2 ^2 Ah 2 ~ I #,», n

2 + -— - —- -„2 — T

= n

i " niT " U8)

-> Ah G

2

Since the term -n_ —- is small with respect to the other terms of

2 r 2 r

(28) it may be neglected, and thus:

n Ah B

Z 9

2

n

2 + ~— " n

2 T

= n

i " ni T (20)

9

2

9

2 n Ah

or

_n

2 T2-12, = " n

i T " ~T~ + (n

i " n

2

)

' {30)
-28-

n r

n

z

If one now divides both sides of equation (30) by n and assumes

n l n

2 = n -n and — 1 (30) may be arranged to yield

1 2 n2

e

2

£ e^ t ^-h

-2(n

1 -n2

) . (31)

Writing (31) in terms of N units,

0. (mr) = 9^ + i^- x 10

6

-2(N -N ) (32)

2 1 r

.

12

if 9 is in milliradians .

Generalizing (32) for the kth and the (k+l)th layers,

yy 2 < rk+ r

rk

) 6 Vmr>= Vmr)+

r Xl ° - 2{N

kNk+ l

}

' 133)

Also from the geometry shown in figure 2, a useful relationship

for t can be obtained. Tangent lines drawn at A. and B will be respec

tively perpendicular to r and r , since r and r describe spheres of

refractive indices n, and n concentric with o. Therefore,

angle AEC = angle AOB = f* ;

also, in triangle AEC

angle ACE = 180° - angle CAE - angle AEC = 180° -9 -0 . (34)

But from triangle DCB

angle ACE = angle DCB = 180° - t -9 (35)
-29-

Since angle DBC and 9 are vertical angles, (34) and (35) are equal.

Thus

180° - t - 8 = 180° -G - <|> ,

'1|2

= *

+ le r 92»' (36)

Now since <j> i*1 radians = d/a, where d is distance along the earth's

surface:

t = f

+ (6 - 9 ), (37)

1, Z a 1 2

or the bending of a ray between any two layers is given in terms of the

distance, d, along the earth's surface from the transmitter (or receiver),

the earth's radius, a, and the elevation angles 8 and 8 (in radians) at the

beginning and end of the layer.

If one considers figure la, Snell's law in polar coordinates can

be obtained from the more familiar form of Snell's law:

at layer n + An:

n cos (9 + At) = (n + An) cos 9 (38) or n(cos 9 cos At - sin 9 At) = n cos 9 + An cos 9 or n(cos 9 - At sin 9) = n cos 9 + An cos 9

. „ ^ An

- At sin 9 = cos 9 n An

or At = - cot 9 , (39) n which, in the limit, becomes the integrand of (2).

Now from (37)

and since

and

A9 = Ac}) - At, (49)

As = rAcj),

. Q ~ As rA<J> cot 9

_ Ar Ar

then

a a Ar A$ = — cot 9 . r
-30-

Then (40) becomes

AG = -- cot 6 + —cot . r n

(41)

If one considers the limiting case of (41), then

Km

Ar —

An —

AG = — cot G + — cot G

r n

dG = -- cot + — cot G

r n

becomes

(42)

because the layer n + An shrinks toward n as r + Ar shrinks toward r.

Integrating (42)

G

S

2 r

2

n

2

tan 6 dG = J

\ '^r + f — J n

1 n

1

yields

in | sec 6 1 j = Hn\ r | + £n|n|

n.

or :n cos

1 cos G.

. n n

2

r

2

n

i rl

or n

i r

i COS 9

l =n

2

r

2

C ° S V (43)
31-

APPENDIXII

TABLES OF REFRACTION VARIABLES FOR

THE EXPONENTIAL REFERENCE ATMOSPHERE

The following table of estimated maximum errors should serve

as a guide to the accuracy o£ the tables.

Errors in elevation angle, 9;

4 milliradians

4 mr. < 6 < 100 mr.

o >100 mr.

± 0. 00005 mr.

± 0.000005 mr.

± 0.00004 mr.

nearly

independent

of N . s

Errors in r, c (in milliradians):

N = 450 404.8 377.2 344.5

s •

e =o ±o.ooi

o

313 252.9

0.00065 0.0005 0.0004 0.0003 0.0002

200

0.000

=1 0.0003 0.00015 0.0001 0.00008 0.00006 0.00005 0.000

o0=3° 0.00004 0.000025 0.00002 0.000017 0.000015 0.000013 0.000

o

Errors in R , R, R , or Ah (in meters):

o e

N =

s 450 404.8 377.2 344. 5 313 252.9 200

= ±

o 5.0 2. 7 1.8 1.2 0.8 0.65 0.6

>1°

o 0.4 0.3 0.25 0.2 0. 17 0.15 0. 14

Assume that error in AR or AR is ± 0.5% or i 0.1 meters, whichever

e

is larger.
N •AN

-32-

Table A-l

•dN

ZOO 22.3317700 0. 118399435 23.6798870 1. 17769275

210 23.6125966 0. 119280212 25.0488444 1. 18991401

220 24.9668845 0. 120458179 26.5007993 1.20315637

230 26.3988468 0. 121916361 28.0407631 1.21752719

240 27.9129385 0.123642065 29.6740955 1.23314913

250 29.5138701 0. 125626129 31.4065323 1.25016295

260 31.2066224 0. 127862319 33.2442030 1.26873080

270 32.9964614 0. 130346887 35.1936594 1.28904048

280 34.8889558 0. 133078254 37.2619112 1.31131073

290 36.8899932 0. 136056720 39.4564487 1.33579768

300 39.0057990 0. 139284287 41.7852861 1.36280330

310 41.2429556 0. 142764507 44.2569972 1.39268608

320 43.6084233 0. 146502381 46.8807620 1.42587494

330 46. 1095611 0. 150504269 49.6664087 1.46288731

340 48.7541501 0. 154777865 52.6244741 1.50435338

350 51.5504184 0. 159332141 55.7662495 1.55104840

360 54.5070651 0. 164177379 59.1038565 1.60393724

370 57.6332884 0. 169325150 62.6503054 1.66423593

380 60.9388149 0. 174788368 66.4195799 1.73349938

390 64.4339281 0. 180581312 70.4267116 1.81374807

400 68.1295015 0. 186719722 74.6878887 1.90765687

410 72.0370324 0. 193220834 79.2205420 2.01884302

420 76. 1686780 0.200103517 84.0434770 2. 15232187

430 80.5372922 0.207388355 89. 1769927 2.31525447

440 85.1564647 0.215097782 94.6430240 2.51823286

450 90.0405683 0.223256247 100.4653113 2.77761532
-33-

Table A-2

AN N dN

o

20 180. 226277 .117626108 21. 1993155 1,,15617524

22 197. 316142 . 118216356 23. 3259953 1 .17457412

24 212. 917967 . 119594076 25. 4637276 1,.19366808

26 227. 270255 . 121491305 27, 6113599 1 ,21348565

28 240. 558398 . 123746115 29. 7681671 1.,23406110

30 252. 929362 . 126255291 31. 9336703 1,,25543336

32 264. 501627 . 128950180 34. 1075325 1. 27764560

34 275. 372099 . 131783550 36. 2895127 1.,30074523

36 285. 621054 . 134721962 38. 4794288 1.,32478398

38 295. 315731 . 137741207 40. 6771452 1, 34981825

40 304. 513148 . 140823306 42. 8825481 1. 37590934

42 313. 261483 . 143955014 45. 0955611 1. 40312414

44 321. 602888 . 147125889 47. 3161107 1, 43153539

46 329. 573439 . 150328075 49. 5441405 1. 46122250

48 337. 204713 . 153555418 51. 7796106 1. 49227226

50 344. 524418 . 156803056 54. 0224815 1. 52477960

52 351. 557000 . 160067149 56. 2727266 1. 55884863

54 358. 324138 . 163344614 58. 5303179 1. 59459364

56 364. 845143 . 166633002 60. 7952415 1. 63214058

58 371. 137293 . 169930326 63. 0674811 1. 67162830

60 377. 216108 . 173234984 65. 3470266 1. 71321044

62 383. 095581 . 176545680 67. 6338699 1. 75705732

64 388. 788373 . 179861358 69. 9280046 1. 80335830
-34-

Table A -2

(Continued)

AN N dN

66 394.305974 . 183181171 72. 2294300 1.85232456

68 399.658845 . 186504431 74. 5381454 1.90419225

70 404.856538 . 189830583 76. 8541525 1.95922635

72 409.907798 . 193159183 79. 1774555 2.01772514

74 414.820650 . 196489873 81. 5080567 2.08002556

76 419.602477 .199822385 83. 8459677 2.14650999

78 424.260086 .203156494 86, 1911914 2.21761358

80 428.799768 .206492043 88. 5437400 2.29383429

82 433.227348 .209828917 90. 9036251 2.37574437

84 437.548229 .213167031 93. 2708570 2.46400458

86 441.767432 .216506335 95. 6454475 2.55938222

88 445.889634 .219846812 98. 0274147 2.66277367

90 449.919193 .223188453 100. 4167688 2.77523207

92 453.860184 .226531281 102, 8135290 2.89800399

94 457.716416 .229875327 105. 2177108 3.03257531

96 461.491458 .233220637 107. 6293319 3.18073184

98 465.188659 .236567271 110. 0484114 3.34463902

100 468.811163 .239915290 112. 4749663 3.52694820
-35-

Table A-3

N •AN dN

1.0 0.0 0.0 0.0 0.0

1.2 217.689023 24.6471681 0.120160519 26.1576259

1.3 275.037959 33.9367000 0. 131692114 36.2203304

1.4 312.297111 41.7747176 0. 143600133 44.8459068

1.5 339.003316 48.4839018 0.154339490 52.3215987

1.6 359.298283 54.2941700 0. 163827653 58.8629945

1.7 375.341242 59.3759008 0. 172203063 64.6349115

1.8 388.391792 63.8586055 0. 179626805 69.7655765

1.9 399.243407 67.8426334 0. 186242834 74.3562234

2.0 408.424907 71.4070090 0. 192172034 78.4878451

2.1 416.304322 74.6148487 0. 197514185 82.2260091

2.2 423.146728 77.5171828 0.202351472 85.6243634

2.3 429.148472 80. 1557288 0.206751820 88.7272277

2.4 434.458411 82.5649192 0.210771674 91.5715264

2.5 439.191718 84.7734613 0.214458304 94. 1883108

2.6 443.438906 86.8054237 0.217851443 96.6038056

2.7 447.272272 88.6811886 0.220984823 98.8403838

2.8 450.750273 90.4181120 0.223887193 100.9172131

00 523.299600 135.5109159 0.299693586 156.829534

4/3 289.036274 36.6922523 0. 135758874 39.2392391
-36-

Tables of Coefficients, Standard Errors of Estimate, and

Correlation Coefficients for use in the Statistical Method

Table I, h - h s = 0. 1km

o r b a S.E.

0.0 0.2665 0.0479 -8.7011 6.7277

1.0 0.2785 0.0257 -4. 1217 3.4363

2.0 0.2881 0.0162 -2.3732 2.0960

5.0 0.3048 0.0073 -0.9181 0.8792

10.0 0. 1915 0.0053 -0.6085 1.0551

20.0 0.2070 0.0025 -0.2639 0.4555

52.4 0.2100 0.0009 -0.0973 0. 1688

100.0 0.2105 0.0005 -0.0507 0.0879

200.0 0.2105 0.0002 -0.0250 0.0435

400.0 0.2107 0.0001 -0.0120 0.0208

900.0 0.2108 0.00004 -0.0040 0.0070

Table II, h - h s = . 2 km

o r b

a S.E.

0.0 0.2849 0.05801 -10.4261 7.5726

1.0 0.2979 0.0348 -5. 6431 4.3330

2.0 0.3104 0.0239 -3.5707 2.8357

5.0 0.3415 0.0117 -1.5287 1.2512

10.0 0.2306 0.0073 -0.7436 1. 1990

20.0 0.2550 0.0035 -0.3184 0.5122

52.4 0. 2604 0.0013 -0. 1162 0. 1890

100.0 0.2610 0.0007 -0.0603 0.0983

200.0 0.2613 0.0003 -0.0299 0.0486

400.0 0.2613 0.0002 -0.0143 0.0233

900.0 0.2604 0.00005 -0.0047 0.0078
-37-

r

Table III, h - h. s =0.5km

Qo b a S.E.

0.0 0.3615 0.0769 -14.6443 7.6170

1.0 0.3997 0.0510 -9.0567 4.4954

2.0 0.4369 0.0384 -6.5408 3.0395

5.0 0.5205 0.0228 -3.6605 1.4376

10.0 0.3933 0.0140 -1.9055 1.2733

20.0 0.4563 0.0071 -0.8926 0.5365

52.4 0.4731 0.0027 -0.3308 0.1966

100.0 0.4753 0.0014 -0.1721 0.1022

200.0 0.4760 0.0007 -0.0851 0.0505

400.0 0.4761 0.0003 -0.0408 0.0242

900.0 0.4764 0.0001 -0.0137 0.0081

Table IV, h-hs = 1.0 km

S.E.

0.0 0.3936 0.0840 -15.1802 7.6151

1.0 0.4620 0.0607 -10.3739 4.5217

2.0 0.5238 0.04918 -8. 2066 3.1040

5.0 0.6348 0.0337 -5.4816 1.5931

10.0 0.5718 0.0224 -3.2378 1.2574

20.0 0.6598 0.0124 -1.6959 0.5531

52.4 0.6823 0.0049 -0.6495 0. 2071

100.0 0.6851 0.0026 -0.3388 0.1080

200.0 0.6859 0.0013 -0.1676 0.0534

400.0 0.6860 0.0006 -0.0803 0.0256

900.0 0.6864 0.0002 -0.0270 0.0086
-38-

Table V, h - hs = 2. Okm

©o r b a S.E.

0.0 0.4524 0.0985 -17.7584 7.5391

1.0 0.5490 0.0752 -12.9451 4.4420

2.0 0.6316 0.0636 -10.7566 3.0277

5.0 0.7707 0.0475 -7.8969 1.5234

10.0 0.7634 0.0345 -5.3712 1.1421

20.0 0.8515 0.02111 -3.1571 0.5086

52.4 0.8668 0.0089 -1.2770 0.2003

100.0 0.8679 0.0047 -0.6705 0.1057

200.0 0.8681 0.0023 -0.3323 0.0524

400.0 0.8682 0.0011 -0.1593 0.0252

900.0 0.8684 0.0004 -0.0535 0.0084

r

Table VI, h - h s -5.0 km

Qo b a S.E.

0.0 0.4962 0.1115 -19.1704 7.5676

JLO 0.6101 0.0881 -14.3543 4.4401

2.0 0.7030 0.0764 -12.1589 3.0001

5.0 0.8504 0.0601 -9.2514 1.4422

10.0 0.8674 0.0464 -6.6445 1.0420

20.0 0.9484 0.0308 -4.0706 0.4028

52.4 0.9674 0.0139 -1.6236 0.1426

100.0 0.9695 0.0075 -0.8348 0.0739

200.0 0.9701 0.0037 -0.4098 0.0365

400.0 0.9702 0.0018 -0.1960 0.0175

900.0 0.9703 0.0006 -0.0658 0.0059
-39

Table VII, h -hs

"= 10. km

e r b a S.E.

0.0 0.5099 0.1149 -18.5627 7.5227

1.0 0.6290 0,0915 -13.7469 4. 3895

2.0 0.7250 0.0799 -11.5514 2.9443

5.0 0.8734 0.0635 -8.6434 1. 3733

10.0 0.8950 0.0498 -6.0729 0.9713

20.0 0.9723 0.0338 -3.5012 0. 3179

52.4 0.9907 0.0157 -1.1441 0.0844

100.0 0.9927 0.0085 -0.5084 0.0406

200.0 0.9931 0.0043 -0.2310 0.0197

400.0 0.9932 0.0020 -0,1078 0.0094

900.0 0.9932 0.0007 -0.0359 0.0032

Table VIII, h - h = 20. km

S.E.

0.0 0.5155 0.1165 -17.9573

1.0 0.6367 0.0931 -13.1413

2.0 0.7336 0.0814 -10.9463

5.0 0.8815 0.0651 -8.0397

10.0 0.9028 0.0514 -5.4747

20.0 0.9785 0.0353 -2.9228

52.4 0.9968 0.0169 -0.6738

100.0 0.9984 0.0093 -0.1802

200.0 0.9986 0.0047 -0.0467

400.0 0.9986 0.0023 -0.0161

900.0 0.9986 0.0008 -0.0048

7.5131

4. 3763

2.9281

1. 3521

0.9573

0.2909

0.0535

0. 0203

0.0096

0.004b

0.0016
•40

Table IX, h - hs = 70 . km

eG r b a S.E.

0.0 0.5174 0.1170 -17.9071 7.5113

1.0 0.6391 0.0936 -13.0912 4.3738

2.0 0.7361 0.0820 -10.8960 2.9251

5.0 0.8837 0.0656 -7.9895 1. 3481

10.0 0.9051 0.0519 -5.4209 0.9539

20.0 0.9797 0.0358 -2.8696 0. 2862

52.4 0.9979 0.0173 -0.6246 0.0445

100.0 0.9997 0.0096 -0.1402 0.0095

200.0 1.0000 0.0048 -0.0212 0.0013

400.0 1.0000 0.0024 -0.0027 0.0002

900.0 1.0000 0.0008 -0.0002 0.0001
-41.

XW<H

U

0) X!Oh

01OS

<*00oIax

d) oo

IN

coa>

i-H s•H>do

•iHuu

4)

CD©

O

CD CD

CD CD O

CD CD

-H

h

CD CD O t- o

CD CD

O CO "itf 00

m o rg F-H o vO "* o O I-H <* 00 r^- rg co r- —4 CO r~ rg O o o o —i

<* oo

m CO ~4 sO co

O o o o o o o i-H CO in nO t^ O o o o o o o o o o o o

•* o <# o o o .—1 o 00 -* 00 ^ o o rg ^ <J^ ^ o i-H r- r- r- i-H 00 00 00 00 oo o CO 00 cO o m o

r-H 1—1 i—<

-H .-H rg rg rg * r- rg o^

vO vO v£> vO vD vO -o nO nO o r^ a^

OJ rg eg rg rg rg rg rg rg rg rg rg

in

O

m 00 o r~ CO r^ o o co rg * O rg <* 00 vD vD o r^ rg m o o O rg ^ 00 i—i i-H >£> ^r nO r^ m o o o o o rg •tf r~ m rg r«- o o o o o o o o o ^H rg rg rg

in r—1

r~ «* r- rg nX3 sO 00 in vO r~- oo i-H 00 r-H vO rg sD o m o^ O 00

CO * tT ^ 00 v£> 00 CO rg r^ o "*

<M rg rg rg rg co •^ r- «* -T rg

m

in in in in in in

m

m sO rO in

-H o^ r- rg o p» o m rt* rg ^H o r- r- m o r- co o »o >o r-4 r- r> ^H o —4 CO r- m nO o •^ m sO vO vO o o O o i-H CO vO rg CO Ol 00 o o o o o o o o "H rg CO co -r o

00o

i-H •—i r- CO o CO vO -V co o

rg * r- in 00 00 CO r^ rg 00

rg <tf 00 •-H ** o 00 rg "f ~H o o o o o rj <* t*- r- ~H rg o

CO CO CO CO CO co CO cO * vO 00

-H

m r- CO in CO i-H 00 o ^ ^ co CO co O •tf ^H 00 r-* o >-H rg rg rg -o rg tf r—1

rg

—1

^H CO co <NJ o v© vD

O

O I-H rg <* o m "* CT^ o vD 00 o o o O o o —H

rg co

mm

m

cN] CO in >* ^r rg

O in <*• rg 00 rr CO vO * in 00 O o ir> 00 nO 'NJ o r-l rg vD rg CO rg rg co m CO > ^r o o O —i rg m o m 00

co in vO

-^ vO CO I—1

INJ o r- rg CO rg co o ~H

m co vO

•"** o vO i-H r- in in rg rr •~* -o

m

rg (3s CO CO i-H o ^r

m o rm rx£> vO o m rg >* ~H ^H 00 o o vO

t—I

rg ^ r-

^H 00 nO sO

—.

rg o ^H

—<

»H rg t*y

m nO vO rO r- r- v£> 00 O tf m rg 00 o •—4 o rg o o r- ^ in ^ o CO rg r—i ^H o m r~ in

t—i

i-H 00

m r~ i—i

CO rg r- ^ 00 •—> O^ IT) f- rg CO vJD 4" co 'VJ rf r- vO

1—1

rg cO in r- •-H n£) CO r~ CO ^o •^

-H •-H rg rr\

m r- ^H in m o o .—<

CO o CO o ^ o o ^H

m -H in o o CO sO -H <—t ^H rg CO CO vO 00 o O^ r» CO f- o sD 00 cO ~H o ~H o> o ^f I—»

00 ~H 00 00 vO r^- co ^

tN] * nD o rg o Is- r- CO ^ ^H -,^ •—i

rg rg CO in v£> r~- r- o nD f- rg vD o 1—1 o CO r^ 00 v£> vO 00 nO —i *

I-H 00 o CO in 00

O o^

mm *-t in

co o in vD o vO

O o N© co 00 ^H vD CO nD —i CO in ^ rg rg 't r^ vO

r—1

CvJ CO

m r- i-H vO CO r^ CO o *

-H i-H rg CO

m r~ *-H

r-H CS] in o o o O o O o o o o O o ^H rg m o o o o o o o O o o o o r—> rg m o o

rg o
-42- 1 r—r

fM "St* s 00

m Is- 00 fM 00 <* <* ro

fM

J- —1 fM

m 1—1

ro Is- i—i

O s- r—i Is- 00 o O o —4

C^

m l—l w—t ro rm fM o o o o o l—l fM M< vO 00 s

|| o o o o o o o O

O o o O

-<* s- ro vO ro -* fM 00 fM o r- o

CD o o fM -* s ro r- in r-- fM s in

CD co co 00 00 co

O fM r-- fM s ro oo

I—1

1—4 i—4

r—4 I—1

fM fM fM M* sO fM oo vO •X) \D N0 ^o O vO vO nO vO Is- s- fM fM fM fM fM fM (M fM fM fM (M fM

<tf r—1

fM ro <# CO s s- vO ro M* •* M<

cu • l- >£> fM

O o s —i o fM s- Is- fM ro

fM o —i

ro vO 1—1 s in ro

m Is- O s m o o O o 1—1 fM in O o s- vO Is0< ii o o O o o O o r—4 fM fM ro ro

COO£

"tf CO s- 00 vO CO Is- 00 00

O i—i

m

<D

CO o t- s- ro •** fM

m o ro vO sO

CD ro ^ -t

m 00 in Is- o 00 fM ro 00

+>

a

A < (M fM fM fM fM ro M< Is- ro "tf fM "^ in in in m m m in m vO I ^ s- s- in

r—4

fM

O

rO

h nO fM co o ro 00 o i—4 •—i ro M< fM vO o r—1

fM in Is- 00 s 00 v£> 00

M4 ^

fM 1—1

fM

m o o o ro r- -* fM vO Isr—4 o o o •—i fM "* s o r—4 rO O fM „ II • • H m a n a • • a t B o o o o o o o o #—i

ro ^ in in 1 CD CD (M ro 00 ^

fM s in 00

O fM CO in

B- * co o 1—I

fM vO o ro 00 o ^ o o fM "tf 00 o s- nO (M in 00 oo

(D • • • R • a o a

a o o o o o fM ro Is- Is- o i—1

00

l—l X s- ro

m ro ro ro ro ro rO <* vO 00 tf » i—i fMID O h r- o "* vO © fM r- in ^ fM s- o fM r—4

ro "* o^ Is- s 00 ^ "* vO ro CO

II fO nO in s vO ro s l—l r-4 —I fM * m . „ o o I—1

fM in fM o ro ro Is- m Is- <: 45 II o o O O o —i fM ro m vO Is- Is- H £ O

CD "tf Is- r- fM ro i—i r- s s- l—l ro fM

<U CD fM * o oo

m ro o 00 fM 00 Is- m

45 1—4 fM nO 1—1 fM

O 00 Is- oo in Is- r—i d o O o r—1

fM

m 00 ^ Is- ro vO vO

•-t fM ro in r~- ^

03 1—1

i—i

ro CO fM in in -H l—l o ro o fM vO

43 p—i h fM o r—4 f—< ro o s l—l s- 00 nO Isrt fM vO 00 tf m in 1—1 s o "tf vO 00

•rH

r-4 fM ro sO o m

m vO s- r—4 in ro

m • - o • o a a B a a >

II

OO o r—4 i—1

fM ro -* Is- 00 s s

•—( s- in fM r—1

vO

O in 00 o <£> vO

d Is- ^ Is- O

m vO xO o o^ in fM i—i

1-4 CD CD co «tf v£» i—i

i-H fM s r- ^ nO i—i 00 i • • a • o a o i o o o • o

1—4 fM ro

m Is- l—l in fM n£) fM vOm

Ki- i-H r—1

(M ro

m Is- -*

rn

i—4 *H

cu vO vO

O ro * <* M< fM 00 s- r—4 vO

ri ro o —4 •—< o ro

m s oo r- vO Is- O h O Is- O Is- s s vO ro m s p—1

ro

<tfm o fM Is- r- 00 fM ro Is- vO 00 • « a a » a a a a • o

II o o O r—4 l—l fM ro

m Is- 00 s s1—1

sO >£> ro l—l fM s ro <tf r—1 o fM CO ro ro O CO fM fM oo oo <tf fM i—4 o m fM in o o fM O vO M^ vO i—I

00 1 CD CD a B • • , u B o • • • * !I 1—4

fM ro in Is- 1—1

l—l

mi—i

fM

fM ro

fMm rm

M<

i—i

I—1

fM

m o o O o o O o o o o O o i—i fM m o o O o o o i • a o • a B a a a o o i s*-' o O o o o o I—i fM in O o o 1 4-> i—i

fM Is- X
-43-

W< u

<o r£

a,

COO

^3

p-in

rooi x o00

(M

0) r£

co

0)

r-H s•H>dO

•H

-t-> CJ

rrj u

MH

0)

rg

CD CD

rg

in

II

O

CD CD

ro

II

O

CD CD O

«D CD

—i r- O

CD CD O h O

CD CD

mi—i oo Orgoo

rO

p-8

in

•—i o

00

rgo

cB

p~o

00

ro

r—1 O

in

rg

rgOrg

in o

inoprgmo

rggoO o o o o o o o o o o o r-H o

00

00o

00

r-H CMCO

roCO

oo

00co

—-1 rgo

00

CM

rg

-HPro

CO

F-H oo

oo go

rg

rg

rg rg rg (M rg

rgrg

rggD

rg g3

(M

go

go

rg

rg

P-rg

00Org

If) Is-oo o•>* t—i o

ro

p-rOO

r\l <*

p-o

i-H pr-H 00

in

ro

in

00 oo

vO

rg

ga

r-oo

p-r-m

rg

rg

r- goo o o O o o O o r-H rg ro ^ ^

rgoo

CO

If) o p- CO

Lf) o

00

roCO sOo

v£) 00

ro CO vO

gs

r-H m

f\1 If) If) if) rg

Lf)

r\l If) rom

m in

ro ro f-H

mr—r)

r-H ro

•—i O

F-H SOo O

If) o Oo(M

r-H

rLf) rgvO

(M

-hvO

i—«

rgm

rg

ino

CO

r-H 00 ooo-H

F-H oooo

ino00

r-H o o o o O o ^h rg rO in m o

COo

00

p-o r—1 ooro

00Is- rrgO

rg00

rg

ro

rg

r- roO

rg

ro

r-H rggo o

ro o oro

O

ro oro

H

ro ro

rO ro o

g3

r-H 00

00

r-H

t—i

ono

if) rp-o

•—i o

i—i

r-H O

ro

00o

vO

00

rorom

*—* •—t O

OroO oo

ro

in

ronO

-H

r-H o

rgrgin

rg o o O O o F-H rg * vO 00 o o

p-

•—t ro

(SI

ro

Pif)

Pr-H -H

ro

r—4

00o om

r—1

rg gao

00org

mrgo o o O •—1

rg <tf 00

rg ro

rgm

go in

r—4

p-p00rgvO

pLf) in

00

P-o

•—> ro

rgpLf) ro

fM oro

inm

rgmrg

in

ro

00

ro

ro

r-H in

roO

ing? o o o »—t

r-H ro -* vO 00 o

*-H

r-H—4

r—i

CO

00

00

ro

if) rif)

Oif) ro

inO

00

rom

f-H ro

00

rg

rgo

P* r-00

in

r-H (M ro rf vO o

f-H inr-H rg

rg

mro

rgin

in

rm

ro

rg

.—i if)

ro

rgpri—i or-H

rgNO

r\J

ino

ro

Lf)

roin

r-00

•—* ro00

in

O

-H

ro

ro

oo org

rg

O r-H r-H rg vO

O

O r-H r-H

<* o rg r~ gD 00 gs 00 in ro r-H ^

ro r- ro

m r~ o o m rg r-H ^ CO

lO r-H <* 00 00 o in i-* 00 o gD in • • o • a • « a r-H rg ro ^ vO o in rg m rg m in

r-H r-H rg ro Lf) p-

—

r-H rg

m o o o o O o o o o o o o i—4 rg m o o o o o o • • * • u o a o o o o o O f-H rg in o

r-H o

rg o
44-

X!W< u

Cv> aX

CO

r-H oI

ft o

CO

t—I

oo

II

0) -d

[0CD r—I

•tHu>do

•rHo

Sh m

<U

eg

CD

egin

CD o

CD CD O H

00

O

CD <D

CD CD

—i r- o

CD CD O t- O

CD CDX

Oooo

co oo o o o o

r-H (\J

egO

no co —i m o o

ooo

eO in m o r-H f\J IT)

^t< o o Tf O ^

CO 00 IT) 00 o —< o o o o o o o o o o r-H r-H o

CO

00o

00 o

eg

00 00

CO00

00

i—i

r-H O

oo

og

og

r- gD

gD

r-H

og

r-H in

rr—1

oggj

g?

eg CsJ On] eg

Os] og

og

og

og

og gD

og

gD

og

og

r~

og

00

og

oooo

r-H r-

-ho

CvJ o

r—1 m

ooo

oo

r-H

ogo—<

CO

og

co

CO

CO

—(00

gDvO gD

ogCO

CO

00 o o o o o o o r—i og CO ^ *

1—1

oo

COo Is- r-m

oo

00

-h

cO om

g3 rr—( m

00oogD

egin00m

m

Oxl in CO

in

COm

m

gD

in

COgDCO r-H

mi-H o

in

i—i o o(MO

mo o

00

•—<

i—»

rsj O

mogo

ogOog

CO

r—<

CO

O

in

CO

CO

m

og

CO00

00m

gD oCOoo o O o o O o .—i

og *#

m gD gD

COo

in

r-o

gD

oo

•—i

Osl roo

r-H

fCO00

CO

i—i

in

ogo

CO

CO

m

CO o m o

CO oco o

CO oCO o

oo

r-H CO

CO

CO CO vO

m

r-H oo

00i—i o o0000o

COco

.-H eg

(N)

00mo

00

ogo o

vOoo o

00

CO 00

COo

r-og

00

r-oo

r-H

CO

CO00

CO o o O o o i-H CO 't r- O or-H Or-H

r-H

r-H

i—l

r—1

CsJ eg

mmmo

ego o

-sO gD

gD

00

coo

r-om

og

m o o o r—<

eg <tf 00 COog ooogin in

m

r-H

1—i

CO

00

COm oooo

eg

oo

oom

og

OnI oog

i-H googo

CO

CO

Omo

CO

00

og

CO

og

00

rg

r-O

CO

COO

r-H

O o o —* og co in r- or-H og

r-H og

r-H COi—( oo

00

r-H

CO

<*o<*

eg

00

CO00

r- o

g^

r—

H

og o

CO

CO

mmO

CO

r-H og

r-H eg CO * vO o—4

m og in

CO

m

in

m

•—*

i-H o oO

inoo or—<

oo

Osl O00

mcox>

m

r-H

r-H

gDovO

CO

og

00

in

00

CO

ogCOm

O —i

gD r~ —I -^

—I CO

r-H og

co

O

l> o mo

nO o o r-H og

ogingDom

in

co

in

coo

.-H og oo oO

gD

r-H COm

COoo

og

r-moo

—I og

CO

ogo

co

—i in

in Is- oo o o r-i og

oo

i—i mooo

l>
45. >i—i Xw<H

0)

COO&<

XiCOoi

a,X

CO m

ro

CDm

CO I—I 3

•i-t >

C!o

•Hu

Cm*

O

© CD

rgin

O

CD CD

rO

CD CD O h o

<-h

H

O

CD CD O rCD CDX ooo oo 8

i—i o oo

rgo Is-oroO

vOo

00

inin

ro oo oc—1 o rg

00s vO

rg o o o o o o o O o •—1 -^ f—t mo

00 Is-o

oo

r—4

00

00

ro00 Isr-00

in 00o—\

ro

rg 00

I—1 o

00

in

roino

r—1

in

IM rg CM rg NO

rg rg

rg

rg

rg

rg rg rg

rgIsrg

00s- rg

COo1—I o O

rgo

CO

r—4 mo

*—4 rgo

r—1

1—(

rgorg

O

m

ro o

rg

rrg

m

mrgIsrg

roIs-—* o o o O 3 o O f-H ro ^ m in o

00

ro en

00

in

r-mm

in

rgm

ro

rgro o

rg

<*

in

rg

rOs- Is-

—(

ro s rg

rgin in

rgin

rgm

rom in in

rg rg s- m o

00

t—1 o s- m3 o

ooo m

>>

r- omro o

rg

oo o

00m o

rgO

00 oi-H

•-H

•-* rOmo

vO o

insOIs- o o o o O o r-H rg m vO r- r- <*mo o .—1 Is-~H

rg

ro 00

i—< o

rornroro

mm

r- sOo IsrgvO o

in

rg om

CD

<-o o

OO

O

ro o

ro ro ro

ro ro

m

m o

00

oo

i—4 Isroino

00 NOo

f—t

rgrg o1—i m

O

rg

rro

ro

*-*

m ro

rgos- 00

00

ino o o O O o rg ro m 00 o i—4

i—1 rg

-H

O

r—1

1*o

rg

roO

in

m

00omo00

sO

ro

in rg00

rg o

o^

00

in

00

m—4

in

sOm o O o o *—* •* Is- m

rg

m

m

mIsin

—4

rg

rgvO

vO

r—4

ro

rg

r-oo00

roin

00

rogD o

vO ooo

in rg

•-H in

rg

in

m

rg o C i— —'

rg rf o o rg * m~H

in

»—4 Isro

rg

rogD

ro

mm

in

cr- rgm

vO omo Om

vOo00

m

oo

00org

-h rg ro <* vO o

rg ro om Isin

i—i

rgmIsrgmo

rg00Isrg

ro

00Is- rg

ro

rgo o

ro

in

m

r~ rgO

rgsOsin

»-h rO

_"5 r—1

-* rg ro in r~ o rg

1—«

in

m<—4 sO

»—4

in

CO omo

i-H rg

m

rg

m

rg

ro

rg vDrgO

ro

inm

0000 s-

<cf< 00

morg i—t

rg ro -* n£)

Or—t rg ro om Isin

i-H «—»

rgo

ino o

r—1 O

rg o

in oo oo Oo oo oo oo o o o O o o ^-1 rg m o o

rg oIs-
46-

u

<D A

Oh

toO

+-> g

fM

co

fXX >

<u

X (M

r~

w rhh m

CQ ii < _^ H rCl £

d) A

coD

i—i

crj >o

•H

-t-> o

Sh

cu

fM

II

CD CD

fM*

IT)

CDO

CO

CD CD

CD o

CT> cd

rl

h

CD°CD O

CD CDX

M<oo oo

i—i

fM

~Ho

fMfMo or-o

•—I —1

CO

fM

fM

in

•-H

i—l fM00

m—I omCO

fM

CO o o o o o o O o o .—i >—i

^h

COo

oo o

oo

i—i

CO

CO

00 00

m r—1

i—t unm

fMo

o^

fM

CT^

m

in

M<

fM fM fM fM fM

fMvO

fM

fMvO

fM

COvO

fM vO

fM .—ifM

oo

fM

in

fM

i-H o O

»MO

O

fMvOOCO

fM

i—i

fMM<

fM

•—i

fMm

00mfM

i—i

fMoo

fM

r-00

COo

in

fM

fM

CO fM

inO o O O O O o —<

fM CO in in sO

co CO

in

COm

fM

i—i rfMm

fM

CO

—1 o fM

fM

co

mfM

M1 00o

fM

IT) (Min

fMm

fMin

fMm

COm

m in

fM fM o

in

—I Isi—i (MO

T+4 COM<O

i-H cooi—i

i—i fM

CO

fM

fM~Ho mo

CO

•—A o

co00m

oo

1—1 m

r>- vO

M*

fMO

00

O o o OO !—<

—I CO v£> r~ 00 oo

i—i

COO

I—l o

coin

O

CO

fM

•—1 vO

CO

COin

fMo

M<

•—* vOo

fM

fM

fMoo

-H

COo o

CO o

CO oCO

O

CO

O

CO —1CO

CO

CO CO

m com o

oo

co

i—i Omo

fMO(M

r—1

i—i

oo

i-H CO

COo

fM o

00

i-H

CO

IM

•—I

rfMOfMO

00o

inCO

COo om

M<o ooo o O o O i-H fM •^ r- oi—i

ro

r—i CO

1—i i-H

O

fMoo

1—1

COoo

00

00o o i-H

i—l

O

in

COo m

fM

oo

COo

CO

fM

CO

O o o o —I CO r- fM

fM

m

CO

m

m

r- in

M<

i—i

in

00

m

in

fMoo

m

COmm o1—I o

COo

00

fM om

co

CO00

fMCO

'M00

(M

—I

00vO

i—i fMO

O

CO o o i—i

fM CO

m 00 •—I in r- oo 00

COCO

i—i

i—i

(MOr—1

COo

i-H

r-4

00ino

r—< Oi—1

vO

co

—I CO CO

.—1 fM CO -tf vO

O co

i-H o

fM COCO om f—1

r—1

O^

co

fM

CO

ino

in

i-H oo

coo

i—l o

in

fM oo

rCOo

fM

i—i O

fM

ri—i

i—i

00m

CO o

fM

fM fM

<*

O i—<

fM CO M< sO o fM sO 00 (T* o

IT) co

mi—i

—I

COo

COofM COofM

i-H reg

fM

1—1 CO 00

CO

,—i

O

CO

1—1 00

^H i—i

CO ^ vO o CO

i—i o

fM COCO om i-H

i—l o

fMO

mo or—1 o

fM o

in oo O oo oo oo oo o O o o o o 1—1 fM m o ofM o
47- >Xw<

a)

Oh

(0On

4300o00oI O

fM

O

CD CD O

a> cd

ro

II

CD CD O h

—> II 4£z

CO 43

TO

<D #—i

•s

•H>CoH»-> CJ m

<D

CD CD O

CD CD O

CD CD

O

moo

IN

i-H o

ro00

fMO

1—1 vOmo

ro

COo

vO

fM

in

fMO

CO

in

rfM

vOO o o o O o o O o O ^H 1—«

^-i

fMO

00

mo

00

i-H co o

roCO o

00 in fM ro

CT- fMro

r-4 00

in

ro

fM (N N0(N fM IN IN

IN

fM

fMvO

fM

ro

IN

vO

IN

1—1

rfM

00

fM

1—1 o

INOrgo

fMo r—1 o

vO

00

fM

fM

fM

rOfMro ro

m 00o00 o

in

Oo—4 o o o O O O »-H fM -* in vO nO

in

r-m

rO

M^ in

r-l roin

i—i OO

00

<*mr-o

rom o

<Nin

fMin

fMmfMm

IN

in

roin

roin

mm

—4 IN oo

roin

in

ino oo

ino

00sO

(N

CO—<

in

fM

moi-H

00fM

IN

00

•—< m

m

oo

00

ro ro o o o o o i-H IN ro o 00 o crr-<NO

ino

CO

r—1

<N o in

rO 1>-

in

-*1 00

ino

00 oo

d

ro

d

rO o

ro

O o

ro

—

Hro

(N

ro inro

<*<*

m i—i

fMo

•—< m•—<

00

rOrro

ino

CO

00oro

ro

i—i

r-ofM

ro ro

i—«

ro vDro00O

m

ro

ino o o o o i—i

ro in 00 fM—< in

m o

00o o

i-H oro o

CO i—i m

00

inin oinm

m

ro

r- fM ro rO o o O O i—i

ro vJD fM

<*ro om

m

rOm

00m

CO

1—1 oo

in

rooo

in00m00

rITI

i—i ro

rorO

fMooro

O

fM

00

00

O

<*

fM

fMO

fMIs00

ro

COOO o 1—1

-* fM *tf rO

r—t ro 00

—1 IN (N

fM

fM in

i—i mo o

CO

<*

mo

00 oofM ro

fM fMro

fM inort-H IN m rf in O rO fM

ro

CT- 4" ro

i-H r-o

<n o

ri—i

CO

I—» o

ro

fMro

<NvO

IN

fM *fM

ro

INm oroO

•—1 00

00o

IN

i—ii—i fM ro

m 00 i—i o^i—i

fM fM

fM fMfM 00 O

00

fMO

00

i—«

rm

in—i

i—i

00

(N

m

rO

CTfM

ri—i i—i fM "5f m o ro

—1

fMro CO

O

fMO

ino oi-H o

IN om

Oo oo Oo oo oo oo o O o o o o —H

IN

m o—H o

(N o
-48-

u

0)

a,

COoI<A

egco

(NJ egOi > X X <D O W in J ^

CQ

n <! ^-» H

-C£

0) ^3

CO

<u i—I

•i-i u>

•rH

4-> O

n$

Jh VM

<0

r- in vO i—1

CO CO h- o eg o- r~- o

co r- 00 r- CO rm * O .—i

oo rO o 1—1

co r- r- CO o eg o^ * go

t*- o o o o o —> CO sO 1—t ^ gD g?

t- o o o o o o o o 1—1 -H ^H ,—i eg eg ^ o -* CO m o m 1-1 in ON CO o o 1—1 eg <# i-M "* vO O i—i CO -H

II oo 00 00 00 00 o^ o

eg

CO

eg

in

CO

i—i g3 eg

oo nD vO vO vO vO vO O NO vO vO r- ON

CD CD INI eg eg eg eg eg eg eg eg eg eg eg

1—1 eg IN] CO

O r- r- ON 00 —1 eg oo o 00 in On •* i—i

00 1—1 <* gD

m vO

1—1

CO o oo r- o on eg in r- eg r~ o o o i—i

CO on vO

O <tf ON in

m ^. t- o o o o o o -H CO in sO r^ r»

fVl

IT) 1—1 (Nl *# o CO r^ i—l o o 1—1 o m » oo —1

r*- oo ^ esj 00 v£) 1—1

m i—i

II co CO <# ^ in ON vO 1—1 oo eg o sO

eg' eg CsJ eg eg evj CO in o —1 o co in in in m in in m in vO r» on in CD CD F—1

-«* r- r—t o -* i—i o on ^^ CO —4 co

en vO sO o o 00 "#

—i

vO ON CO

m

m vO nO eg in in o sO ^H ^ o NO

O

CO

h

O o ^-1 CO sO in o^ o r- r- •* * o o o O o i-H eg m 00 o r^ —-I

II ,—i

•-H —1 o oo in i—i ON ^ eg i—i

eg i—i

00 eg

i—i CO o o 00 i—l

m o^ sO <* 00 *

CD o o o r-i CO o —< vO r—1

00 1—1

m CD o o o o o o 1—1

eg ^ CO sO ON t~- CO CO CO CO CO CO CO CO ^ m r~- c 1 O CO "* eg 00 r- eg eg 00 CO m ON o o fM sO eg ^ ON 00 "* O CO CO

O

h o o O vO vO o <tf -* CO

m 00 ^ •—i i-H 1—1 ^ On co eg eg CO 00 co o r—1 o o o O ~H -<* r- 1—1

vO O^ o o

II 1—<

F-^ <-H eg eg N0 co —1 O eg CO O 00 r- o ^H m o m i—1

00 vO r—t eg o^ 00 o I—l r~- eg CD CD o r—1

(NJ

m

—'

r- eg o -o CO ON 00 o O o o 1—1 eg m oeg

eg O^ CO .—1 * O o r- eg 00 O o vC- "* 00 vO m r- 00 o in vO 00 m co eg -£> [N- r—t CnJ 00 o o ^ eg cO eg CO r~ CO

H CO ^ r- CO vO o ^ *f ^ o r- 00 O i—<

eg -v vO 1—1 in o v£- ON o^ ON

II 1—1 —I eg eg eg eg eg o 00 o 00 00 o o ^H CO r~ CO m O

CD sO o 00 CO <* CO i—i

inm ON o 00

CD

•* CO

m in ON o^ ^o ^ o eg CO ^

1—11—1

eg co -* r~ 1—1 ,—i

i—i

CO

00 CO

1—1 ^ eg o o in co 00 m o oo CO (—1

in * o r^ eg •—I ^o gD vO o in go

CO vO o o r- 00 00 o o i—i in r-*

O

h 00 vO eg on eg vO o ^^ 1—1 r- ^ m

i—i (M •*

m 00 eg r~- eg 00 o ,_; -i

II

i—< .—i

eg eg CO CO CO

•sh g3 o <*

m r- r- eg r- r- o eg vO O 00 ON <# vO vO eg co 00 o^ 00 o o m CO CO 00 00 m "tf o eg eg •* CD CD o „ o „ „ —i —<

eg CO * r~ i—i in- i—i 1—1

CO oo CO ^-*

i—i £

•—I

eg

m

O o o o o o o o o x o O o i—i

eg

m o o o o o o o • • o • 4 •

+->° o O o o o O i-H eg

m o o o 1—1 eg r- E
-49-

TABLE XVIII

Initial N Gradients, AN , in the C.R.P.L.

Exponential Reference Atmosphere

Range of AN (N units/km) N(h)

200 exp(-0. 1184h)

252.9 exp(-0

.

1262h)

289 exp(-0. 1257h)

313 exp(-0. 1428h)

344.5 exp(-0. 1568h)

377. 2 exp(-0. I732h)

404.9 exp(-0. 1898h)

450 exp(-0. 2232h)

* Note height, h, is in kilometers.

AN < 27. 55

27. 55 < AN < 35. 33

35. 33 < AN

e < 47. 13

42. 13 < AN

€ < 49. 52

49. 52 < AN < 59. 68

59. 68 < AN < 71. 10

71. 10 < AN

e < 88. 65

88, 65 < AN

€

<.
-50-

iDfOO^r^Or^O^ffMu-ifV]

^rOLn^OOOOOOOO

+

CD+

cd"

CD

fM

i—

i

O +

II

<

CD

CD vD U O«—

«

o X

fi M

A

X

t-i —

i

t-H oj 4. X 3 ^

w

Xi

-1 (NJ <

o H a) ,* a, £

a fM

Rj W r—

1

T5 +

o y,

,£3 x\

-«-> >—** 0) r\J 2

CO v ti 32

i—

<

sJ 4dow

fM

^5

fM

O Xl

-"loo^o^^^i^inoh-vo roLnr\jir>sOco^r^-oooooo o-^oofOfMOsooo^oo

NiflNOOOOONNOOO^O v£>moO'-H mr^coaO'-*\£>oo

vOfMmo^o^r-cO'-^(MOfvj

mvO'f^iOfOO^inocO'-* [^-fMsOonsOmoOLnfMrOOrO

^•vOO^fMoo^r^oo-^r^-

r

s-ir>fy">sOfv">00r>-fM h- O n£> f"> sO 1—< fM —*\rl \r\ot~-—*\r~oo

-—I ^ vO -^ fM fM fM •— nO ft

omrO'-ir^fM^fvDor-unvX)

OOOOOOOOO^vOO

O —« -"^ fM —i O O —« O ^ fM ooooooooor-oo

ooofo^unLnLOvDr^r^-a^O

r^r^r-r-r^-t^-r^r-r>-r~-r--(--oo rOr'lflflflfir'lflr'lClMflrO

sO vO vO vO O sO vO O vD v£) \£) sD xO

U

6

CO

fM

-*" fM

+

Kl

w

^5 CJ oomoinooooouioo

omror— sOrnfMLnfMr^moom OvOrooOor-r^-tntnonrMOOO ^ rO on fM —• —i

.(^-<r-ir-4—i

-~M ooooooooot^oo

OrOOOflOf^OfM—H^OvOOO

oof^-^Lnu-iLnor^-r^-cr-o

^ I O^fMrOT^iLrisDr^COOO—'fM
-51-

X u

cdx

CD

rsj

t^-rMcoo^cOvor-cOvOcoo

rMrvi-^cosoo^o^ —< in c\j

CO C\J <sf< —iOOOOOOOO

r^-OvOui^fMcoc^rNjco-^o ^ o o s r- •-•r^oo'MvOOxD

ooox>r*-oor~<MON ^-, oO'-<o

o co —< o co m r-- oo —iror^'-'

s

COo

o

X!X

W

m

<

ao

< CDuo

co

•r-l +->

as" a

aoU

a) r—

<

a

£

CO

T3X

4-

CD

(M

4-X

CD

sO

Xi

4-

r\J

(M

X

X

4-

,—

V

.—

.

1—

I

(15 ^ 4- +J £

^ a sr* ?

£

4- x J*!g

XI *

fM i

—i Nnr--r\]mvOOOON c\jT^ioco

sOm-^mcor-^m —• —'vOO oOsOcomTfOr^vOcxj^rNjco

in-| OOoNoo>^1 —' in oo ^

r^-inoOvOonoor«-f\]r-OvOoo

O —i (M —

i

O s vO o —• -i sD ^

m iTi O

on vO O

p- —> r~- o ocr^ o^ co rf r~-

rsl —i vD co

OrOrO—'r^-iNl^xDNOr^-LOvO

r- v£> O CO *f co ro — m <NJ OOOOOOOOO^sDO OONMvOOMOONvDf^OON vOfMMiniriTf'NJNOOOooO'^*

K3

w o ^ (M —i O o ^r <Ni

(fl

•rH

i—l X

o

CO

rt 6

ooooooooor-oo ^in^OTfoo-^moo —• vO r- oooro^iriLnirixDr-r-oo

r^r^-r-r^r^-r^-r-r-r-r^r^r-oo

X c OOmOmOOOOOmOO

Oirionr— vOf^i'MirifM^'inoOiri ONOcoonor^r^minoorMoao ^fOfON-i-i-i-H-H-'-

J4

o o o o

-tf m v£) -^f O CO O O CO O O co -^

cr* O O O O r» O O

o 7f in oo -< sO nco s c\) —i vD vO 00

in m vO f— r- cr o

o —'cvjco-sfirivor-oooo—• c\j
-52-

vOOM HH vOOO\OfnvD^OO

o o ooooooooo

X

w

CO

CD

om

C!o

• iH +>

4-»

a,aou

cu I—

I

a

to.o

CD+ I cd"

M

CD

a (M

CD o

o

I

H

V

f\J

.—

1

en + +> M ri Z 3

i Z

£<^0 ?1

mroooo'-^orMfMvOfMmo

fMrOOOO(M(M'<^-<^vOr^--H^

ino^o-^fMo^ooor^o

ro *-i o •—toroo-Hr^-r-rvjvD

00-—i M O 00 in iM -1 sO CO h- N

r~unpOsOt>ooor^<^ r^-OsDoo

v£> ^h (M --h m m o

O O vO O oo vD o —* ~* O ^t1 M fM

r*- —i r- o o

o o (^o ^ rN vX) (> CO Tt* oo rn —i m N

OOOOOOOOO^fOO OONNvDONOOM^iNOON vONNinin^NvOMCOO't

in

II

J4

^ fM O h O h o ^ M

COd• H

J*

u

CO

-* -^a ooooooooor-oo

OOOOO-^ U-iiDuDsDr^r-OO

sD vD ^D sO nO ^£> nO vO vO ^D vO vO vO O O in O in OOOOOmOO

omr^r^ororomfM^mooin OvOflciONi^ininfONoco

•* a ooooooooot^-oo ^ in vO ^ oo^moo-HsDt^ O oO O O CO CD CO O (M •—< O vO 00 ooooTf Lnminvnr^r^-oo

OHiMflTfUlvOI^OOOOHN
53-

5 ~

- s

mrOvOmOOO-NO'-^vOmO^ r-n^-irO'-HOOOOOO'^O

OOOOOOOOOOOO

X

X

W

m

u8

00

CDU

om

d

O

•i-i +-> rd

+-> d

BoO

O

CD^

CD

r\J

CD

sO

I

r\]

+ U

CD S

i i

(0 M z

z

M(Mnfl^ON^Iv-'-| 00N

omoooinvO<M',3

< i^minr^

r^'-^0'-< t^os ^^ooNoos OLr) oO'-HOvOoo-^tnorgcno

—irvJcororo-^T^^munvOr^

vO O nO vO vO vO vO nO nO nO vD vD oooor^or^^ooo^iNm omaoinoor^inr^^mfviin

m tji ^ ^ m ^ fl

h- o r- o oo ^ —i m r- —' -^ o Is- o oo oo s o o o o

vO ^O nD vO vD vD vO

in h oo (> n oo n m ^ vO <*i oo vO H n h m m

o o vO o m vO

r-H r-l ^O ^JH f\]

o m m —tr-iNJ^NOor^invO

OOOOOOOOOrfxOO 00<M(Nv0O<NJ00<N)v£>rs-00<NJ vONNinm^NvOXooo-t

^ o o v£) o

o O oo O ^H o <NJ PO m

o o o O —

r^ r^ t^ r^ r^

r» rt< — in 00 m

f- fNJ r^ o vO ro o r- -h r- o o^ o o o oo ^ r^

^J (M ^H vD rn

s

O H ^ (\J H O O —i O <^ fv]

CO "d

^5

r—

d

o

CO

ooooooooor^oo ^•mvO'^'OO-^inoo^Hsor^

roos omorndNrvJ'^NOv00o

ooooo^ininmvOi^r^-oo

t^t^r^r^-t^r^^-r-r^r^r-t^-oo

vDvDvOvOOOsOvOvOsOv£)sOvO

4*5 d

2 d

oomomoooooinoo

Oinfi^^^NmN^incoin OvDfnonor-^r^LnmroiNJooo

-^ fO M N -H _« ^ _,^h_h^h

^ ooooooooor-oo ^mv0^oN o^mcoH Vr)> oonc-omooooiN^OvOoo

oooo-^minmvor^t^o^o

o—'fMfn-^insOt^oooo^fM
-54-

TABLE XXIII

Departures Method Sample Computation for 9 =0 mr,

0.340

0.950

3.060

4.340

5.090

5.300

5.940

6.250

7. 180

7.617

9.660

10.870

N(h)

400.0

365.0

333.5

237.0

196.5

173.0

172.0

155.0

152.0

134.0

125.5

98.0

85.0

A(N,h)

400.0

397.8

419.5

459.7

475.7

478.5

484

485

490

493

493

495.9

495.2

AA

- 2.2

+21.7

+40.0

+ 16.0

+ 2.8

5.6

1. 5 4.9

2.9

0. 1 2.6

0.7

9+9

k k+1

6.388

11.230

22.683

28. 195

31.383

32.

34.

35.

39.

40.

47.

50.

191

654

811

123

581

236

645

0.3131

0. 1135

0.0590

0.0393

0.0336

0.0315

0299

0284

0267

0251

0228

0204

ideparture i"(h)

term

+0.689

-2.463

-2.359

-0.629

-0.094

-0. 176

-0.045

-0. 139

-0.077

+0. 003

-0.059

+0.014

N 1,2

30.776 25.441

5. 335

TABLE XXIV

De partures Method Sample Computation for 9 = lOmr.

400.0 400.0 10. 000

o

0.340 365.0 397.8 - 2.2 11.873 0.0914 +0.201

0.950 333. 5 419.5 +21.7 15.038 0.0743 -1.613

3.060 237.0 459.7 +40.0 24.789 0. 0502 -2.008

4.340 196.5 475.7 + 16 .0 29.947 0. 0365 -0. 585

5.090 173.0 478. 5 + 2.8 32.938 0.0318 -0.089

5.300 172.0 484. 1 + 5.6 33.712 0.0300 -0. 168

5.940 155.0 485.6 + 1.5 36.069 0.0287 -0. 043

6.250 152.0 490. 5 + 4.9 37. 181 0.0273 -0. 134

7. 180 134.0 493.4 + 2.9 40.383 0. 0258 -0.075

7.617 125. 5 493.3 - 0. 1 41.801 0.0243 +0.002

9.660 98.0 495.9 + 2.6 48.282 0.0222 -0.058

10.870 85.0 495.2 - 0.7 51.630 0. 0200 +0.014 19.

I

= -4.556

19.414 14.858
-55-

TABLE XXV

Departures Method Sample Computation for 9 = 52.4mr (3 )

N(h) A(No> h) AA 9

Q + Q

departure T(h) N

t

k k+1 term s

s

400.0 400.0 52.4

0.340 365.0 397.8 - 2.2 52.750 0.0190 +0.042

0.950 333.5 419.5 +21.7 53.550 0.0188 -0.408

3.060 237.0 459.7 +40.0 57.061 0.0181 -0.724

4.340 196.5 475.7 +16.0 59.567 0.0171 -0.274

5.090 173.0 478.5 +2.8 61.045 0.0166 -0.046

5.300 172.0 484.1 + 5.6 61.477 0.0163 -0.091

5.940 155.0 485.6 + 1.5 62.792 0.0161 -0.024

6.250 152.0 490.5 +4.9 63.433 0.0158 -0.078

7.180 134.0 493.4 +2.9 65.368 0.0155 -0.045

7.617 125.5 493.3 - 0.1 66.277 0.0152 +0.002

9.660 98.0 495.9 +2.6 70.512 0.0146 -0.038

10.870 85.0 495.2 - 0.7 72.920 0.0139 +0.010 7.024 5.350

I

= -1.674

TABLE XXVI

Departures Method Sample Computation for 9 =261.8 mr (15")

F ci

o

400.0 400.0 261.8

0.340 365.0 397.8 - 2.2 261.880 0.0038 +0.008

0.950 333.5 419.5 +21.7 262.035 0.0038 -0.083

3.060 237.0 459.7 +40.0 262.758 0.0038 -0.152

4.340 196.5 475.7 +16.0 263.308 0.0038 -0.061

5.090 173.0 478.5 +2.8 263.633 0.0038 -0.011

5.300 172.0 484.1 + 5.6 263.733 0.0038 -0.021

5.940 155.0 485.6 +1.5 264.034 0.0038 -0.006

6.250 152.0 490.5 +4.9 264.184 0.0038 -0.019

7.180 134.0 493.4 +2.9 264.647 0.0038 -0.011

7.617 125.5 493.3 - 0.1 264.872 0.0038 +0.000

9.660 98.0 495.9 +2.6 265.934 0.0038 -0.010

10.870 85.0 495.2 - 0.7 266.592 0.0038 +0.003 1.506 1.143

1

0.363
HIT' NT OF CO^

H. Hodges, Secretary

INA1 BUREAU OF STANDARDS

\ \ . Vstin -rector

THE NATIONAL BUREAU OF STANDARDS

rtie scope of activities of the National Bureau of Standards at its major laboratories in Washington, O.C., and

Boulder, Colorado, is suggested in the following listing of the divisions and sections engaged in technical work.

In general, each section carries out specialized research, development, and engineering in the field indicated by

its title. A brief description of the activities, and of the resultant publications, appears on the inside of the

front cover.

WASHINGTON, D.C.

Electricity. Resistance and Reactance. Electrochemistry. Electrical Instruments. Magnetic Measurements.

Dielectrics. Metrology. Photometry and Colorimetry. Refractometry. Photographic Research. Length. Engineering Metrology.

Mass and Scale. Volumetry and Densimetry.

Heat. Temperature Physics. Heat Measurements. Cryogenic Physics. Equation of State. Statistical Physics.

Radiation Physics. X-ray. Radioactivity. Radiation Theory. High Energy Radiation. Radiological Equipment.

Nucleonic Instrumentation. Neutron Physics.

Analytical and Inorganic Chemistry. Pure Substances. Spectrochemistry. Solution Chemistry. Standard Refer- ence Materials. Applied Analytical Research.

Mechanics. Sound. Pressure and Vacuum. Fluid Mechanics. Engineering Mechanics. Rheology. Combustion

Controls. Organic and Fibrous Materials. Rubber. Textiles. Paper. Leather. Testing and Specifications. Polymer Structure. Plastics. Dental Research.

Metallurgy. Thermal Metallurgy. Chemical Metallurgy. Mechanical Metallurgy. Corrosion. Metal Physics. Electrolysis and Metal Deposition.

Mineral Products. Engineering Ceramics. Glass. Refractories. Enameled Metals. Crystal Growth. Physical

Properties. Constitution and Microstructure.

Fire Research. Mechanical Systems. Organic Building Materials. Inorganic Building Materials. Computation. Statistical Engineering. Mathematical Physics. OpDigital Circuitry. Digital Systems. Analog Systems.Building Research. Structural Engineering.

Codes and Safety Standards. Heat Transfer. Applied Mathematics. Numerical Analysis,

erations Research. Data Processing Systems. Components and Techniques.

Applications Engineering.

Atomic Physics. Spectroscopy. Infrared Spectroscopy. Solid State Physics. Electron Physics. Atomic Physics.

Instrumentation. Engineering Electronics. Electron Devices. Electronic Instrumentation. Mechanical Instru- ments. Basic Instrumentation. Physical Chemistry. Thermochemistry. Surface Chemistry. Organic Chemistry. Molecular Spectroscopy. Molecular Kinetics. Mass Spectrometry.

Office of Weights and Measures. BOULDER, COLO.

Cryogenic Engineering. Cryogenic Equipment. Cryogenic Processes. Properties of Materials. Cryogenic Technical Services. Ionosphere Research and Propagation. Low Frequency and Very Low Frequency Research. Ionosphere Research.

Prediction Services. Sun-Earth Relationships. Field Engineering. Radio Warning Services.

Radio Propagation Engineering. Data Reduction Instrumentation. Radio Noise. Tropospheric Measurements.

Tropcspheric Analysis. Propagation-Terrain Effects. Radio-Meteorology. Lower Atmosphere Physics.

Radio Standards. High Frequency Electrical Standards. Radio Broadcast Service. Radio and Microwave Materials. Atomic Frequency and Time Interval Standards. Electronic Calibration Center. Millimeter-Wave Research. rowave Circuit Standards. Radio Systems. High Frequency and Very High Frequency Research. Modulation Research. Antenna Research.

Navigation Systems.

Upper Atmosphere and Space Physics. Upper Atmosphere and Plasma Physics. Ionosphere and Exosphere

5 utter. Airglow and Aurora. Ionospheric Radio Astronomy.