Ooer1967

~A7

"MP06A

an
0o

nee

4r

0f
1025

0=ZBCSM.01BRTR

WHITE SANDS MISSILE RANGE, NCE

EU~

IC

.....
....

KROKSCMMN
I~f

MEXICO

A SIX-DEGREE-OF-FREEDOM DIGITAL COMPUTER PROGRAM

FOR TRAJECTORY SIMULATION
by
LOUIS D. DUNCAN
and
BERNARD F. ENGEBOS

October 1967

ECOM-5158
DA TASK 1V0145O1R;3A-1O

ATMIOSPHERIC SCIENCES LABORATORY

WHITE SANDS MISSILE RANGE, NEW VEXIC0

Distribution of this
report is unlimited.

ABSTRACT
A documentation of a six-degree-of-freedom model for digital
simulation of tlh- trajectory of an unguided, fin-stabilized, windsensitive rocket is presented. This model was developed by the Atmospheric Sciences Laboratory, White Sands Missile Range, New Mexico,
to study both theoretical and empirical performance characteristics
of unguided rockets.
The basic equations of motion and their mathematical formulation for this model are presented without derivation.
A general flow chart, a listing of the program, a list of the
principal flads used, and a listing of a typical irput data deck
are incl-led.

iii

CONTENTS
Page
ABSTRACT---------------------------------------------

iii

INTRODUCTION-------------------------------------------

COORDINATE SYSTEMS AND TRANSFORMATIONS ------------------

TWE MAIN PROGRM LRBol-----------------------------------3

INITIALIZATION -----------------------------------------

MHE TWVO
SUBMONITORS ------------------------------------

11

THlE EVAL ROUTINE --------------------------------------

16

THlE INPUT SYSTEM --------------------------------------

23

THlE OUTPUT SYST31I-------------------------------------

31

THE INTEGRATION ROUTINES -------------------------------

36

LISTING OF PROGRAM INSTRUCTIONS------------------------42

INTRODUCTION

In recent years many theories and computer programs have been

developed to simulate missile trajectories. These models range from
the extremely simple to the very complex, the defrce of complexity
usually depending upon the specific simulation requirement placed upon
the developer.
This report describes a six-degree-of-freedom computer program
developed by the Atmospheric Sciences Laboratory, White Sands Hissile
Range, New Mexico, to simulate the trajectory of an unguided, finstabilized, wind-sensitive rocket. The simulation theory upon which
the program is based is presented in reference 1. The references
may be consulted for derivations of the equations and coordinate transformations; however, the coordinate systems used and the principal
equations which are evaluated by the program are included for completeness.
The system is programmed in Fortran IV language. It consists
of a main program (monitor), two subprograms which serve as submonitors for specific simulation options, and a group of subroutines,
each designed for a specific task.
An attempt has been made in this report to present a reasonably
complete program documentation without boring the reader with trivia.
With this goal in mind most of the routines are documented in four
parts: (1) a statement of the purpose of the routine and the equations to be evaluated; (2) definition of the principal flads; (3)
a macro subroutine; and (4) a listing of the program instructions.
It will be observed that in some of the minor routines (2) and/or
(3) have been deleted.
COORDINATE SYSTEMS AND TRANSFORMATIONS
Three right-hand coordinate systems are used in the program.
These are: (1) The launcher system (denoted X', Y', Z') which has
its origin at the launcher and rotates with the earth. The positive
X'-axis points east; the positive Y'-axis points north; and the positive Z'-axis along the outward normal of the earth. The X'-Y' plane
is tangent to the earth at the launcher. (2) The inertial system
(X, Y, Z) has its origin at the center of the earth. The system is
oriented so that the X-Y plane lies in the earth's equatorial plane

Preceding Page Blank

with the positive Y-axis pointing initially through the longitude of

the launcher. The Z-ais is coincident with the earth's axis of rotation and is positive toward the North Pole. This system ooes not
rotate with the earth. (3) The body system (x, y, z) has its origin
at the center of gravity of the rocket. The x-axis coincides with the
longitudinal axis of the rocket and is positive toward the nose.
The position of the y- and z-axes is determined by the rocket's motion. The initial positions of these axes are defined as follows.
Let 0 be the angle between the x-axis and the positive Z'-axis measured from the Z'-axis. The y-axis lies in the X'-Y' plane and is
positive in the direction of positive 6. It is easy to see that the
angle between the positive Y'- and y-axes is the launch azimuth measured clockwise frui, the Y'-axis.
The transformation between any pair of these coordinate systems
depends upon, among other things, the assumed earth model. The program considers the earth to be an oblate spheroid with an equatorial
radius of 20,926,428 feet and an eccentricity of .00672267. The transformation from the launcher system to the inertial system is easy to
obtain from the geometry of the problem. This transformation is, in
matrix form,

f.Cos

xl
X

sin Ag sin wt

wt

iig

- sin wt

- sin A

cos

cos wt

Co

g sin W

cos X

cos w

[xi
Y

'I

sin

where A is the geodetic latitude of the launcher, w is the earth's

rotatioR speed and t is time after launch. The derivation of the
preceding transformation, except for the assumption of a spherical
earth, is presented in reference 1.
The transformation from the body system to the inertial system
is obtained by integrating the derivatives of the elements of the
transformation matrix -- the direction cosines. This transformation
is denoted by

Y
L

n,

x
Y

m2

n2

t3

m3

n 3

It i, shown in [1] that the required derivatives are

m.

= pni

n.1

rm.

qn.

rZi
Pmi

i = 1, 2, 3

14i

where p, q, r are the x, y, z components of the rotation of the body

system with respect to the inertial system.
It is shown in reference 1 that the initial conditions for these
direction cosines are

-1

0
00

m.,

n2

n3

sin

sin a
0

cos a
0
-sin a
0

12

sin X

cos X

sin e0 cos a
0
0

Cos X

sin X

cos 0

sin a cos
0

-sin e

where 00 is the launcher elevation angle and a 0 is the launcher azimuth

angle measured clockwise from north.

TiL !'MIN PROGivAM LRBM
LRBM is used to monitor the other routines, to initialize the
trajectory simulation, and to monitor the type of trajectory desired.
INITIALIZATION
(this is the

elevation of the x-axis angle and is measured from the :'-axis); and

3,

cos a cos n
0
0

Initially the x-axis and Z'-axis form an angle e

the projection of the x-axis in the X'Y' plane forms an angle a with
the Yl-axis (this is the azimuth angle, of the x-axis and is measured
clockwise from the Y'-axis). Since the y-axis lies initially in the
X'Y' plane, it forms an angle a + 90 with the Y'-axis. (See Figure
1). Thus initially,

sine cosa

sine sin

-sina

y Cosa

LI

-sno

Cosa cose

sina.cosO
0

-sine

IY1
IZ'
I00II&

and
=

Y=

Y' sinX

Z =

Y' cosA

x -Z U

+ Z' cosA

+ Z' sinX

wX

- t. 2 U + wY

Zu
3 0
3 Uo0
t

where A is the geodetic latitude of the launcher, U is initial x

0
g
component of velocity in body system, and w is the angular rate of

the earth's rotation.

The initial conditions for Zi. mis Ili, i a 1, 2, 3, are given

by

Z6Cs
3

L12

3J

in sia
in 0 si 0

0i~cs 0I
-sina

] cosa

00

00

s
sina cose

Cosa cosO

cose

-1

-sinA

IL

-sine

0
cOSA
g

-cOSA

sinA

FIGURF 1
IUIAL POSITiOtI OF THE BODY SYSTEM

NTYPE is an indicator, re ad in the input routine, that detemines

the type of trajectory to be simulated (1, a regular trajectory; 2,
an anguiar displacement table; 3, iterative type trajectory; 4, the
balAistic factox option; or 5, a parameter variability type trajectory). LRB:4 us., NTYPE to monitor the program and ottains the desired
type of trajectory simulation.
The equations of motion of an unguided rocket are numerically
integrated by a Runge-Kutta integration technique. This technique
is discussed in the section, the Integration Routines. LRBM sets
up the initial integration interval for each phase and monitors the
integration henceforth.
Several ar-ays used by all routines are defined. They are the
Y, DY, and ROFF arrays. All three arrays have double subsc.ipts.
The second subscript, J, denotes the step in the Runge-Kutta integration. The Y array then has the following definitions:
Y(l,J) - time of the trajectory simulation.

Y(2,J) - x component of velocity vector in body s).teri.

Y(3,J)

Y(4,J)

-z

y component of velocity vector in body system.

component

f ve> city vecto- in body system.

Y(5,J) - x component of the rotation of body system with respect to the inertial system.
Y(6,J)

y component of the rotation of the body system with

respect to the inertial system.

Y(7,J) - z component of the rotation of 'he body system with

respect to the inertial system.
Y(8,J), Y(9,J), Y(lO,J) - direction cosines of the body X-axis
with respect to the inertial axes.
Y(l,J), Y(12,J), Y(1l.,J) - direction cosines of the body
Y-axis with respect to the inertial axes.

Y(14,J), Y(15,J), Y(16,J)

direction cosint of the body

Z-axis with respect to the inertial axes.

Y(17,J)

component of the inertial position.

Y(18,J) - Y

component of the inertial position.

Y(19,J) - Z

component of the inertial position.

Each DY(I,J), I equals 2 through 19, is the derivative of the

corresponding Y(I,J) entry.
Each ROFFI(I,J), I equals 2 through 19, is the round-off error
associated with the corresponding Y(IJ) entry due to the Runge-Kutta
integration. The principal flads are:
AL

- Azimuth angle of rocket on launcher.

BIll

- Current height for cutoff of wind profile for ballistic

factor option.

CLAT

- Cosine of geocentric latitude of launcher.

CLATG

- Cosine of geodetic latitude of launcher.

CNTR

- Number of the previous trajectory.

DELC

- Cross unit wind effect (m/mph)

DELH

- Head unit wind effect (m/mph)

DELT

Tail unit wind effect (i/mph)

DONE

Indicates if integration should be continued (1) or

not (2).

DYSTOR(22)

Storage of DY array for future use.

ENF

- Indicates whether the winds in the subroutine, AOPTUN,

have been iterated for (1) or n3t (0).

EPSQ

- Square of the eccentricity of the earth.

GCONI, GNU

- Constants used in evaluation of gravity for geodetic

earth.

- Integration interval (sec),

liP(SO)

- Heights in wind table (ft).

I!,I2,13,14

- Last entry used in atr-,pheric, wind, Mach, and time

tables.

IBFEND

- Indicates if all wind profiles of ballistic factor option have been simulated (1) or not (0).

ITOUT

- Indicates print out only at end of each phase (1).

Ili

- Indicates whether to call ballistic factor option (2)

or not (1).
IOUT

- Indicates print out (2) or nit (1), end of phase (3),

or impact (4).

IWIND

- Entry of wind currently being used in ballistic factor

option.

X and Y entry in table for angular displacement table.

IXWD, IY1WD

3l

- Phase currently beirn

J2

- Phase at which to pick up trajectory for ballistic factor

option.

NBST

- Phase end of which booster drops off.

NFAZE

- Total nmber of phases for trajectory.

NOGOOD

- Indicates whether integration is accepted (1) or not

used.

(2).
NTYPE

- Indicates whether a regular trajectory (1), iterative

type trajectory (2), an angular displacement table
(3), the ballistic factor option (4), or a parameter
variability trajectory (5) is to be simulated.

NUF

PI

- Print interval (sec).

Indicates if the height of the simulated trajectory

is less than or equal to height of impact area (1)
or not (0).

PSL

- Atmospheric pressure at sea level (lbs).

RANGE

- Distance to desired impact point. (ft).

REQ

- Radius of earth at equator (ft).

RO

- Magnitude of radius vector of oblate earth along radius vector to rocket,

ROA

- Sum of RO and hcight above sea level of impact area.

ROFSTR(22)

- Storage of ROFF

SLAT

- Sine of geocentric latitude of launcher.

SLATG

- Sine of geodetic latitude of launcher.

TBO

- Time, end of current phase.

TI1

- Tilt angle of rocket on launcher.

TIME

- Current time of simulated trajectory (sec).

TIMEO

- Time rocket leaves launcher (sec).

TO

- Time of end of last phase.

TOhTIL

- Tower tilt effect (ft/ra.ian).

UO)VO,WO

- Initial body coordinates of aerodynamic velocity.

- Angular rate of rotation of earth.

IVR

- Theoretical range displacement due to wind

nXS,

IVYS

array for future use.

(ft).

- E-W and N-S components of wind in launcher coordinate

system (mph).

XDIF, YDIF

- Theoretical E-V and N-S components of wind displacement of rocket (ft).

XLAT

- re!ocentric latitude of launcher.

XLATG

- Geodetic latitude of launcher.

XNO

- Total number of trajectories.

LREM

IWITMIZE
1

FOR, ITERATION POR

IRSTIT AMOXEMATION
LLANCHER
FPoICHICIL
FIRST
SEM14GS
WEDID)
(17(17

ccenoz
INITIALIZING

MDNT.TOR
INTEGRATION

MMDNTTOI&t
oUrPvr

91(D
OF PRASZ
YSS

PICK UP
SECT PFASR
11pust

YES

BOOSTER
TO GRD

TORS Y AND
ARRAYS

NO

Nf

WACX?
's
SOOST
TO GROUND
YES
RK ALL T
AND DY

NO

LAST
ES
mmim

XUWA, YUW

- Cross and range Uniit wind effects (m/mph).

XWANT, YWANT - E-W and N-S components of desired impact point (ft).
ZL

- Height

f impact area above sea level (ft).

THE TWO SUBMONITORS

For convenience of the user, techniques are available for simulating a series of trajectori.es from one computer submission. The
first of these techniques simply consists of submitting a number of
consecutive data decks - one for each trajectory. Two other techniques are available for computing a series of trajectories of a specialized nature. These computations are controlled by the submonitors, BALFAC and AOPTUN.
SUBROUTINE BALFAC (INTA)
This subroutine is used to monitor the calculations of ballistic
factors (wind weighting factors) and unit wind effects. A complete
trajectory is simulated for each entry in the wind table (see section
on Input). Each card in the wind table has an entry of wind and height.
The trajectory corresponding to a given entry in the wind table is
computed for a wind profile having the wind value of the entry when
the simulated altitude is below the specified height and a value of
zero when the simulated altitude is above this height. The first entry in the wind table should have zero winds, i.e.,a no-wind trajectory is required.
The X and Y components of impact for each entry are subtracted
from the no-wind impact point to obtain a range displacement. The
range displacement of each trajectory is then divided by the range
displacement of the final trajectory to obtain a ratio or percentage
of the total range displacement. Each ratio is then subtrac'ted from
the previous ratio to give the Lallistic factor for that layer.
The unit wind effect is obtained by dividing the last range displacement by the magnitude of the wind.
flow.

Three indicators are used by the subroutine to control the logic

T!,cse are INTA, IWIND and IBFEND.

112

f.

The routine is entered after the completion of cach integration

step.
INTA determines the cntry point into thc subroutine. During
tie simulation of a given trajectory the value of INTA is 1. ',Oien
impact is detected, 1,y LR['!, INTA is set equal to 2 and
UALFAC
is
entered to store tLC iimpact point and cstal'lish the initial conditions
for the next simulation.
IMIND is the indicator to deterrinc which wind table entry to use
for a given trajectory.
IIi-lNLD is an indicator whica tells the routine if the simulations
have been corplcted. This determinatioi is made in I.,C'.
If the
sin.ulations have been completed, ITFEXD is set equal to 1 and
UALFAC
is entered to c.nIpute the ballistic factors and unit v.'ind effect.
S[11VOIJTINL AOPTUN (K)
The simulation model contains a technique for deternining the
launcher settings necessary to compensate for a given wind profile
so that the computed impact will be in a desired area. The launcher
settings arc determined by an iterative procedure uhich is discussed
in reference 3.
This subroutine is used to monitor the iteration procedure. In
addition it is used to monitor the computation of a nomogram of launcher
angles versus ballistic winds. (See reference 4Q) These latter calculations are referred to as angular displacement trajectories.
K is an indicator.used to determine the entry point to the subroutiWe and has admissible values of 1 and 3. When K is equal to
1 the initial approximation of the launcher setting is made and the
wind profile is set up. The subroutine is entered with a value of
K = 3 at impact and a check is made to determine if the simulated
impact is within a required tolerance of the desired impact. If it
is not, a new set of launcher angles is determined and the iteration
is continued. If the tolerance is met and an iterative trajectory
is being computed, a return is generated. If a nomogram is being
generated, the launcher setting is stored in tabular form and the
next wind profile is set up unless the nomogram has been corpleted;
whereupon a return is generated.
The first approximation for the
launcher settings for this profile is made and the new trajectory
is begun. The principal flads for the two submonitors are:
Al.

- Azimuth angle of roclet.

ALPHA (11,11)- Storage of azimuth angles in tabular form for angular

displacement option.

12

EF(50)

- Ballistic factors of the rocket.

BFACT

- Ballistic factor.

ILEI

- Height of current entry in wind table (feet).

DLLC

- Cross unit wind effect (m/mph).

DELhI

- Head unit wind effect (r/mph).

DELR

- Ballistic wind displacement of current trajectory

(miles).

DELT

D)IFX

- L-1-, ballistic wind displacement of current trajectory

(miles).

DIFY

Tail unit wind effect (m/mph).

N-S ballistic wind displacement of current trajectory

(miles).

- Indicator, set equal to one when angular displacement

table finished.
iP(50)

- Heights used in the wind table (feet).

liT

- Height of rocl,ct above oblate earth (feet).

IBFEND

- Indicates whether or not the entire wind table has

been used.

111

- Indicates whether to use wind or not.

MIVIND

- Entry in wind table currently being used.

MiF

- Number of entries in ballistic factor table.

'TYPE

Indicates type of trajectory to be simulated.

Nl,

Number of entries in wind table.

PA\N(W

Sxne as DELR.

RATIO

Percentage of wind effect.

13

BALFAC(ITj

BALLISTIC FACTORS

II~'I

ARLL

ITR

UNI

AREF

YES

PODFI1E

WIN'?

EFEC

AOPTUN (K)

IS K:1,3?

YE

IS 27NYP1.

ISTBLEESRTR

S NTPE--2211IS

IPAC

SLMX

- E-W component of ballistic wind.

SUMlY

- N-S component of ballistic wind.

TESTVA

- Range difference between simulated and desired

impact.

R11

- Tilt angle of rocket.

TIETA (11,11)

TOIVrL

- Tower tilt effect (ft/rad).

UNIT

- Unit wind effect (m/mph).

WIND

- Magnitude of wind (mph).

Storage of tilt angles in tabular form for angular

displacement option.

wxP(5o), WYP(SO) - Entries of E-W and N-S components of wind.

WXS - WYS

- Value of E-W and N-S components of wind.

XDIF - YDIF

- E-W and N-S components of difference between simulated and desired impact.

XCI ANG - YC!ANG

Increments of E-W and N-S component of wind for

angular displacement option.

XFIRST

Initial and terminal values of E-W compenent of

wind for angular displacement option.

XLAST

XtJW - YUW

- Cross and range component of unit wind effect.

XhANT - YWANT

- E-W and N-S components of desired impact.

XYZ(l)

XYZ(2)

E-W and N-S components of simulated impact.

THE EVAL ROUTINE

The subroutine EVAL is used to evaluate the rocket equations

of motion and to compute the derivatives of the direction cosine matrix which specifies the transformation from the body coordinate system to the earth-centered inertial system. The derivation of the
equations can be found in reference 2.

16

Let
represent the velocity of the rocket and let i represent
the rotation of the x, y, z system with respect to the X, Y, Z system. The x, y, z components of V and
are designated by u, v, w and
p, q, r, respectively.
The equations of motion

derived in (1] are

mu

m(rv - qw) + F_

mv

m(pw - ru) + F

mw

m(qu-I

pv) + F

Iyyq

(Izz

zzr

(13

)qr - I

+ L

yy q + M

Ixx) pr -

-I yy) pq - Izzr + N

where m is the rocket mass, I

Iy, i are the moments of inertia,

XXI -zz
Fx, F and F are the body components of the applied forces and L, M,
y
z
N are the body components of the moments. The forces and moments are
computed from the equations

-C

-C

y
;

*~g

qS
I
x

Fx

sin

mgx
q' S -mg

ny

Ct

[Cm

na

sin a*q' S -mg

6 + C1
sin a* + Cm
a

m(C sin 8
aq

q
Cm

(-v-']q'
a

Sd

( 2vJ

17

q '

m(l-cg)2 q

~-g2r

where
q

1/2 Pv2 is the dynamic pressure,

S is the reference area,

d is the reference length, and
va is the magnitude of the aerodynamic velocity.
m(l-cg)q and m(l-cg)r are jet damping terms.
The thrust, T, is computed from the
formula
T - Ts.l. * Ae (Ps.I. " Pa)
where
T

is sea level thrust

Ps.

is sea level pressure

Pa is atmospheric pressure
and

A is the area of the exit nozzle.

e

The body components of the acceleration due to gravity are designated gxP gy, and g." These are obtained from the inertial components of gravity which are assumed to be

gx

gR

gz

R2 v(1 - 3 cos 2 X
eq

y Me X
R5
3

1-

gy Me Y
-

1 - 3

(1 - 3 cos 2 X)/R

2/R
2

.5eq
e
R-

[l

3 R2

(1 -

eq

18

cos 2 X)/R2

where y

'-

2
1.40775 x 1016 ft 3/sec

g
v=
A

.273 x 10-'

is geodetic latitude.

Wind enters the equation in the terms sin a, sin a* and v a

precise formulation is presented in [1]. The wind components and the

atmospheric pressure and atmospheric density, p, are obtained from
a table using the subroutine TABL. The wind components in the body
system are obtained by applying the transformations described in the
section "Coordinate Systems and Transformations."

The stability derivatives C, C

, Ct,
C , C and C , which
cx
a
)p a
q
are defined in the section on Input, are also obtained by a table
lookup. The routine has an option of either obtaining Cm from the

table or computing this value from

cMCa

-cn 1C -cPI

The derivatives of the direction cosines are computed from the

formulas
1,rM.
n
i = 1,2,3
i.
i. - qni
,,

m. = pmi - rli

i = 1,2,3

qI.

i = 1,2,3

n. =

pm.

The inertial component of the accelerations are obtained by applying

the body-to-inertial transformation.
The only input parameter to the routine which is not transmitted through Common is the integer J which appears in the calling sequence. This parameter indicates the step in the Runge-Kutta integration routine at which EVAL 4as called. The principal flads are:
AE

- Exit area of motor nozzle (sq. ft.).

CG

Center of gravity (feet from nose).

19

CHEK

CLA7

- Cosine of geocentric latitude.

CMA

- Center of pressure (feet from nose).

CMQ

- Damping moment coefficient.

CNA

- Normal force coefficient.

CP

- Center of pressure (feet from nose).

CX

- Axial force coefficient.

DAMP

- Jet damping term.

DMIKOMK (3)

- Derivtives of moments of inertia.

DSQ

- Reference area of aerodynamics.

EPSQ

- Square of the eccentricity of the earth.

FLOW

- Mass fl,,: rate,

PM

FX

- Moment of inertia abotit x-axis.

FY

- Moment of inertia about y-axis.

FZ

- Moment of inertia about z-axis.

- Gravity.

GCONI
CNUJ

- Constants used in evaluation of gravity for geodetic

earth.

IIRO

- Ma,,itude of radius vector to rncket.

1IT

- Height of rocket above oblate earth.

II!

- Indicates whether to call ballistic factor option (2)

or not (1).

IWIND

- Entry of wind currently being used in ballistic factor

option.

Magnitude of radius vector to rocket (feet).

ass.

20

if burning (1) or coasting (2) phase.

KAN

.-Indicates

MAL

NTYPE

- Indicates what .ype of trajectory being simulated.

- x-component of the rotation of the body system with

respect to the inertial system.

PA

- Atmospheric pressure.

PSL

- Atmospheric pressure at sea -evel.

- y-component of rotation of body system with respect

inertial system.

QP

- Dynamic pressure.

- z-component of rotation of body system with respect

to inertial system.

REFL

- Reference length for aerodynamics.

REQ

- Radius of earth at equator.

RHO

- Atmospheric density.

RM

- Mach number.

RS

- Radius of oblate spheroid along radius vector to rocket.

SA

- Sine a*.

SB

- Sine S.

SLAT

- Sine of geocentric latitude of launcher.

SLATG

- Sine of gewdetic latitude of launcher.

TIMEO

- Time rocket leaves end of launcher.

TSL

- Thrust at sea level.

UP, VP, WP

- Body coordinate system components of aerodynamic velocity.

VA

- Atmospheric velocity.

Indicates if CP (1) or CM (2) is input.

21

EVAL (J)

NORMALIZE
DIRECTION
COSINES

INTERPOLATE FOR
WINDS AND COMPUTE
COMPONENTS IN
BODY SYSTEM

COMPLrrATION
L

OF

CHER TO ECI

LROTATION MATRIX

COMPUTE TOTAL
FOCE|

ARODYNAIC

KND MOMENT

COEFFICIENTS

-- ----COMPIUTE
ACCELERATIONS

IN BODY SYSTEM

22

VS

Speed of sound.

Angular rotation of the earth.

1XOT (1-3)

- Body components of p, q and r.

WDOT (4-6)

- Body components of p, q, r.

IVX, 1WY, WZ

- Components of wind in inertial coordinate system.

WXP (50), WYP(50)- Components of wind in launcher coordinate system.

XLT

Length of rocket (feet).

THE INPUT SYSTE01

The data required for a simulation is read into the computer

at the beginning of the program. These data are stored in memory
for later use. The data required for a given phase (a specific portion of the trajectory) are extracted from this storage area and placed
into the working storage area at the beginning of the phase.
SUBROUTINE

Pintrj

This subroutine is used to read in all input data and to store

the data phasewise. All card formats are in columns of eight with
a few minor exceptions. The first six columns of all data cards can
be used -or identification purposes. The first card of the input
deck, containing the number of distinct trajectories in the job, is
read in the monitor routine, LRB11, in columns 17-24.
Description of the Input Deck

The second input card determines the type of trajectory and an

indicator to determine if the regular printout is desired or not.
If column 9 contains a:
1 - a regular trajectory is calculated.
2 - an iteration for launcher angles is run.
3 -an

angular displacement table is computed.

the ballistic factors and unit wind cffect are calculated.

a dispersion analysis is run.

23

If also on the second card, column 17 contains a 1, printouts

occur only at the end of each phase.
The input for cards 3 through 6 is as follows:

Card
3

Column

Parameter

Description (Units)

1-72

IDEN

Alphanumeric information

used for identification.

4

9-16

TO

Time (sec) rocket leaves

end of launcher.

17-24

ZL

Mean sea level altitude at

end of launcher (ft).

25-32

ZIK

Hean sea level altitude of

impact area (ft).

33-40

XLAT

Geocentric latitude of launcher

(degrees).

41-48

XLONG

Geocentric longitude of launcher

(radians).

49-50

NFAZE

Total number of phases (right

justified).

57-58

NBST

Phase, end of which, to pick

up booster (right justified).

9-16

UO

x component of velocity in
body system at end of launcher
(ft/sec).

17-24

VO

y component of velocity in
body system at end of launcher
(ft/sec).

25-32

11O

z component of velocity in
body system at end of launcher
(ft/sec).

33

Indicator; if equals 1, input CP; if equals 2, input

CA.

MAL

24

Card
6

Column

Parameter

Description (Units)

9-16

Tilt angle (degrees).

17-24

AL

Azimuth angle (degrees).

25-32

XWANT

East-west component of desired impact (ft).

33-40

YWANT

North-south component of
desired impact (ft).

41-48

DELTA

Iteration tolerance.

49-54

TOWTIL

T,,er tilt effect (ft/radian).

The remainder of the input consists mainly of tables. All tables

use one card for each value of the independent parameter (height,
Mach number, or time).
The first
card of each table contains a right
justified integer in columns 7 and 8 which indicates the number of
cards in the table.
There are five types of tables:
1.
Column

Atmospheric Table

Parameter

Description (Units)

HH

Height above sea level (ft).

17-24

RHO

Density (slugs/ft 3).

25-32

VS

Speed of sound (ft/sec).

33-40

PA

Pressure (lbs/sq ft).

9-16

2.
9-16

Wind Table

liP

Height above sea level (ft).

17-24

WXP

East-west component of wind (mph).

25-32

WYP

North-south component of wind (mph).

25

3.

Mach Table

Parameter

Description (Units)

9-16

PM

Mach number.

17-24

CX

Axial force coefficient.

25-32

CNA

Normal force coefficient (per radian).

33-40

CMA

Restoring moment coefficient or center

of pressure (ft from nose).

41-48

OIQ

Pitch damping moment coefficient.

49-56

CLD

Roll driving coefficient.

56-64

CLP

Roll damping coefficient.

Column

4. Time Table
9-16

Time csec,

17-24

FX

Roll moment of inertia (slug/ft 2).

25-32

FY

Pitch moment of inertia (slug/ft 2).

33-40

TSL

Thrust (lbs).

41-48

Cc

Center of gravity (ft from nose).

49-56

FM

Mass (slugs).
S.

9-16

BF

Ballistic Factor Table

Ballistic factor.

For each phase there are two Phase Cards containing in:
1-8

TRO

Time phase terminates (sec).

9-16

AE

Area of motor exit nozzle (sq ft ).

17-24

DSQ

Reference area used in aerodynamic calculations (sq ft).

25-32

REF.

Reference length used in aerodynamic calculations (ft).

26

Column

Parameter

Description (Units)

33-40

XLENT1I

Length of rocket (ft)6

41-48

EPTINY

49-56

EPBIG)

Integration tolerances for Runge-KuttaGill numerical integration.

57-64

FINC

65-72

T-II

Maximum integration interval (sec).

73-80

P1

Print interval (sec).

1-8

TFP

Time of first print (sec).

KAN

Indicator; if equals 1, input Mach and

time tables; if equals 2, Mach table
only.

Depending on the type of trajectory desired, several other cards

may be necessary. These cards include:
Unit Wind Effects Card
Parameter

Description (Units)

9-16

DELC

Cross unit wind effect (m/mph).

17-24

DELT

Tail unit wind effect (m/mph).

25-32

DELH

Head unit wind effect (m/mph).

Column

Angular Displacement Wind Card

9-16

XFIRST

The first E-W wind used (mph).

17-24

XLAST

The last E-W wind used (mph).

25-32

XCHAYG

Increment for E-W wind change (mph).

33-40

YFIRST

The first N-S wind used (mph).

41-48

YLAST

The last N-S wind used (mph).

49-56

YCHANG

Increment for N-S wind change (mph).

27

The input for a regular trajectory is as follows:

Cards 1-6
Atmospheric Table
Wind Table
Repeated for each phase as
needed

Phase Cards
Mach Table
Time Table (if KAN is 1)

The input for the interative type run is as follows:

Cards 1-6
Atmospheric Table
Wind Table
ballistic Factor Table
Unit Wind Effect Card

Phase Cards

Repeated for each phase as

needed

Mlach Table
Time Table (if KAN is 1)
The input for the angular displacement table is as follows:
Cards 1-6
Atmospheric Table
Angular Displa enent Wind Card
Unit Wind Effects Card
Repeated for each phase as
needed

Phase Cards
Mach Table
Time Table (if KAN is 1

28

rI
The input for the ballistic factor option is the same as for the
regular trajectory. The first card of the wind table should contain
a zero wind. The remainder of the cards of the wind table should
contain winds of constant magnitude and direction.
The input for the dispersion analysis is the same as for the
regular trajectory except after cach 'Iacl. and time talle, another
card follows. This card contains a percentage of change in decimal
form plus one in the same field as the parameter to be changed. That
is, if 1.1 appears in columns 17-24 on the card following the :1ach
table, the axial force coefficient is increased ten percent. If the
parameter is not to be changed, a 1. must appear in the same field
as the parameter itself.
If one also desires to bring the booster to impact, another phase
must be added. This phase follows tle last regular phase and consists
of:
Phase Cards

Repeated as required.

Mach Table
Time Table (if KAN is 1)
On card 2, the value of NFAZE, should include the booster phases.
SUBROUTINE

Phasin (J)

This subroutine selects the data required for the Jth phase and
places these data in working storage areas. Recall that all data
input is at the beginning of the trajectory. The principal flads
for the input system:
AE

- Area of exit nozzle (sq ft).

AER(10)

- Storage of AE for all phases.

CG(40)

- Entries of center of gravity in time tables.

CB(40,10)

- Storage of CG for all phases.

CMA(20)

- Entries of center of pressure in Mach table.4

CMAB(20,10)

- Storage of CtA for all phases.

6q(20)

- Entries of coefficient of damping moment inMach table.

CMB(20,10)

--Storage of CMQ for all phases.

29

Entries of coefficient of normal force in

CNA(20)

ach table.

CNAB(20,10)

- Storage of CNA for all phases.

CX(20)

- Lntries of coefficient of axial force in Mach table.

CXB(20,10)

- Storage of CX for all phases.

DSQ

- Reference area used in aerodynamic calculations.

DSQB(10)

- Storage nf DSQ for all phases.

EPBIG3
LPTINY)

- Integration tolerances for Runge-Kutta-Gill numerical

integration.

EPBIGB(10)j
'PTINB (10 j

- Storage of EPBIG and EPTINY for all phases.

FM(40)

- Entries of mass in time tables.

FMB(40,10)

- Storage of FM for all phases.

FX(40)

- Entries of pitch moment of inertia about x-axis in

time table.

FXB(40,10"

- Storage of

rX for all phases.

Entries of pitch moment of inertia about y-axi5

time table.

FY(40j

in

FYB(40,10)

Storage of FY for all phases.

IOPITlB(10)

Storage of IOPTUN for all phrses.

KAN

Indicator, if equal to two, only MaL> table is necessary.

KANB(10)

- Storage of KAN for all phases.

Ni

- Number of values in Mach table.

NTIB(O)

- Storage of NM for all phases.

NT

- Number of value3 in tine table.

NTB(1O)

- Storage of NT for all phases.

PI

- Print interval.

30

PIB(lO)

- Storage of PI for all phases.

PM(20)

- Entries of Mach number in Mach table.

PMB(20,10)

- Storage of PM for all phases.

T(40)

- Entries of time in time table.

TB(40,10)

- Storage of T for all phases.

TBO

- Time at end of phase.

TBOB(10)

- Storage of TBO for all phases.

TFP

- Time of first print to be used in phase.

TFPB(lO)

- Storage of TFP for all phases.

TMII

- Maximum integration interval.

ThIIB(1O)

Storage of maximum integration interval.

TSL(40)

Entries of thrust in time tables.

TSLB(40,10) - Storage of TSL for all phases.

XLT(IO)

- Entries of length of rocket.

XLNTh

- Length of the rocket.

TIHE OUTPUT SYSTM

The output from the simulation is a printed record of the simulation at discrete preselected times. The output is referenced to
the launcher coordinate system. Three subroutines are used in the
output system. These are: (1) CHKOUT is used to determine when output is required. (2) TI2L transforms the data from the Inertial
System to the Launcher System. (3) XOUT prepares the magnetic tape
for the printef. XOUT and CIKOUT are called by LRBM; TI2L is called
by XOUT.
SUBROUTINE CHKOUT
Output from the trajectory is required at the following times:
(1) specified intervals within a phase; (2) end of each phase; and,
(3) impact. The subroutine CHIKOUT is called at the end of each integration interval to see if any
of these conditions have been met.

31

N1JP

IS TIMEzTF~o

32

An indicator IOUT is set equal to 2, 3, or 4, respectively, if condition (1), (2) or (3) is satisfied and set equal to 1 otherwise.
The logic flow of LRBM is, of cuurse, affected by the value of IOUT.
In performing the checks it is assumed that impact will not occur
before a simulated time of 20 seconds.
The integration interval is also controlled by CIIKOUT. This
control is limited to control necessary to preclude "jumping over"
(1) or (2).
SUBROUTINE TI2L (KWIND, JF).
This subroutine is used to convert values of position, velocity,
and acceleration from the earth-centered inertial coordinate system
to the launcher coordinate system for use in the subroutine, XOUT.
The total aerodynamic velocity and the rzimuth and elevation angles
of the velocity vector are also calculated.
The position, velocity, and acceleration in the launcher system
are obtained from the following equations:

RL

T(R1 -1ix0)

t
RL

"t
T(R I -

x (1ixi21)

2 xRL)

whereL is the position vector in the launcher system,

4I is the position vector in the inextial system,
4.

W iF the earth's rotation vector,

T is the linear transformation from the inc-.ial to the

launcher coordinate system, and
4.

is the radius vector from the inertial to the launcher

0system.

KWIND is an indicator used in the calling sequence t, determine

if the accelerations should be computed in the launcher system or

33

IT2L (KWIN9,JF

MARIX FROM ECI

TO LAU~NCHER

FCOMPUTE

TRANSLATION
VECTOR

0HUTE POSITION
IN LAUNCHER
SYSTEM

OKPUTE RELATIVE
MOTIONS

AND ACLERATION

SYST2

(RETURN

34

M4

not. Only when KIND equals one are these accelerations computed.
The principal flads of the output system are:
All

Azimuth angle of velocity vector (deg).

B(l)-B(2)-B(3)

Launcher coordinates of aerodynamic acceleration

(ft/sec/sec).

CHEK

- Magnitude of radius vector to rocket (ft).

CLAT

CLATG

- Cosine of geocentric and geodetic latitude.

D(l)-D(2)-D(3)

Launcher coordinates of aerodynamic velocity (ft/sec).

DIF

Difference between the time of the simulated trajectory and the time of the beginning of the phase

(sec).

EMPT

- Sum of the radius of the oblate spheroid along radius vector to the rocket and the height of the
impact area above sea level (ft).

ENDPII

EPSQ

= Square of the eccentricity of the earth.

Integration interval (sec).

IDEN

Alphanumeric information used for indentification.

IJK

Current line of print out on page.

NUF

Set equal to 1 if CIIEK is less than or equal to

EMPT.

P1

Print interval (sec).

PO

If PO equals zero, print out occurs.

RO

Radius of earth (ft).

RS

Radius vector of oblate spheroid along radius vector to rocket (ft),

Difference between current time of simulated trajectory and time of the end of the phase (sec).

5,'3S

iihi

SLAT

- Sine of geocentric latitude of launcher.

SLATG

- Sine of geodetic latitude of launcher.

TBO

- Time of end of current phase (see).

TFP

- Time of first print out for current phase (see)..

T1!

- Elevation angle of velocity vector (deg).

TIMEO

- Time (see) rocket leaves end of launcher.

TO

- Time of beginning of phase (sec).

- Aerodynamic velocity in launcher coordinate system.

- Omega, rate of angular rotation of earth.

WT

- Distance earth rotates in time, T.

WS, %YS

- E-W and N-S components of wind (mph).

XL, YL, ZL

- Translation vectors from earth-centered inertial coordinate system to launcher coordinate syster.

XYZ (3)

- Launcher coordinates of position of rocket (ft).

YD (6)

- Relative motion velocities of inertial system with

respect to launcher system.

ZIM

- Impact height above sea level of impact area (ft).

THE INTEGRATION ROUTINES

The trajectory is determined by integrating IS differential equations of motion. Three are of second order, while the other 12 are
first order. Iterated integration is used for the second order equations. The integration is monitored by LRBM through the subroutine
TNTEG.
The subroutine INTEG monitors the mechanics of the numerical
integration through the routines RKG and RKERR2.
The SUBROUTINE RKG (II,TIME) is used to integrate the equations
of motion and uses the fourth order Runge-Kutta-Gill numerical integration scheme. RKG calls upon EVAL for evaluation of the equations

36

SAVE INITIAL
CONDITIONS

CALL
RK(HflE

NOGOOi) -:

ZCAL
RKERR2 (H)

NO
S NOGOOD

DONEj

NO

IS H

C.
O1'?
)
NO

INTVGRAT ON

LNMA

3,

E2

aRMvvN

of motion. IIis the integration interval currently being used.

is the current time in the trajectory simulation.

TIflE

The F!JRROUTINE RKERR(l) is used to check the validity of the

integration. This check is performed by comparing the value, say Yis
obtained for an interval of length 11with the value, say YIII obtained
!y performing two integrations using the value 1I1/2.
follo;ing ratio
(Y1I

Specifically the

YI)/YII is compared with two tolerance values c1 and c2 .

(This comparison is made separately for each of the integrated values.)

If (YII - Yd/YII is greater than

E1

for any of the parameters,

iiis replaced by 11/2 and the integration is redone. If all values of

f
Y 1 )/Y1 1 < E then not only is the integration accepted but the
integration interval for the next step is doubled.
If at any time during the simulation I becomes less than .0001,
a message is sent to LRB! using the indicator DONE to terminate the
simulation. The indicator NOGOOD is used to tell whether the integration is acceptable or not.
SUBROUTINE 'Mtrxl (Al, A2, NC2, PRO)
This subroutine is used to multiply two matrices. Al, A2, NC2,
and PRO are variables used in the calling statement. Al is a 3 x 3
matrix. A2 is either a 3xl, 3x2, or 3x3 matrix depending on whether
NC2 is equal to one, two, or three. The product of these two matrices,
Al and A2, is then stored in the matrix, PPO.
SUBROUTINE TABL (LOOKUP, J)
This subroutine is used to interpolate linearly for values from
the atmosphere, wind, :lach, and time tables. The parameter LOOKIIP
is used to designate which table is to be used.
To conserve computer time the routine is programmed to "remember"
the arguments used in the previous entry to the table. The "search"
is either forward or backward depending upon the present value of the
independent variable.

38

Suppose the independent variable lies outside the arguments of the

table. For values from the ,Mlach or atmospfiere tables linear extrapolation is used; for the wind tables the value is set to zero; and,
for the time tables the last value of the table is chosen.
The interpolation routine for the time table also provides for
the calculations of the derivatives of the moments of inertia, and
of the mass flow rate. The principal flads are:
CD

CG (40)

- Entries of center of gravity in time table.

CL

- Interpolated value of normal force coefficient.

CLD

- Roll driving moment.

CLP

- Roll damping moment.

CMA(20)

- Entries of center of pressure in Mach table.

CHQ(20)

- Entries of coefficient of pitch damping moment in

Hach table.

CMT

- Interpolated value of coefficient of pitch damping moment.

CNA(20)

- Entries of normal force coefficient in M1ach table.

CND(40)

- Entries of roll driving coefficient.

CNP(40)

- Entries of roll damping coefficient.

CON

- Constant used for interpolation.

CP

- Interpolated value of center of preijre.

CX(20)

- Entries of coefficient of drag in Mach table.

DENS

- Interpolated value of atmospheric density.

SDMKOMK(3)

Interpolated value of axial force coefficient.

- Derivatives of moments of inertia.

F (40)

- Entries of mass in time table.

FX(40)

- Entries of moment of inertia about x-axis in time

table.

39

FY(40)

- Entries of moment of inertia abput y-axis in time table.

11H(44)

- Entries of height in atmospheric table.

IP(SO)

- Entries of height in wind table.

Irr

- l1eighIt of rocket above oblate earth.

11,12,13,14 - Entry last used in the various tables.

LOOKUP

- Indicates which table to interpolate.

- Number of entries in the atmospheric table.

I'M

NT

- Number of entries in the time table.

NW

- Number of entries in the wind table.

PA(44)

- Entries of pressure in atmospheric table.

PM(20)

- Entries of Mach number in Mach table.

PRES

RIO(44)

- Entries of atmospheric density in atmospheric table.

SR

SPD

- Interpolated value of speed of sound.

T(40)

- Entries of time in time table.

IINE(6)

TSL(40)

- Entries of thrust in tine table.

WXP(S0)

- Entries of E-W components of wind.

WXS

- Interpolated value of E-W components of wind.

WYP(SO)

- Entries of N-S components of wind.

WYS

- Interpolated value of N-S components of wind.

Number of entries in the Mach table.

Interpolated value of pressure.

Mach number.

Interpolated values of moment of inertia about the xand y-axes, thrust, center of gravity, mass, and center
of gravity of the propellant.

40

TABL (LOOKUP ,J)

INTrRPOLATE

INTERPOLATE

INTERPOLATE

F'OR ATMOSPHERIC

FOR WIND

FRPRMTR

PAWME ELS

VRSU MAH

fA

41

IINTERPOLATE
O

AAEES

VESUST:M

LISTING OF PROCRM
SATENO
1.EXECLJTE
11E3JOt
%I9FTC LRHM
c
C

INSTRUJCTIONS

3a...0.777717.6.o'JMP
113J08
GO.LO6.1CMAP
LIST.kEF~f)ECK(

IN ALL COM~vAN STATE~i-NT

.2=Kq2
4
k RG
TAL

030020
000.30
00040
OQ0bO

J=JWIVL.A=AOiTuN.=FACT.=CVAL.C=CHK.1=INT
PA' N.PPINJ
=O)T
L=TI4LP

C
C0MMON/DAZ

IXWD.IYWD.',ITHP.COTHP.ENF.THP.ALP

C
COMM0N/0AP/TH.AL.XANT#YANT.DELC.)LTOELHl.T0WTIL
C

COMMON/DAEP/NTYPE
c
COMMON/OAETX/ W~,Y

C
COMMON/D8/

YSTOI.22).UYSTOPC22).RJFSTW(221. IUFEN0

C
C0MMON/DFC/IWIND. 1,.8HZ
C
C0MM40N/3f3TP/NW.HP(50)
C
COMMO0N/DE/WEGPSL .GCO-NI .GNU
C
COMMON/OEL/ SLAT.CLAT.ZLATG.CLATG.TIMLO.W
C
COMMON/OCP/ TO
C
COMMON/DC/ NOF. lOUT
C
CONIMON/DCF/PI .TfO
C
C0OM0ON/flCl/

C
COMMON/OI/ DONE
C
CON9MON/DIX/ TIME
C
IP2/ NOGOOD
C
LOMMON/VT/ 11*12.13.14
COMMON/

C
COMMON/OP, UO.VO.*0.ZL.XLAT.NFAZE.NbST.IFOUT
C
COMMON/OF/ TMIl
C
COMMON/OX/

1.3K

C
COMMON/DL/ PO
C
COMMON/DEC/EPSO
C

4200540

00070
oooao
00090
00100
00110
00120
00130
00140
00150
00160
00170
00180
00190
00200
00210
00220
00230
00240
00250
00260
002!70
00280
00290
00300
00.310
00320
00330
00340
00353
00360
00370
00380
00390
00400
00410
00420
00430
00440
00450
00480
00470
00480
00490
00500
00510
00520
00530

DIMENSION AIC(3.3) .BIC(3*3)

C
C
C

INITIALIZE
CNTRO0.
READC5.600 IXN0
600 FORMA'rCI6XF8.O)
101 CALL PINTRJ
DATA REO.PSL*rWEPSO.GCON1,GNU/20926428..2116...72921E-04..0067227
1.1I.40775E 16..273E-03/
TH-TH*.01745329
AL=AL*.01745329
XLAT=XLAT *.01745329
SLAT=SINCXLAT)
CLAT=COS (XLAT 1
EPSO=EPSO**2
TEN=EPSO* SLAT *CLAT
FIVE=I .- EPSO*(CLAT**2)
DLTA=ATAN2 (TENoFIVE)
XLATG-XLAT+OLTA
SLATG-SI NCXLATGI
CLATG=COS (XLATG)
RO=REO/(1.+EPSU*SLAT*SLAT/(I.-EPSO) )**.5
ROA=RO.ZL
B1CC I n--I.O
BIC(2*1 10.0
BIC(3.1 1=0.0
8ICCI .21=0.0
BIC(2.21=-SLATG
B IC (3.2 =CLATG
BICCX .31-..
BICI2.3)=CLATG
B1Cc 3,3) =SLATG
TI MEO=TO

C
C
C

SET UP FIRST APPROXIMATION FOR

ITERATION FOR LAUNCHER SETTINGS

IF (NTYPE.EO.2) GO TO Jos
IF (NTYPE.NE.3) GO TO I
IYWD. 1
IXWD= I
rH-P.TH
ALP=AL
105 RANGEsSORT CXWANT*XWANT+YWANT*YWANT)
SI THP-XWANT/RANGE
COTH-P-YWANT/RANGE
CALL AOPTUN(12
WR=WVS*COTH.P+WXS*S ITHP
IF (WR.GE.0.) GO TO 102
DETH-OEL T

00550
00560
00570
00580
00590
00600
00610
00620
006J0
00640
00650
00660
00670
00680
00690
00700
00710
00720
00730
00740
00750
00760
00770
00780
00790
00800

00a10
00820
00830
00840
00850
00860
00870
00880
00890
00900
00910
00920
00930
00940
00950
00960
00970
00980
00990
01000
01010
01020
01030
01040

o01050

GO TO 103

01060
01070

102 DETHxDELH
103 YUV=OETH*COTHP-ELC*SITH~P

43

XUT;!ETHl*SI THP4OELC*COTHP
X0IFWXS*XUW*528 0
YOIFNWYS*YUW*528 *
STCAwSIN(TN )*COSIAL)
COSTsCOS(ITN)
THI-.ATAN2(STCA.COST)
STSAwSIN(TH )*SIN(AL)
TH2=ATAN2(STSA.COST)
THnt THI-YOIF/TOWTIL
TI42oTH2-XD IF/TOWT IL
STI42uStN(TH2)
CTHImCOS(THI)
CTHZ.COS (TH2)
STMImSINITNI I
STCT21-STHZ*CTHI
CTS2)nCTH2*STNI
AL*ATAN2(STCT2loCTS21)
COASTuCOS (AL)
STOCT1-STHI/CTHI
BSTORmASS
(STOCT 1/COAST)
THeATAN (BSTOR)
C
C
C
C
C

CNTR.WHICH TRAJECTORY
CONTINUE

IN--

NUF-IsIMPACT--IJK

INITIALIZING

100 NUF-0
I JK=O
11=1
I2.1
13-1
4.1
Y(j .5)-TIMEO
STN.SIN(TN)
CTHwCOS (TN)
SALoSIN4AL)
CALNCOS(AL)
Y( 17.5).0O
V (18 *5) .POA*CLAT
Y( 19*5)=ROA*SLAT
A IC (I* 1) =ST*SAL
AIC(2.1 )uSTH*CAL
AIC(3. )-CTN
AIC(I'2)=CAL
A! C(2.*2).-SAL
AIC(3.21.0*0
AIC(I.3)-SAL*CTH
AIC(2.3)-CAL*CTH
AIC (3e3) -- STH
CALL MTRXL (BIC*AICa3.Y(8t5))
DY(117.5) =Y(895) *UO-W*Y( 18.5)
OY(IS.5)=Y(9.5)*UOW*Y(17.5)
DY( 19.5)-Y(1095)*UO

44

IS LINE OF OUTPUT ON PA

01080
01090
01100
01110
01120
011 30
01140
01150
01160
01170
01 IBO
01 190
01200
01210
01220
01230
01240
01250
01260
01270
01280
01290
01300
01310
01320
01330
01340
01350
01360
01370
01380
01390
01400
01440
01450
01460
01470
01480
01490
01500
01510
01520
01530
01540
01550
01560
01570
01580
01590
01600

Y(2.5)=Y(BS)DY(17.5)ey(995)*DYC18*5)eY(IOo51*OY119.S]

~Y(3-5)=Y(Ill5)*OYllT,5)+Y(1295)*DY(18*5)eY(1395)*DY419.51
Y(4951=YlI45)*OY(1T,5)+YiI5#5)*DY(18,51+YlI69Sl*OY119*51
Y(5,5)=0=
Y(6.5)0O.
Y(7.5)=W
Y(20.51:DY(17*S)
Y(21.5)=DY(18o51

2000

C
C

Y(22.5)=DY119.5)
DO 3 1=2e22
ROFFI(195)=O.O

01640
01650

CONTINUE
CALL PHASIN(I)
DO 2000 I=1.22
YSTOR(i)=Y(I.5)
DYSTOR(1)=DY(I.5)
ROFSTR(I)=ROFFI(I.5)

01660
01670
01680
01690
01700
01710

CONTINUE
Ji=
IF (NTYPE.NE.4) GO TO

01720
01730
01740

IwIND

WHAT

ENTRY

AT

IN

BALL

FAC

SUBROUTINE

(HEIGHT TABLE)

C
IWINOsI
J2.1
JI IS THE PHASE
J2 IS THE PHASE

C
C
It
C
C

WE ARE CURRENTLY WORKING IN

(SF) TO PICK UP TRAJECTORY AT FOR 8F PROGRAM

BHI=HP(IWINO)

IH-2

IN BALLISTIC FACTOR

ROUTINE

01850
01860
01870
01880
01890

IH-I
CALL EVAL(S)
CALL XOUT(5.0)

C
C

MONITOR THE

INTEGRATION

01900

C
I H=AMINI(PI.
50 IF (IH.NE.21
2

01910
01920
01930

TMII)
GO TO 990

CALL BALFAC(I)
IF(IH.EO.I) TMII*TMII*10.

01940
0194S
01950

iJ2-Jl
990 TIME=Y(I.5)
C
C
C

NOGOOOO
991
5

INTEGRATION

NOGOO0O
CALL INTEG
IF (DONE.EO.2.) GO

01960
INTERVAL

ACCEPTED

MONITOR

OUTPUT

SET H=H/2.

02040

C
C

=1

01970
01980
01990
02000
02010
02020
02030

TO 9

C
C

01770
01780
01790
01800
01810
01820
01830
Ole40

C
10

01750
01760

02050
IOUT-1

NO

OUTPUT--IOUT=29 PRINT OUTPUT

4S

02060

C IOUT=3. END OF PHASE-IOUT=4. IMPACT

C
CALL C'KOUT
GO TO 150.67,8)1.IOUT
9 CALL XOUT(51)
GO TO 8
C
C
IFOUT-I*PRINT OUT ONLY AT END OF PHASE
C
GO TO 5
6 IF (IFOUT*EO.)
CALL XOUT(5.O)
GO TO 50
7 JI-Jl+l
CALL XOUT(501
C
C
NFAZE IS THE NUMBER OF PHASES USED
(TOTAL)
C
NBST
PHASE END OF WHICH BOOSTER IS TO BE PICKED AT
C
IF (JI.GT.NFAZEI GO TO 1100
IF (NOST.E *(JI-II)GO TO 20
25 CALL PHASIN(JI)
CALL EVAL(5)
CALL XOUT(SoOl
GO TO
C
C
C

C
C

STORE Y AND DY ARRAYS FOR BOOSTER TO GROUND

20 DO 160 K=-,22
YSTOR(K)-Y(K.5)
360 CONTINUE
00 161 L=2,22
DYSTOR(LI=OY(L .)
ROFSTR(L)-RCFFI CL.15
161 CONTINUE
GO TO 25
8 CALL XOUT(5#OI
IF (NTYPE.EO.41 GO TO 12
IF (NTYPE.EO.31 GO 70 13
IF (NTYPE.EO.2) GO TO 13
IF (NBST*EO.O) GO T3 11
UNSTORE Y.DY*T ETC
FOA BOOSTER
DO 151 K*I22
Y(K.51sYSTOR(K)
DO 152 Lw2*22
DY(L95)=DYSTOR(L)
ROFFI (L.5)-ROFSTR(..)
152 CONTINUE
JI=Jl !
TO=Y( .55
CALL PHASIN(JI)
151

CALL

02070
02080
02090
02100
02110
02120
02130
02140
02150
02160
02170
02180
02190
02200
02210
02220
02230
02240
02250
02260
02270
02280
02290
02300
02310
02320
02330
02340
02350
02360
02370
02380
02390
02400
02410
02420
02430
02440
02450
02460
02470
02480
02490
02500
02510
02520
02530
02540
02550
02560
02570
02580

XOUT(5-O)

GO TO 1

02590
46

13 CALL AOPTUN(3)
C
C
C

ENPwl

02600
02610
02620
02630
02640
02650
02660

FINISHED

IF (ENF.EQ.I.) GO TO
GO TO 100

11

C
C

IN

IS

INDICATOR TO CHECK

IF WE ARE PAST HEIGHT OF CURRENT BF

12 CALL BALFAC(2)
JI =J2

'02670
02690
02700

C
IF CIWIND.GT.NWIIBFENDmi
IF (IBFEND.EO.I) GO TO 1010
CALL PHASIN(JI)
TMI I=TMI 1/10.
HwAMINI (PIsTMIII
TO-TIMEO
CALL CHKOUT
11-1
12-1
13-1
14-1
IH=2
BHI=HPC IWINDI
CALL EVAL(5)
CALL XOUT(5tOl
GO TO 990
1010 CALL BALFACtI)
1100 CNTR-CNTR+I.
IF(CNTP.LT.XNOI GO TO 101
RETURN
END
SIBFTC LOPSUN
LISToREF*DECK
SUBROUTINE AOPTUNCK 1
C
C
IN ALL COMMON STATEMENTS D.DRIVER.A-A0PTUN.B=8FACT.EsEVAL.CCHiKIuNT
C
R2-RKERP2.R-RKG.T=TABL.F.PHASINP.PINTRJ.X.XOUTAND L-T I2LP
C
COMMON/OA/ IXWO. IYWD.SITHP.COTHP.ENF.THP.ALP
C
C0MMON/0AP/TH.AL.XWANT.YWANT.DELC.oELT.DELHT0WT IL
C
COMMON/DAEP/NTYPE
C
COMMON/OAETX/ WXS.WYS
C
COMMON/AP/ XFIRST.'FIRST.XCHANG.YCHANG.XLAST.YLAST.N8F,8F(50)
IDELTA
C
COMMON/ABETP/ WXP(50) .WYP(50)

47

02730
02740
02750
02760
02761
02762
02763
02764
02770
02780
02790
02800
02810
02820
02830
02840
02850
02860
02870
02880
02890
02900
02910
02920
02930
02940
02950
02960
02970
02980
02990
03000
03010
103020
03030
03~040
03050
03060
03070
03080

03090
COMMON/ABXL/

XYZ(3)

c
DIMENSION
THETA(IIllIg)ALPHA(Ito1I)
GO TO (1.2*3).K
1 IF (NTYPE.EO.21 GO TO 6
C
C
c

SET UP

INITIAL WIND

PROFILE

WXS.XFIRST
WYS-YFIRST
RETURN
c
C

SET UP

SUBSEQUENT

WINO PROFILE

C
2

IF (WXS.EQ.XLAST) GO
IXWDNIXWD+1
WXSwWXS+XCHANG
GO TO 21
10 IF (WYS.GE.YLASTI GO
RYWD-IYWD+1
IXWO-I
WXS.XFIRST
WYS-WYS+YCHANG
GO TO 21
C
C

OBTAIN BALLISTIC

TO

10

TO

WINOS FOR

ITERATIVE

C
60

61

C
C

03360
03370
03380
03390
03400

SUMYOSUMY+WYPCI)*BF(I)
CONTINUE
WXSuSUMX

03410
03420
03430

WYS.SUMY
RETURN

03440
03450

OUT

03460
03470

ANGLES

C
C
C

03230
03240
03230
032u0
03270
03280
03290
03300
03310
03320
03330
03340
03350

SUMX-O.
SUMYO.
DO 61 I=I*NBF
SUMX-SUMX+WXPII)*BF(I)

PRINT
11

03100
0-110
03120
03130
03140
03150
03160
0317C
03180
03190
03200
03210
03220

DO 75 MI-IvIYWD
DO 75 NI=I*IXWD
THETA(NIMI)=THETA(NI

MI)*57.29578

03480
03490
03500
03510

ALPHA(NI.Mi) ALPHA(NIMI)*57.29578
75 CONTINUE
00 71 M1ItlYWO
WRITE (6.70) (THETA(I.M).ALPHA(I.M).I=I.IXWOI
70 FORMAT (IHO.2X.11(F8.3.2X)/(11C2X.F8.3)))

03520
03530
03540
03550
03560

71

03570
03580
03590
03600

CONTINUE

ENF-i

FINISHED TABLE

12 ENFEI.

03610

48

RETURN
C

03620
03630
03640
03650
03660
03670
03680
03690
03700

C
CHECK IF IMPACT IS WITHIN TOLERANCE
c
.1 XOIF.XYZCI )-XWANT
YOIFXYZ (2)-YWANT
TESTVAOSQRT (ZWIF*XOIFI.YOIF*YDIF)
IF CTESTvA.LE.DELTA) GO TO 20

C
(~ITERATIVE

TECHNIQUE FOR COMPUTING LAUNCHER SETTINGS

30 STCA-SIN(TH)*COS(AL)
COST-COS THI
THI-ATAN2(STCA*COSTI
STSA-SINE TH)*SIN(AL)
TH2-ATAN2(STSA*COST)
THI .TNI-YDIF/TOWTIL
TH2T2-X IF/TOWT IL

03730
03740
03750
0:760
03770
03780
03790

STH2=SINITH2)

03800

CTH2*COS(TH2P
STHZ-SINITHI)
STCT21-STN2*CTHI
CTS2I =CTH2*STHI
AL-ATAN2(STCT21 9CTS21)
COAST.COS (AL)

03820
03830
0:3840
03850
03860
03870

STOCTI-STHI/.CTHI

03880

BSTOR=ABS

03890

(STOCTI/COAST)

TH-ATANCSSTOPI
RETURN
20 IF (NTYPE.E0.2) GO TO 12
C
C
C

STORE LAUNCHER SETTING FOR ANGULAR OUTPUT

IN TABLE

THETAC 1XWD. IYWD)OTH

ALPHA(IXWOIYWDIwAL
GO TO 2
C
C
C

DETERM4INE APPROXIMATE
21

03710

03720

THEORETICAL DISPLACEMENT FOR NEXT TRAJECTORY

IF (WXS.EQOO.AN0.WYS-E0.O.l GO TO 23
IF (IYW0.EQ.1) GO TO 22

THWtHETACIXWD.IYWO-1)
ALsALPHA( IXWO. IYWO-I 1
IF(COTHP.GE.O.) GO TO 41
DETHxOELT
GO TO 42
23 THETA(IXWD.IYWOIUTHP
ALPHA( IXWD. IYWO)SALP
GO TO 2
41 DETH=DELH
42 YUWODELT*COTHP-OELC*S ITHP
YD IF= YCHANG*YUW*5280 *

0390-3
03910
03920
03930
03940
03950
03960
03970
03980
03990
04000
04010
04020
04030
04040

04050o
04060
04070
04080
04090
04100
04110
04120
04130
04140

49

XDIF-0.
GO TO 30
22 THmTHETAC IX~l)-I *IYWDI047
ALuALPHA( IXWl)-1 IYWD)
I~tSITHP.GE.O.) GO TO 51
DETH-DELT
GO TO 52
131 DETH.OELH
52 XUW-DELC*COTHP+DETH*S ITP
YOIFnO.
XD IF
VXCHANG*XUW*5280.
GO TO 30
END
SIBFTC 8FACT
LIST*REF*OECK
SUBROUjTI NE BALFAC (I NTA)
COMKV)N Y(22*5:.DY(22*5)*ROFFIC22s5)
c
C
IN ALL COMMON STATEMENTS ODRIVERAXAOPTUN.BUBFACTExEVAL.CCH(,IUINT
C
R2-RKERR2.RRKGTrAL.FnPHASINPPINTRJ.XXOUT.AND L&T12LP
C
COMM0N/08/ YSTOR(22i .DYSTOR (22 I ROFSTR (22), ISFENO
C
COMM0N/DSE/IWIND, 11HB
C
COMMON/DBTP/NW. NP 50)
C
COMMON/ABETP/ WXP(50 .WV~l501
C
COMMON/ABXL/ XYZ(3)
C
COMMON/SET/ HT
C
COMMON/B 1R2/ YSAVE(22,O0YSAVE122) .ROFFS(22)
C
COMMON/SPX/
IDEN

C
BALFAC MONITORS BALLISTIC FACTOR PROGRAM
IBFEND=I WE FINISHED LAST TRAJECTORY
IF (IBFEND.EQ*1) GO TO 102
IF (INTA.EO.2) GC TO 15
C
C
C

CHECK( HEIGHT OF

SIMULATED TRAJECTORY VERSUS HEIGHT OF

IF 4HT.GE.BHII
GO TO 101
C
C
C

04180
04190
04200
04210
04220
04230
04240
04250
04260
04270
04280
04290
04300
04310
04320
04330
04340
04350
04360
04370
04380
04390
04400
04410
04420
04430
04440
04450
04660
04410
04480
04490

04500
DIMENSION XIMPAT(50I*YIMPAT(5O)

C
C
c
C

04150
04160

STORE

GO TO

I8

Y AND Dy ARRAYS FOR NEXT

TRAJECTORY

180 00 103 1s1.22

103 YSTOP(I:-YSAVEII)

so

CURRENT WIND PR

04510
04520
04530
04540
04550
04560
04570
04580
04590
04600
04610
04620
04630
04640
04650
04680
04670

104 MP=.22
DYSTOR(M P DYSAVE(141
ROFSTR N I-OFFS (MI
104 CONTINUE040
I 1.1
101 RETURN

0*E80

0O

04690
04710
04720
04730

04740

CALCULATE BALLISTIC FACTORS

AND UNIT WIND EFFECT

j.750

04760
102 WRITE (6.1101 IDEN
DIFXuXIMPAT(NW)-XIMPATII)
DIFY'YIMPAT(NWI-YIMPAT(1)
DI FX= IFX/528' *
O IFY-DIFY/.5280.

04770
04780
04790
04800
04810

ELRSORT (0IFX**2+D IFY**2 I

0420

WIND.SORT(WXP(NW)**2+WYP(NW)**2)

04830

UN ITOELR/W IND
DEN0MO *
00 105
IINW
DIFX=XIMPAT(I )-XIMPAT(l)
DIFvwYIMPAT(I)IYIMPAT(1 1

04840
04850
04860
04870
04880

DIFX-01FXo5280*
DIFY-OIFY/5280.
RANGE-SORT(OIFX**2+DIFY**2)
RATI OuRANGE/OELR
8FACT 'RAT I -DENOM
WRITE (6.111) HP(I)eXIMPAT(I),YIMPAT(I).DIFXDliFY.RANGE.RATIO#
IBFACT
DENOMwRAT 10
105 CONTINUE
WRITE (69120) UNIT
120 FORMAT CIHO/IHOI X17HUNIT WIND EFFECT=FIO.6)
RETURN
C
C
C
C
C

STORE CURRENT
IWINO

IMPACT AND UNSTORE

04890
04900
04910
04920
04930
04940
04950

04960

Y AND DY ARRAYS FOR NEXT TRAJECTORY

IS THE ENTRY OF WINO TABLE CURRENTLY USING

ISO

IWIND-IWIND+1
XIMPAT(IWIND-I IXYZ(I I
YIMPAT(IWINO-1I.-XYZ42)
D0 151 Kul.22
Y(K*S)-YSTORiK1
151 CONTINUE
00 152 Lu2v22
DY (L#5)wOYSTOR(L 1
ROFFI (L.i)-ROFSTRCLI
152 CONTINUE
GO TO 101
110 FORMAT (IHI/lH012A6/lHO/IHOIOX,6HHEIGHT.13X.IHX,13X.IHY.14X,
IOELTA X.8X97HOELTA Y.9XSHRANGE.IOX.SHRATlOI1X.4H8.F./INOI
III FORMAT (1HO5XFIO.2*7Fl5*4)
END
p5

7o4

04970
04980
04990
05000
05010
05020
05030
05040
05050
05060
05070
05080
05090
05100
05110
05120
0513o
54
05150
05260
05170
05180
05190
05200

SI6FTC

C
C
C
C

EVALU
LIST.PEF.OECC
SUBROUTINE EVAL (J)
COMMON Y(22.5).OY(22.5).Q0FI

22.5)

IN ALL COMMON STATEMENTS DuDRIVER.A-AOPTUN.8V8ACT.EEVAL.CCH(.1.INT

P2-KRR2.RRKG.TTABL.FPHASIN.PPINTRJ.XXOUT.AND L-712LP

c
C
C
C

C
C
c
C

COMMON/DAEP/NTYPE
X.Y
0NNDEX
OMNDEX
X-Y
C0MNSEIIO
OMNDEIIDI-H
COMND/CO.S.CN.N
OMND/E.5.CN.N
C0MNDC
COMNDL
COMNAcT/~P5).Y(O
OMNAEP
MT50
C0MN/3T
COMMN/qE/
C0K.R54,o
COMN/C
HKP
OMNEZ
A
COMNCS
OMNEP
A
0MNDCE5
OMNnCES
COMNcT~H.SP.XCACACQF.YTLC.MR.M0,
COMNEL0OVSP*C.NACCLD.F.Y..CLPM*MDMOK()

05280
050
05300
0HaN5310
053)0
05330
05340
05TCL350T.LAOTIE.
05370TSLT*LAGTIE.
05370
3X(0oY(0
5380
HT054Z00
05420
054450
054403
05460o
05460
0537.
05500#CL
05500

C
COMNEF/A
CCMNE~

05210
05220
05230
05240
05250
05260
05?70

REL

COMMN/EP/KAN05510
FIN
E 00.XL

DIMENSION A 13.3) ,F(3).E(3).AM(3*31,WS(3)

LUUIVALFNCL cEll)*WXI.(E(2loWYI*(E(3)*WZ)
DIMENSION WOOT(6)
EOUI'IALENC
(WDOTI1).P).IVDOT(23.Q).(WDOT(3).RI
L
C NORMALIZE DIRECTION COSINES
C
DO 1
J)P-~1+2
AYS(JPT(Y(I.J)4e2+Y(I+I.J)*2+Y(JP.J)**21
DO 1 K1I.JP
1 Y1K.J)-YtK.J)/AY
C
C INTERPOLATEC POPW INOSi AND COMPUTE COMPONENTS IN BODY SYSTEM

~.1.

CHEKF-S83T(Y(17.J)**24YCI8.J)**2+Y( I9.J)**2)
RS2962.(I+PQ((9.)CE)*/l-PU)*5
HRO- SORT(II7.J)*42+YI l8.J)**2+Y( l9.J)**2Z
HT -HPD-P',

S2

05540
05550
b7

05580
05590
05800
05610
05620
05630
05640
05650
05660
05670
05680
05690
05700
05710
05720

CALL

05730

TABL(19J)

IF CNTYPEoEO.31 GO TO 12
IF (NTYPE.NE.4) GO0 T0 11
174 IH.NE.21 GO TO 1305760
IF (NT.GT.8N23 GO T1305765
WXS-WXP(IWINO)057
WYS=WYP( IWINO)
GO TO 12
13 WXS.O.050

05740
05750

WYS-O*
1A2.J
CAL
GO
J
12 CALL T)LWXS
12VS (2) .WS
VS (31.0.
VT.V(Y(I J)TIoO

05810
05820
05830
05850
05850
05870

C
C COMPUTATION OF LAUNCHERTO ECI
-O(T
A(.)

05780
05790

ROTATION MATRIX

05880
05890
05900
05920

A(2.1)= -CS(WT)
A(3-1)- 0.0
A(3*21= CLATG
A(3.3)- SLATG
A(1*2)- -Ai3o3)*A(2*1)
A(2.2)- A(3*31 *A(1911
A(3.2)*A(2.1 I
A(1 .31=
A(2e3)2 -A(3v2)*A(I~lI
CALL MTRXL (AWSI.F)
F(l)u-F(I)*1.466667-W*Y(18@Jl
F(2)--F42*1.466667+V*Y(179J?
F(31.-F (3)*1 .466667
WXINT.F(1 I
WYINT-F (21
WZ INT-F43)
00 2 11,o3
1+7oJ)
AMU , I IY
AM(29I)-Y(1+10*JI
2 AM(3.1 )aV(1413*J)064
CALL MTRXL (AM.F.1.EI
OYC 17.J=-Y(20*J)
OY(18.J)-Y(21 ,J)
DY(C19.J).Y(22*J)
CALL MTPXL(AM.oY(17*J1.1,Y(2.J))
UPuY(29JI-WX
VPaY(3.J)-WZ
VP=Y (3 *J -WY
VA.SORT (UP**2+VP**2+WPO*Z)
SS-VP/VA
SA-VP/VA
OP. 4.54H0*VA 1*VA
RMwVA/VS
CALL TABL(3#J)

S3

05930
05940
05950
05960
05970
05980
05990
06000

06010
06020
06030
06050

06090
06100
06110
06120
06130
06140

CALL MTRXL (AM*Y(5.J).1,*ID0T(j))

DVxREFL*.S/VA
PO3P*DV

200

CD-CNP*PO
oqxPD*S8
CY-+CNA*SB-CD*SA
CZ-CNA*SA+CD*SB
GO TO(200920092
CALL TABL(4*J)

C
C COMPUTE
C
C

MAL=I

1.200-200)-J

TOTAL AERODYNAMIC FORCE

INPUT

CENTER

OF PRESSURE

AND MOMENT
9 MAL=;

COEFFICIENTS

'NPUT

RESTORING MOM.ENT COEFFIC

C
201 GO TO (3.4)oMAL
3 CPxCMA
C'A--CNA*ABS(ICG-CP)I/REFL
4
CMY-CMA*SA+CM04:)*DV
CMX=CL0*FI1NC+ I,-*REFL*.*5/VA)I*CLP
CMZ.CMO*R*DV-CMA*SB
CTaTSL+AE* IPSL-PA)
C
C COMPUTE ACCELERATIONS IN BO')Y SYSTEM
C
IF IKAN.EQ.2) CT=0.
C0S2L= (Y(17,J)*Y(17-J)+Y( 18,JI*Y(18.j)-Y(19..)*Y( 19,J))/(H-RO**2)
G=GCONI/HRO**3*(I1.-3.*REQ*REU*GNU*(1 .- 3.*COS2L)/HRO**2)
CALL MTRXL (AMoY(17.J)1
FQA=QP*DSO
FOR.FOA*REFL

F)

FZ-FY

-F(1)*G+(CT-CX*F-uA)/FM
0V12$J)-F12)*G-ICY*FQA)/FM
0Y(3oJ)=
-F13 )*G- (CZ*FQA I/FM
DV (49.J I
WD0T(4)=((FY-~FZI*Q*R+CMX*FOR)/FX
WD0T(5)I( FZ-FXI*P*R+CMY*FQR)/FY
WDOT(6).=I(FX-FYI*P*a+ILMZ*FaR)/FZ
IF(KAN.EQ.2) GO TO 41
wDOrI4)=wDOTI4)-DMKOMK(1 )*P
WDOT151 =WDOT (5) -UMKOMK (2) *Q
DMKOMKI2)*R
WDOT(6Is*DOT(6)
DAMP=FL0W*(Xi-T-CGI**2
WOOT (5) :WDOT 15 +DAMP*Q/FY
WDOT (6)=WDOT16 )+DAMP*R/FZ
41 CALL MTRXL(YI8.J)tWDOT(4)o1 DVI5,J))
DO 5 KA=8 9 10
KD-KA+3
KEsKA6
DV IKA .J I R*Y IKj. J -Q*Y( KE.J*
DVI KD, J I P*Y KE J I-R*Y (KA. JI
5 DY(KE.J)=G*Y(KAiJ)-P*Y(KD.J)
CALL MTRXL (Y (8,JI *DY s2v*JI* DY (209 J))
RETURN

06i50
06160
06170
06180
06190
06200
06210
06220
06230
06240
06250
06260
06270
06280
06290
06300
06-310
06330
06340
06350
06.(o
06370
06380
06390
06400
06410
06450
06460
06470
06480
06490
J6500

06540
06550
06560
06570

501

54

06580
06590
06600
06610
06620
06630
06700

06710

END
STE- IT
SIBFTC

MTR
LIST*REF*OECK
SUBROUTINE MTRXLcA1,A2*NC2*ORO)
DIMENSION Al (3oS)9A2(3o3) sPRO(3.3)
DO 1 I113
DO I J=lNC2
PRO(I,J)- 0.0
DO I K=1.3
1 PROClI,)
=-PRO(Ioj)+AICIvK)*A2(KvJ)
RETURN
END
SIBFTC CHK
LlSTsREF9DECK
StBP'<UT INE CHKOUT
CC,' -N Y(22*5)9DY(2295)oROFF1(22s5)
C
C
IN ALL COMMON STATEMENTS DsDRIVERsA=AOPTUNB=BFACTE=EVALC=CH-K,1
C R2=RKERR2iR=RKGT=TABL.F=PHASINP=PINTRJX=XOUT.ANO L=TI2LP
c
Z;O;:MON/DCP/

06890

TO

06900
06910
06920
06930
06940
06950
06960
06970
06980
06990
07000
07010
07020
07030
07040

C
OMMON/DC/ NUF' lOUT
C
COMMON/DCF/P! .TBO
C
COMMON/DC1/ H
C
COMMON/EC/

CHEK*RS

COMMON/CF/

TFP

COMMON/PC/

ZIM

C
C
C
C OMMON,/DEC/EPSO
C
C10UT1.o
C
C

IOUT=39

NO OUTPUT--1OUT=29

TO

07060

PRINT OUTPUT

END OF PHASE--IOUT 49

NUF=O
IF(YI.1o5).EQ*TFP)G0

INT

06720
06730
06740
06750
06760
06770
06780
Q6790
06800
06810
06820
06830
06840
06850
06860
06870
06880

07070
07080
07090
07100

IMPACT

CHEK=S IRT(Y(17.5)**2+Y(18t5)**2+Y(19,5)**2)

07110

RS=20926428./(l.+EPSQ*(Y(1995)/CHEK)**2/(1.-EPSQ))**.5)

07120

EMPT=RS+Z TM
IF (CHEK.LL.EMPT )NUFxl
IF(NUF.EO.1.AND.Y(1.5).GT.20.)GO
ENDPH-ASS(TBO-Y( 1.5))
IF(ENDPH.LT..OO1jGO TO 3

07130
07149
07150
07&60
07170

DiFnABS(Y(I .5)-TO)
PO=AMOD(DIF*PI)
lF(ABS(PO-PI).LT..0O1)GO
IF(PO*LT*.01)GO TO 4
POSR 1-P0

TO 4

TO

07180
07190
07200
07210

07220

HuAMINI(H9PO4ENDPH)
lOUT&I
GO TO 1
4 IOUT&2
H*AMINI(H.PI ENDPH)
GO TO 1
3
OUT=3
TOmTBO
GO TO 1
2 IOUT=4

07230
07240
07250
07260
07270
07280
07290
07300
07310
07320

GO TO 1

0733G

5 TO=TFP
ENDPH=TBO-TFP

07340
07350

GO TO 4

07360

1 RETURN
END
SIBFTC INTEGR
LIST.REF#DECK
SUBROUTINE INTEG
COMMON Y(22s5).DY(22v5)4ROFFI (22.5)
C
C
IN ALL COMMON STATEMENTS DaDRIVER,A=AOPTUNB=BFACTE=EVAL.C=CHK.I=INT
C
R2=RKERR2,R=RKG.T=TABLF=PHASIN.P=PINTRJX=XOuT.AND L=TI2LP

07370
w7380
07390
07400
07410
07420
07430
07440

07450
COMMON/DC1/ H

C
COMMON/DI/ DONE
C
COMMON/DIX/ TIME
C
COMMON/OIR2/ NOGOOD
C
COMMON/BIR2/

YSAVE(22)9OYSAVE(22),ROFFS(22)

C
COMMON/IR2/ INO#YHALF(22)
C
C INTEG IS THE INTEGRATION MONITOR WHIC CONTROLS
C AND THE EkROR CHECK
C
SAVE THE INITIAL CONDITIONS
DO 1 1=1%22
YSAVE(I)=Y(1.5)

THE RUNGA KUTT

DYSAVE(I)=OY(i'5)
1 ROFFS(I)=ROFFI(i.5)
CALL RKGIHoTIME)
3 NOGOODzO
CALL RKERR2(H)
IFIIND.EQ19)GO TO 24
IF(NOGOODOEO.O) GO TO 23
C IF NOGOOD EQUALS 0 THE INTEGRATION IS ACCEPTABLE. IF
C EQUAL TO 19 DT IS JOUBLED FOR THE NEXT PASS.
IF(H.LT..0001)
H=H/2*
Y(II)zTIME
DO 2

GO TO 25

INTEGRAT

IN ADDITION

IND 1

07460
07470
07480
07490
07500
07510
07520
07530
07540
07550
07560
07570
07580
07590
07600
07610
07620
07630
07640
07650
07660
07670
07680
07690
07700
07710
07720
07730
07740

1=2.22

07750

56

2 Y(I,5)=YHALF(l)
GO TO 3
25 DONE'e2.O
C DONE IS AN INDICATOR TO TELL IF THE SIMULATION SHOULD BE CONTINUED 19Y
RETURN
24 H=2.O*H
23 DONE=l.O
RETURN
END
SIBFTC RKER2
LIST9REFvDECK
SUBROUTINE R1KERR2 (H)
COMMON Y(22v5)9DVC22s5)9ROFFl (22.5)
C
C
IN ALL COMMON STATEMENTS DSDRIVERAAOPTUN.BB5FACTE=EVAL.CaCKliUNT
C R2=RKERR2,RZRKGT=TABL.F=PHASINP=PINTRJX=XOUT.ANO L=T 2LP
C
COMMON/B1R2/ VSAVE(22)9DYSAVE(22) .ROFFS(22)
C
COMMON/1R2/ INOVHALF(22)
C
COMMON/FR2/EPT INYt*EPB IG
C
COMMON/D IR/ NOGOOD
C
DIMENSION YI(22)oTRUNC(22).YI 1(22)
VI (l)=V(I .5)
DO 1 1=2.22
V 1(1) V(1,5)
I CONTINUE
00 7 K=2922
Y(K,5)=YSAVE(K)
DV(K95)=DYSAVE(K)
ROFFI (K*5)uROFFS(K)
7 CONTINUE
INDml
HALF=H/2 *
CALL RKG (HALF*V(tol))
DO 10 KmIs22
1D VHALF(K)=Y(K*5)
CALL RKG (HALFoY(I95))
VII
1~V(
.5)08160
DO 3 11.2922
VII (II )VuY1.5)
3 CONTINUE
DO 4 1=5922
IF(ABS(Vli(I)-YI(i)).LT..OOOOOI) GO TO02
TRUNCCIl'.02222222
*ABS(cVII(I)-VI(I))/Vmi()
IF(TRUNCCI).GT.EPBIG)GO TO 5
IF(TRUNC(I ).LTEPTINY)INDOIND+I
GO TO 4
2 IND*IND+I
4 CONTINUE
GO TO 6

S7

07760
07770
07780
07790
07800
07810
07820
07530
07840
07850
07860
07870
07880
07890
07900
07910
07920
07930
07940
07950
07960
07970
07980
07990
08000
08010
08020
08030
0804G
08050
08060
08070
08080
08090
08100
08110
08120
08130
08140
08150
08170
08180
08190
08210
08220
08230
08240
08250
08260
08270
08280

5 NOGOOD=1
6 RETURN
END
LIET*REF*DECK
SIBFTC PKGA
SUBROUTINE RKG CH*.TIME)
COMMON Y(22,5)qDY(22*5),POFF1 (22q5)
DIMENSION AC4).8(4) ,C(4)
DATA AC ),8(1 )sC(1 ),A(2),8C2),CC2).AC3)83)C3h*AC4hdC34)9CC4)/
6 66
666.2.
1.5.2...5..2928932.1...292893291. 7071068.1..1.70710683..16b
2. .5/
Ja)
99 GO TO (I92o3*4)*J
1 Y I1.) =T IME
DO 11 1=2.22
Y (1.1 )=Y(1.5)
DY(I.! )=DVCI'5)
ROFF1 CI * 1 =ROFFI (1.5)
It CONTINUE
GO TO 5
1=1-4/2.
2 Y(I921=Y(1I
GO TO 6
3 Y(I,3)=Y(1*2)
G0 TO 8
4 Y( 1 4)=YC1 ,3)eH/2*
6 CALL EVAL (J)
5 DO 50 1=5.22
YC1.J+1 )nY(1.J)+H*CACJ)4(DYC1.J)-BCJ)*R0FFI CIoJ)I')
ROFF1CI.J+1I) a R0FF1(1,J)+3*A(J)*CDY(I.J)-B(J)*ROFFI.IJ)))
I-CC J )*QY CIOJ
50 CONTINUE
IF(J.EQ*4)GO TO 52
GO TO 53
52 JTzJ+1
S() .5) VII .4)
CALL EVALCJT)
GO TO 100
53 J-J+l
GO TO 99
100 RETURN
END
IST*REFsDECK
$18FTC TAB
SUBROUT INE TAB..(LOOKUP j)
COMMON Y(22*5)o0Y(22s5)sROFFI (22.5)
C
C
C

D0RIVER.A=A0PTUN.8=BFACT.EEVALCCHK,1=INT
IN ALL COMMON STATEMENTS
P2= KERW2.R=RGTTABL.F-PHASIN.P=PINTRJ.XmXOUT.AND L=T 2LP

C
COMMON/PT/ N9HH-(44) ,RH-4(44) .VSC44)*PA(44)
C
COMMON/DAETX/ WXSvlWYS
C
C0MMON/DB3TP/NW*HP(50)
C

58

08290
08300
08310
08320
08330
08340
08350
08360
08370
08380
08390
08400
08410
08420
08430
08440
08450
08460
08470
08480
08490
08500
08510
08520
08530
08b!00
08560
08570
08580
08590
08600
08610
08620
08630
08640
08650
08660
08670
08680
08690
08700
08710
08720
U87J0
083740
08750
08760
08770
08 780
08790
08800
013810

COMMON/DT/

1912,13.14

C
CON MON/ABETP/ WXP (50) ,WYP (50)
C
COMMON/BET/ HT
COMMON /FT/
PM(20), CX(20)hCNA(20) .CMAC2O).CMG(20).T(40).Fx(40) *
IFY(40)sTSL(40),CGC4O).FM(40OhNMNT.CND(20),CNP(2O)
COMMON/ET/DENSSPD.PRES,CD.CL.CP.CMTTME(5) .RMDMKOMK (3) .FLOW

08820
08830
08840
08, , 50
08860
08870

086

I .CLD9CLP
C
C
COMMON/EFT/KAN
C
C
C

I1.12913914
GO

C
C
C

TO

IS WHAT LAST ENTRY FROM THE ATMOS*WINDvMACH AND TIME TAB

CI.293*4)oLOOKUP

INTERPRET

FOR ATMOSPHERIC PARAMETERS

1 IF(HH(1tJ.LT.HT) GO TO 5
7 IF(HH(II-1).LT.HT) GO To6
IF (I1.EO.2) GO TO 6
1I I1-1
GO TO 7
5 IF (II.EQ.N) GO 'TO 6
11=11+1
IF (Il.LT.N) GO TO 1
6 CON=(HH(1I1-1 )-HT )/(HHI I1-I )-HH( 11))
DENS3RHO(II-1H-CRHO(il-I1-RHO(Il))*CON
SPD=VS( 11-1 )-(VS( Il-I)-VS( 11 )*CON
PRESaPA(11-1U-(PA(11-1)-PA(liIH*CON
RETURN
C
C
C

INTERPRET FOR WINO

2
14
12

IF (NW.EQ*.)10; TO 15
IF (HPCI2).LT.HT) GO TO 10
IF (HP(.2-11*LT.HT) GO TO 11
12=12-1
IF (12*EO*I) GO TO 17
GO TO 12
10 IF (12.EQ*NW) GO TO 16
12W12+)
11

08900
08920
08930
08940
08950
08960
08980
08990
09000
09010
09020
09030
09040
09050
09060
09070
09080
09090
0'0100
09110
09120
09130
09140
09150
09160
09170
09180
09190
09200
09210
09220
09230
09240

GO TO 14

09259

CON=(HP(12-1IhHT)/(HP(12-1l-HP(12fl

09260

WXSu=WXP(12-1 )-(wxPC 22-I )-W.:(P( 12) )*CON

0927o

WYSUWYP( 12-1 )-(WYP( 22-I )-W.P(12) )*CON

GO TO 13

09280
09290

17 12=2

09295

16 WXSmO.
WYS=O *

09.300
09310

GO TO

13

09320
59

:1w

15 WXSBWXP(I)
WYSSWYPII)
13 RETURN
C
C

INTERPRET FOR PARAMETERS

VERSUS MACH

09330
09340
09350
09360
09370

C098
3
22

IF

21

C
C

4
32

30

31

34

09390

(PM(13)*LT*RM)GO TO 20

IF (PM(13-1).LT.RM) GO
IF (13*EO.2) GO TO 21
13x13-1
GO TO 22
20 IF (I3.EO.NM) GO TO 21

TO 21

I3F34
I(13*LT*NM)GO TO 3
CONu(PM(13-1-RM/(PM(13-1U-PM(13))
CO.CX( 13-' )-(CX( 13-1 )-CX( 13) )*CON
CLUCNA(13-1)-(CNA(13-'I)-C:NA(13n)*CON
CP.CMA( 13-1 )-(CMA( 13-1 )-CMA(13) )*CON
CMT.CMO( 13-1 -(CMO( 13-1 )CMO( 13) )*CON
CLOCNO(13-I )-(CND( 13-1 )-CND( 13) )*CON
CLPzCNP 13-1 )-(CNP( 13-1 )-CNP. 13) )*CON
RE TURN
.NTERPRET FOR PARAMETERS VERSUS TIME
C
IF (NT.EO.I) GO TO 40
IF (T 14).LT*Y(I*JI) GO TO 30
IF (T(14-1)*LT.Y(19J)l GO TO 31
IF (14.EO.2) GO TO 31
l4w14-I
GO TO 32
IF (14.EQO.NT) GO TO 33
14=14+1
IF (14.LT.NT)GO TO 32
CON.(T( 14-I )-Y(1I J) )/(T( 14-I )-T( 14))
tME( I )FX( 14-I P(FX( 14-I )-FX( 14) )*CON
TME(2)xFY(14-1)-(FY(14-1)-FY(14))*CON
TME(3)aTSL(I4-I) (TSL(14-1)-TSI.(14))*CON
TME(4 )aCG( 14-1 )-(CG( 14-I )-CG( 14) )*CON
TNE(5).FM( 14-I )-(FM( 14-1 )-FM( 14) )*CON
DMKOM4K(I)S(FXU14 I)-FX(14))/(T(14-1)-T(14U)/FX
OMKGMK(2)=(FYU14 l)-FY(14))/(T(14-I)-T(14)),'FY
OMKOMK (3) uDMOM( 2)
FLOW. (FM( 14-I )-FM( 14) 1/ITI14-I )-T( 14))
GO TO 13

09400
09410
09420
09430
09440
0945G
09460
09470
09480
09490
09500
09510

09520
09530
09540
09550
09560

095 7 U
09580
Q9590
09600
OiI
09620
09630
09640
09650
09660
09610
09080
C4690
09700
09720
09730
09740
09750

33 TME(1,.FX(NT)

09760

TME (2)-FYi(NT)
TME (3) sTS4 (NT)
TME(4 )-CG(NT)
TME (5)aFIM (NT)
IF (KAN.EO.Z2
GO
GO TO 34
40 TME(I)-FX(1)

09770
09780
09790
09800
09820
09830
098'40

TO 201

60

TME(2!=FY(1!
TME(3)=TSL(1)
TME(4)2CG11)
TME(5)=FM(1!
201 DO 200 121.3
200 DMI(0MK(1!=0.
RETURN
END
SISFTC FAZIN
LIST*REFvDECK
SUEROUTINE PHASIN(JI
C
C
IN ALL COMMON STATEMENTS D.ORIVERA=AOPTUN.B=FACTEEVALC=Cli(.
C
R2=RKERR2.R2RKG.T2TABL.F=PHASINPoPINTRJ.XZXOUTANO L=TIZLP

siNT

C
COMMON/DCF/P I *TBO
C
COMMON/DF/

TMI!

C
COMMON/CF/ TFP
C
COMMON/FRF2/EPT W''* EPB I
C
COMMON /FT,'
PMC2O!,CX(20),CNA(20!.r MA(2O!,CMQ(2O!,T(40!.FX(4O!.
IFY(40)*TSL(40).CG(40)vFMC4O!.
NMeNT.CND(20!.CNP(20)
C
COMMON /PF/
TBOB(10),AEB(I0.#DSQB(10OhREFLB(1O!.TMIIB(10).
IPIBI 10! gEPTINS(I !,EPBIGB( IO!,TFPB(1O!.KANS(IO!.
NMB( 10!
2.PMB(20. 10),CXB(20,1O!.CNAB(2O.1O!.CMAB(20,10).CM0B('O,910),
3NTB(1O! .TB(40.10) .FXB(40.14J).FYB14Q. 10) .TSL.B(40. 10) .Ct.o40,10! *
4FM8(4O.20)*XLT(1 !,CNDB(20.10).CNP5i2QtO.1!,INCB(IO)
C

09990
10000
10010
10020
10030
10040
10050
10060
10070
10080
10100
10110
10120
10130
10140
10160
10170
10180

COMMON/EFT/KAN
C
COMMON/EF/

09850
09860
('9670
09880
09900
09910
09920
09930
09940
09950
09960
09970
09980

REFLiAEgOSOXLNTHPINC

10200
10210
10220
10230
10240
10250

TSCaTBOB 14!
AE&AEB (4!
osomosoB (4!
REFL2REFLP (J!
IKANaKANH I4
TMI IuTMI 18(J)
P1 uP 8(J)
EPS IGuEP51GS J)
EPT IN~EPT INB (4!
71 NC=F INCB (4!
TFPuTFPB (4!
XLNTM.XLT (4!
NMuNMB (J!
DO 200 IsIgNN
PM II! PM5 CIq !
CX (H SCXB (1,4!
CNA (I ) wCFNA8( I o J!
CMA( I )CMAS(
1.4!

10260
10270

10280
10290
10300
'40315
10320
10330
10340
10350
103,0

10370
61

CMO I )UCMOB4 I.jiI

CNB(I l*CNBII$Jl
200 CONTINUE
IF (KAN.EO.2) GO TO 202
201 NT.Nre(j1
D0 203 Is1.NT
Tit)-*TB (1 *J)
FXC(I) -FXB CI.*J)
FYI!) -FYB CI.*J)
TSLCI -TSLBCl.J)
CG CII-CGB (I *.Ji
FM (II
FMB E.lJI
203 CONTINUE
202 RETURN
END
SIOFTC PINTR
LIST.REF.DEC(
SUBROUTINE PINTRJ
C
1
C
IN ALL COMMON STATEMENTS D.DR1VER.AsA0PTUN.JaBFACTEuEVAL.CCI~e 1N
L
R2=RKERR2.R=RK.TTABLF.PHASNP=PINTRJ.X-XUT.AND
L-TI2LP
C
COMMON/PT/ N.NH(44).Rl0(44) .V5(44).PA(44)
C
C0MMON/DAP/TH.AL.XWANT.YWANT.0ELC.OELT.OELH.T0WTlL
C
COMMON/Z;AEP/NTYPE
C
COMMON/OBTP/NW.*HPi50)
C
COMMON/DCP/ TO
C
C0MMON/A'
XFtRST.YFIRST.XCHANG.YCHANG.XLAST.YLAST.NBF.aF(50).
IDELTA
C
COMMON/DP/ UO.V0.WO.ZL.XLAT.NFAZE.N1ST.IF0UT
C
COMMON/ABETP/ WXP(50) .WYP(50)
C
COMMON/EP/
MAL
C
COMMON/PC#' ZIM
C
COMMON/89X/
IDEN
C
COMMON /PF/
TBOB(10).AEB(I0IDS08(I0).REFLBIO).TMIIB(10).
NMa(IO)
IPIBE IO).EPTINB(1 ),EP8IGB( 10).TFP8( 10) eKAN8(IO)9
2.PMB(20910).CXBC2O. I)O.CNAB(2C,10).CMABC20. IO).CMCIB(20,10).
3NT8(11.840.1O).FXB(40.10)FY3140. 10I.TSLB(40.101.CGBI4O.10).
4FM814O.20).XLT( 1 1.CN0dI2O. 10) ,CNPB(2O.10).FINC8CIO)
C
DMACH(6).DTIME(5IlDEN(I. j.CNO(20).CNPC2OI
DIMENSION
t'CMA(20).CMQ(20).T(40).
PM(20hCX(20).CN'4!1,
UIMENSION
.
IFX(40).FY(40).TSL(40),CGI40).F-l

10380
i0390
.0400
01
10420
10430
10440
10450
10460
10470
10480
10500
10510
10520
10530
10540
10550
10!D60
10570
10580
10590
10600
10610
10620
10630
10640
10650
10660
10670
10680
10690
10700
10710
10720
10730
10740
1075o
10760
10770
10780
10790
10800
10810
10820
10830
10840
10850
10870
10890

10910

C NTYPE -TYPE OF OUTPUT DESIRED

!-1.REGULAR TRAjl*(u2oINTERATESI
C
f=3.ANGULAR OUTPUT).(-4.BALL FACTORS).(85PAAMETER VAR1A8ILITY*0)
C IFOUT=I CONJTINUE PRINTING AT ENO OF PHASE-NO SKIP TO NEXT PAGE
C
I READIS. IUOO)NT'ePE, IFOUT
READ(5#150OO1DEN
WPITE (6.*499)1 DEN
C
C
NFAZE IS THE TOTAL OF THE PHASES-- NBST IS PHASE TO PICK BOOSTER UP AT
READ (5.2000) TO#ZL.ZIM.XLAT.XLONG.NFAZE.NBST
WR~ITE (6.500) TO.ZLZIMXLAT.XLONG.NFAZE.NBST
J=1
C
C MAL-I
C

INPUT

CP

MAL=2

INPUT CMA

10920
10930

10940
10950
10960
10970
10980
10990
11000
11020
21030
11040
11050
11060

READ(5s7G0O )U-VO9WOMAL

11070

WRITEi6.5IO)UO.V0.WO.MAL
READ(5.3600)TH.AL.XWAN!.YWANT.0ELTA.TOWTIL
WP1TE(6.52O)TH.AL.XWANT.YWANT.0ELTA.TOWTIL
REAO(5.2500)N,(HH(I),RHO(I).VS(I).PA(UIa1lx.N)

11090
21100
11110
11120

IF(*ZTYPE.NE.3)G0 TO 2
READ(5.3UOOIXFIRST.XLASTsXCHANG.YFIRST.YLAST.YCHANG
WPITE(6.540)XFIRST.XLASTqXCHANGYFIRST.YLAST*YCHANG
GO TO 12
2 5REAO(5.3500)NW.2HP(I).WXP(I).WYP(1).1.1.NW)
IF(NTYPE.NE.2)GO TO 4
REA0(5.4OOO)NBF.I8F(I)o.1.lBF)
WRITE(6.560)NBF9.BFU).=1.NBF)
12 READ(5s3Q)DELC.DELTDELI
WRITE (6.570 )DELC.OELT*OELH
c
C TFP-TIME OF FIRST PRINT
KANz1.MACH AND TIME TA8LES
KAN=2vMACH ONLY
C IOPTUNw30 TO GROUND
C
4 READ (5.4500) TBO.AE.DSO.REFLsX..ENTH.EPTINYEPaIGFINC.TMII.PI.TFP
2 ,KAN
WRITE (6.580) TBO.AE.OSGREFL.XLENTH.EPTINY.EPBIG.FINC.TMI I PI .TFP
I ,KAN
READ (595000) NM.(PMII),CX(I).CNA(l),CMA(I).CMQtI).CNO(I,.CNPII ).Z
12 *NM1
WRITE (6.590) NM. (PM( I).CX( 2.CNA(I ).CMA( 2).CMO(I ).CND( I )CNP(I 2.1
1.1 .NMI
tF(NTYPE.NE.5)GO TO 7
READ(b.5500) (OMACH( I)o 1*1.6

11140
11150
11160
11170
11160
11200
11210
11220
11230
11240
11250
11260
11270
11280

I
.

11330

WRITE(6.600)(OMACH(2).t1,6)I
DO 8 KsINM
CX(W)=CX(K)*DMACH(l)
CNA(K )=CNA (K)*DMACH(2)
CMA (K )=CMA (K )*DMACHC3)

11360t
21370
11380

6311390

If

CM (K ) CMQ (K I*OMACN (4 2
CNO(K )-CNDCKI*ONACN(51
* CNP(K,CNPCK*P4ACH(6)
7 T8OBIJI-TSO
AE8 (J IAE
OSOs (j 2 S0
REFLA(J ) REFL
TMI !R(J).TMi 1
P1 ( J 2 P
EPTINB(J.-EPTINY
EPOIGO(J.zEPRIG
FINcntj).FtqC
TPIJ
ITFP
KANBa J )-KAN
XLT 1!) XLENTH
NMS(J)I-NM
00 10 l.1.N4
PMB (I ,J 1PM CI)
CXB II JloCX( 1
CNAB (1. J)2CNA (2)
CMA6 I .JI.CMA( II
CMOB( I.JI-CMO( I
CNDSB( I .J I cNo( I)
%0 CNPB(IojlaCNPII)
IF(KAN.NE.lI)GO TO 6
REAO(5.6000INT9(T(l).FXCIIFY(12,T3L(1 2.CG(II.FMCII.
WRITE(6.610)NT,(TII).FX(I),FY(1).TSL(I 2.CG(I1,FM(12.
IF(NTYPE.NE.S)GO TO 9

15
9

20
6

499
500
510
520

WR1TE(6,620)(DTIME(I2.1=1,5)
D0 15 K1.*NT
FX(KI.FX(K)*DTIME(1I
FY(K)AFY(K)*DTIME(2)
TSL(K 2-T-SL(K)*OTIME(3)
CG(K)ICG(K IDTIME(4)
FM(K1.FM(Kl*0TlME(5)
CONTINUE
NTO(J)=NT112
00 20 I=1.NT
TB(1 .J)-T(I
FX8( I.Ji.FXC II
FY8 I .J)vFY(1)
TSLB(1.JI3TSLiI)
CG8( I .JIUCG( I)
FM8 (I 'J)IFM (12
CONTINUE
J-J+I
IF(J.LE.NFAZEIGO TO 4
RETURN
FOPMAT (IHI.I OX. 1246)
FORMAT (11X.F11.2.2FI1.2.2F11.5.7X.14,7X,14)
FOPMAT(IIX.3F11. .7X*14)
FGPMAT(IIX.4F11.2.FlI.O.FIl.1)
64

11410
11420
11430
11440
11450
11460
11470
11480
11490
11500
11515
11520
11530
11540
11550
11560
11570

11,.NT)
Izl.NT)

11590
11600
11610
11620

11650
11660
11670
11680
11690
11700

11730
11740
11750
11760
11770
11780
11790
11810
11820
11630
11840
11850
11860
11870

530
540
550
560
570
580
5'J

F04RMAT5.!42XF'9.0.2X.Fl1.9.2F~l.3/(11X.F9.O.2X.Fl1.92F11.3))
FCORMATl11Xe6F11. )

FORMAT(5X,14,2X.3FI1.2/11lX.3FI1.2))
FORMAT(5X. 14,2XFI I.6/( I X.F i .6))
FOPMAT(llX*3F11.5)
FORMAT (1liXF11.2.2FlI1,7.4FI1 14 9 FI 1.8.2FI1 .4/51 1.4.11)
FORMAT (SX. 14.2X.FlI .3.Fli.5.F11.2.FI 1.3,, ,1.2,FI1,3,FIl.3/Cl1XFI
11 .3,F1 1 *,FI I,2.F1 1.3.Fi I *F1 1 3.F1 1.311
600 FORMAT (I1X.6F!1.2)
610 FORMAT(5X.14.2X.Fl.2.F11.3tFl1.1.F11.2,F11.4,F11.3
620
1000
1500
2000
2500
3000
3500
4000
4500
5000
5500
6000
6500
7000

FORMAT(lI1X6FlI.3)
FORMAT(8Xo2(11-7X) 1
FORMAT(12A6)
FORMAT (8Xt5F8*0%2(I2.6X))
FORM-T(6Xgl2*4F8.O/(8X*4F8.0))
FORMAT(8X96FB*Q)
FORMAT(6X.1243F8.0/(8Xv3F8.O))
FORMAT(6Xsl2,F8. )
FORMAT (10F8.0/F8.0. (11.7x))
FORMAT (6X.I2*7F8.Q/(6Xs7F8.0))
FORMAT (16X96F8.
FORMAT(6X. I246F8.0/(sX.6F8.0))
FORMAT(16Xo6F8.0)
FORMAT(BX*3F8.0*.i1
END
.AST.REF.DECK
SIBFTC XOUTP
SUBROUTINE XOUT {JF.IX)
COMMON Y(22,5)*DY(22.5).R0OF1 (22v5)
C
C
C

IN ALL COMMON STATEMENTS D=0RIVLR.A=AOPTLN.B=FACTEEVAL.C=C1K,1=1NT

KGTTABL.F=PHASIN.PPINTRJ.X=XOUTAND L=T I2LP
R2=RKERR2,R
TH9AH.DC3) .EU3) ,V

C
COMMON/DAETX/ WXS.WYS
C
TIME

COMMON/DIX/
C
COMMON/OX/

IJK

C
COMMON/AE3XL/

XYZ (31

COMMON/BPX/

bDEN

I 1900
11910
11920

119
11990
12000
12-110
12020
12030
12040
12050
12060

12110
12120
12130
1214G
12150
12160
12170
121830
12190
12200
12210
12220
12230
12240
12250
12260
12270
12280
12290
12300
12310
12320

C
COMMON/XL/

11880
11890

C
C
IDEN(12)
DIMENSION
FORMAT(23Xo3F15.2)
FORMAT(22X.ZFI5..2x,2F15.3)
FORMATC(/9XF15.3.FIJ.1.2F15.1 .2X.3F15.3)
FORMAT(IH0O,17X.4HTIME. 12X. 1HX. 13X. 1HY,12X..H-Z. 7X,3HI--. IOX.5HTHET
IA,1IX.5HALPHA/33Xe2HXD,12X.2HNDO,1X,2HZD,16X,3HWYS.1OX.IdV/33X.3HX
200.I1Xo3HYODl0X93HZDD)
699 FORMAT (///4OX40HINTEGRATION INTERVAL IS LESS THAN .00001)
695
696
697
698

65

12330
12340
12350
12360
12370
12380
12390

IF(IX.EOQ.1)GO TO 11
CALL T12L'1.S)
I JKM-MODCI JK*13)
T IME-V (1 *JF)
IFC1JKM.GT.O)GO TO 10
WRITE (6,700110EN
700 FORMAT (IH112A6)
WRITE (6,698)
10 WRITE(6,697)TIME.XYZ(1 ).XYZ(2).XYZ(3)1.AJXS.TH,AH
WRITE(6.696)D(1h0D!2),D(3).WYS.V
WRITEC66Q5)BI)9Bf2)*8(3)
GO TO 12
11 WRITE (61699)
12 IJK=IJK+l
RE TURN
END
LISTqREFDECK
SIBFTC T12LP
SUBROUTINE TI2L(KWIND-JF)
COMMON Y(22s5)9DY(22.5),R0FF1 (22.5)
C
C
C

IN ALL COMMON STATEMENTS D=DRIVERA=AOPTUNB=8FACT.E=EVAL.C CHK,1=1NT

P23RKERP2-R=RKG, T=TABL.F=PHASIN.P PINTRJ. x=xour .AND L=T ILP

C
COMMON/XL/

T.HeAH*D(3)*B63),V

C
COMMON/OEL/

SLATCLATSLATGCLATG.TIMEO.w

C
COMMON/OL/ R0
C
COMMON/ABXL-/

XYZ (3)

C
DIMENSION A(393)oYD(6)
C
C IF KWIND-It THE SECOND DERIVATIVES wiLL- BE COMPUTED
C COMPUTATION OF THE ROTATION MATRIX FROM THE LCI TO THE LAUNCHER SYSTEM
C W=ROTAlION OF EARTHS1RSRAOIUS TO LAQNCHER SYSTEMSCLAT-COS(GEOC LATIS
SLATG-SIN((.EOO) LAT) ARE ALL U
CLATG=COS(GEOU LAT)
C SLAT-SIN(GEOC LAT)
C IN DRIVER
C
WT-W*(Y(1,JF)-TIME0)
A(191I= -COS(WT)
A(192)= -SIN(WT)
0.0
A(I.3)
A(293)= CLATG
A(393)= SLATG
AC2. )=-A(3*3)*A4I.2)
A(391)= A(2.31*A(I.2)
A(3s2,)-A(2*3)*A(ItI)
A(292)= A(393)*A(I1
C
VECTOR
C COMPUTATION OF TRANSLATION

12400
12410
12420
!2430
12440
12450
12460
12470
12480
12490
12500
12510
12520
12530
125-40
I 2550
12560
12570
12580
12590
12b00
12610
12620
12630
12640
12650
12660
12670
12680
12690
12700
12710
12720
12730
12740
!2750O
12760
12770
12780
12790
12800
12810
12820
12830
12840
12850
12860
12870
12880
12890
12900
12910
12920

C
XLzRO*CLAT*A( 1.2)

66

12930
12940
12950
12960
12970

YL=-RO*CLAT*A(I,I
ZL-RO*SLAT
C
C COMPUTATION OF POSITION
C

IN LAUNCHER

SYSTEM

12980
12990
13000
13010
13020
13030
13040
13050

8(1 )=Y(17,JF)-XL
B(2)=Y(!8,JF)-YL
B(3)=Y(19,JF)-ZL
CALL MTRXL (A*,.IvXYZ)
C
C RELATIVE MOTION COMPUTATIONS
C
YD(t)=DY(17.JF)+W*Y(18,JF)
YD(21=DY(18JF -W*Y(17JF)
YD(3)=DY(19.JF)
IF(KWINO.NE.I)GC TO 10
YD(4)=DY(209JF1+2.*W*YD(2)+W*W*Y(17.JFI
YD(5)=DY(21,.JF)-2.*W*YD(I)+W*W*Y(18,JF)
YD(6) DY(22.JF)
C
C COMPUTATION OF VELOCITY AND ACCELERATION

IN LAUNCHER SYSTEM

C
CALL MTRXL(AYDO4),*IB)
10 CALL MTRXL(A*YDI),ID)
IF(KWINO.NE .1)GO YO 20
V=SQPT (D(1)**2+D(2)**2+O(3)**2)
TH=57,2957795*ARSIN(DC3)/V)
2U AH=57.2957795*ATAN2(D(I).D(2))
RETURN
END
SDATA
XNO

Is

AEROBEE NASA 4.195 GAGI

4000.
4143.
TO-Z-LATI*
283.
0.
VEL
2.90
348.
.00210161100.99
ATMOS 404001.
.00207931099.04
4501,
ATMOS
00204811097.09
5001.
ATMOS
ATMOS
6002.
.00198681093.18
ATMOS
7002.
900192681089.25
*0018683108b.31
8003.
AIMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATM05
ATMOS

13060
13070
13080
13090
13100
13110
13120
13130
13140
13150
13160
13170
13180
13190
13200
13210
13220
13230

10005.
12007.
14009.
100'2o
19017.
22023,
25030.
29040.
33052.
37066.

.0017553107739

.00164761069.4
*00154 41061.36
.00144741053o25
1040.97
.00131
.00118271028.55
*00106511015*9a
.0009225998*96
.000795 981.65
.0006759968.08

32.40392!.8B58
11
1821,09
1787.51
1754.34
1689.76
1627-N03
1566,21
144871
1341.01
1236.69
1142.77
1010.27
890.17
782.48
655.2
545,.24
450.8

67

ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
.7MOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS
ATMOS

41081.
45097.
50120.
55145o
60173.
65203.
70236.
75271.
80308.
85348.
90390.
100482.
120695.
140946.
160222.
180550.
200917.
220303.
240748.

90005577968.08
o00046 1968.08
.00037 4968.08
.0002845968.08
.0002237968.08
e000176 968.08
*0001384971.06
*0001063974,44
o0000856977.82
.0000665981.19
.0000515984,54
o0000314991.21
.00001231022.01
.00000521057.52
.00000241082.02
.00000121071.74
.00000 61046.84
.00000 3999.33
.00000 1945.08

260720. 0.
280216. 0*
299233. 0.
506000. 0.
1000000.0.
264196.
3.0
4325.
2.0
4500.
2.5
4700.
2.5
5000.
5400.
5800.
6250.
6750.
7250.
7750.
8500.
11500.
16500.
21500.
26500.
31500.
36500.
41500.
46500.
91500.
9900n.
69000.
79000.
89000.

2.5
MACH

1.0
1.0
1.5
0.0
1.5
5.0
5.5
6.0
7.0
0.0
-4.0
6.5
9.5
7,0
0.0
-9.5
-8.5
10.0
4.0
10.0
14.0
0.0
99000.
.296871 1.562t
1 60.0

.389

889.
884.
892.8
922.7
922.7
6,5
5.5
8.5
8.0

3-1*95
306.9
241.34
189.77
149.24
117.36
92.14
72.89
57o46
45.44
35.99
22.68
9.28
4,03
1.92
.88
.4
.17
.067
.02
.008
.001
0.
0.

6.5
6.5
7.0
3.5
2.5
3.
2.
3.5
2.5
18.0
21.5
18.0
13.5
12.0
16.0
10.5
10.0
1.5
.0
1.5
2.5
13.0
1.25

31.2142

.0001

23.31

26.15

-1300.
68

.001

1.

1.

.389
0389
.400
.557

21.50
22.10
23o09
24.87

26.5125 -1355.
26.68125-1425.
26o8125 -1510*
-1650.
27.35

2.0
I
1 40.
1.
1
2.
1
2.5
1
.296875

.558
14.2
13.65
13.18
12.8
1.5625

18.45
4350.
3980.
3620.
3425.
1.25

27.35
20050.
20050.
20050.
20050.
28.2142

217.5
.75
2
1.
2
1.15
2
2.0
2
2.5
2
3.
2
3.5
2
4.
2
4.5
2
5.0
2
5.5
2
6.0
2
6.5
2
7.0
2
7.5
2
8.
2
2102.5
5o
2

.28
.293
.445
.657
.438
.3825
.338
.302
.271
.244
.220
.206
.194
.185
.181
.180
.176
11.2
10.9

16.0428
17.1887
19.3659
19.8243
12.3071
10.1127
8.9095
8.1360
7.3625
6.7895
6.3885
6.0447
5.8155
5.5577
5.3572

-1300.
24.075
-1380.
24.275
-1520.
25.2
25.6375 -1650.
23.7875 -1148.
2Z.6875 -977.
21.8375 -912.
21.0875 -860.
20.33125-828o
-805.
19.6
18.9125 -786.
-772.
18.4
-762.
17.925
17.4375 -754.
17.1875 -749.

MACH
MACH
MACtMACH

1
I
1
1

MACH
TIME
TIME
TIME
TIME
51.8
MACH
MACH
MACH
MACH
MACH
MACH
MACH
MACH
MACH
MACH
MACH
MACH
MACH
MACH
MACH
MACH
MACH
TIME
TIME
TIME
TIME
TIME
TIME
TIME
TIME

2
2
2
2
2
2

TIME 2
TIME 2
10(0

.25
.5
.75
1.0

-1650.
21.5
21.042
20.5
20.25
.0001

5.1566
4.9274
1750.
1740.

17.0
16.85
4100.
4100.

-747.
-745.
17.8333349.91629
17.9083348.30503

4100.
4100.
4100.
4100.
4100.
4100.
4100.
4100.

18.0

.001

10.
15.
20.
25.
30.
40.

0o.35
9.78
9.2
8.64
8.07
6.93

50.
51.8

5.8
5.76

1725.
1710.
1690.
1660.
1620.
1460.
1120.
1000.

1.5625

1.25

28.2142

.475
.414
.365
.325
.289
.259
.235

12.3071
10.1127
8.9095
8,1360
7.3625
6.7895
6.3885

23.7875 -1148.
22.6875 -977.
21.8375 -912.
21.0875 -d60.

5.5
6.
6.5
7.
7.5
8.

.216
.200
.190
.186
.185
.184

6.0447
5.8155
5.5577
5,3572
5.1566
4.9274

3
3
3
3
3
3

l.

45.08253

17.9583341.86003
17.8416738.63752
17.6666735.41502
17.3916732.19252
16.4833325.74751
14.4666719.30250
18.14240
14.15

2
MACH 3132.
2.5
MACH 3
3.
MACH 3
3.5
MACH 3
4.
MACH 3
4.5
MACH 3
5.
MACH 3
MACH
MACH
MACH
MACH
MACH
MACH

69.72185
66.16309
62.28208
60.34158
.001

20.33125-828.
-805.
19.6
18.9125 -786.
18.4
17.925
17.4375
17.1875
17.0
16.85
69

-772.
-762.
-754o
-749.
-747.
-745.

.01

2.

4.

REFERENCES

1. Duncan, Louis D., and Ronald J. Ensey, "Six Degree of Freedom

Digital Simulation Model for Unguided Fin-Stabilized Rockets,"
Atmospiteric Sciences Laboratory, White Sands Missile Range,
Ntw Mexico, 1964.
2. Ralston, A:thony and Herbert S. Wilf, Mathematical Methodq for
Digital Computers, John Wiley and Sons, New York, 1960.
3. Duncan, Louis D.,and Bernard F. Engebos, "Techniques for Computing Launcher Settings for Unguided Rockets," Atmospheric
Sciences Laboratory, White Sands Missile Range, New MIexico,
1966.
4. Engebos, Bernard F.,and Louis D. Duncan, "A Nomogram for Field
Determination of Launciier Angles for Unguided Rockets," Atmospheric Sciences Laboratory, White Sands Missile Range,
New Mexico, 1966.

70

ATMOSPHERC SCIENCES RESEARCH PAPERS

1. Webb, W.L., "Development of Droplet Size Distributions in the Atmosphere," June
1954.
2. Hansen, F. V., and H. Rachele, "Wind Str, cure Analysis and Forecasting Methods
for .Iockets," June 1954.
3. Webb, W. L., "Net Electrification of Water Droplets at the Earth's Surface," J. Meteorol., December 1954.
4. Mitchell, R., "The Determination of Non-Ballistic Projectile Trajectories," March
1955.
5. Webb, W. L., and A. McPike, "Sound Ranging Technique for Determining the Trajectory of Supersonic Missiles," #1, March 1955.
6. Mitchell, R., and W. L. Webb, "Electromagnetic Radiation through the Atmosphere," #1, April 1955.
7. Webb, W. L., A. McPike, and H. Thompson, "Sound Ranging Technique for Determining the Trajectory of Supersonic Missiles," #2, July 1955.
8. Barichivich, A., "Meteorological Effects on the Refractive Index and Curvature of
Microwaves in the Atmosphere," August 1955.
9. Webb, W. L., A. McPike and H. Thompson, "Sound Ranging Technique for Determining the Trajectory of Supersonic Missiles," #:3, September 1955.
10. Mitchell, R., "Notes on the Theory of Longitudinal Wave Motion in the Atmosphere," February 1956.
11. Webb, W. L., "Particulate Counts in Natural Clouds," J. Meteorol., April 1956.
12. Webb, W. L., "Wind Effect on the Aerobee," #1, May 1956.
13. Rachele, H., and L. Anderson, "Wind Effect on the Aerobee," #2, August 19,56.
14. Beyers, N., "Electromagnetic Radiation through the Atmosphere," #2, January 1957.
15. Hansen, F. V., "Wind Effect on the Aerobee," #3, January 1957.
16. Kershner, J., and H. Bear, "Wind Effect on the Aerobee," #4, January 1957.
17. Hoidale, G., "Electromagnetic Radiation through the Atmosphere," #3, February
1957.
18. Querfeld, C. W., "The Index of Refraction of the Atmosphere for 2.2 Micron Radiation," March 1957.
19. White, Lloyd, "Wind Effect on the Aerobee," #5, March 1957.
20. Kershner, J. G., "Development of a Method for Forecasting Component Ballistic
Wind," August 1957.
21. Layton, Ivan, "Atmospheric Particle Size Distribution," December 1957.
22. Rachele, Henry and W. H. Hatch, "Wind Effect on the Aerobee," #6, February
1958.
23. Beyers, N. J., "Eromagnetic Radiation through the Atmosphere," #4, March
1958.
24. Prosser, Shirley J., "Electromagnetic Radiation through the Atmosphere," 45,
April 1958.
25. Armendariz, M., and P. H. Taft, "Double Theodolite Ballistic Wind Computations,"
June 1958.
26. Jenkins, K. R. and W. L. Webb, "Rocket Wind Measurements," June 1958.
27. Jenkins, K. R., "Measurement of High Altitude Winds with Loki," July 1958.
28. Hoidale, G., "Electromagnetic Propagation through the Atnmsphere,'" #6, February 1959.
29. McLardie. M., R. Helvey, and L. Traylor, "Low-Level Wind Profile Prediction Techniques," --1, June 1959.
30. Lamberth, Roy, "Gustiness at White Sands Missile Range." #1, May 1959.
31. Beyers, N. J., B. Hinds, and G. Hoidale, "Electromagnetic Propagation through the
Atmosphere," #7, June 1959.
32. Jieyers, N. J., "Radar Refraction at Low Elevation Angles (U)," Proceedings of the
Army Science Conference, June 1959.
33. White, L., 0. W. Thiele and P. H. Taft, "Summary of Ballistic and Meteorological
Support During IGY Operations at Fort Churchill, Canada," August
1959.
34. Hainline, D. A., "Drag Cord-Aerovane Equation Analysis for Computer Application,"
August 1959.
35. Hoidale. G B., "Slope-Valley Wind at WSMR," October 1959.
36. Webb, W. L, and K. R. Jenkins, "Hiogh Altitude Wind Measurznents," J. MeteoroL, 16, 5, October 1959.

37. White, Lloyd, "Wind Effect on the Aerobee," #9, October 1959.
38. Webb, W. L., J. W. Coffman, and G. Q. Clark, "A High Altitude Acoustic Sensing
System," December 1959.
39. Webb, W. L, and K. R. Jenkins, "Application of Meteorological Rocket Systems,"
J. Geophys. Res., 64, 11, Novemnber 1959.
40. Duncan, Louis, "Wind Effect on the Aerobee," #10, February 1960.
41. Helvey, R. A., "Low-Level Wind Profile Prediction Techniques," #2, February 1960.
42. Webb, W. L, and K. R. Jenkins, "Rocket Sounding of High-Altitude Parameters,"
Proc. GM Rel. Symp., Dept. of Defense, February 1960.
43. Armendariz, M., and H. H. Monahan, "A Comparison Between the Double Theodolite and Single-Theodolite Wind Measuring Systems," April 1960.
44. Jenkins, K. R., and P. H. Taft, "Weather Elements in the Tularosa Basin," July 1960.
45. Beyers, N. J., "Preliminary Radar Performance Data on Passive Rocket-Borne Wind
Sensors," IRE TRANS, MIL ELECT, MIL-4, 2-3, April-July 1960.
46. Webb, W. L, and K. R. Jenkins, "Speed of Sound in the Stratosphere," June 1960.
47. Webb, W. L, K. R. Jenkins, and G. Q. Clark, "Rocket Sounding of High Atmosphere Meteorological Parameters," IRE Trans. Mil. Elect., MIL-4, 2-3,
April-July 1960.
48. Helvey, R. A., "Low-Level Wind Profile Prediction Techniques," #3, September
1960.
49. Beyers, N. J., and 0. W. Thiele, "Meteorological Wind Sensors," August 1960.
50. Armijo, Larry, "Determinaticn of Trajectories Using Range Data from Three Noncolinear Radar Stations," September 1960.
51. Carnes, Patsy Sue, "temperature Variations in the First 200 Feet of the Atmosohere in an Arid Region," July 1961,
52. Springer, H. ., and R. 0. Olsen, "Launch Noise Distribution of Nike-Zeus Missiles." July 1961.
53. Thkeo. 0. W., "Density and Pressure Profiles Derived from Meteorological Rocket
Measurements," September 1961.
54. Eqamond. M. and A. B. Gray, "Accuracy of Missile Sound Ranging," November
1961.
55. Lamberth, R. L and D. R. Veith, "Variability of Surface Wind in Short Distances,"
#1, October 1961.
56. Swanson, R. N., "Low-Level Wind Measurements for Ballistic Missile Application,"
January 1962.
57. Lamberth. R. L and J. H. Grace, "Gustiness at White Sands Mis.ile Range," #2,
January 1962.
K8. Swanson, R. N. and M. M. Hoidale, "Low-Level Wind Profile Prediction Techniques," #4, January 1962.
59. Rachele, Henry, "Surface Wind Model for Unguided Rockets Using Spectrum and
Cron Spectrum Techniques," January 1962.
60. Rachel., Henry, "Sound Propagation through a Windy Atmosphere," #2, February 1962.
61. Webb, W. L, and K. R. Je-kins, "Sonic Structure of the Mesoephere," J. Acous.
Soc. Amer., 34, 2. February 1962.
62. Tourin, M. H. and M. M. Hoidale, "Low-Level Turbulence Characteristics at Vhite
Sands Missile Range," April 1962.
63. Miers, Bruce T., "Mesospheric Wind Reversal over White Sands Missile Range,"
March 1962.
64. Fisher, E., R. Lee and H. Rachele, "Meteorological Effects on an Acoustic Wave
within a Sound Ranging Array," May 1962.
65. Walter, E. L., "Six Variable Ballistic Model for a Rocket," June I'A2.
66. Webb, W. L, -:. et&Jed Acoustic Structure Above the Tropopt ,," J. Appied Meteoro/, 1, 2, June 1962.
67. Jenkins, K. R., "Empirical Comparisons of Meteorological Rocket Wind Sensors," J.
Appl. Meteor., June 1962.
68. Lamberth, Roy, "Wind Variability Estimates as a Function of Sampling Interval,"
July 1962.
69. Pachele, Henry, "Surface Wind Sampling Periods for Unguided Rocket Impact Prediction," July 1962.
70. Traylor, Larry, "Coriolis Effects on the Aerobee-Hi Sounding Rocket," August 1962.
71. McCoy, J., and G. Q. Clark, "Meteorol ,gical Rocket '1hermometry," August 1962.
72. Rachele, Henry, "Aaal-Time Prel rich Impact Prediction System," August 1962,

7,.

Beyers, N. J., 0. W. Thiele, and N. K. Wagner, "Performance Characteristics of

Meteorlogical Rocket Wind and Temperature 7ensors," October 1962.
74. Coffman, J., and R. Price, "Some Errors Associated with Acoustical Wind Measurements through a Layer," October 1962.
75. Armendariz, M., E. Fisher, and J. Serna, "Wind Shear in the Jet Stream at WSMR," November 1962.
76. Armendariz, M., F. Hansen, and S. Carnes, "Wind Variability and its Effect on Rocket Impact Prediction," January 1963.
77. Querfeld, C., and Wayne Yunker, "Pure Rotational Spectrum of Water Vapor, I:
Table of Line Parameters," February 1963.
78. Webb, W. L., "Acoustic Componeni of Turbulence," J. Applied Meteorol., 2, 2,
April 1963.
79. Beyers, N. and L. Engbergq, "Seasonal Variability in the Upper Atmosphere," May
1963.
80. Williamson, L. E., "Atmospheric Acoustic Structure of the Sub-polar Fall," May 1963.
81. Laimberth, Roy and D. Veith, "Upper Wind Correlations in Southwestern United
States," June 1963.
82. Sandlin, E., "An analysis of Wind Shear Differences as Measi' by AN/FPS-16
Radar and AN/GMD-1B Rawinsonde," August 1963J.
83. Diamond, M. and R. P. Lee, "Statistical Data on Atmospheric Design Properties
Above 30 km," August 1963.
84 Thiele, 0. W., "Mesospheric Density Variability Based on Recent Meteorological
Rocket Measurements," J. Applied Meteorol., 2, 5,October 1963.
85. Diamond, M., and 0. Essenwanger, "Statistical Data on Atmospheric Design Properties to 30 kin," Astro. Aero. Engr., December 1963.
86. Hansen, F. V., "Turbulence Characteristics of the First 62 Meters of the Atmosphere," December 1963.
87. Murris, J. E., and B. T. Miers, "Circulation Disturbances Between 25 and 70 kilomaters Associated with the Su..den Warming of 1963," J. of Geophys.
es., January 1964.
88. Thie!e, 0. W., "Some Observed Short Term and Diurnal Variations of Stratospheric Density Above 30 km," January 1964.
89. Sandlin, R. E., Jr. and E. Armijo. "An Analysis )f AN/FPS-16 Radar and AN/
GMD-1B Rawinsonde Data Differences," January 1964.
90. Miers, B. T., and N. J. Beyers, "Rocketsonde Wind and Temnerature Measurements Between 30 and 70 km for Selected Stations," J. Applied Meteorol., February 1964.
91. Webb, W. L., "The Dynamic Stratosphere," Astronautics and Aerospace Engineering, March 1964.
92. Low, R. D. H.. "Acoustic Measurements of Wind through a Layer," March 1964.
93. Diamond. M., "Croas Wind Effect on Sound Propagation," J. Applied Meteorol.,
April 1964.
94. Lee, R. P., "Acoustic Ray Tracing," April 1964.
95. Reynolds, R. D., "Investigation of the Effect of Lapoe Rate on Balloon Ascnt
,
May 1964.
96. Webb, W. L., "Scale of Stratospheric Detail Structure," Space Reseach V, May
Iiro4.
97. Barber, T. L., "Proposed X-Ray-Infrared Method for Identification of Atmospheric Mineral Dust," June 1964.
98. Thiele, 0. W., "Ballistic Procedures for Ungunided Rocket Studies of Nuclear Environments kU)," Proceedings of the Army Science Conference, June 1964.
99. Horn. J. D., and E. J. Trawle, "Orogra hic Effects on Wind Variability," July 1964.
100. Hoidale, G.. C. Querfeld, T. Hall, and W. Mirele, "Spectral Trananissivity of the
Earth's Atmoophere in the 250 to 500 Wave Number Interval," #1,
September 1964.
101. Dunck,, L D., R. Ensey, and B. Engebos, "Athena Launch Angle Determination,"
,September 1964.
102. Thiele, 0. W., 'Feasibility Experiment for Measuring Atmospheric Density Through
the Altitude Range of 60 to
KM Over White Sans Mi
Rae ,
R
October 1964.
103. Duncan, Li D., and R. Ensey, "Six-Degroe-of-Freedom Digital Simulation Model for
Unguided, Fin- Stabilired Rockets," Novsmber 1964.

104.

Hoidale, G., C. Querfeld, T. Hall, and R. Mireles, "Spectral Transmissivity of the

Earth's Atmosphere in the 2.50 to 500 Wave Number Interval," #2,
November 1964.
105. Webb, W. L., "Stratospheric Solar Response," J. Atmos. Sci November 1964.
106. McCoy, J. and G. Clark, "Rocketsonde Measurement of Stratospheric Temperature,"
December 1964.
107. Farone, W. A., "Electromagnetic Scattering '-nm Radially Inhomogeneous Spheres
as Applied to the Problem of Clear Atmospht-e Radar Echoes," December 1964.
108. Farone, W. A., "The Effect of the Solid Angle of Illumination or Observation on the
Color Spectra of 'White Light' Scattered by Cylinders," January 1965.
109. Williamson, L E., "Seasonal and Regional Characteristics of Acoustic Atmospheres,"
J. Geophys. Res., January 1965.
110. Armendaiz, H., "Ballistic Wind Variability at Green River, Utah," January 1965.
111. Low, R. D. H., "Sound Speed Variability Due to Atmospheric Composition," January 1965.
112. Querfeld, C. I., 'Mie Atmospheric Optics," J. Opt. Soc. Amer., January 1965.
113. Coffman, J., "A Measurement of the Effect of Atmospheric Turbulence on the Coherent Properties of a Sound Wave," Jar",tary 1965.
114. Rachele, H., and D. Veith, "Surface Wind Sampling for Unguided Rocket impact
Prediction," January 1965.
115. Ballard, H., and M. Izquierdo, "Reduction of Microphone Wind Noise by the Generation of a Proper Turbulent Flow," February 1965.
116. Mireles, &, "An Algorithm for Computing Half Widths of Overlapping Lines on Experimental Spectra," February 1965.
117. Richart, H., "Inaccuracies of the Single-Theodolite Wind Meisuring System in Ballistic Application," Febrary 1965.
118. D'Arcy, M., "Theoretical and Practical btudy of Aerobee-150 Ballistics," March
1965.
119. McCoy, J., "Improved Method for the Reduction of Rocketsonde Temperature Data," March 1965.
120. Mireles, R-, "Uniqueness Theorem in Inverse Electromagnetic Cylindrical Scatterin ," April 1965.
121. Coffmnan, J., The Focusing of Sound Propagating Vertically in a Horizuntally Stratified Medium," April 1965.
122. Farone, W. A., and C. Querfeld, "Electromagnetic Scattering from an Infinite Circular Cylinder at Oblique Incidence," April 1965.
123. Rachele, H., "Sound Propagation through a Windy Atmosphere," April 1965.
124. Miera, B., "Uppr Stratospheric Circulation over Aicension Island,' April 1965.
125. Rider, L., and M. Armendariz, " A Comparison of Pibal and Tower Wind Measurements," April 1965.
126. Hoidale, G. B., "Meteorological Conditions Allowing a Rare Observation of 24 Micron Solar Radiation Near Sea Level," Meteorol. Magazine, May 1965.
127. Beyers, N. J., and P. T. Miers, "Diurnal Temperature Change in the Atnosphere
Between 30 and 60 km over White Sands Missile Range," '. Atmos.
Sci, May 1965.
128. Querfeld. C., and W. A. Farone. "Tables of the Mie Forward Lobe," May 1965.
129. Farone, W. A., Generalization of Rayleigh-Gans Scatteiing from Radially Inhomogenous Spheres," J. Opt. Soc. Amer ,June 1965.
130. Diamond. M., "Note on Mepob,gheric Winds Above White Sands Mis;le Range," J.
Applied Meteorcl., June 1965.

Clark, G. Q.. and J. G. McCoy. "Measurement of Stratoap'.1eric Ternperature," J.

Applied Meteorol.. June 1965.
132. Hall, T., C. Hoidale, R. Mireles, and C. Querfeld. "Spectral Transmissivity of the
Earth's Atmosphere inthe 250 to 500 Wave Number Ii.terval," .3,
July 1965.
133. McCoy, J., and C. Tate, '"he Delta.T Meteorological Rocket Payload," June 1964.
134. Hern, J. D., "Obetacle Influence in a Wind Tunnel," July 1965.
135. McCoy, J., "An AC Probe for the Measurement of Electron Density and Collision
Frequency in the Lower Ionesphere," July 1965.
135. Miera, B. T., M. D. Kays, 0. W. Thiele and F. M. Newby, "Investigation of Shc'Term Variations of Several Atmospheric Parameters Above 30 KM,"
July 1965.
131.

137.
138.
139.
140.
141.
142.
143.
144.
145.
146.
147.
148.
149.
150.
151.
152.
153.
154.
155.
156.
157.
158.
159.
160.
161.
162.
163,
164.
165.
166.
167.
168.
169.

Sema, J., "An Acoustic Ray Tracing Method for Digital Computation," September
1965.
Webb, W. L., "Morphology of Nocti'ucent Clouds," J. Geophys. Res., 70, 18, 44634475, September 196.5.
Kay3, M., and R. A. Craig, "On the Order of Magnitude of Large-Scale Vertical Motions in the Upper Stratosphere," J. Geophys. Res., 70, 18, 4453-4462,
September 1965.
Rider, L., "Low-Level ,Jet at White Sands Missile Range," Septemoer 1965.
Lam erth, R. L., R. Reynolds, and Morton Wurtele, "The Mountain Lee Wave at
White Sands Missile Range," Bull. Amer. Meteorol. Soc., 46, 10, October 1965.
Reynolds, R. and R. L. Lamberth, "Ambient Temperature Measurements from Radiosondes Flown on Constant-Level Balloons," October 1965.
McCluney, E., "Theor2tical Trajectory Performance of tke Five-Inch Gut. Probe
System," October 1965.
Pena, R. and M. Diamond, "Atmospheric Sound Propagation near the Earth's Surface." October 1965.
Mason, J. B., "A Study of the Feasibility of Using Radar Chaff For Stratospheric
Temperature Measurements." November 1965.
Diamond, M., and R. P. Lee, "Long-Range Atmospheric Sound Propagation,' J.
Geophys. Res., 70, 22. Novemer 1965.
Lamberth, R. L., "On the Measurement of Dust Devil Parameters," November 1965.
Hansen, F. V., and P. S. Hansen, "Formation of an Intermal Boundary over Heterogeneous Terrain," November 1965.
Webb, W. L., "Mechanics of Stratospheric Seasonal Reversals," November 1965.
U. S. Army Electronics R & D Activity, "U. S. Army Participation in the Meteorological Rocket Network," January 1966.
Rider, L. J., and M. Armendariz, "Low.Level Jet Winds at Green River, Utah," February 1966.
Webb, W. L., "Diurnal Variations in the Stratospheric Circulation," February 1966.
Beyers, N. J., B. T. Mie:s, and R. J. Reed, "Diurnal Tidal Motions near the Stratopause During 48 Hours at WSMR," February 1966.
Webo, W. L., "The Stratospie,'ic Tidal Jet," February 1966.
Hall, J. T., "Focal Properties , Plane Grating in a Convergent Beam," February
1966.
Duncan, L. D., and Henry Rachele, "Real-Time Meteorological System for Firing of
Unguided Rockets," Februarv 1966.
Kays, M. D., "A Note on the Comparison of Rocket and Estimated Geostrophic Winds
at the 10-mob Level," J. Appl. Mete ., Februar" 1966.
Rider, L., and M. Armendariz, "A Comparison of Pibal and Tower Wind Measurements,' ,V. Appl. Meteor., 5, February 1966.
Duncan .,D., "Coordinate Transformations in TlRajectory Simulations," February
19b6.
Williamson, L. E., "Gun-Launched Vertical Probes at White Sands Missile Range,"
February 1966.
Randhawa, J. S., Ozone Measuremcnts with Rocket-Borne OZonesondes," March
1966.
Armendariz, Manuel, and Laurence J. Rider, "Wind Shear for Small ThickntS Layers," March 1966.
Low, R. D. H., "Continuous Determination of the Average Sound Velocity over an
Arbitrary Path," March 1966.
Hansen, Frank V., "Richardson Number Tables for the Surface Boundary Layer,"
March 1966.
Cochran, V. C., E. M. D'Arcy, and Florencio Ramirez. "Digital Computer Program
for Five-Degce-of-Freedom Trajectory," March 1966.
Thiele, 0. W, and N. J. Beyers, "Comparison of Rocketsonde and Radic*onde Temp..
eraturee anl a Verification of Computc. Rocketsonde Pressure and Density." April '966.
Thiele, 0. W.. )b.erved Diurnal Oscillations of Pressure mid Density in the Upper
Stratosphere and Lower Mesotphere," April 1966,
Kays, M. D., and R. A. Craig, "On the Order of Magnitude of Large-Scale Vertical
Motions in the Upper Stratosphere," J. Geoph'N.Res.. April 1966.
Hansen. F. V., "The Richa:dson Number in the Planetary Boundary Layer," May
1966.

170.
171.
172,
173.
174.
175.
176.
177.
178.
179.
180.
181.
182.
183.
184.
185.
186.
187.
188.

189.
190.
191.
192.
193.
194.
195.
196.
197.
198.
199.
200.
201.
202.

Mew"The Measurement
Ballard, H. N.,
June 1966. of Temperature in the Stratosphere and
spheHe,"
Hansen, Frank V., "The Ratio of the Exchange Coefficients for Heat and Momertum
in a Homogeneous, i hermally Stratified Atmosphere," June 1966.
Hansen, Frank V., "Comparison of Nine Profile Models for the Diabatic Boundary
Layer," June 1966.
Rachele, Henry, "A Sound-Ranging Technique for Locating Supersnic Misiles,"
May 1966.
Farrt, _, W. A., and C. W. Querfeld, "Electromagnetic Scattering from Inhomogeneous
Infinite Cylinders at Oblique Incidence," J. Opt. Soc. Amer. 56, 4, 476480, April 1966.
Mireles, Ramon, "Determination of Parameters in Absorption Spectra by Numerical
Minimization Techniques," J. Opt. Soc. Amer. 56, 5,644-647, May 1966.
Reynolds, R., and R. L. Lamberth, "Ambient Temperatuiii Measurements from Radiosondes Flown on Constant-Level Balloons," J. Appl. Meteoro,., 5, 3,
304-307. June 1266.
Hall, James T., "Focal Properties of a Plane Grating in a Convergent Beam," Appl.
Opt., 5, 1051, June 1966
Rider, Laurence J., "Low-Level Jet at White Sands Mizsle Range," J. Appl. Meteorol., 5, 3, 283-28', June 1966.
McCiuney, Eugene, "Projectile Dispersion as Cauped by Barrel Displacement in the
5-inch Gun Probe System," July 1966.
Armendariz, Manuel, and Laurence J. Rider, "Wind Shear Calculations for Small
Shear L-vors," June 1966.
Lamberth, Roy L., anci .%nuel Armendariz, "Upper Wind Correlations in the Central Rocky Mountains," June 1966.
Hansen, Frank V., and Virgil D. Lang, "The Wind Regime in the First 62 Meters of
the Atmosphere," June 1966.
Randhawa, Jagir S., "Rocket-Borne Ozonesonde," July 1966.
Rachele, Henry, and L. D. Duncan, "The Desirabilit-- of Using a Fast Sampling Rate
for Computing Vind Velocity from Pilk,-Balloon Data," July 19%.
Hinds, B. D., and R. G. Pappas, "A Comparison of T)ree Methods for the Correction of ,adar Zlevation Aryle Refraction Errors," August 1966.
Riedmuller, G. F., and T. L. Barber, "A Mineral Transition in Atmospheric Dust
Transport," August 1966.
Hall, J. T., C. W. Querfeid, anA G. B Hoidale. "Spectral Traismissivitv of the
Earth's Atmosphere in tle 250 to 500 Wave Number Interval," Part
IV (Final), July 1'5 6.
Duncan, L. D. and B. F. Engebo,

'Techniques for Computing Launcher Settings

for Unguided Rockets," September 1966.

Duncan, L. D., "Basic Considerations in the Development of an Unguided Rocket
Trajectory Simulation Model," S'r+ember 1966.
Miller, Walter B., "Consideration of S)me Proulems in Curve Fitting," September
1966.
Cormak, J. E., and J. D. Horn, "Th, Tower Shadow Effect." August 1966.
Webb, W. L., "Stratospheric Circulation Response to a Solar Eclipse," October 1966.
Kemnedy, Bruce, "Muzzle Velocity Measurement," October 1966.
Traylor, Larry E., "A Refinement Technique for Unguided Rocket Drag Coefficients," Oct-ber 1966
Nusbaum, Henry, "A Reagent for the Simultaneous Microscope Determination of
Quartz and Halides," October 1966.
Kays, Marvin and R. 0. Olsen, "Improved Rocketsonde Parrchute-derived Wind
Profiles," October i966.
Engebos, Bernard F. and Duncan, LAuis D., "A Iomogram for Field Detexmination of Launcher Angles for Unguided Rockets," October 1966.
Webb, W, L., "Midlatitude Clouds in the Upper Atmosphere," November 1966.
Hansen, Frank V., "The Lateral Intensity of Turbulence as a Function of Stability,"
November 1966.
Rider, L. J. and M. Atnendariz, "Differences of Tower and Pibal Wind Prfiles."
November 1966.
Lee, Robert P., "A Comparison of Eight Mathematical Modes for Atmospheric
Acoustical Ray Tracing," November 1966.
Low, R. D. H., et al., "Acoustical and Meteorological Dc Report SOTRAN I and
II," November 1966.

203.
204.
205.
206.
207.
2WS,
209.
210.
211.
212.
213.
214.
215.
216.
217.
218.
219.
220.
221.
222.
223.
224.
225.
226.
227.
228.
229.
230.
231.
232.
23.3.
234.

Hunt, J. A. vnd J. D. Horn, "Drag Plate Bala-Am," December 1W6.

Anmendariz, M., and H. Rachele. "Determixwtion of a Representative Wind Profile
from Balloon Data," December 1966.
Hansen, Frank V.. 'The Aerodynamnic Roughness of the Compiox Terrain of White
Sands Missie Raaga," January 1967.
Morris, James E., "Wind Measurementa in the Subpolar Mesopause Region," January 1967.
Hall, James T., "Attenuation of Millimeter Wavelength Radiation by Gaseou
Water," January 19A7'
Thiele, 0. W.. and N. J. Beyers, "Upper Atmosphere Pressure Measurementa With
thermal Conductivity Gauges," January 1967.
Armendariz, M., and H. Rachele, "Determination of a R p e'ntativ'e Wind Profile
from Balloon Data," January 1967.
Hansen, F. V., "The Aerodynamic Roughnw's of the Complex Terrain of White Sands
Missile Range, New Mexico," January 1967.
D'Arcy, Edward M., "Some Applications of Wind to Unguided Rocket Impact Prediction," -Marh 1967.
Kennedy, Bruce, "Operation Manual for Stratosphere Temperature Sonde," March
1967.
Hoidale, G. R., S. M. Smith, A. J. Blanco, and T. L4. Barber, "A Study of Atincephez.; Dust," March 1967.
Longvear, J. Qi., "An Algorithm for Obtaining Solutions to lApace's nthad Equations," March 1967.
Rider, L J., "A Comparison of Pibal with Raob and Rawin. Wind Measurements,"
April 1967.
Breeland, A. H,, and R. S. Bonner, "Results of Tests Involving Hemispherical Vind
Screens in the Redaction of Wind Noise," April 1967.
Webb, Willis L., and Max C. Bolen, 'The D-rngion Fair-Weathxer Electric Field,"
April 1967.
Xubinski, Stanley F., "A Comparative Evaluation of the Autornz..c Tracking PilotBallooni Wind Measuring System," April 1967.
Miller, Walter B., and Henry Rachele, "On Nonpararnetrc Testing of the Nature of
Certain Time 5S*'s, April 1967.
Hansen, Frank V., "Spacial and Temporal Distribution of the Gradient Richardson
Number in the Surface and Planetary layers," May 1967.
Randhawa, Jagir S., "Diurnal Variation of Ozone at High Altitudes," May 1967.
Ballard, Har-old N., "A Review of Seven Papers Concerning the Measurement of
Temperature in the Stratosphere and Mesosphere," May 1967,
William, Ben H ., "Synoptic Analyses of the Upper Stratospheric Circulation During the late Winter Storm Period of 1966," May 1967.
Horn, J. D., and J. A. Hunt, "System Design for the Atmospheric Sciences Office
Wind Resexarh Facility," May 1967.
Miller, Walter B., and Henry Rachele, "Dynamic Evaluation of Radar and Photo
Tracking Systems' " May 1967.
Bonner, Robert S., and Ralph H. Rohwer, "Acoustical and Meteorological Data Report - SOTRiAN III and IV," May 1967.
Rider, L. J., "On Time Variability of Wind &t White Sands Missile Range, New Mexico," June 1967.
Randhawa, Jagir S., "Mesospheric Ozone Measurements During a Solar Eclipse,"
June 1967.
Beyers, N. J., and B. T. Miers, "A Tidal Experiment in the Equatorial Stratosphere
over Ascension Island (8S) ", June 1967.
Miller, W. B., and H. Rachele, "On the Behavior of Derivative Processes," June 1967
Walters. Randall K., "Numerical Integration Methods for Ballistic Rocket Trajectory Simulation Programs," June 1967.
Hansen, F~rank V., "A Diabatic Surface Boundary layer Model," July 1967.
Butler, Ra'ph L., and James K. Hall, "Comparison of Two Wind Measuring Systems with the Contraves Photo-Theodolite," July 1967.
Webb, Willis L., "The Source of Atmospheric Electrification," June 1967.

Hinds, B. D., "Radar Tracking Anomalies over an Arid Interior Basin," August 1967.
Christian, Larry 0., "Radar Cross Sections for Totally Reflecting Spheres," August
1967.
237, D'Arcy, Edward M., "Theoretical Dispersion Analysis of the Aerobee 350," August
1967.
238. Anon., "Technical Data Package '-r Rocket-Borne Temperature Sensor," August
1967.
239. Glass, Roy I., Roy L. Lamberth, and Ralph D. Reynolds, "A High Resolution Continuous Pressure Sensor Modification for Radiosondes," August 1967.
240. Low, Richard D. H., "Acoustic Measurement of Supersaturation in a Warm Cloud,"
August 1967.
241. Rubio, Roberto, and Harold N. Ballard, "Time Response and Aerodynamic Heating
of Atmospheric Temperature Sensing Elements," August 1967.
242. Seagraves, Mary Ann B., "Theoretical Performance Characteristics and Wind Effects
for the Aerobee 150," August 1967.
243. Duncan, Louis Dean, "Channel Capacity and Coding," August 1967.
244. Dunaway, G. L., and Mary Ann B. Seagraves, "Launcher Settings Versus Jack Settings for Aerobee 150 Launchers - Launch Complex 35, White Sands
Misvile Range, New Mexico," August 1967.
245. Hinds, R. D., "Radar Tracking Anomalies over an Arid Interior Basin," August
1967

235.
236.

DISTRIBUTION LI
ID#/CYS

304/2

Divisioi. (ESTI)
Electronic Systems
J, G. Hanscom Field
Bedford, Massachusetts 01730

305/1

AFCRL (CREW)
L. G. Hanscom Field
Bedford, Mamchusetts 01730

307/1

HQS, Air Weather Service

DEPARTMENT OF DEFENSE
101/20" Defense Documentation Center
ATTN: DDC-IRS
Cameron Station (Bldg. 5)
Alexandria, Virginia 22314
Ofice of Asst Secretary of Defense
102/1

ATTN: AWSAE/SIPD

102/1

(Research and Engineering)

ATTN: Technical Library, RM 3E1065
Washington, D. C. 20301
103/1

Director, Defense Atomic 6apport

AgencyDocument Library Brau-ch
A'Mr:
Washington,D. C. 20301

308/1

309/1

s
seSIl
Att AW
Scott Air Force Base, Illinois 62225
U. S. Air Force Security Service
ATTN: ESD
San Antonio, Texas 78241
Air Proving Ground CTR (PGBPS-12)
ATTN: PGAPI
Eglin Air Force Base, Florida 32542

DEPARTMENT OF THE NAVY

310/1

Headquarters
Research & Technology Divi.ion
ATTN: RTITC
Bolling AFB, D. C. 20332

200/1

Chief Of Naval Research

ATTN: Code 427
Department of the Navy
Washington, D. C. 20325

201/1

Naval Ships Systems Command

ATTN. Code 6312 (Technical Library)
Main Navy Building, Room 1528
Washington, D. C. 20325

400/2

Director
U. S. Naval Research Laboratory
ATTN: Code 2027
Washington, D. C. 20390

Chief of Research ard Development

Department of the Army
Washington, D. C. 20315

402/1

Commanding Officer and Director

U. S. Navy Electronice Laboratory
ATTN: Library
San Deigo, California 92101

CommandingGeneral
U.S. Army Materiel Command
ATTN: AMCRD.RV-A
Washington, D. C, 20315

405/1

Commanding General
U. S. Army Missile Command
ATTN: AMSML.RRA
Redstone Arsenal, Alabama 350

406/3

Redstone Scientific Information Center

ATTN: Chief, Document Section
U. S. Army Missile Command
Redstone Arsenal, Alabama 35809

410/2

Cemmanding Officer
Aberdeen Proving Ground
ATTN: Technical Library, Bldg. 313
Aberdeen Proving Ground, Maryland
21005

411/2

Headquarters
S. Army Supply and Maintenance.
Haqatr
Command
AToN: AMSSM-MR-M
AMSSM-MM
Washington, D. C. 20315

206/2

207/1

208/1

209/1

Commander
U. S. Naval Ordnance Labora~ory
ATTN: Technical Library
White Oak, Silver Spring,
Maryland 20910
AFSC STLO (RTSND)
Naval Air Development Center
Johnsville. Warminster, PA. 18974

210/1

Officer of Naval Weathei' Service

(Code 80)
Navy Yard (Bldg. 200)
Washington
Washington. D. C. 20390

211/1

Inn ChrgeU.
ChargeOffcer
Officer 211/

DEPARTMENT OF THE ARMY

412

U. S. Navy Weather Research Facility

Building R-48, U. S. Naval Air Station
Norfolk, Virginia 2W,11
DEPARTMENT OF THE AIR FORCE

302/1

303/1

Headvarters
Strr-egic Air Commari
ArrN: DOCE
Offutt Air Force Base, Nebraska 68113
Group (SEPIR)
Systems Engineering
Wright-Patterson Air Force Base
Ohio 45433

Increase to 50 copies if releasable to

CFSTI. See para 5e (4). ECOMR 70-31, for
types of reports not to be sent to DDC.

413/1

Commanding General
U. S. Army Combat Development
Command
ATTN: CDCMR-E
Virginia 22060

Fort Belvoir,

414/1

Chief of Research and Development

ATTN: CRD/M
Department of the Army
Washington. D. C. 20310

ULASS I FIED
socuvity Classification

DOCUMENT CONTROL DATA.- R & D

(S.'..ity .0.60ticati,
.1 rtil, body .1 .6tw'.I and indaming anotaIi- ...

be e.flISd Oeh.'

I. OR IGINA TING A CTOIVIT Y O'Colpolso 04Audo)

CLASSIFICATiON

UNCLASSIFIED______

U. S. Army Electronics Comand

Fort Monmouth, New Jersey

ispmclassiied)

the 0".,011

ast. RIPORIT SECURITY

07703

;b

RU

SREPORT TITLE

A SIX-DEGREE-OF-FREEDOM DIGITAL COMPU)TER PROGRAM FOR

TRAJECTORY SIMULATION
4. O9SC"IPTSVE Nora$ES"f'.;

"U TNORtsI

(Pitel mo..

.1~P,.t
ad I.Cte..

dae.)

WOIdds SIiil lst ne,)

Louis D. Duncan
Bernard F. Engebos
S.

ORPORT

OATE

7m, TOTAL NO. OF PAGSS

October 1967
Se. CONTRACT OR GRANT NO.

6.

No. oil RaeS

go. ORigONATOR'l REPORT NUNIVRRISI

ECOM-5158

PROJECT NO0.

,.

70

TASK NO. 1V014S01B53A-l0

S.OTHER REPORT

NOS) (AW *dI.'.eR"

lkt*ey

6~.0
001001

16. DIST'!*JTION STATEMENT

Distribution of this report is unlimited.

11ISUPPLEMENTARY

NOTES

12. SPONSORIN4 MILITARY

ACTIVITY

Atmospheric Sciences Laboratory

U. S. Army Electronics Coaffand

White Sands Missile Range. New Mexico

________________________

IS. AGATRACT

A documentation of a six-degree-of-freedom model for digital

simulation of the trajectory imf an unguided, fin-stabilized, windsensitive rocket is presented. This model was developed by the Atmospheric Sciences Laboratory, White Sands Missile Range, New Mexico,
to study both theoretical and empirical performance characttzistics oF
unguided rockets. The basic equatiort of motion and their mathematical
formulation for this model are presented without derivation. A genf
flow chart, a listing of the program, a list of the principal flads
used, and a listing of a typical i".Lut data deck are included.

DDr"*V e'4 3

::PLAC*0 ".
0.276

V0,11

FORM IS7,
ARMY

Vag-

1 JAN".

RICSE IS

IJI&LASSI FIED

UN~CLASSIFIED
IsLINK

____

Kgy WOROS-PIOLE

WT

LINK 9
ROLIE

WT

1,, Trajectory Siula'tions

F 2. Computer Programs
S. Unguide-1 Rockets
4. Equations of Motion

UN~CLASSIFIED
Sacwity Clamaifcellso

LINK C
OE