NASA
User's
Report
4705
Manual
for Space N.C.
Contractor
Debris
Surfaces
(SD_SURF)
Elfer
Contract Prepared
for Marshall
Space
NAS8-38856 Flight
February
Center
1996
NASA
User's
Report
4705
Manual
for Space N.C. Elfer Lockheed Martin
Contractor
Debris Marietta
Manned
Surfaces Space
(SD_SURF)
Systems
National Aeronautics and Space Administration Marshall Space Flight Center • MSFC, Alabama 35812
Prepared
for Marshall Space Flight Center under Contract NAS8-38856
February
1996
FOREWORD
The
SD_SURF
computer
contract
NAS8-38856
contract
title
Impacts." Jennifer 1990
was
The
September,
NASA
"Structural
Technical
Robinson. and
from
programs
The
Marshall
Damage Monitors
code
1992.
and
and
Updates
user's Space
Prediction were
Joel
user's are
manual included
iii
guide Flight and
Williamsen, were through
were
prepared
Center.
The
Analysis
for
Greg created
study
Hypervelocity
Olsen, between
October,
under
1995
and October,
ACKNOWLEDGEMENTS
The author wishes to acknowledge
the support and assistance of the
following people:
From
Martin Marietta Manned
Space Systems
Robert Meibaum John
From
Magyari
NASA
Marshall Space Flight Center
Joel WiUiamsen Greg Olsen Jennifer Robinson ScottHill
TABLE
FOREWORD
..................................................................................... iii
ACKNOWLEDGEMENTS TABLE
OF CONTENTS
APPENDICES FIGURES
OF CONTENTS
.................................................................. v
.................................................................... vii
................................................................................... viii ......................................................................................... x
1
Intzx_uction
........................................................................
1
2
BACKGROUND
......................................................................
2
2.1
Environment
.......................................................................
2
2.2
Ballistic
2.3
Probability
3
SD_SURF
Analysis
4
SD_SURF
- FORTRAN
4.1
Platform Selection
5
SD_SURF
6
Installation
.......................................................................
6.1
Installation
- SD_SURF
6.1.1.
Select
One
6.1.2.
Install
the SuperTab
6.2
Installation
- SD_SURF
VAX FORTRAN
6.3
Installation
- SD_SURF
EXCEL
7
Performing
a FORTRAN
7.1
Running
7.1.1.
Finder
7.1.2.
Input
7.2
BUMPERII
GEOMETRY
7.3
BUMPERII
RESPONSE
7.4
A SURF
7.5
P_SURF
7.6
R PLOT5
Limit
Surface Analysis
..........................................................
2
..............................................................
3
Approach
...............................................
5
.............................................
7
Version
.
eo*eoeeieeeeeeeeeeeoee.eeoeeeeeoeoeeJeeeeeeele"
- EXCEL
of the
&0 Version
Open
..............................................
Macintosh
Application
Options
Files
9 11
Applications
Disk
Universal
Applications
7
*Deeleoloeooaeee
...................
12
.............................
12
....................................
12
.............................
12
..........................................
Analysis
12
.......................................
on the Macintosh
13
...............................
13
Box .............................................. ..................
from a BATCH.COM
13
File ...........................................
14
.................................................... .................... •
ooeeeo.oeoeeee,eeoeoe
•
De---e•
eeeee.o•eeeeeeeo•e,.eeeeeee,
•
•
.eeoeeee.ao
14 ee
iooeoe.,
e.o
J.
oo.leoloeeooeoeeoee
ie
............................................................................ • .oeoee
............................................. •
,e
cote.
vii
eeee
eeeel
ooo
14
•e
14 1,5
•
•
15
8
Performing
8.1
SD_SURF
8.2
Ballistic
8.3
AREAMAKER
9
Probabmty
9.1
Effective
9.2
Penetration
10
Problem
10.1
Macintosh
10.2
EXCEL
Macros
and Templates
10.2.1. EXCEL
Macros
- Error
10.2.2. EXCEL 10.2.3. EXCEL
- Errors - Updating
10.2.3. EXCEL
- Hints
11
References
an EXCEL Macro
Analysis
/ PNP_Template
..........................................
16
.......................................
Limit .................................................................... MACRO
St_lies Area
and
Area
Template
19 ........................
..............................................................
21
....................................................................
Anslysis Hesolu_on
21
..........................................................
22
............................................................
Applications
23
......................................................
on Opening
23
............................................
on Open SD_SURF R_PLOT5
23
or AREA_MAKER
or A/P_SURF
Output...
Links ..................................................... for Custom
Analyses
.....................................
.........................................................................
viii
19
23 24 24 25 26
APPENDICES
P_SURF
Source Code
B.
A_SURF
Source Code
C.
R_PLOT5
D.
Limits to Language
E.
SD_SURF
F.
AREA_MAKER
Macro
G.
SD_FUNCTION
Macro
H.
Typical EXCEL
I.
Typical FORTRAN
A.
Source Code Systems FORTRAN
Macro
Output Input and Output.
ix
for the Macintosh
FIGURES
Fig.
1.
Impact
flux versus space debris diameter ....................................... 28
Fig.
2.
Angular
and velocitydistributionof debris flux............................... 29
Fig. 3.
Penetration
mechanisms
....................................................... 30
Fig. 4.
Ballistic Limit
Fig.
5.
BUMPERII
Fig.
6.
SSF Model for BUMPERII-GEOMETRY
analysis ............................ 33
Fig.
7.
SD_SURF
- FORTRAN
Modules ............................ 34
Fig.
8.
SD_SURF
- EXCEL
Fig.
9.
SD-SURF-AREA_MAKER
Surface ......................................................... 31
Modules, Input and Output ......................................... 32
Fig. 10.
SD-SURF-AREA_MAKER
Fig. 11.
AREA_MAKER
and BUMPERII and SD_SURF Macro Macro
- FORTRAN
Modules .................. 35
Available Geometries ................... 36 Dialog Box .................................. 37
analysis of a plate edge on to x ............................... 38
(The surface normal isin the y axis directionon Fig. 9.) a) The projected areas in each threat direction. b) The probabilitydistribution(as in Fig. 2.) c) The effectivearea at each velocityand obliquity. Fig. 12.
A_SURF
analysis of the same plate in Fig. 11 ................................. 39
(45 Threats used in GEOMETRY) Fig.
13.
Area Analysis of a sphere ........................................................ 40
Fig.
14.
P_SURF
analysis ofthe fiatplatein Fig. 12 ..................................... 41
Fig.
15.
P_SURF
analysis of a SSF
module ................................................ 42
(1995 exposure environment).
I INTRODUCTION The SD_SURF
code takes a differentapproach than that used by other codes
such as BUMPER
[1],BUMPERII
[2],or Space Debris Vulnerability (SDV) [3].
All of the codes treata space vehicle as a faceted geometry. The space debris environment
is considered to be a seriesof threats from differentdirections.
Each directionhas a corresponding impact velocity. The other codes calculate the probabilityofno penetration foreach facetbased on the exposure area and the penetration resistance (ballistic limit)to each threat'simpact velocityand obliquity. This output tellsthe designer which areas are most vulnerable. However, the output does not give any information to help selectthe most efficient
shield
designed narrow
design
for a given
to cover
all possibilities,
the number
of variables
area.
While
additional
parametric
studies
information
to be optimized
can be
can help
for a given
a designer
penetration
mechanism. The exposed
new
approach,
area
on the spacecraf_
be generated or the
from
code
geometry small
can
used
read
the
analyst
FORTRAN
most
in Section
A_SURF,
module) analyst
If different
SHIELD,
must
review
repeated
table
or cylinder)
allows and
can
a complex stored
as a
or the VAX.
of this
manual
at home
with
Installation
the
is
is not critical,
so the
here.
multiple
files with
the same
name
on the
i
I
PIDs
element
PIDs
must
be analyzed
select
the model!
ID ranges
in all of the selected
are summed
know
This
will feel right
the Macintosh
range
must
The
sphere,
calculations)
BUMPERII
in each
the analyst
an
file will be deleted!
and multiple
facets
(plane,
BUMPERII.
the VAX, ff you save
the earlier
Unlike,
geometry
from
summarizes
or optimization.
6. A complete are
first
and obliquity.
self shielding
with
on both
Unlike
Macintosh
NOTE:
studies
warnings
WARNING:
PID.
trade
who is familiar
significant
The
(including
applications
described
output
code,
of velocity
of a simple
GEOM
to be processed
SD_SURF
in a table
a description
file for further
The
in the
only those
to only separately, elements
are lumped
ranges
one area
by--]
are ignored. array
regardless
(eg. a window with
together
the desired
along PID.
of a The
2 BACKGROUND Previous approaches willbe explained first,since SD_SURF work.
Space debris codes probabilitycodes BUMPERII
expands on that
(and itspredecessor
BUIVIPER), and Space Debris Vulnerability (SDV) analyze a space vehicle as a faceted geometry. 1-3 These codes calculatethe probabilityof no penetration for each facetbased on the exposure area and the penetration resistance (ballistic limit)to each threat'simpact velocityand obliquity,as described in the following sections. This output tellsthe designer which areas are most vulnerable.
2.1 Envlronment The space debris environment diameter, d, or larger,dependant in the environment
isdefined in terms of a flux of particlesof on the year of interest(due to assumed
growth
as well as solar flux)and the spacecraftaltitude. 4 Figure 1
shows a flux versus d curve fortypicalparameters of interest. The space debris environment
may
be modeled as a seriesof threats from
discretedirections. For low earth orbit(LEO), space debris may
be assumed
to
existin circularorbits. This assumption fixesthe orbitalvelocity.Debris cannot intercept a spacecraf_ from more than approximately 10° above or below a plane tangent to the localEarth normal, otherwise the debris would enter the Earth's atmosphere in LEO
and be removed
as a threat. Therefore, the relativeimpact velocity
is determined by the orbitalvelocity,Vo, and the intersectionangle, g},of
the two orbits.The impact velocity,Vi, is: Vi = 2 Vo. cos (18v_-0)^^ Figure 2 shows the fractionof the totalfluxcoming from angles relativeto the directionof flight.The relativeimpact velocityforthe intersectionof 388 km orbitsis also shown When
on the plot.
the spacecraftattitudeisfixedrelativeto the earth,the orientationof
each faceton the surface willdetermine the most probable impact velocitiesand obliquities.
2.2
Ballistic
Limit
The spectrum
Surface of debris
shield lead to a variety 3. Figure
4 Klustrates
multi-watt
shield.
surface will
sizes, velocities,
of penetration a ballistic
A projectile
penetrate
limit
mechanisms.
These
which may impact are illustrated
surface for hypervelocity
diameter
the shield.
and obliquities
impact
at a velocity and obliquity
A diameter
2
in Fig. on a
above the
below the surface will
not
a
penetrate the shield.Changes mechanism
in shield parameters affecteach penetration
differently. Therefore, itis important for the designer to know
what penetration mechanism
has the greatest effecton the overallprobabilityof
no penetration.
2.3 Probability Analysis The probabilityof no penetration (PNP) from each directionand for each element isbased on the Poisson distributionfor zero events: ( nthreats PNPel-
expI-
_(Ni.Ai).t i---I
where (with consistent units) Ni
= flux which penetrates from each threat direction,i. -- 4-fi.Nr(di)
Nr
= flux on a randomly
tumbling plate of diameter di or larger.(As defined
in the specifications.) di
= diameter to penetrate at the velocityand obliquityof the iththreat. = fractionof flux from threat direction
Ai
=
projected area of the facetin the flux direction.
t
=
exposure time.
The totalPNP
is determined by the product of the PNP
foreach element.
nelements PNPtotal
=
I'_PNPj j=l
Figure 5 shows the BUMPERH
modules and theirinput and output as they
calculate PNP. BUMPERII
startswith a SuperTab
output filefiniteelement model of the
spacecraft. Figure 6 shows a model of Space Station Freedom.6 GEOMETRY
module of BUMPERH
The
calculatesthe projected area of the
elements exposed to each threat directionbased. A significantpart of this calculationis intercomponent
shadowing.
This can be a very time consuming
process for a large model. The RESPONSE
module creates a ballistic limit surface from a menu
of user
selectedpenetration equations. The ballistic limitfor each shieldof interestis
3
stored
in a matrix
binary
form in the computer.
binary
file and
and °60 obliquity The SHIELD requested
for every
0.25
km/s and 5 ° obliquity.
Another
puts out a formatted
file with
the
code, ballistic
RPLOT, limit
stored reads
in the
at 0 °, 15 °, 30 °, 45 °,
for 2D plots. module
calculates
by the analyst.
SHIELD
plot probability
BUMPERII
This is also
contours
the PNP
for any range
also has an option
on the original
4
geometry
of element
to create
model.
numbers
a SuperTab
file to
3 SD
SURF
ANALYSIS
To design penetration Debris
the most
or damage
Surfaces
effective
each
the
analyst
must
is predominant.
computer
with
by the probability
PNP(V,_)
shield,
mechanism
(SD_SURF)
The flux associated weighted
APPROACH
It is the
programs
point
to provide
on the ballistic
of an impact
at that
know
particular
goal this
limit
which of the
information.
surface
velocity
Space
and
can be obliquity.
= exp[-N(d).f(V).A(V,_).t]
where A(V,_) = totalprojected area ofthe spacecraftthat willbe impacted from a debris particleat an obliquity,_,at velocityV. f(V) = the fractionof the totalflux at velocityV. N(d) = the flux associatedwith the diameter d that just penetrates at V and ft.
approach isto store the elements of area in an array in small
The SD_SURF
increments ofvelocityand obliquity.The totalPNP
-tot° = e
_,
_(N(dij)-A(Vi,_j)) i=l j=l
There is a differencein the PNP
isthen given by:
1
calculatedfor a unit area at a singlevelocity
and obliquityversus distributingthe area over two bracketing velocitiesand two bracketing obliquities.This isdue to the non-linear relationshipbetween
flux
and diameter. On the other hand, the analysis of a curved surface in BUMPERII
is more accurate than SD_SURF
only ifthe angle subtended by the
facetsis smaller than the fivedegree increments used on the RESPONSE
and
AREA_SURFACE
be
tables. SD_SIYRF
used as confidentlyas BUMPERII
overallprobabilitycalculationsmay
for models that have coarser increments than
5° facets and 90 threat directions,regardless of the ballisticlimitsurface. Finer models are prohibitivelytime consuming
and will not necessarily produce a
differentresult. Like any computer model that treatsa continuous process as a discreteor finiteelement, there is a chance of introducing errors. Of course, BUMPERII and SD_SURF
do not require the same levelof debugging as a finiteelement or
hydrocode model.
One sources of potentialerror is in the shadowing
and area
calculation,
which
a good quick affect
BUMPERII
way to determine
the effective
ballistic
limit
average
impact
area.
surface
introduced
fiat facet Each
and
velocity
and
curvature
However, limit
surface
is possible
shape
were
information
caught
small limit
into
models
and
the
is potentially
a discreet
surface
number
as if it occurs
the magnitude
a small into one
of threats. at one velocity
of which
depends
surface. increments
or sharp
to the
of using
debris.
of a curved
area
error,
obliquity
surface
to judge
by the model
area
is
is how well the
method
there
option
and geometry
for space
varying
exposed
peaks,
In other
matters
geometry
inappropriate
shadowing
of error
The old meteoroid
distribution
and
troughs
limit
source
all of the exposed
key areas.
of the ballistic
provides
angular
of the ballistic
What
potential
is smoothly
is a relatively
has deep
and
environment
all of the
facet treats
to miss
environment
if the discretized
surface
if the velocity
undersampled.
The partial
is certainly
the debris
This
well.
is interrogated.
limit
and
quite
second
by lumping
obliquity.
on the
The
velocity
If the ballistic error
does
terms,
then
are large, a larger
the ballistic
and the ballistic
error
limit
is possible.
surface
analyst
is whether
it affects
has a direct
impact
on the fidelity
needed
whether
the
whether
they
6
to sample cusps
it.
in the ballistic
will influence
can be
the result. of the
The SD_SURF limit
the PNP.
It
output surface
The
4 SD SURFThe Fig.
interrelationship
7.
SD_SURF
GEOMETRY from
FORTRAN
source
acts
FORTRAN
modules
as a post-processor
output.
P_SURF
codes
The create
of SD_SURF
is shown
of BUMPERII-RESPONSE
It provides
additional
ASURF
are required
and
are in Appendicies
ASURF
module
the exposed
area
than
lump
uses
the lever
over
the four nearest
equal
the area
information
in
and
not readily
obtainable
uses
flux.area.time should
be compatible
final
with
Appendix
C. It is used
formatted
files.
increments
rather
RESPONSE.
Commas
4.1
Platform
Selection
SD_SURF
and
computers. on the
Macintosh
The
to Language Macintosh
meteoroid RESPONSE
analyses module.
area,
for one facet
The sum
the output,
based
graphics
of the
file and module.
The
map
file.
by the EXCEL
as well
of this
The
source
listing
is given
and
increments
and
Systems
FORTRAN
version
3.0 was
used
minimal
changes
FORTRAN cannot
due to memory
are given
handle
limitations. meteoroid
7
finely
analysis
with
Macintosh for compilation
FORTRAN
in Appendix
resolved
However,
by
by EXCEL.
on VAX
the original
files to text
used
to ease import
in
5 degree
compiled
Systems
plots
paper.
file is set up at 0.5 km/s
from
file which
of the contour
output
5 degree
and
which
as a text
BUMPERII-RESPONSE
and
files,
the
is generated
Examples
section
as delimiters
have
and
is
formatted
output
been
Debris
areas
a formatted
to calculate
packages.
0.25 km/s
applications
one threat,
exposed
and RESPONSE
platform,
is R_PLOT5.
formatted the
and
or it can be read
contour
in the next
BUMPERII
with
A_SURF
by the P_SURF
any FORTRAN
are used
Language
Limitations
A text
to translate
than
obliquity,
section.
module
The text
and
in the A_SURF
in the examples FORTRAN
velocity
in BUMPERII-SHIELD
for sophisticated
will be shown The
reads
to
Rather
an unformatted
check
in the next
array.
output
obliquity.
obliquities.
both
to manually
(NAT)
be used
and
file can be read
flux routines
binary and
the projected
creates
same
The
of velocity
at the nearest
module
module
the
as a function
by BUM:PERII.
P_SURF
an analysis.
B respectively.
reported
as described
The
matrix
binary
to perform
the BUMPERII-GEOMETRY
velocities
file can be used
modules
reads
to distribute
A_SURF
The unformatted
A and
of one facet
rule
to the area The
may
of the
BUMPERII. Only
text
VERSION
this
code.6 D.
models
or
does not
145 threats
in
affect
the
BUMPERII
can be used on models with 2100 elements. If partialshading is
used in BUMPERH,
additional elements are created so the totalnumber
of
elements must stay within the allocation. The Macintosh
SD_SURF
to share data with EXCEL.
and BUMPERII
networked may
However,
with Macintosh
and IBM
be run on a mainframe
computer
many
applicationsmake mainframe
compatible PCs.
and the text filesmay
and used in an EXCEL
itconvenient
computers are The FORTRAN
programs
be transferred to a personal
analysis or forimport into any available
charting package. Language
Systems FORTRAN
loops,as allowed by VAX
did not support jumping
FORTRAN.
changes to the originalBUMPERII
However,
into IF-THEN
this requires only 3 minor
code, and ithas been requested that these
changes be incorporated in future releases of BUMPERII. Language
Systems
FORTRAN
and Absoft FORTRAN
variables. A REAL*8
with a REAL*8
dummy
argument
mismatch
must be matched
listof the calledsubroutine,
willbe transferred. The same holds true for integer
variables as well. FORTRAN-LINT Corpation, (Palo Alto,CA
and the subroutine
variable in the main program
variable in the argument
otherwise wrong numbers
Furthermore,
for the Macintosh require
that the variable sizematch between the callingprogram dummy
or DO
by Information Processing Techniques
(415)-494-7500)provides a means
of checking that this
is detected,since itisnot identifiedby either Macintosh
compiler.
8
5 SD
SURF
The final
- EXCEL
the EXCEL
product
3.0 VERSION
version
offers
is not as fast
or as
an alternative =turn
key"
to the FORTRAN
as a FORTRAN
However,
it has
the advantages
of a spreadsheet.
checking
is very
easy
there
is easy
access
The
structure
of the EXCEL
version
is shown
PNP
calculation
and
is the PNP_Template.
version.
The
application.
Customization
and
error
to graphing.
There
in Fig.
8. The backbone
are several
different
of the
areas
on the
worksheet: • Ballistic
Limit
five degrees imported
surface,
of obliquity.
from
and
• Area
for each
macro
that
Surface,
of 0.5 km/s
Limit
Template
and
or
5.)
year,
solar
flux
level
(explicit
or
limit
surface.
(This
altitude.
calculation
function
including
in increments
on a Ballistic
via RPLOT
definition
calculated),
to penetrate
(It is created
RESPONSE
• Environment
• Flux
diameter
diameter
is defined
A(V,_),
in the ballistic
on the function
created
using
macro
Area_Maker
is a
worksheet.)
Macro,
or imported
from
A_SURF. • Flux
• Area
• Time,
is used
to calculate
Function
macros
limits
or flux
files
and
name but
operate
for later
_. (The
input
of named
and
values.
summation
arrays
Hardcoding
from
provided
the
of these
may
limit
cells
and
make
names
ballistic
provide area
be customized
would
general
to calculate
macros
ballistic
names
by using
are used
Command
Any of the templates
use.
the flexibility
V and
as subroutines
for appropiate
Template.
for each
the PNP.)
the pasting
the PNP_
N.A.T,
saved
to
by any
for a new
was deemed
of
templates
and
it easier
control
user,
to be more
important. The
Area
Area_Maker menu.
each
Macro.
The
entered
maybe The
standard
in customized
velocity
velocities), angles
compared
from be used
in SHIELD, However,
the
created
analyst
geometries dialog
increment.
SD_SURF must
Surface
selects are
boxes
This
on the Area_Template
shown shown
for EXCEL
geometry in Fig.
in Fig
is effectively
to the 45 threat
direction
the
default
10. Each
the
from
specific facet
a pull-down
geometry
is analyzed
(at equally
in BUMPERH
is at
spaced
(at equally
spaced
of flight). lacks
some
of the features
for shadowing
analysis
in GEOMETRY,
or the extensive
iterations
required
the GEOMETRY
desired
9. The
64 threats
using
and RESPONSE 9
of BUMPERII. multiyear
to run PEN4
output
BUMPERII flux
averaging
in RESPONSE.
may be imported
via the
FORTRAN
A_SURF
be programmed
into
and R_PLOT5 the EXCEL
programs.
macros
analysis time.
10
with
Multiyear
flux
a corresponding
calculations increase
can in
6 INSTALLATION
6.1 Installation 6.1.1.
- SD_SURF
Select
One
The selection processor,
of which
co-processor,
option,
and
The
should
the output
supplied
handle
options.
are smaller
than
(and
SD
SURF
of elements
ITH
=
145
number
of threats
IPFUNCS
=
31
number
of PIDS
=
700
ITH
=
145
IPFUNCS
=
12
options
do not have
a Language
customized
application.
for array Systems
file processing
spaces
and whether
was
version
CONTOUR
1.3
will not
(good
for meteoroids)
13M_Small
IELM
out excess
memory
v. 1.6)
number
generation)
(Virtual
BUMPERII
apply:
2100
If Batch
2.2
and a math
BUMPERII
option.
=
If different
a either
for the largest
all other
BUMPERII
require of ROM
IELM
MacBumperHv
names
1.3 applications
from any
limits
on the machine
RAM.
programs
MacBumperIIvl3M_Large
surface
to run depends
8 Meg ROM is recommended
The SD_SURF
following
Options
(MacBumperIIv13M_Large)
and 5 Meg ROM
will handle
run.
version
Therefore
not tested.)
Disk
set of applications
or 6 megabytes
coprocessor.
Application
and available
The MacBUMPERII megabytes
Applications
of the
Macintosh
sizes,
FORTRAN
is desired
create
compiler,
file without
a word processor!
English
etc. are required,
on the Macintosh
a batch.corn
" " using
no coprocessor,
or Metric
This
please
request
(useful
and a
for RESPONSE
comment
lines.
Then
will
avoid
some
errors
Then
use
is read properly.
you
option
strip in file 8 to
read in the file. All of the files the floppy
should
be copied
to a new folder
on your
hard
on top of the hard disk will put all of the files in a folder
disk.
11
disk.
Dragging
on the hard
6.1.2.
Install
the
Four SuperTab EXCEL
SuperTab
Universal
Universal
Files
Files were stuffed and placed on the "SD-Surf
I Stuffed UNIs" disk. The fileswere placed in a self-unstuffingarchive
using Stufl_tTM Deluxe by Aladdin Systems, Inc. The filesare: *
MB17-ALLUNI
•
MB17-CR1.UNI
•
MB6-CR1 .UNI
•
PLATE.UNI
The
first
coordinates
three
files were
in the last
distributed
with
file may be edited
BI.TMPERII.
by hand
The plate
for a plate
edge
of any size and
orientation. Open or "Double files
will be unstuffed.
BUMPERH
They
application.
approximately
Installation
(Stuffed
SuperTab.Uni
Files.sit) and the
should be placed in the same folder as your
When the files
1,800 K on your hard
OK to overwrite
6.2
click" the application
are unstuffed
disk.
they
If PLATE.UNI
require is already
there,
it is
it.
- SD_Surf
VAX
FORTRAN
The source code isprovided on both Macintosh and IBM
compatible PC
formatted disks. The source code and applicationsshould be transferred to the VAX
and compiled. The SD_SURF
BUMPERII
programs
version 1.5 work with any
version 1.2a output. The SD_SITRF
any BUMPER]]
programs
version 1.3 output. The BUMPER]]
version 1.6 work with
versions 1.2aM and 1.3M
(Martin Marietta Modified) contains the features necessary to compile on the Macintosh. memory PRV
Version 1.2aM also fixesone error in BITMPERII
allocationfor the variable IDG.
Version 1.3M fixesan error in function
(with negligibleimpact on overallPNP). Note that on the IBM
compatible PC formatted disk,the filenamed
Solar_Flux.Dat was renamed The SD_SURF BUMPER]]
due to PC naming
output files(interms of number
- SD_SURF
The EXCEL
files are provided
disks.
of elements, threats,and PID
be adjusted in the COMMON*.BLK
Installation
formatted
restrictions.
filesshould be compiled with large enough variables to open
cases). These may
6.3
regarding the
The files
files.
EXCEL
should
on both be copied
disk.
12
Macintosh to a folder
and IBM
compatible
or directory
on your
PC hard
7 PERFORMING
A FORTRAN
To perform BUMPERII should
an analysis,
should
be run.
ANALYSIS
the GEOMETRY
be run
as described
R_PLOT5
should
and
below.
be run
RESPONSE
Then
modules,
A_SURF
to plot RESPONSE
and
data
of
P_SURF
or prepare
data
for EXCEL.
7.1 Running
Applications
The applications or in System work
I
7. This
using
Macintosh
are compiled
to run in the background
slows down
to be performed
30 PIDs
on the
while
the calculations
the analysis
RESPONSE
is very
is being
time
the
hold down the
standard
command
When
Macintosh the code
the application, results
stops,
or they
should
may
be viewed
using
but it allows
other
performed.
The calculation
of
with
the same
button
and type
name
on the
a period.
This
window
to review
a uniformly
may be printed
or print
spaced
with
font
directly
any text
such
from
editor.
as Monaco
The
or
Courier. Double
clicking
application
that
application.
on a text
will open
If using
file may not open
text
EXCEL,
files,
and
open
using
it directly.
then
open
the
the
comma
Instead
start
files from delimited
any
that option
for TEXT
files. There The
are two
FINDER
possible
additional
open
features
on the Macintosh
box may be accessed
for "batch"
processing.
These
with
a "?" and
features
are
version input
of BUMPERII.
from
described
a text
in the
file is
following
section.
7.1.1.
Finder
Open
On the VAX, when may
respond
retained
with
in SD
a "?" for a print
SURF
On the Macintosh, The
analyst
to jump
may
down
Box BUMPERII
then
modules
scroll
requests
a filename
from
the
out of files in the directory.
analyst,
This
feature
FINDER
open
he was
for the VAX.
a "?" response
brings
up the normal
to the file of choice
the list. 13
is
a process.
in the
be saved
somewhat,
files
"Apple"
to stop
the results
Multi finder
consuming.
WARNING: Unlike the if you save Macintosh, the earlier file will be deleted_
To stop an analysis,
under
or type the first
letters
box.
of the file
7.1.2.
Input
from
The analyst BUMPERII,
a BATCH.COM
may create a file using the existing
and
the BUMPERII
then
use
initial
in BATCH.COM
continues
from there.
instead 7.2
output
The
7.3
files.
feature
is particularly
The BATCH.COM
subroutine
to support a SHIELD
in
file. The initial
text
and the analysis
useful
for generating
file may be edited
is run as normal
subroutine
analysis
to create a certain
using
a text editor
in BUMPERH.
number
for an analysis
wall
thickness
is unaffected.
is perfectly of PIDS
acceptable.
as required
A RESPONSE However,
by SHIELD.
of a specific range of elements.
is to create a series of shield for parametric rear
option
process.
of the RESPONSE
is required
and/or
eighth
RESPONSE
The operation
shield
the new
option in
GEOMETRY
BUMPERII
necessary
by using
by BUMPERII
the BATCH.COM
GEOMETRY
analysis
file for input
are ignored
This
of repeating
BUMPERII
that
BATCH.COM
options list: Read from a BATCH.COM
responses
RESPONSE
File
and/or
analyses
it is not Only one
Another
(eg. step through
option
bumper
spacing).
7.4 A_SURF Run
A_SURF.
A typicalinput session is shown are identicalto BUMPERII
NOTE: A_SURF, facets
know
PIDs must
dement
ID ranges
in all of the selected
are
summed
must
be analyzed only those
to only one area separately,
elements
with
ranges array
(eg. a window the desired
options
user's manual.
are lumped are ignored. regardless along PIT).
together
charting
file may be reviewed
with
a text
package.
14
editor,
of PID. a module)
The analyst
or used by EXCEL
by
The
the model!
The text other
PIDs
multiple
range
select
I. The environment
and described in the BUMPERII
SHIELD,
and multiple
in each
different analyst
Unlike,
in Appendix
or any
If the must
7.5
P_SURF Run
P_SURF.
RESPONSE
output
A typical
input
are identical
based
use the A SURF The
The
output
charting
The
maximum
and plot
file may
any
It contains
single
is used
for input,
to the PID
screen
as is a
environment user's
as shown.
or all of the PIDs
a text
the text
options manual.
The
analyst
may
in the RESPONSE
studies. editor,
based
is included
from
The
in the BUMPERII
with
both
used
I.
for parametric
of all PIDs
influence
in Appendix
is output
be reviewed
be remembered
a greater
file (.ASB)
described
is useful
subroutine
It should can have
is shown
with
case
description
of RESREAD
7.6
output
package.
values.
session
carpet
latter
binary
file (.RSP).
to BUMPERII
The text
output.
The A_SURF
or used
carpet
in the
plot
output
by EXCEL and the
or other
calculated
due to the structure
BUMPERII-SHIELD.
that
a group
on the
of cells
overall
PNP
with than
moderate a single
values
of NAT
cell with
NAT.
R_PLOT5 Any or all of the PIDs in a RESPONSE
format
by R_PLOT5.
one or all of the PIDs output
Run R_PLOT5. for output.
due to the structure
output
Select
the Response
The description
of RESREAD
SHIELD.
15
file may be converted
subroutine
output
of all PIDs used
to text
file and
is included
select in the
from BUMPERII-
8 PERFORMING
EXCEL
AN
The EXCEL
analysis
If A_SURF
and
may be performed
R_PLOT5
a VAX or other
computer
Macintosh
import
may
There
are minor
ANALYSIS
text files
on which
differences
The
graphics
in the AREA_Template
PC.
A text
macros.
down
The
the
shorter
PC names
are to open the Macros
and
let them
so that
the names
commands
to change
and
were
open
first
adjusted
the templates
(SD
can be recorded
to the macro
sheets
so that
or save,
EXCEL
in the
macros.
not convert
SURF
(or guide
files.
Windows
would
from
The
application
on the Macintosh
guidelines
impo_
can be compiled.
Macintosh
is included.
or a PC.
the PC must
the Macintosh
description
AREA_MAKER) process)
the the codes
between
and
a Macintosh
are to be used,
from a VAX or use
spreadsheets
General
on either
to the
and
you through
sheets.
the
Use the
the current
names
pull
can be
recorded. Watch
the message
The Macros given
are documented
in Appendix
The analyst reviewed •
box at the lower
in Appendices
be familiar
EXCEL
manuals:
Command
and
Function
Command
Macros
would
E, F, and
G.
Typical
output
is
H. should
in the
left for instructions.
by "Run"
under
EXCEL.
These
features
should
be
Macros.
can
do manually
with
automatically
in EXCEL.
"Macro."
perform
The macro
Most
almost
any
commands
are available
function
may
installed
you
be accessed
in pull
down
menus. Function
Macros
Variables
may
be passed
to the function
variables
may
be returned.
Function
templates
to perform
of the list
generated
like
arguments
for the function
determined
by looking
"Window"
files
any worksheet
- The
on Excel
get to the inner
when
macros
one or more
macros
sheet.
16
and
at the end are available
is open.
in the name,
They
The or may
are identified
be
by the
are received.
and" UnHide..."
respectively.
appear
sheet
commands
can be used to keep your work area
workings
by the worksheet
command,
the macro
they
and
or FORTRAN.
are used
may be abbreviated
in the order
"Hide..."
macro
Function
at the macro
function
in BASIC
Function
by the Paste
into
Hiding
a sub-routine
calculations.
for pasting
ARGUMENTS •
perform
under neat
If after you make
or to let you
a change
you
want the sheet to open up hidden the next time you run it,just change a cell(e.g.add and deletea space),hide the sheet (without saving it),quit, and say yes when
itasks ifyou want to save changes. Itwillbe hidden
the next time you open it. •
Changing
filelinks (especiallyuseful for charts or function macros) is
under the FILE *
menu.
Automatic/Manual
recalculation- how
to change (Options - Calculation)
and avoid on saving or printing (Apple-period stops recalculationon saving). •
Auto_open/close macros.
These run automatically at the open or close
of a macro sheet. Pull down
menus
and opening dialog boxes are added
(or deleted) with these macros. •
Excel Startup Folder - Filesplace in thisfolder(inthe System folder)will automatically be opened at the startof an EXCEL SD_Function_Macros done the SD_SURF
may
be moved
Auto_Open
Macro
session. The
here for general usage. If this is should be changed
so itwillnot
try to look forthe file. •
Open/Save
File Options -Review Open comma
delimited textfiles(CSV =
comma separated variables)and Save.As Options. See Problems Section.
17
8.1
SD_SURF
Macro
/ PNP_Template
Open the SD_SU'P_
Macro
SD_FUNCTION_MACRO sheet.
sheet.
The auto_open
prompts
are in the
PNP_Template
macro
Message
has been
included
analyst
the next
time.
down
menu
SD command Open
An auto_open
If it fails
asks window
is provided.
surface
A pull
first.
the analyst
the analyst at lower
If a particular
with
Opens
R_PLOT5
left of the screen)
these
to open the
a PNP_Template.
limit
by a new name
Function (This is displayed
R_PIot5 Output
is requested
to open
ballistic
it may be saved
is installed
macro will try to open the
(NB:
A
or effective
area
and
by the
opened
commands:
in the Message
output to Paste to PNP/Ftux
Bar.) Template
Open BL Template
Keeps track of which file to use as Ballistic Limit template
using macros.
Ballistic
Copy Ballistic Limit from
Template
Limit to PNP
Active BL Template
to PNP/Flux
.
Open Area Maker Macro
Use Area Maker
Macro to open A_Surf
output or create
new geometries.
.
Open PNP Template
Keeps track of which file to use as PNP/FLUX
template.
Save PNP Template
Keeps track of which file to use as PNP/FLUX
template.
Set PNP/Flux
Keeps track of which file to use as PNP/FLUX
template.
Template
Close SD Surf Macro
These
Closes Macro and deletes SD menu.
are the instructions
included
on the
PNP
Template:
QUICK INSTRUCTIONS: BALLISTIC LIMIT Calculations Use Ballistic Limit Worksheet & BALLISTIC LIMIT TO PNP (in Pull Down Menu) or... Run RESPONSE and R_PLOT5 and Usa OPEN R_PLO'r5 OUTPUT (In Pull Down Menu) GEOMETRY Analysis Use Area Maker Macro to make a new table or open GEOMETRY/A_SURF output ENVIRONMENT Enter in Q15-Q21 this worksheet (or Velocity distribution on Area_Template) or... Change SD Rux Function MACRO (and Quick Flux macro if necessary) or... Change flux formulas D79-V110 (Fill right and down) Use the PULL Down Menu on Right end of Menu Bar and watch Message Bar at bottom.
Printing of the resultsis set up fora Macintosh LaserWriter. Other printers must be formatted by the analyst. The print area contains some pages which are only filledunder certain conditions (eg.pasting from the Area Template). These may
be deleted from the Print Area ifdesired.
18
8.2
Balllstic
Limit
The Ballistic Whipple
Lizmt Template
Bumper
and Multi-shock
SD_Fu_ctdon_Macro for an interpreted "Dreadbom'd"
may be used to create equations
because the PEN4 spreadsheet.
new equations
Note that BL_Paste
a surface.
are include
routine
The Ballistic
m the
requires
Limlt
too much computation
Template
or custom modify a RESPONSE
copies specific cell ranges rather
wish to modify the template
OPJy the JSC
can be used to / R_PLOT5
text file.
than named areas if you
or create a new template.
The SD pull down menu can be used to open the template
or the R_PLOT5
text file and paste the results to the PN'P_'I_M_PI._TE.
8.3
AREA_MAKER
MACRO
Open the Area the analyst A pull
Maker
and
Macro
Area first.
Template An auto_ open macro will
then
request
to open the Area Template. down
menu
is installed
with
Areas Command
Status Bar Text
Clear
Clears
Arrays
these commands:
Area_Array
& Description_Array
Rectangle
Adds a Rectangle
Disk
Adds a Disk to Area..Array
Cone
Adds a Cone to Area_Array
Cylinder
Adds a Cylinder
Sphere
Adds a Sphere
Whole Sphere
Adds a complete sphere to Area_Array.
on Area
Template
to Area_Array
to Area_Array to Area_Array Faster than Sphere!
.
Open Template
Opens a file to be used as the Area Template
Save Template
Saves Template.
Set Template
Identifies
active document
AreaS to PNP
Transfers
Effective Area
Open A_SURF file
Opens
A_Surf/Output
Close AreaS
Closes
Area Maker Macro.
The pull down
menus
Identifies new name as the Area Template as the Area Template to PNP Template.
and puts in 0.5 km/s increments.
are used to create an area array or open an A_SURF
textfile.Figures 9 and 10 shows the geometries and the Dialog Boxes used to describe each geometry. The analyst is advised to take advantage of symmetry
19
and the Area facets.
Multiplier
(Eg. a cylinder
one fourth
the
distribution
need
Macro.
next to the velocity and set the width
need
The
array
The cells functions,
to stop
then
With
If a different
which
still
need
zero cell width.
when
the macro.)
20
can be modelled one axis
the
macro
is desired,
on the
to be normalized, Select
adjacent
by the last warning
are cells
cells.
to the PNP_Template.
"Revert"
as
of rotation,
inclination
to see the function
be copied
symmetric
and the velocity
to a function
open it will be replaced
saved,
calculating
direction
of 4.
refer
to standard
may
time
on the template
but have
is already
(If it has not been
"Cancel"
the flight
only put the value
of all of them area
extra
be modelled.)
distribution
PNP_Template
template. select
effective
spending
its axis along
will re-calculate.
AREA_MAKER
desired
than
and an area_multiplier
cylinder
the analyst
The
with
of a cylinder
only one half then
rather
If the saved
box comes
up
9 PROBABILITY 9.1
Effective
STUDIES
Area
The A_SURF
program
and the Area_Template
calculate the effective
exposed area, f(V).A(V,_),at each velocityand obliquity. Figure 11 illustratesthe analysis ofa fiatplate that is oriented edge on to the directionofflight.The firstpart ofthe analysis isthe calculationof the projected area, A(V,_), relativeto each impact velocitydirection.Figure 11 (b)shows the probability,f(V),associatedwith each impact velocity.Figure 13(c)shows the finalresult,f(V).A(V,_),afternlultiplyingthe projected areas by the relative probability. A_SURF
reveals the coarseness, or granularity,in the spacecraft model and
debris threat in the GEOMETRY
analysis. Solving the firstproblem (a plate
edge on in Fig. 11) using BUMPERII directionsin BUMPERII
produces Fig. 12. The default of45 threat
gives only 22 velocitiesdue to symmetry.
There are
now gaps along the velocityaxis. The _waves" on the surface are an artifactof the coarseness of the modelling. This does not imply the overallmodel isin error,but rather itshows how BUMPERII
and SD_SURF
are querining the
ballisticlimit surface. Ifthe "wave" spacing is small compared the ballisticlimit surface,then the overallPNP
to changes in
calculationis correct.
Since the distributionsare not smooth, the analyst must recognize that adjacent cellswith moderately high impact rates can be more significantthan a single cellwith the maximum The A_SURF model.
impact rate.
output can be used to double check the originalSuperTab
If some elements were entered with normals in the wrong
there may
direction
be unexplained gaps in the model.
The sphere is an easy shape to analyze sinceitlooks the same from any direction. (That is why itisa separate option in the AREA_Maker projected area from any directionis shown
in Fig. 13. Also shown
macro.) The is what it
would look likeifmodelled using facetsthat cover 15 degrees of curvature. The granularity, or waviness is obvious. The sphere is also a good representation of the surface area of any spacecraft which is not Earth oriented. It willappear to be randomly tumbling to the debris flux and average out to the obliqueimpacts on a sphere with the same surface area.
21
9.2 Penetration
Analysis
Figure 14 shows the P_SURF
analysis of the effectivearea in Fig. 13. This is
an example of the textbased contour plot. The ballistic limitwas the RESPONSE output for a 0.050 inch bumper, 4 inch standofF,MLI, and a 0.125 inch 2219 aluminum
rear wall, using the regression equation and default analysis of
Wilkinson
momentum
failure.
Figure 15 is an illustrationof the velocitiesand obliquitiesforwhich most penetrating impacts could occur on one early concept for a space stationmodule. (The same
RESPONSE
ballisticlimitsurface is used as in the previous
example.) It can be noted that BUMPERII 99.88305%, while P_SURF
analyzed the PNP
calculateditas 99.88475%.
identical,but as mentioned
forone year as
The effectivearea was
previously,partitioningthe area to discretevelocities
and obliquitieswill affectthe result,just as assuming
a curved surface is
represented by a flatfacet.The probabilityofpenetration (POP = 1 - PNP) was 0.11695%
forBLTMPERTI
to 0.11525% forP_SLrRF.
the two is 1.5% of the POP.
The percent change between
This differenceisnegligible.
22
10
PROBLEM
RESOLUTION
This
is meant
section
to help trouble
shoot
any errors
encountered
in an
analysis.
10.1 Macintosh
Applications
There are only limited options ifa FORTRAN (Language
Systems FORTRAN
compiler isnot available.
version 3.0 was used to compile the
applications.) RAM
requirements were suggested by the Language
linker. Ifany unusual errors are encountered, the RAM increased using the "Get Info" command
Systems FORTRAN allocationmay
under "File"in the "FINDER."
applicationmust be closedfor thisto work.) Ifavailablememory these techniques may •
be (The
is a problem
be used to increase availablememory:
Run using the FINDER
only. Use System
software 6.X, and turn oR"
Multifinder. •
Remove
non-essential software from the System folder.
•
Use Virtual memory.
The software has worked
Connectix to set memory
to 8 Megabytes
well using VIRTUAL
on a machine
by
with 5 Megabytes
installed(running System Software 6.07). (However, at the relativelylow price of RAM,
a hardware
upgrade should be considered.) It has not
been tested using System 7 virtualmemory. All data filesshould be kept in the same folderas the applications.Use of the Finder open box by responding with a "?"may be able to use filesoutside the applicationfolder,but this has not been tested.
10.2 EXCEL
Macros
If a Macro command,
and
Templates
halts you may
unhide the macro, using the '_vVindow"
and see what operation itwas attempting to perform.
Display values
to see cellresults. On the Macintosh the "Apple -'" willtogglebetween
formulas
and values. (NB: "" is a singlebackquote at the upper leftof the keyboard under the "-" mark.
10.2.1.
EXCEL
If an error command saved
It isnot an apostrophe.)
may
when
occurs
Macros
on opening
be at fault.
quitting
- Error
SD_SURF
If either
from EXCEL,
on Open
or AREA_MAKER
or AREA_MAKER,
document the document
23
SD...SURF
is hidden
and
the any
HIDE()
changes
will open as hidden.
are
Then
when running the HIDE() command, it will either hide an open sheet, or an error will occur. To avoid this problem: •
Save
changes
normally unless
to the macro
and
hide
itself.
it is unbidden.
(To avoid
being
with
it unbidden.
Do not save changes You should
asked
The macro
if you want
not have
will then
open
when
quitting
EXCEL
to save
changes
normally.
to save changes
go to the Auto_close
macro and replace the "x" in the "x SAVE.AS(,0)"
line with an equal
sign. This willstop you from being asked ifyou want to save changes when •
quittingEXCEL.)
Disable, clear or delete the HIDE() document
command
as hidden. (Save changes when
macro is hidden.
Make
in Auto_Open. quittingEXCEL
sure the "SAVE.AS0"
Save the when
command
the
in Auto_Close
is disabled.) •
Work
around the error on open by continuing rather than halting.
•
Work
around by halting,unhiding the macro, and the run Auto_open
using the MACRO
10.2.2.
EXCEL-
When
opening
Errors text files,
tab is encountered. in the arrays. and use the
spaced the
command.
that
of FORTRAN
font such
entire
sheet
Data
should
and
command
or the
Smart
or A/P_SURF
everything
opening
Output
in the first col-m-
the Desktop
delimited.
text files
commas
or Finder.
until
aRer every Start
commands
button
in SD_SURF
will be enhanced Use Styles
with
the use of a uniformly
to redefine
Normal,
or select
the font.
may
be broken
into
many
Parse
in the Flat File Macro.
using
EXCEL's
Consult
Parse
the EXCEL
Manuals.
10.2.3.
EXCEL-
If files
are moved
function
macros
function macros.
are.
Updsting from
Links
different
field
EXCEL
the file use the Text options The macro
a
automatically.
or Courier.
change
in one column
from
do this
as Monaco
R_PLOT5
files are set to write
Before
it is comma
AREA_MAKER Display
puts
EXCEL
The FORTRAN
Open
menu.
on Opening
Do not open directly
to tell EXCEL and
pull down
locations,
It may be necessary
EXCEL
to change
This is explained in the EXCEL
24
can lose track links
of where
to the current
Manuals.
10.2.4.
EXCEL
The EXCEL Paste_Special monitor
Custom
to speed ballistic
author depends
limit
commands.
Table
range
command or year)
if input
should
be consulted
commands
if the screen
in a selected
(eg. altitude
manual
Analyses for customization
are particularly
useful.
and white (in the control panel of a Macintosh)
time
Paste_Special
The
for Custom
and the Table
ECHO(FALSE)
values
user's
to black
up response
- Hints
and
and output
is updated
or redrawn
often.
hints.
Setting
a color
can also speed
(Several
macros
set
up recalculation.) surfaces
may
This
can be used
of cells.
generated
to add,
See the EXCEL
can be used record
be easily
to parametrically
output
(eg. PNP).
are on different
on a table.
25
sheets,
multiply Manual vary A macro
using
the
Copy
or replace
and
the
for details. input
parameters
is available
or if the output
already
from
the
11
REFERENCES
(1)
Coronado,
A. et al.: "Space
Damage
Control,"
Center, (2)
Graves,
Station
Contract
Integrated
NAS 8-36426,
Wall
Design
NASA-Marshall
and Penetration Space
Flight
1987. R.; and Smiley,
XD683-99402-1
J.: User's Guide for Design Analysis Code BUMPER/I,
on Contract NASS-50000
(3)
Elfer,
N.; et al. Martin
(4)
Space
Station
Marietta
Program
Natural
Report
(1991).
IR&D
M-01S,
Environment
unpublished Definition
research,
1987.
for Design,
NASA
SSP 30425. (5)
Elfer,
N.; and Rajendran,
Wierzbicki,
N. Jones
York, p. 41-78, (6)
"Language Corporation,
A. M.: "Space
Eds. Structural
Debris
Failure,
Protection," John
Wiley
in T. & Sons,
New
1989.
Systems
FORTRAN
441 Carlisle
Drive,
version Herndon,
26
3.0," Language VA, 22070-4802,
Systems (703)-478-0181.
FIGURES
Fig.
1.
Impact
flux
Fig.
2.
Angular
and velocity
Fig.
3.
Penetration
Fig.
4.
Ballistic
Fig.
5.
BUMPERII
Fig.
6.
SSF
Fig.
7.
SD_SURF
- FORTRAN
Fig.
8.
SD
- EXCEL
Fig.
9.
SD-SURF-AREA_MAKER
Fig.
10.
SD-SURF-AREA_MAKER
Fig.
11.
AREA_MAKER
a)
Limit
Surface
12.
surface
of debris
normal
effective analysis
(45 Threats
and Output
BUMPERII
SD
SURF
.........................................
32
analysis
............................
33
Modules
............................
34
Modules
..................
35
...................
36
- FORTRAN Available
Geometries
Macro
Dialog
Box ..................................
37
of a plate
edge
on to x ...............................
38
in each
direction
threat
distribution
on Fig.
9.)
direction.
(as in Fig. 2.)
at each velocity
of the same
plate
and
obliquity.
in Fig.
11 .................................
39
used in GEOMETRY)
Fig.
13.
Area
Fig.
14.
P_SURF
analysis
of the fiat plate in Fig.
Fig. 15.
P_SURF
analysis
of a SSF module
(1995
29
31
is in the y axis
area
flux ...............................
..................................................................
Macro
areas
28
30
and and
.......................................
..............................................................
Input
analysis
The projected
A_SURF
diameter
for BUMPERII-GEOMETRY
b) The probability
Fig.
debris
distribution
Modules,
SURF
c) The
space
mechanisms
Model
(The
versus
Analysis
exposure
of a sphere
.............................................................
environment).
27
12 .....................................
.................................................
40 41 42
Orbital Debris Flux
lO0
Inc_inllJon. 28.5" I
0.01 Flux
Nom_ Oe_Qkm (W_mout
[ImpacLs °'°°0' /y/m^2] "__--'_......,._
Avoidance)
1E-o6 % I I
1E-Oa
I"_-'[ I I
1E-lO 0.001
0.01
0.1 Minimum
Fig.
I.
Impact
flux
versus
space
28
I Debris
Diameter
debris
10 [cm]
diameter
Wire lOO*/, Collision Avoidance over 10 cm
lO0
0 120
608 11
150
30 1 4
Velocity
[kin/s]
16
180
Angle from Direction of Flight
210
330
Fraction
240
of Total
300 270
Fig. 2.
Angular and velocitydistributionof debris flux
29
Flux
Impact ORIGINAL PARTICLE
O
Process CRITICALITY
EJ_C'rA
•
BUMPERFRAGMENT PARTICLE
I VAPORIZE
INTERMEDIATE SHIELD DEFEAT FRAGMENTS AND REDUCE VELOCITY • .)
.
.. /
VAPOR
CLOUD
• . . • • " '
_...,,.. _JT_
B.LaE
Fig.3.
.-
7" T'_
# ...PALL
AND
-_/PARTICZ'il 5
REAR WALL . ABSORB MOMENTUM RESIST FRACTURE
_T.Ro.o.
_
CRATER
,ND,VIDUAL CRATERS
Penetration mechanisms 30
AND
Penelradon Mechanism
S;ngle Projectile Penevates Intact.
Q
1.8 Q
Proiect_le Shaners on Bumper and F_gmenLS Penevme Rear Wail
Q
Rear Wall Fails from M,el;s, Momentum,Proiect_te Va_onzes. or Fragments am C_ughl by Thermal Blanket.
Diameter to Penatra|e (cm)
Q
o-- l 0 Velocity (_/sec)
I 16
Fig.
4.
Ballistic
Limit
Surf'ace 31
Rear Wall PeneCated F_ments from Bum.,.,.,.,.,.,.,.,.,._er Alter Projecdle Riccc_els from Bumper.
BUMPERH • Supertab • NASTRAN translator
•
• Space Debris • Meteoroids
Space Debris
J
RPLOT
BUMPERH
Modules,
SHIELD • Space Debris • Meteoroids
• Meteoroids
Fig.
5.
Input 32
and Output
Fig. 6.
SSF Model for BUMPERII-GEOMETRY 33
analysis
Space
Debris
Surfaces-
FORTRAN
Version
Ask;,. FUmF_onTablo
• Super'cab • NASTRAN Model Generation translator
• Space Debris • ._eceoroids RESPONSE
• Space Debris
,li Lookup Tables [binary] r
"Tab Universal File • Nodes Elements
vel. and obl. (FORTRAN vers.) • Space debrisonly • BUMPERII flux subroutines
GEOM * Space Debris , Meteoroids
||
Key:
] BUMPZml
Output
sv SURF
SD_SURF
1
- FORTRAN
and BUMPERII 34
Modules
R_PLOT5 (5dc_ccs)
Space
Debris
Surfaces
- EXCEL
Version
SD_Surf EXCEL vers. (0.5 km/s by 5 degree increments) Command Macros Function Macros • Create Pull Down Menus • Ballistic Limit calculations • Dialog boxes and messages * Flux calculations • Open/Save files and templates Templates • Control order of calculation * Predefined areas, calculations • Cut and paste from templates and formate Optional Input from FORTRAN programs:
• FunctionMacro JBL-Macro • Paste usingSD_SLrRF Template _PNP/FIux
R_PLOT5
Obliquity Templat_
Output BL-RPLOT Template • Open and Paste using SD_SURF Macro
,
V[ Ballistic
Limit
° lFlux (d,environment1 i J " Functionmacro
A_SURF
AREA Template • Aroa_Maker Macro fills in template • Vel.Distrib. function
Output
AREA-Maker Macro • Directfrom textfile usingArea_Maker
Fig. 8.
SD_SUP_
-EXCEL
and SD_SURF
35
[
t
(
• Flux.Area.time .....
- FORTRAN
[
'
Modules
,....' .-......:.,::::::_"
AXES
ROTATION
CYL I NDER
RECTANG_
01SK
SPHERE
CONE
'
I / .o=-_s,art-
.,oo.,..e
.no...,s;=,.. s,.,:_ ,._ :.. cr._ ,,..)[
I Length=L° I I,_(al0ng x) l '_
Nega¢lve numoers may 0e useo. {Eg.-90 to 90 rcr -y sic- of _'yllr_Ger.) Cone z_o _lin(:ar are not symmetric. 0" to 5" ooesnot ._isogestate ;/-'5" r.o 180" The Olrterenc.= _nthe Start aria stoo angles must :e evenly Givlsl01eoy the :no'emetiC
Fig. 9.
SD-SURF-AREA_MAKER
Macro
36
Available Geometries
|[ ]_
RECTANGLE _' Pitch
[+_90"]
Yaw [-+180"] Area
=El
Multiplier
C cancel
CYLINDER
Radius
[m]
Length
[m]
Area
start [I
Multiplier
-L,_
(no closure)
Pitch
[_*90 ° ] [+_180"]
finish
a ngle _-'_
Yaw
facet
angle
[Cancel)
(no closure)
Radius
aft
Radius
fore
11
finish
angle _'7"_
Length
[m]
I1
facet
angle
start
Multiplier
----_I
angle
Cone
Area
I
angle
Units:
[m]
_ Pitch
[_+90°I
Yaw
[_+180 ° ]
_-_ [deg]
[ Cancel
]
|F'1_SPHERE Radius Area
[m]
[rill
Multiplier
start
finlshLat. start
facet
Fig.
10.
angle
Lat.
_'_
Pitch
I-.90 ° ]
_
Yaw
[+180 ° ]
Long.
finish
Long.
SD-SURF-AREA_MAKER
[ Cancel
Macro
37
Dialog
)
Box
0.45. 0.4 -_ 0.5 0.35
At.[..,, o._z___L]p_W.L_V.:._. _0.1 !
_
g.O
" _.".,'.... e:
'0
_
so" '_25' a
oo,5i 0.04
Velocity Distribution
I
Obliquity [degreesl
='_.
[kmle]
_.
-_" _n td-)
j Ilelin
J=" O=5o= i -
•
0.02 0.015
"
" i
eiaglDmalal •
;_
Velocity
Z
.
.
t
0.01
0.005
j"
i
0 ....
almm
0
2
4
6
8
10
12
14
2.5E-2
V_oc_y [k_sl
1.5E-2 .! Aree (sq ml
i t .E-2 t 5.E-
.
'
75"
Obliquity m Velocity
Fig.
11.
,,_J_A__R ('['he a)
surface
analysis normal
The projected
of a plate
in each
-
[degrees} _i
o -
edge on to z
is ]n the y a_s
areas
[km/a]
o
thzeat
direction
on Fig.
direction.
b) The probability distribution (asin Fig.2.) c) The effective area at each velocity and obliquity. 38
9.)
2.50E-02-
2.00E-02
1.50E-02Area
[sq m] 1.00E-025.00E-03-
75 50 25 0 Velocity
Fig. 12.
A_SURF
[kin/s]
_
analysis of the same
(45 Threats
used in GEOM:ETRY)
89
--
plate in Fig. 11
Obliquity [degrees]
12
i0 m diameter sphere I threat direction
I0
j\
Area Ira'2]
78 m^2
-
_
total
projected
: 5 ° facets
5 4
0 0
I0
20
JO
40
Obliquity
Fig.
13.
Area
Analysis
50
Angle [degree._
of a sphere
4O
60
70
80
90
area
RESPONSE PID : 1 A SURF FILE: PLATE PNP(%)= 99.99709 CONTOURS .12345 at
RESPONSE FILE: ONE EDGE.ASB Total Flux x Area equal increments from
RESPONSE.
RSP
ON
IMPACT
VELOCITY
x Time (NAT) 0 to max NAT
- 0.29084E-04 - 0.20709E-05
km/s
Obl
1
2
3
4
5
6
7
8
9
10
1!
12
13
14
15
16
Deg
I
I
I
I
I
I
I
I
I
I
I
I
I
I
i
I
0
,..,.,,e,ee,e,ee,eeo.,,e,,oe,ee.,e,o,,oe,,oe,oe,e,oo,e,,,oo.,.,.
5
,,.,.eeo.ee.eo,.,..,,e.oe.eeeee,,ee,ee,.e..ee.ee.,e,e.,.,,,....°
0
,..,....,,e°o,,e,e..ee,ee,oeee,e,oe,ee.eeo.,,,.e...e.e....o..,,.
5
.,.o..,..,..oe,eoe.e,.eee,o.,,e,,oeo.e,,eo,..e.e....o....,,...,o
70 65 60 55 50 45 40 35 30 25 20 15 10 5
....................................................... 111 ...... ..................................................... 11111 ...... ................................................... 2111 ......... ................................................ 14.11 ........... ........................................... 11.5111 .............. ........................................ 21.31.11 ................ .................................... 14..31 ...................... ................................. 32.12 .......................... .......................... 11.15..11 ............................. ...................... 21..51.11 ................................. .................. 21..31 ........................................ ............. 12...11 ............................................ ..... 11..11..11 ................................................. .11..11..11 ..................................................... ,11,.ee,..eee,.e,e,eo,,.e,eee,.e,e,,.,,,.,..e,.eee,.ee,,e.,...,-
Fig. 14.
P_SURF
analysis ofthe fiatplatein Fig. 12
41
1 P_: 11 _ A SI_ r FILE: MB17-AZ_.A_ 1_L_(%)99.88475 _
FILE: 30 0_._ Flux x _
03RIT313L_ .12345 a/: _
x _
Oh1
1
2
3
4
5
6
7
8
9
Dec]
I
I
I
I
I
l
I
Z
Z
90
. e ..°.
o,...*o..o.o.ooeeooeeoe
85
o.°.°..,oe,.o.,,oo,
80
...°°o..o....o..........°ee..o°o.o..*...,.e....o
75
..............................
_
- 0.11532Z-02
Lnc:L-mem_ J_-cm 0 to max _
10
- 0.26100E-04
1_
12
Z
Z
1"
e..oo°.,°....ooo,.,
3.3 14 I
I
15
16
Z
I
o o.o,..oo....o..
o.o..°.°o.I*o°oe.eeeo.e....°..e.......o..°,.. .......,,,°,..o°
1 ......
1..i
.....
1..1.1
............
70 ..................................................... 1.1 ........ 65 ........................................i..i..2..2.!.212........ 60 ......................1. •.2.• .1..I.•.1.....I..2..2.51/2......... 55 ..........................1 .............1..21.3.15.41........... 50 .................................11..3..2..11.51.3.............. 45 ..........................i.• .2..1......31.41.1 ................. 40 .....................................5..3 ..........i ............ 35 ..............................2.,11..i...........i .............. 30
°°.
....
°..
°,,°°.,..o.,o°,o..,
....
°..,°°.°,°o,,,.°
,o°,.,.,°o°°.,°
25 .....................................i ...........I."............ 20 ..................1 .....................1 ....................... 15 ...........................................i..1..1.1............ 10.,
• • ..°°.,
• • • • • .,°,
°,,°,..o.,
..,..eo.o,.,
,.°°°_..,_°.,o,,e.
°.,
.°.
5
,,,,.°,o,..°.°°..o...,o,,.°°,,,..°,,.°,.,°°o°°.°.o._....°.°°
,,°°
0
,,,,.,,.oo,.o°,,.°....°.°eo.,
,,°.
,°,°.°.°,,..°o,.°,o°°.,°..°,,.,
3.00E-052.50E-052.00E-05Area
[sq
m]
1.50E-05 1.00E-055.00E-06 O.OOE+O0.
80 40 Obliquity [degrees]
Fig.
15.
P_SIIEF (1995
analysis
exposure
of a SSF module environment).
42
SD_SURF User's Manual
Appendix
A.
P_SURF Source Code
A-i
P_SURF
Usting
Listing
from Language Systems FORTRAN (Version 3.0 Tue, Nov 19, 1991) Sat, Sep 12, 1992 1:27 PM Options OFF: A BKG=e CASE CCD CCX CRAY DYNE EXTENDED F77 I2 LINEFEED MC68020 MC68881NOIMPLICIT OV R S SANE SYM T72 TRACE W X Z Options ON: ANSI C L SAVEALL U VAX 0001 0002 0003 0004 0005 0006 0@@7 0008 8009 @@10 0011 0012 @@13 0014 0015 0016 @017 @018 0019 0020 @021 0022 @@23 @@24 0025 0026 0027 0028 0029 0030 0031 0032 0033 0034 0035 0036 0037 0038 0039 0040 0041 0042 0043 0044 0045 0046 0047 0048
MC68848
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C P_SURF VER 1.6 8/23/92 C C C C MARTIN MARIEl-rA C C MANNED SPACE SYSTEMS C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C C P_SURF VER 1.6 will compute the Probability of No Penetration (PNP) by C space debris for a designated area on a spacecraft. P_SURF calculates C the flux (N) which penetrates the spacecraft multiplied by the exposed C area CA) and the exposure time (T) as a function of velocity & obliquity. C One data point on the surface represents the sum of alI projected C areas that can be hit by a particle at a certain velocity and obliquity, C multiplied by the fraction of the total flux that will cause a C penetration and the exposure time. C C P_SURF VER 1.6 works with BUMPERII Version 1.3 C C The code requires two files generated by other code as input. One C output file is from the A_SURF code. This file contains a selected C exposed area of a spacecraft, summarized in a matrix as a function of C velocity and obliquity. The other required file is the output file of C the RESPONSE portion of the BUMPERII code. This fiIe contains the C ballistic limit (minimum diameter to penetrate) as a function of velocity C and obliquity. Multiple shield designs may be included in the RESPONSE C output, and the PNP calculation may be performed for a specific shield C or for each shield in turn. C C The RESREAD and FLUX subroutines are taken directly from BUMPERII version 1.2a except for the COMMONPS.BLK instead of COMMON2.BLK Other modules were modeled after BUMPER for continuity. BUMPER was developed under the NASA contract 'Integrated Walt Design Guide and Penetration Control Plan' by M.A.Wright & A.R.Coronado. Note that peaks or waves in the area, flux or probability surfaces may be artifacts produced by granularity in the spacecraft model or threat models used in the GEOMETRY portion of BUMPER. Surface contours reflect the way BUMPER interogates the ballistic limit surface created by RESPONSE. P_SURF code was developed under the NASA contract Prediction and Analysis for HyperveIocity Impacts direction of N. Elfer.
A-1
'StructuraI Damage Study' under the
P_SURF
0049 0050 0051 0052 0053 0054 0055 0056 0057 0058 0059 0060 0061 0062 0063 0064 0065 0066 0067 0068 0069 0070 0071 0072 0073 0074 0075 0076 0077 0078 0079 0080 0081 0082 0083 0084 0085 0086 0087 0088 0089 0090 0091 0092 0093 e_94 0095 0096 0097 0098 0099 0100 0101 0102
C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C
Listing
Version 1.S corrects an error in the PNP calculation (found by Ben Hayashida). The FLUX from BUMPERII version 1.2a and 1.3 returns the Debris flux for the old environment, but the flux times the exposure time for the new environment. The was not recognized in version 1.4. Version
1.6
reads
Include
module
BUMPERII
COMMONPS
ver.
1.3
variable
Response
fires
list
alt = operating altitude , km asfite = the output Area Surface filenome binc = impact angle (beta) increment , deg conf = text description of wall configuration diam = critical diameter , cm etime = spacecraft exposure time , years flx = number of impacts per projected area per year of diameter D or larger inclin = orbital inclination, degrees idens = debris density, 1- constant density, Z-size function ienv = environment type, 1- 3SC 2_1&60_, 2- 7/90 memo it - current threat case itype = analysis type , i- debris, Z-meteoroids nb= number of angtes in the response array nc - number of wall configurations in the response array nee - the total number of exposed elements summed nr = number of element ranges to sum over nt = number of threat cases nv = number of velocities in the response array ptd = the property id associated with a11 elements of the ranges psfile = the Probabitity (Flux Area Time) Surface filename rsfile = the Response Surface filenome sflevel - solar flux level units = english or metric vr = impact (relative) velocity , km/sec vinc - impact (relative) velocity increment
, km/sec
Arrays area = array containing the vatue of the surface area for each element, sq-meters areas - the area surface containing the summed area fractions for each velocity and obliquity for oli elements in the specified element id ranges. (vr,beta) exposed = fist of the number of exposed elements for each threat angle fluxs = array containing flux corresponding to the diameters in the response surface (for each velocity and obliquity) geometry 1 array containing the values of the cosine of the impact angle for each exposed element for each threat angle. id = array containing the values of the element and property id for each etement 1- id 2- pid natmox = maximum
Flux*Area*Time
A-2
on one nots surface,
[impacts]
P_SURF
0103 0104 0105 0106 0107 0108 0109 0110 0111 0112 0113 0114 0115 0116 0117 0118 0119 0120 0121 012Z 0123 0124 0125 0126 0127 0128 0129 0130 0131 0132 0133 0134 0135 0136 0137 0138 0139 0140 0141 0142 0143 0144 0145 0200 0201 0202 0203 0204 0205 0206 0207 0208 0209 0210
C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C
Listing
natsflux*area*time surface as a function of (velocity,obliquity, pid), [impacts] her - array containing the range number for each element pids = PID (see scalar) number to process point = array of the element numbers corresponding to the elements in the geometry array. range = array containing the starting and ending elment id for each range to sum over 1-starting id 2- ending id response = array containing the values of the critical diameter as a function of impact angle and velocity. (vr,beta,pid) standm - shield stand-off, cm shden - shield density, g/cc shthkm - shield thickness, cm tnattotal flux * area * time for each PID and the areas array vwden - vessel wall density, g/cc vwthkm - vessel wall thickness, cm
Main
Program
Variable
List
Scalers answer - user input areae - the area times ob - Obliquity for the
LOGICAL
the threat probability current threat/element.
FIRST
C CHARACTER*80
ANSWER
C INTEGER*2 C C C
REAL*4
PROB
INCLUDE C C C
IC
'COMMONPS.BLK'
Initialize
the Velocity
increment
and number
of velocities
VINC-Q.25 NV=68 C C C
Initialize
the Obliquity
increment
BINC-5.0 NB-19
IBATCOM
= 0
A-3
and number
of angles.
P_SURF
0211 0212 0213 0214 0215 0216 0217 0218 0219 0220 0221 0222 0223 0224 0225 0226 0227 0228 0229 0230 0231 0232 0233 0234 0235 0236 0237 0238 0239 0240 0241 0242 0243 0244 0245 0246 0247 0248 0249
ezse 0251 0252 0253 0254 0255 0256 0Z57 0258 0259 0260 0261 0262 0263 0264
C C C
Write
header
Listing
to screen and read in orbital
parameters
CALL PSINPUT C C C
Read in
the
A_SURF output
file
CALL ASREAD C C C
Calculate
the
total
effective
DO 10 I2-1,NB DO 10 II-I,NV taeff - taeff 10 CONTINUE C C C C C
Read
exposure
+ AREASCII,12)
in the Solar flux data
IF ( ISoI.EQ.I.OR.ISoI.EQ.2 C C C C C
area.
)CALL
Read in the RESPONSE output rite. the RESREAD subroutine in BUMPER.
SOLREAD
This
is
identical
to
CALL RESREAD C C C
Verify
the
Response
file
has the
same increments
CBINC.NE.5.0 .OR. VINC.NE.O.25 ) THEN WRITE C6,*)'RESPONSE FILE HAS DIFFERENT STOP ELSE CONTINUE ENDIF
as the
Area_Surface
IF
C C C
Check
array
size
and
FORMAT THAN AREA_SURF!'
set to A_SURF size
IF
CNV.LT.68 .OR. NB.LT.19 ) THEN WRITE (6,*)'RESPONSE FILE IS SMALLER THAN AREA_SURF!' WRITE (6,*) NV,NB STOP ELSE NV-68 NB-19 ENDIF C C C C C
Determine
the
If
of
number IF
RESPONSE PIDs cases
CNC.EQ.1) PIDS(1)=I
CNC) is
to
process.
only
one then
THEN
A-4
proceed.
P_SURF
0265 0266 8267 8268 0269 0270 0271 0272 0273 0274 0275 0276 0277 0278 0279 0280 0281 0282 0283 0284 0285 0286 0287 0288 0289 0290 0291 0292 0293 0294 0295 0296 0297 0298 0299 0300 0301 0302 0303 0304 0305 0306 0307 0308 0309 0310 0311 0312 0313 0314 0315 0316 0317 0318
WRITE ( 6,20 ) FORMAT (/lX,'The
20 C C C C
one case
Listing
in
the
RESPONSE file
will
For multiple PIDs select one or all. If only one, NC is Write number of PIDs and first PID in A_SURF to screen. ELSE WRITE ( 6,25 ) 25 FORMAT (/lX,'The WRITE ( 6,26 ) 26 FORMAT (/1X,'The
be used' set
NC Number of PIDs in the RESPONSE file PID first PID processed by A_SURF was
to
is
1.
',I4) ',I4)
C 30 1 2
WRITE ( 6,30 ) FORMAT (//lX,'Enter to use the A_SURF PID. ' /IX,'Enter the PID number to use o specific PID.', /lX,'Enter
to use all PIDs.')
C 35
4@
45
READ ( 5,35 ) ANSWER FORMAT (A) IF ( ANSWER(I:I).EQ.' ' ) THEN PIDS(1) = PID NC=I ELSE IF (ANSWER(I:I).EQ.'A' .OR. ANSWER(I:I).EQ.'o') DO 4@ II=I,NC PIDS(I1) = II CONTINUE ELSE READ ( ANSWER(I:80),45 )PIDS(1) FORMAT ( BN,I2 ) NC=I ENDIF
THEN
C ENDIF C C C C C C
Calculate the Flux surface the Response surface
using
the critical
diameters
from
DO 120 I3=I,NC NATMAX(13) = 0 TNAT(I3) = 0 DO 100 I2=I,NB DO 100 II=I,NV DIAM=RESPONSE(II,12,PIDS(I3)) CALL FLUX FLUXS(II,I2,I3) = FLX C C C C
Calculate the FLUX x AREA x TIME surface (NOTE that FLUX returned FLUX x TIME for the new environment) IF (IEnv.EQ.I) THEN NATS(II,12,13) = FLUXS(II,IZ,13)*AREAS(II,12)*ETIME ELSE
A-5
)
P_SURF
0319 0320 0321 0322 0323 0324 0325 0326 0327 0328 0329 8330 0331 0332 0333 0334 0335 0336 0337 0338 0339 0348 0341 0342 0343 8344 0345 0346 0347 0348 0349 0350 0351 0352 0353 0354 0355 0356 0357 0358 0359 0360 0361 0362 0363 .0364 0365 0366 0367 0368 0369 0378 0371 %372
NATS(II,I2,I3) END IF C C C
Calcutate
the total TNAT(13)
C C C
Listing
= FLUXS(II,I2,I3)*AREAS(II,I2)
FLUX x AREA
x TIME
= TNAT(13)+NATS(II,I2,13)
Find the max NAT for the PID number IF (NATMAX(I3).LT.NATS(II,I2,I3)) NATMAX(13) = NATS(II,12,13) ELSE CONTINUE ENDIF
THEN
CONTINUE
I_W_
Calculate
PNP for Ranges
PNP(I3)
in Exposure
Surface
= (DEXP(-TNAT(I3)))*I_.DO
120 CONTINUE C C
Print out the carpet CALL
C C C C
CARPETPLOT
Print out the flux x area x time surface in a comma delimited format to be read by spreadsheets. CALL
C C C C C C
plot
Close
NATTEXT
summary
file
CLOSE ( UNIT-I@,STATUS-'KEEP' ) WRITE( 6,60003 ) PSFILE FORMAT( /' The PNP Surface file is complete.'/ ' filename: ',A ) C C C
Finished
END C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C SUBROUTINE CARPETPLOT C
A-6
P_SURF Usting 0373 0374 0375 0376 0377 0378 0379 0380 8381 0382 0383 0384 0385 0386 0387 0388 0389
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C This subroutine takes the NAT surface and writes out a text C based carpet plot to the output file C C ARRAY LIST C C CARPET = LINEAR ARRAY IN VELOCITY OF SYMBOLS FOR CARPET PLOT C NATINC - LINEAR ARRAY OF INCREMENTS TO BREAK CARPET PLOT C C C C CHARACTER*I CARPET(70) REAL*8 NATINC(6) INCLUDE 'COI4N)NPS.BLK' DO 2500
0445 0446 0447 0448 0449 0450 0451 0452 8453 0454 0455 0456 0457 0458 0459 0460 0461 0462 0463 0464 0465 0466 0467 0468 0469 0470 0471 0472 0473 0474 0475 0476 0477 0478 0479 0480
Calculate
1990 C C Write C
2000 i
2005
2010 1
2020 1
13-1,NC the increments
in the carpet plot
DO 1990 3-1,6 NATINCC3)=3*NATMAXCI3)/6.DO CONTINUE out
header
information
to
screen
and
to
file
WRITE (10,2000) PIDSCI3), RSFILE WRITE (6,2@@@) PIDS(I3), RSFILE FORMAT(IHI,/,IX,'RESPONSE PID: ',13,5X, 'RESPONSE FILE: ',A) WRITE (10,2005) ASFILE WRITE (6,2@@5) ASFILE FORMAT(IX,'A_SURF FILE: ',A) WRITE (10,2010) PNP(I3),TNAT(I3) WRITE (6,2010) PNP(I3),TNAT(I3) FORMATClX,'PNP(%)',FIO.S,SX, 'Total Flux x Area x Time (NAT) =',E12.5) WRITE (10,2020) NATMAX(13) WRITE (6,2020) NATMAX(13) FORMAT(IX,'CONTOURS .12345 at equal increments ' 0 to max NAT -',E12.5)
C
2030
2040
2050 C
WRITE WRITE FORMAT WRITE WRITE FORMAT WRITE WRITE FORMAT
( 10,2030 ) ( 6,2030 ) ( /SX,19X,'IMPACT VELOCITY ( 10,2040 ) (I,I=i,16,1) ( 6,2040 ) (I,I-1,16,1) ( lX,'ObI',IX,1614 ) ( 10,2050 ) ( 6,2050 ) ( IX,'Deg',IX,16(' I'))
A-7
km/s',/)
from',
P_SURF
0481 0482 0483 0484 0485 0486 0487 0488 0489 0490 0491 0492 0493 0494 0495 0496 0497 0498 0499 0500 0501 0502 0503 0504 0505 0506 0507 0508 0509 0510 0511 0512 0513 0514 0515 0516 0517 0518 0519 0520 0521 0522 0523 0524 0525 0526 0527 0528 0583 0584 0585 0586 0587 0588
C C
Calculate
the carpet
Listing
plot
DO 2500 I2-1,NB DO 2380 II-I,NV C C C
To plot with
max obliquity
at the top, 14 is substituted
for 12
14-20-12 C
2300 C C Write C C
2400
IF (NATS(II,I4,I3) .LE. NATINC(1)) CARPET(I1)-'. ' ELSE IF (NAT5(II,I4,I3) .LE. CARPET(I1)='I' ELSE IF (NATS(II,I4,I3) .LE. CARPET(I1)=' 2' ELSE IF (NATS(II,I4,I3) .LE. CARPET(I1)=' 3' ELSE IF (NATS(II,I4,I3) .LE. CARPET(I1)-'4' ELSE IF (NATS(II,I4,I3) .LE. CARPET(I1)=' 5' ENDIF CONTINUE out
the
contour
THEN NATINC(2))
THEN
NATINC(3))
THEN
NATINC(4))
THEN
NATINC(5))
THEN
NATINC(6))
THEN
marks
Write (6,2400) (14-1)*5,(CARPET(II),II=I,64) Write (10,2400) (I4-1)*5,(CARPET(I1),I1-1,64) FORMAT ( IX,13,1X,64AI )
C 2500 C
CONTINUE RETURN
C END C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C SUBROUTINE NATTEXT C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C Write NAT array for 5 degree and 0.25 km/sec increments C C INCLUDE 'COMMONPS.BLK' REAL*4 NATK(19) C C C C DO 700 I=I,NC
A-8
P_SURF
0589 0590 0591 0592 0593 0594 0595 0596 0597 0598 0599 0600 0601 0602 0603 0604 0605 0606 0607 0608 0609 0610 0611 0612 0613 0614 0615 0616 0617 0618 0619 0620 0621 0622 0623 0624 0625 0626 0627 0628 0629 0630 0631 0632 0633 0634 0635 0636 0637 0638 0639 0640 0641 0642
Listing
WRITE ( 10,600 ) PIDS(I) FORMAT ( 1H1,/,IX'RESPONSE PID: ',I3,/ ) WRITE ( 10,630 ) PIDS(1),(J,J-0,90,5) FORMAT (IX,16.2,19(',',I12.2) )
600 630
DO 690 K-I,NV DO 648 J-l,19 NATK(J)=NATS(K,J,I) CONTINUE VR-K*VINC WRITE(10,650) VR,(NATK(J),3-1,19) FORMAT (F6.2,19(',',E12.4) ) CONTINUE CONTINUE
640
650 690 700 C
RETURN C END C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C SUBROUTINE RESREAD C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
Resreod reads in the output from the RESPONSE code. consists of the critical diameter data as a function id, impact angle, and impact velocity.
note:
Variable
for variables contained listing for definition.
in the common
list
answer = character string represnting user itf = analysis type for rfile rfile = response output filenume ienvr = environment for response file
C C
block
input
CHARACTER LENGTH*2 CHARACTER*80 ANSWER,RFILE,Form CHARACTER*46 A46 CHARACTER*15 B15A, B15B CHARACTER*8 C8A, C8B CHARACTER*2 D2 REAL*4 STND INTEGER*4 ITF, ITA, IC, ICT, ICB, IPF !!! CHANGES TO BE COMPATIBLE WITH BUMPERII vet1.3 INTEGER*2 INTEGER*2
ITF, ITA, IENVR
IC,
ICT,
ICB,
C
A-9
IPF,
IPFUNC3
This output of property
refer
to the main
P_SURF
0643 0698 0699 0700 0701 0702 0703 0704 0705 0706 0707 0708 0709 0710 0711 0712 0713 0714 0715 0716 0717 0718 0719 0720 0721 0722 0723 0724 0725 0726 0727 0728 0729 0730 0731 0732 0733 0734 0735 0736 0737 0738 0739 0740 0741 0742 0743 0744 0745 0746 0747 0748 0749 0750
Listing
INCLUDE 'COMMONPS.BLK' IF (IBOTHS.EQ.2) GOTO 60 C C C
Read tn
the
RESPONSE output
filenome
, set
default
to
resp.dat
IF (INDEX(ROOTFILE,' ') .EQ. O)ROOTFILE='STATION.' ANSWER-ROOTFILE(I:INDEX(ROOTFILE,'.'))//'RSP' 30T - INDEX( ROOTFILE, ' ' ) WRITE (LENGTH, '(I2)' )30T+3 FORM-'(/IX,"RESPONSE OUTPUT FILENAME (=",A'//LENGTH//
•
V
,
WW)
>
I!
,$)'
10 WRITE ( 6,FORM )ANSWER READ ( 5,30 ) RFILE 30 FORMAT CA) C C C
C C C
!!t!
THIS OPEN FOR THE MAC WILL GIVE THE NORMAL FINDER DIALOG BOX. THE DIRLIST METHOD IS SKIPPED IF (RFILE(I:I).EQ.'?') THEN OPEN ( UNIT=Z3,FILE=*,STATUS='OLD',FORM='UNFORMATTED',ERR=48 INQUIRE(UNIT-23,NAME=RFILE) GOTO 60 END IF
!!!!
)
END OF MAC OPEN IF (RFILECI:I).EQ.'?') THEN CALL DIRLIST GOTO 10 END IF IF ( RFILE(I:4).EQ.' ' ) THEN RFILE-ANSWER ELSE ROOTFILE = RFite(I:INDEX(RFiIe,'.')) ENDIF
C IFCIBATCOM.EQ.1) WRITEC13,'(A)') RETURN END IF C C C
THEN RFILE
Open the rite OPEN
( UNIT=23,FILE=RFILE,STATUS='OLD',FORM-'UNFORMA1-FED',ERR-48
C GOTO C C C
Error
60
control
on open
48 WRITE ( 6,50 ) 50 FORMAT ( /IX,'UNABLE GO TO 10 C C
Read
in the analysis
TO OPEN
FILE'
)
type and the number of property
A-IO
cases.
)
P_SURF
0751 9752 9753 @754 0755 0756 0757 0758 0759 0760 @761 0762 0763 0764 0765 0766 0767 0768 0769 0770 0771 0772 @773 0774 0775 8776 0777 0778 0779 0780 0782 @782 0783 0784 0785 0786 0787 8788 0789 9799 0791 0792 0793 0794 9795 0796 0797 0798 0799 0800 0801 0802 0803 0804
C C C C C
C C C
Listing
60 READ (23) ITYPEIN,ITF,IDens,NC I!I NO ERROR CHECKING ON IENVR 60 READ (23) ITYPEIN,ITF,IENVR,IDens,NC WRITE(6,*)'ITYPEIN,ITF,IDens,NC' WRITE(6, *) ITYPEIN,ITF,IDens,NC IF (ITYPEIN.EQ.3.AND.IBOTHS.EQ.I) ITYPE=I IF (IDens.EQ.l) THEN WRITE (6,63) 63 FORMAT (/5X,' Constant density threat') ELSE IF (IDens.EQ.2) THEN WRITE (6,64) 64 FORMAT (/5X,' Variable density threat') END IF Check that the response
file is the correct
analysis
type
IF ( ITF.NE.ITYPE ) THEN IF ( ITYPE.EQ.1 ) THEN WRITE ( 6,70 ) 70 FORMAT ( /IX,'DEBRIS ANALYSIS SPECIFIED IN GEOMETRY FILE ', 1 'BUT RESPONSE FILE IS FOR METEOROIDS ') ELSE WRITE ( 6,80 ) 8Q FORMAT (/IX,'METEOROID ANALYSIS SPECIFIED IN GEOMETRY FILE', 1 ' BUT RESPONSE FILE IS FOR DEBRIS' ) END IF C 90 1
WRITE ( 6,90 ) FORMAT ( /IX,'DO YOU WISH TO CONTINUE '(=NO) > ',$) READ ( 5,30 ) ANSWER
WITH
GEOMETRY
OPTION
C IF ( ANSWER(I:I).EQ.'Y' GO TO 10 ELSE STOP END IF
.OR. ANSWER(I:I).EQ.'y'
C END IF C C C C C C C C C C C C
Read in the impact angle
information
READ (23) NB,BINC WRITE(6, *) 'NB,BINC' WRITE(6,*) NB,BINC Read in the impact
velocity
READ (23) NV,VINC WRITE(6,*) 'IMPACT WRITE(6,*) NV,VINC Initialize
RESPONSE
information
VELOCITY,
VEL INCR.'
to 0.0
A-11
) THEN
',
P_SURF
0805 0886 0887 0808 0809 0810 8811 8812 8813 8814 0815 0816 0817 0818 8819 0820 0821 0822 0823 0824 0825 0826 0827 0828 0829 0838 0831 0832 8833 0834 0835 0836 0837 0838 0839 0840 0841 0842 0843 0844 0845 0846 0847 0848 0849 0850 0851 0852 0853 0854 @855 @856 0857 0858
Listing
C DO 2_ I-1,NC DO 150 J-1,NB DO 1_ K-1,NV RESPONSE (K,J,I)-e. 1(_ CONTINUE 150 CONTINUE 200 CONTINUE C C C C C C C
Read in
the critical
diameter
data
WRITE(6,*) 'NC,NB,NV' WRITE(6,*) NC,NB,NV WRITE(6,*) 'RESPONSE(K,J,I)' Loop thru the property id's DO 4_ I-I,NC
C C
Loop thru the impact angles DO 3_ 3=I,NB
C C
Loop thru the impact velocities DO 250 K-1,NV
C C C Z50 400 C C C
411
Store the critical diameter READ (Z3) RESPONSE(K,J,I) WRITE(6,*) RESPONSE(K,J,I) CONTINUE CONTINUE CONTINUE IF
(INPUTCD.EQ.2)
in response
CALL SETDIAMS
READ ( 23,END-440,ERR-440 ) A46 WRITE(6,*) 'A46' WRITE ( 6,'( //IX,A)' ) A46 WRITE ( 10,'( //IX,A)' ) A46 READ ( 23 ) C8A,ITA,C8B,ICB,UNITS WRITE ( 10,'(A,I4)') ' Threat (I Debris, 2 Meteoroid) WRITE ( 10,'(A,I4)') ' Density (I Constant, 2 Function) WRITE ( 10,'(A,I4)') ' Number of PID Cases WRITE ( 10,'(2A)') ' Units ',UNITS WRITE(6,*) 'C8A,C8B' WRITE(6,*) C8A,C8B WRITE ( 6,'(A,I4)') ' Threat (I Debris, 2 Meteoroid) WRITE ( 6,'(A,I4)') ' Density (i Constant, 2 Function) WRITE ( 6,'(A,I4)') ' Number of PID Cases WRITE ( 6,'(2A)') ' Units ',UNITS DO 420 I-I,ICB READ ( 23 ) ICT,D2,BISA,BISB,IPF,IPFUNC3 WRITE (10,411) I FORMAT( /IX,'PID NUMBER ',I4 )
C IF (ICT.EQ.2) THEN IF ( IPF.EQ.I ) THEN
A-12
',ITA ',IDens ',ICB
',ITA ',IDens ',ICB
P_SURF
0859 0860 0861 0862 0863 0864 0865 0866 0867 0868 0869 0870 0871 0872 0873 0874 0875 0876 0877 0878 0879 0880 0881 0882 0883 0884 0885 8886 0887 0888 0889 0890 0891 0892 0893 0894 0895 0896 0897 0898 0899 0900 0901 0902 0903 0904 0905 0906 0907 0908 0909 0910 0911 0912
485 486 487 488 484 489 490 491 492 493 494 495 496 497
Listing
WRITE (10,485) ELSE IF ( IPF.EQ.2 ) THEN WRITE (10,486) ELSE IF ( IPF.EQ.3 ) THEN WRITE (10,487) ELSE IF ( IPF.EQ.4 ) THEN WRITE (10,488) ELSE IF ( IPF.EQ.5 ) THEN WRITE (10,484) ELSE IF ( IPF.EQ.6 ) THEN WRITE (10,489) ELSE IF ( IPF.EQ.7 ) THEN WRITE (10,490) ELSE IF ( IPF.EQ.8 ) THEN WRITE (10,491) ELSE IF ( IPF.EQ.9 ) THEN WRITE (10,492) ELSE IF ( IPF.EQ.10 ) THEN WRITE (10,493) ELSE IF ( IPF.EQ.11 ) THEN WRITE (10,494) ELSE IF ( IPF.EQ.12 ) THEN WRITE (10,495) ELSE IF ( IPF.EQ.13 ) THEN WRITE (10,496) ELSE IF ( IPF.EQ.14 ) THEN WRITE (10,497) END IF END IF FORMAT ( /lX,'ORIGINAL PENETRATION FUNCTION') FORMAT ( /lX,'PEN4 PENETRATION FUNCTION') FORMAT ( /lX,'REGRESSION PENETRATION FUNCTION') FORMAT ( /lX,'COUR-PALAIS PENETRATION FUNCTION') FORMAT ( /lX,'BOEING INTERP PENETRATION FUNCTION') FORMAT ( /lX,'DEVELOPMENTAL6, USER INPUT') FORMAT ( /lX,'DEVELOPMENTAL7, USER INPUT') FORMAT ( /lX,'DEVELOPMENTAL8, USER INPUT') FORMAT ( /lX,'DEVELOPMENTAL9, USER INPUT') FORMAT ( /lX,'DEVELOPMENTAL10, USER INPUT') FORMAT ( /1X,'DEVELOPMENTAL11, USER INPUT') FORMAT ( /1X,'DEVELOPMENTAL12, USER INPUT') FORMAT ( /lX,'DEVELOPMENTAL13, USER INPUT') FORMAT ( /1X,'DEVELOPMENTAL14, USER INPUT') WRITE ( 10,'( /A )' ) ' Configuration Shield Wa11' WRITE ( 6,* ) 'ICT,D2,B12A,B12B' WRITE ( 6,* ) ICT,D2,BI2A,B12B IF (ICT.EQ.I) CONF = 'Single Plate' IF (ICT.EQ.2) CONF = 'Double Plate' IF (ICT.EQ.3) CONF = 'Multiwa11' WRITE ( 10,'( lX,A,4X,2A )' ) CONF,BI2A,BI2B WRITE ( 6,'( lX,A,4X,2A )' ) CONF,BI2A,B12B READ ( 23 ) ShThk,VWThk,STND,ShDen(I),VWDen(I),ADEN WRITE ( 6,* ) 'ShThk,VWThk,STND,ShDen(I),VWDen(I),ADEN,I' WRITE ( 6,* ) ShThk,VWThk,STND,ShDen(I),VWDen(I),ADEN,I
A-13
P_SURF
0913 0914 0915 0916 0917 0918 09_ 0920 0921 09ZZ 0923 0924 0925 0926 0927 0928 0929 0930 0931 0932 o933 o934 o935 o936 o937 o938 o939 0940 o941 o942 o943 o944 0945 0946 0947 0948 0949 0950 0951 095Z 0953 e954 0955 e956 e957 e958 e959 e96o e961 0962 0963 0964 0965 0966
Usting
IF (ICT.EQ.3) THEN WRITE ( 10,'(A,A,F8.4)')
' Combined ' of All Shields
+
WRITE
( 10,'(A,F8.4)') ' ,STND WRITE ( 6,'(A,A,F8.4)') '
+
C C C C
+
WRITE
( 6,'(A,FS.4)')
+
C +
410 C
4-
GOTO 410 END IF WRITE ( 6,'(A,F8.4)') ' IF CSHTHK.NE.g.o) WRITE ( 10,'(A,F8.4)') WRITE ( 10,'(A,FS.4)') ' WRITE ( 6,'(A,F8.4)') ' IF (ICT.NE.3) THEN IF (SHTHK. NE .g.g.AND. STND. WRITE ( 10,'(A,FS.4)') WRITE ( 6,'(A,F8.4)') ' END IF IF ( Units .EQ. ' ENGLISH ShThkM(1) - ShThk*2.54 VWThkM(I) - VWThk*2.54 ADAR(1)-ADEN/. 0142233 ELSE ShThkM(1) - ShThk VWThk ADAR(I)-ADEN END IF With
or
without
450 IF (IBOTHS.EQ.I)
Shield Thickness ' Vessel Vessel NE .0.8) '
Standoff = ',STND Standoff _ ',STND
' ) THEN
.RSP file
RETURN
Close the file and return ( UNIT-Z3,STATUS-'KEEP'
Write
Rfile to summary
)
file
WRITE ( 10,5_ )RFILE 5_N_ FORMAT(IX,'RESPONSE OUTPUT
FILE = ',A )
!!!! WRITE(RSFILE,'(BN,A)')RFILE
A-14
= ',ShThk
Shield Thickness ,, ',ShThk Wall Thickness = ',V_F[hk Wall Thickness = ',VWThk
30 MLI
CLOSE
= '
Combined Areal Density', ' of All Shields - ',ADEN Total Standoff = '
READ ( Z3 ) A46 WI_ITE ( 10,'( 4X,A)' ) A46 C WRITE ( 6,'( 4X,A)' ) A46 420 CONTINUE GO T0450 440 WRITE ( 10,42 ) 4Z FORMAT ( /ZX,' No Header following
C C C
Standoff
1
V_k_(I)
C
' ,STND
Total
Areal Density', _ ',ADEN
' )
P_SURF Listing 0967 0968 0969 0970 0971 0972 0973 0974 0975 0976 0977 0978 0979 0980 0981 098Z 0983 0984 0985 0986 0987 0988 0989 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1073 1074
RETURN C END C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C Subroutine ASREAD C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC This subroutine versus velocity
opens and reads the table of exposed and obtiquity created by A_SURF.
areas
CHARACTER*80 ANSWER
INCLUDE
'COMMONPS.BLK'
Read in the ASF filename
, set default
to DATA.ASB
10 WRITE ( 6,'(/IX,"Areo_Surfoce Binary Output 1 " :",$)') READ ( 5,'(A)' ) ANSWER
File"
l!!!
THIS OPEN FOR THE MAC WILL GIVE THE NORMAL FINDER DIALOG BOX. THE DIRLIST METHOD IS SKIPPED IF (ANSWER(I:I).EQ.'?') THEN OPEN ( UNIT-2,FILE=*,STATUS='OLD',FORI&='UNFORMATTED', * READONLY,ERR-10 ) INQUIRE(UNIT-Z,NAME=ASFILE) GOTO 40 END IF
!!i!
END OF MAC OPEN
IF ( ANSWER(I:I).EQ.'
' ) ANSWER-'DATA.ASB'
C WRITE(ASFILE,'(BN,A)')ANSWER C C C
Open
*
the file
OPEN ( UNIT=2,FILE=ANSWER,STATUS='OLD',FORM='UNFORMA1-FED' ,READONLY,ERR=IO )
C 40 CONTINUE C C C
Read in the analysis
type and the number
READ (2) ITYPE,NR,PID,AREATOT IF( NR.GT.IRNGS ) THEN
A-15
of ranges
P_SURF Listing 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1899 1100 1101 1102 1103 1104 1105 1106 1107 1108 1189 1110 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 i128
60001
C C C
WRITE( 6,60001 ) FORMAT( /' ---ERROR--* ,' exceeded.' STOP ENDIF
Read in
the
READ (Z) C C C
Read in
Read in
the
the
Read the
ranges
impact
Close
the
CLOSE C C C C C
Write
angle
impact
information
velocity
information
NV,VINC
Area
READ (2) C C C
Ranges was'
NB,BINC
READ (2) C C C
number of
((RANGE(I,3),I=I,Z),3=I,NR)
READ (Z) C C C
The maximum )
Surface
array
((AREAS(I,3),I-1,NV),3=I,NB)
file C UNIT-Z,STATUS-'KEEP'
A_SURF
file to output
)
file
WRITE ( 10,600 ) ASFILE 600 FORMAT ( IX,'A_SURF BINARY
OUTPUT
FILE - ',A )
C C C
Write
6Z1 C C C
the number
of ranges
and the Property
ID.
WRITEC 10,6Z1 ) NR,PID,AREATOT FORMAT( IX,'RANGES-',IZ,' PIDI',I9, 1 ' EFF. AREA (sq.m) -',FL?..5) Write the start and end Element
ID for each range.
DO 625 I-I,NR WRITE ( 10,622 ) I,RANGE(I,I),RANGE(Z,I) 622 FORMAT(IX,'Range ',I2,' START: ',I12,'END: 625 CONTINUE C RETURN C END C C C C C
A-16
',IiZ)
P_SURF Listing 1129 1138 1131 1132 1133 1134 1135 1136 1137 1138 1139 1148 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1208 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C SUBROUTINE PSINPUT C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C PSINPUT writes the program header to the screen and reads in the C summary ouput filename. It also determines the spacecraft exposure C time and operating altitude. C C C C note: for variables contained in the common block refer to the main C listing for definition C C Variable list C C answer = character string representing user input C Psfile = output filenonme C C C INCLUDE 'COPiuK)NPS.BLK' C CHARACTER*2@ BUM'FrM CHARACTER*8@ ANSWER C C C C C
Write
header
to screen
and summary
file
WRITE ( 6,10 ) 10 FORMAT ***************************************** I 'Space Debris SURFace', 2 //IX,5X,'Ver. 1.5 8/23/92',/IX,SX,'for 3
C C C
BUMPERIIvl.2Q',//IX,
'_____'_
Read
in output
filename,
15 WRITE ( 6,20 ) 20 FORMAT ( /1X,'OUTPUT READ ( 5,30 )PSFILE 30 FORMAT (A)
set default
FILENAME
to SDSURF.PS
(CR=SDSURF.PS)>',$)
C IF ( PSFILE(I:I).EQ.' C C C C C
Open
psfile
!!!! PUT CREATOR='XCEL' OPEN *
' ) PSFILE='SDSURF.PS'
OR
'MSWD' IN OPEN
STATEMENTS
( UNIT=IO,FILE=PSFILE,STATUS='NEW',IOSTAT=IER, CREATOR='XCEL',ERR=4@,RECL=256 )
A-17
ON MAC
P_SURF
Listing
A
1237 1238 1239 1248 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277 1278 1279 1280 1281 1282 1283 1284 1285 1286
### FORTRAN - Warning - This feature is an extension to VAX FORTRAN File "macii_p_surf16.f"; Line 912 # C GO TO 70 C C Error controt C 48 IF ( IER.EQ.2013 ) THEN WRITE ( 6,50 ) 50 FORMAT ( /IX,'FILE ALREADY EXISTS OK TO OVERWRITE (CR=YES,$)>') READ ( 5,30 ) ANSWER C IF ( ANSWER,I:I).EQ.'Y' .OR. ANSWER(I:I).EQ.' ' ) THEN OPEN ( UNIT-10,FILE=PSFILE,STATUS='UNKNOWN',IOSTAT=IER, 1 ERR-48) REWIND 10 ELSE GO TO 15 END IF ELSE WRITE ( 6,60 ) 60 FORMAT (/IX,'UNABLE TO OPEN FILE ' ) GO TO 15 END IF C 70 CONTINUE C C WRITE ( 10,75 ) 75 FORMAT *11X******************************,**lX,3** 1 'Space Debris SURFace', 2 //IX,SX,'Ver. 1.5 8/23/92',/IX,SX,'for 3 ******************************** C C C C
Set anatysLs
BUMPERIIvI.Za',//IX,
type to 1 ,debris)
ITYPE-1 C C C C C C
!!!!!!!!!!!FROM Determine
BUHPERII
Environment
Verl.2.a
Definition,
SHIELD
INPUT!!!!!!!!!!!
set default
to I (original)
51 WRITE ( 6,52 ) 52 FORMAT (/IX,'ENVIRONMENT ?',/2X,'I-]SC 2_w_1&6_ 1 '2- 7/90 MEHO',/lX,'ANSWER 1 OR 2 > ',$) C READ ( 5,53 53 FORMAT (A)
) ANSWER
C IF ( ANSWER(I:4).EQ.'
' ) THEN
A-18
',/2X,
P_SURF
1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299 1300 1301 1302 1303 1304 1305 1306 1307 1308 1309 1310 1311 1312 1313 1314 1315 1316 1317 1318 1319 1320 1321 1322 1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1348
IEnv=l ELSE READ ( ANSWER(l:80),54 $4 FORMAT ( BN,I1 ) END IF C C C
Check that
956
input
Listing
)IEnv
was correct
IF ( IEnv. EQ.1 .OR. IEnv.EQ.2 ) THEN CONTINUE ELSE WRITE ( 6,956 ) FORMAT ( /IX,'INCORRECT INPUT' ) GO TO 51 END IF
365
CONTINUE IF ( ITYPE.EQ.I.AND.IENV.EQ.2 ) THEN 370 WRITE ( 6,380 ) 380 FORMAT (/IX,'SOLAR FLUX LEVEL ?',/,2X,'I-NOMINAL 1 '2-MINIHUM',/,2X,'3-CONSTANT',/,IX,'ANSWER READ ( 5,30 ) ANSWER
',/,2X, 1-3 > ',$)
C
90
385 390
C C C
IF ( ANSWER(I:4).EQ.' ' ) THEN ISOL=I ELSE READ ( ANSWER(1:80),90,ERR=370 ) ISOL FORMAT (BN,I4) END IF IF(IBATCOM.EQ.1) WRITE(13,'(A)') ANSWER IF ( ISOL.EQ.3 ) THEN WRITE ( 6,390 ) FORMAT (/IX,'SOLAR FLUX LEVEL (10"'4 Jy) (=70) READ ( 5,30 ) ANSWER IF ( ANSWER(I:4).EQ.' ' ) ANSWER='70.0' READ ( ANSWER(I:80),120,ERR=385 ) SFLEVEL IF(IBATCOM.EQ.1) WRITE(13,'(A)') ANSWER END IF
Check that
input was correct
IF ( ISOL.LT.1 GO TO 370 END IF END IF C C C
Determine
340 350
.OR. ISOL.GT.3
the spacecraft
exposure
) THEN
date,
set default
IF ( ITYPE.EQ.1.AND.IENV.EQ.2 ) THEN WRITE ( 6,350 ) FORMAT ( /lX,'DATE TO BEGIN EXPOSURE 1 (=1995) > ',$) READ ( 5,30 ) ANSWER
C
A-19
to 1995
( 1994-2025
)
> ',$)
P_SURF Listing 1341 1342 1343 1344 1345 1346 1347 1348 1349 1350 1351 1352 1353 1354 1355 1356 1357 1358 1359 136% 1361 1362 1363 1364 1365 1366 1367 1368 1369 1370 1371 1372 1373 1374 1375 1376 137.7 1378 1379 1380 1381 138Z 1383 1384 1385 1386 1387 1388 1389 1390 1391 1392 1393 1394
IF
C ANSWERCI:4).EQ.' ' ) DATE-1995. ELSE ISpot-IndexCANSWER,'.') IF (ISPOT.NE.O) GOTO 179 k-8e iblank-e do while (iblank.eq.O) if (ANSWER(k:k).ne.' iblank=l
then
goto 1110 end if k-k-1 continue end do ANSWER-ANSWERCI:K)//'.' READ ( ANSWER(l:80),180,ERR=34% FORMAT ( BN,D20.3 ) END IF
1110
179 180 C C C
')
THEN
Check
that date
is
within
)
DATE
range
IF
( DATE.LT.1994 .OR. DATE.GT.Z025 ) THEN WRITE ( 6,36% ) 36% FORMAT ( IX,'---ERROR--Date outside of GO TO 34% END IF IFCIBATCOH.EQ.I) WRITE(13,*) DATE END IF IF (IBOTHS.EQ.Z) GOTO 56 C C C
Determine
the spacecraft
105 WRITE ( 6,110 ) 110 FORMAT (/IX,'SPACE
1 READ ( 5,3e
exposure
STATION
time
EXPOSURE
range'
, set default
to 10 years
TIME (YEARS)
> ',$) ) ANSWER
C IF ( ANSWERCI:4).EQ.'
' ) ANSWER='10.0'
C READ ( ANSWERCZ:80),ZZO,ERR-105 12e FORMAT ( BN,DZe.e ) IF (ETIME.LT.O.) GOTO 105 IF(IBATCOW.EQ.1) THEN WRITE(Z3,*) ETIME GOTO 151 END IF
) ETIME
C 56 IF ( IEnv.EQ.1 ) THEN WRITE (10,57) 57 FORMATE' 3SC-Z_I AND 35C-6%_ FLUX EQUATIONS') ELSE WRITE (10,58) 58 FORHAT(' 7/17/90 MEMO FLUX EQUATIONS')
A-20
)
P_SURF
1395 1396 1397 1398 1399 1400 1401 1402 1403 1404 1405 1406 1407 1408 1409 1410 1411 1412 1413 1414 1415 1416 1417 1418 1419 1420 1421 1422 1423 1424 1425 1426 1427 1428 1429 1430 1431 1432 1433 1434 1435 1436 1437 1438 1439 1440 1441 1442 1443 1444 1445 1446 1447 1448
Usting
ENDIF
130 400 410 420 430
140
IF ( ITYPE.EQ.1 ) THEN WRITE ( 10,130 ) FORMAT ( /1)(,'MAN-MADE ORBITAL DEBRIS ANALYSIS') IF ( ISOL.EQ.I ) WRITE ( 10,400 ) FORMAT ( IX,'NOMINAL SOLAR FLUX LEVEL' ) IF ( ISOL.EQ.2 ) WRITE ( 10,410 ) FORMAT ( lX,'MINIMUM SOLAR FLUX LEVEL' ) IF ( ISOL.EQ.3 ) WRITE ( 10,420 ) SFLEVEL FORMAT ( lX,'SOLAR FLUX LEVEL = ',F8.3 ) IF (DATE.NE.O.) WRITE ( 10,430 ) DATE FORMAT ( IX, 'DATE TO BEGIN EXPOSURE = ',F8.3 ) ELSE WRITE ( 10,140 ) FORMAT ( 1)(,'METEOROID ANALYSIS' ) END IF
WRITE ( 10,150 )ETIME 150 FORMAT ( lX,'SPACECRAFT Read
in operating
altitude
EXPOSURE
TIME
, set default
(YEARS)
=',F8.3
)
to 500 km
151 IF (IBOTHS.EQ.2) GOTO 203 IF(IEnv.EQ.1) THEN AltMin = 350. AltMax = 550. AltMinnm=350.*O.53995680 A1tMaxnm=550.*O.53995680 ELSE AltMin = 100. AltMax = 500. AltMinnm=100.*O.53995680 AltMaxnm=500.*O.53995680 END IF C C C
ALT INTERNALLY
IS IN KILOMETERS.
160 WRITE ( 6,170 )AltMin,A1tMax 170 FORMAT(/IX,'OPERATING ALTITUDE('F4.0'-',F4.0'km) i (=388.92) ') WRITE(6,205) 205 FORMAT(' OR ENTER AN "E" OR "e" TO ENTER IN NMILES > READ (5,'(A)')ANSWER IF (ANSWER(I:4).EQ.' ') THEN ALT=388.92DO ALTNM=210.OODO
,,$)
IF (IBATCOM.EQ.1) WRITE ( 13,* ) ALT GOTO 189 END IF IF (ANSWER.EQ.'E'.OR.ANSWER.EQ.'e') THEN IF(IBATCOM.EQ.I) WRITE(13,'(A)') ANSWER(I:IO) WRITE(6,171)altminnm,altmaxnm 171 FORMAT(/IX,'OPERATING ALTITUDE('F4.0'-' ,F4.0'nmiles) i (=210.00) > '1
A-21
P_SURF
1449 1450 1451 1452 1453 1454 1455 1456 1457 1458 1459 1460 1461 1462 1463 1464 1465 1466 1467 1468 1469 1470 1471 147Z 1473 1474 1475 1476 1477 1478 1479 1480 1481 1482 1483 1484 1485 1486 1487 1488 1489 1490 1491 1492 1493 1494 1495 1496 1497 1498 1499 1500 1501 1502
READ C5,' CA) ' )ANSWER IF CANSWERCI:4).EQ.' ALTNH-210 ._DO ELSE READ (ANSWER(1:80),215) END IF IF (IBATCOM.EQ.I) WRITE ALT = ALTNM / 0.53995680 ELSE 215
C
')
Listing
THEN
ALTNM ( 13,* ) ALTNM
READCANSWERCI:80),215) ALT FORMAT( BN, D20.0 ) IF (IBATCOM.EQ.1) WRITE ( 13,* ) ALT END IF WRITE(6,*) 'ALT' ,ALT
Check
that altitude
is within
range
189 IF ( ALT.LT.AItMin .OR. ALT.GT.AltMox ) THEN WRITE ( 6,190 )AltMin,AltMax 190 FORMAT ( IX,'---ERROR--Attitude outside GO TO 168 END IF C C C
Write
altitude
to output
of range
',2F8.3
)
file
203 IF (IBATCOM.EQ.1) GOTO 204 IF ( Units .EQ. ' ENGLISH ' ) THEN IF (ANSWER.NE.'E'.OR.ANSWER.NE.'e') ALTNM=ALT* 8.53995680 WRITE ( 10,202 ) ALTNM 202 FORMAT ( IX,'OPERATING ALTITUDE (nmiles) = ',F8.3 ) WRITE ( 10,2_W_ ) ALT ELSE 2_
WRITE ( 10,2_ ) ALT FORMAT C IX,'OPERATING END IF
ALTITUDE
(km) = ',F8.3 )
C C C C
204 CONTINUE Finished RETURN
C END C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C SUBROUTINE FLUX C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C
Flux calculates
the meteoroid
or debris
A-22
flux for the given
critical
P_SURF
1503 1504 1505 1506 1561 1562 1563 1564 1565 1566 1567 1568 1569 1570 1571 157Z 1573 1574 1575 1576 1577 1578 1579 1580 1581 1582 1583 1584 1585 1586 1587 1588 1589 1590 1591 1592 1593 1594 1595 1596 1597 1598 1599 1654 1655 1656 1657 1658 1659 1660 1661 1662 1663 1664
C C
diameter
based on analysis
INCLUDE
Listing
type.
'COMMONPS.BLK'
C C C
C
IF (IEnv.EQ.I) THEN This flux definition meteoroids CALL Flux20001 ELSE This flux definition CALL Flux790 END IF
uses 3SC-200001
uses the 7/17/90
for debris
revision
and J5C-6000
for
memo to SSP 30425
RETURN C END CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C SUBROUTINE FLUX2_I C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C Flux calculates the meteoroid or debris flux for the given critical C diameter based on analysis type. C C C note: for variables contained in the common block referr to the main C listing for definition C C Variable List C C ddiam = diam in double precision , cm C ge = gravity focusing factor C intercept = intercept of the flux equation C mass = critical meteoroid mass, g C mden = meteoroid density, g/cc C re = earth's radius, km C slope = slope of the flux equation C C INCLUDE 'CONMONPS.BLK' C REAL*8 DDIAM,GE,INTERCEPT,LD,MASS,MDEN,PI,RE,SLOPE C PARAMETER CPI=3.141592653589793238D@) C mden C Set C C MDEN=O.50DO C C Calculate the focusing factor, equation C is from 3SC-30000
A-23
P_SURF
1665 1666 1667 1668 1669 1678 1671 167Z 1673 1674 1675 1676 1677 1678 1679 1688 1681 1682 1683 1684 1685 1686 1687 1688 1689 1698 1691 1692 1693 1694 1695 1696 1697 1698 1699 17_ 1781 1782 1783 1784 1785 1786 1787 1788 1709 1718 1711 1712 1713 1714 1715 1716 1717 1718
Listing
C RE-6378 .SDe GE-e. 568D_.
432De* CRE/CRE+A LT))
C C C
Convert
diam
to
double
precision
DDIAM-DIAM Calculate IF
the
flux
C ITYPE.EQ.1
) THEN
For debris use 3SC-Z_I, use stated equations for less then 1 ¢m , for those greater use third order curve for region up to S cm . The log of the flux varies to D Kesseler of 3SC.
linearly
between
4_
diameters fit of the
and S_
km according
L_OGIe(DDIAM) IF ( DIAM._I'.S.8 )DIAM-5.8 IF C DIAN.LE.5.8 ) THEN IF C DIAM.LT.I.8 ) THEN SLOPEm-e. _18DO*ALT-Z. 8Zl_De ELSE SLOPE--8._ZZDO*ALT-8.1400DO END IF INTERCEPT-+O._36DO*ALT-7.Z6De FLX-18.SDe**CSLOPE*LD+INTERCEP'O ELSE WRITE ( 6,1_ ) FORMAT ( /1X,'DIAMETER IS GREATER THAN S CM LIMIT') STOP END IF ELSE For
meteoroids
use 3SC-3_,
E-86g
',$
)
ASTFILE.EQ.'
' ) ASTFILE-'DATA.AST'
Open the text output ON VAX USE RECL-256
file for GEOREAD to use. TO WRITE ENTIRE RESULTS.
OPEN(UNITffil0,FILE-ASTFILE,
STATUSffi'NEW',RECL=256,ERR=600)
0523 0524 0525 Read
in
element
$
user
C
0522
f
' ) ASBFILE-'DATA.ASB'
OUTPUT
FORMAT(/' Text READ( 5,50002
0516
output
) ASBFILE
)
60003
0513
)
Binary
FORMAT( IF(
C
6,60002
FORMAT(/' READ(
0507
0526
altitude
list
INCLUDE
0452
0512
operating
C
0451
0505 0506
variables
listing
0443
0504
time,
Listing
sum
C
0444
0502
exposure
B - A_SURF
ranges
to
sum
over
B-8
)
)
1.8',
main
Appendix
0527
IC=0 C
0530
WRITE
0531
250
( 6,250
)
FORMAT(/IX,'
ONE
Area
Fraction
0532
1
IX,'
0533
2 3
IX,' INPUT THE STARTING ' EACH RANGE'/IX,' ENTER
0534 0535
270
IC=IC+I
0537
275
WRITE
0538
280
FORMAT
285 290
WRITE FORMAT
0541 0542
( 6,280
30
of
will
be
element
AND ENDING D OR
created
IDs
'/
selected.'/
ELEMENT ID FOR', WHEN DONE')
' IN
THE
TABLE.
')
( 5,30
ELEMENT
ID
:
',$)
) ANSWER
(A)
IF
( ANSWER(I:I).EQ.'
'
.OR.
ANSWER(I:I).EQ.'D'
) GO
TO
500
) GO
TO
305
then
Ending
C READ
0548
( ANSWER(I:80),300,ERR=285
300
FORMAT
0550
305
WRITE
0551
310
FORMAT
( BN, II2
) RANGE(I,
IC)
)
C
0552
( 6,310
READ
)
( IX,'ENDING ( 5,30
ELEMENT
: ',$)
ID
) ANSWER
C
0554 0555
Table
ranges
'RANGE',I4,
( 6,290 ) ( IX,'STARTING
FORMAT
0547
0553
the
C
0545
0549
of
)IC
( /IX,
READ
0543
0546
ALL
C
054O
0544
from
C
0536
0539
Listing
C
0528 0529
B - A_SURF
IF
( ANSWER(I:I).EQ.'
'
.OR.
ANSWER(I:I).EQ.'D'
C
0556
READ
0557
C
0558 0559
C
Check
(ANSWER(I:80),300,ERR=305) that
ending
id
>
RANGE
starting
(2, IC)
id
C
0560
IF
0561 320
0562
( RANGE(I,IC).GT.RANGE(2,
IC)
WRITE FORMAT
Starting
GO
0563 END
0564 0565
C
0566
C
0567
C
Next
0568
IF
if
have
C
0573
C
0574
C
0575
reached
) GO
,'be
Check 500
0576
that
ID
greater
TO
max
IF
values
270 A
maximum
of
',I2,'
ranges
processed.') were
input
( RANGE(I,I).EQ.0.0 WRITE FORMAT
330
.AND.
( 6,330 ) ( IX,'---ERROR---
RANGE(2,1).EQ.0.0 No
Range
to
the
Values
) THEN Input'
IC=0
0579
GO
0580 0581
END
0582
C
0583
C
0584
C
Set
ID')
allowed.
CONTINUE
0577 0578
not
( IC.LT.IRNGS
*
0572
THEN
275
WRITE(6,340) IC FORMAT(IX,'---WARNING---
340
0571
)
IF
Range
0569 0570
TO
( 6,320 ) ( 1X,'---ERROR---
the
TO
270
IF number
of
ranges
equal
B-9
number
read
in
)
will
'
Appendix
0585
Listing
NR=IC-I
0586
C
0587
C
0588
C
0589
C
0590 0591
B - A_SURF
Finished RETURN END
0592
C C
0593
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
0594
C
0595
SUBROUTINE
GEOREAD
0596
C
0597
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
0598
C
0599
C
0600
C
Georead
reads
0601
C
contains
the
0602
C
exposed
0603
C C
0605 0606
C
0607
C
0608
C
0609 0610
C C
0611
C
0612
C
answer
0613
C
gfile
0614
C
itf
0615
C
0616
C
note:
for
Variable
output
file
threat
and
and
their
variables
listing
0617
the
elements
0604
C
in
global
impact
contained
for
from
the
element angles
in
GEOMETRY
code.
as
as
data
the
for
common
well
each
threat
block
refer
This the
file
list
of
case.
to
main
definition
List = -
character
string
geometry
output
analysis
INCLUDE
type
representing
user
input
filename
contained
in
the
'COMMONAS.BLK'
0618 0653
CHARACTER
0654
CHARACTER*I
0655
CHARACTER*20
BUMDTTM
0656 0657
CHARACTER*80 CHARACTER*30
ANSWER, AA
0658
CHARACTER*40
BB
0659
CHARACTER*I2
CC
0660
C
IX(3)
INTEGER*4
0661
REAL*4 C
IF
0665
C
0666
C
0667
C
AREAMAXSF 'IBOTHS,
(IBOTHS.EQ.1)
Read
in
the
0668
IF
0669
ANSWER=ROOTFILE
0670
JOT
0671
WRITE
0672
FORM
0673
)
0674
10
0675 0676
-
INDEX (LENGTH, =
WRITE READ
30
FIRST
GOTO
GEOMETRY
output
FORMAT
' (/IX,
' (I2) ''GEOMETRY
0)
(ROOTFILE, '.'
GEOREAD',IBOTHS
set
the
default
to
station
?.gem
ROOTFILE='STATION.' '
' ) )//'GEM'
)
' )JOT+3 OUTPUT
'',$)' (6,FORM)ANSWER ( 5,30
IN
filename,
.EQ.
(1 : INDEX (ROOTFILE,
LINE
60
(INDEX(ROOTFILE,'.')
>
Form
ITF
WRITE(6,*)
0664
GFILE,
ITF
INTEGER*2
0662 0663
LENGTH*2
) GFILE
(A)
B-IO
FILENAME
(='',A'//LENGTH//',
''
Appendix
0677 0678
C C
0679
C
!!!!
THIS
OPEN
BOX.
0680
IF
OPEN
0682
INQUIRE( IBOTHS=0
0683 0684 C
0687
C
0688
C
!!!!
0689
END IF
0691 0693
IF
OF
ELSE ROOTFILE ENDIF
MAC
OPEN
=
IF (IBATCOM.EQ.
07O0
WRITE RETURN
0701 0702
END
0703
C
0704
C
0705
C
(13,
the
OPEN
GFile(I:INDEX(GFile,'
'))
i)
THEN
' (A) ')
GFILE
file
(UNIT=22,FILE=GFILE,
STATUS-'OLD',FORM-'UNFORMATTED',ERR-40
)
IBOTHS=0
0709 0710
C
0711
C
0712
C
GO
0713
40 50
0714 0715 0716
C
0717
C
0718
C
0719
C
0720
C
0721
C
0722 C
0724
C
TO
60
Error
control
WRITE FORMAT
( 6,50 ) ( /IX,'UNABLE
GO
i0
TO
TO
Read in the analysis number of elements WRITE(*,*) 60 60
0723 0725
READ
' )
number
of
threat
ITYPe,
IEnv,NT,NELM, IEnv,
NT,NELM,
Inclin Inclin,
IEnv,
NT,NELM,
Inclin'
WRITE(6,*)
IType,
IEnv,
NT,NELM,
Inclin
IF
ITYPEIN,
(ITYPEIN.EQ.3.AND.ITYPE.EQ.1) (ITYPEIN.EQ.3.AND.ITYPE.EQ.2) (IType.EQ.l) WRITE FORMAT
(/5X,'
IF
0731
WRITE ELSE
(6,64)
WRITE
(6,65)
(IEnv.EQ.1)
END
IBOTHS-1 IBOTHS-2
THEN
(6,62)
0730
Debris
cases,
G.READ',IBOTHS
IType, ITYPe,
0727 62
IN
ITYPEIN, ITYPEIN,
FILE
WRITE(6,*)'ITYPEIN,
0726 0728
(22)
OPEN
type,the
'IBOTHS
(22)
READ
IF IF
0734
' ) THEN
C
0708
0733
THEN
IF
Open
0706
0732
)
C
0699
0729
DIALOG
10
0696
0707
FINDER
IF
0695
0698
NORMAL
UNIT=22,NAME=GFILE)
( GFILE(I:4).EQ.' GFILE=ANSWER
0694
THE
SKIPPED
60
GOTO END
GIVE IS
THEN
(GFILE(I:I).EQ.'?') CALL DIRLIST
0692
WILL
METHOD
IF
0690
0697
MAC
Listing
( UNIT=22,FILE=*,STATUS='OLD',FORM='UNFORMATTED',ERR=40
GOTO END
0686
THE
DIRLIST
(GFILE(I:I).EQ.'?')
0681
0685
FOR
THE
B - ASURF
Analysis
THEN
IF
B-11
')
AREAMAX
and
the
Appendix
0735
WRITE
0737
63
(6, 63)
FORMAT (/5X, IF (IEnv.EQ.1)
0738 0739
WRITE
0740 0741
64
' Meteoroid THEN
Analysis')
(6, 64)
FORMAT
(/5X,
' JSC-20001&6000
Environment')
ELSE
0742
WRITE
0743
65
0745
END
0746
C
0747
C
0748
C
0749
C IF
0751 0752
(/5X,'
0753 0754
( NT.GT.ITH
WRITE 105
0756
END
the
number
of
elements
are
'NUMBER
( 6,105
)
( lX,'ARRAY
OF
THREATS
SIZE
MUST
IS BE
GREATER
INCREASED
THAN
ALLOWED',I9)
& CODE
RECOMPILED')
IF
C
0758
IF
0759 0760
110
0761 0762 0763
( NELM.GT.IELM
C
0765
C
0766
C
( 6,110 ) ( /IX,'NUMBER
WRITE STOP
( 6,105
OF
ELEMENTS
IS
GREATER
THAN
MAX
ALLOWED')
)
IF
Initialize
the
DO
I-1,NT
0767
) THEN
WRITE FORMAT
END
0764
150
arrays
0768
THREAT
(3, I) -0.0
0769
THREAT
(4, I) -0.0
0770
EXPOSED
(I) -0
0771
DO
J-1,NELM
140
0772
GEOMETRY
0773
ID (1, J) -0
0774
to
0.0
(J, I) -0.0
ID (2, J) -0
0775
POINT CONTINUE
140
0776 0777 0778
C
0779
C
0780
C
150
Read
in
the
DO
0782
READ CONTINUE
(22)
the
0783
175
175
0784
C
0785
C
Read
0786
C
array.
0787
C
in
0788
DO
0789
ElemLoc CONTINUE
180
(J, I) -0
CONTINUE
0781
180
Threat
data
I-1,NT (THREAT(J,I),J-1,4)
element
id,
and
property
I'i,100000 (I)
-
0
0791 0792
and
)NT
(/IX,
FORMAT STOP
threats
) THEN
( 6,100
FORMAT
0755
Memo')
that the number of the maximum allowed
WRITE I00
JSC-7/90
IF
Check than
0750
(6, 65)
FORMAT END IF
0744
0790
Listing
ELSE
0736
0757
B - A_SURF
DO
200
I-1,NELM
B-12
id
storing
them
in
the
ID
less
Appendix B -A_SURF
0793
READ
0794 0795 0796
200
(ID(J,I),J=I,6)
(ID (I, I) ) =
I
C
0797
C
0798
C
Read
0799
DO
0800 0801 0802
(22)
ElemLoc CONTINUE
Listing
C C
0803
250
in
the
element's
250
I=I,NELM
READ
(22)
surface
area
storing
it
in
the
AREA
array.
AREA(I,0)
WRITE
(6, *)
'AREA
WRITE CONTINUE
(6, *)
AREA
(I, 0) , I, NELM' (I, 0) , I, NELM
0804 0805
Read
in
the
geometry
data
for
the
exposed
elements
0806 0807
DO
500
I=I,NT
0808 Read
0809
in
the
threat
case
and
the
number
of
exposed
elements
0810 0811
READ
0812
(22)
IT,EXPOSED(I)
WRITE(6,*)
'IT, EXPOSED(I)
',IT,
EXPOSED(I)
0813 0814
Loop
thru
the
exposed
elements
0815 0816
DO
400
J=I,EXPOSED(I)
0817 0818
Read
0819
cosine
in
the of
element
the
number
(storing
in
the
POINT
angle
(storing
in
the
GEOMETRY
impact
array), array).
0820 READ
0821 0822 0823
C
0824
C
0825 0826
C
0827 0828
C C
0829
C
0830
C
0831
C
0832
C
0833
C
0834
C
400
CONTINUE
500
CONTINUE
!!!!!
575
(J, I) ,GEOMETRY
LIBSDATE_TIME(BUMDTTM) ONLY RECORDED FOR
CALL
TIME(BUMDTTM) gfile
WRITE FORMAT
to
summary
MAC
( 10,600
FORMAT
(J, I) ,Area
(J, I)
VERSION
file
(10, 575) BUMDTTM ('I', 'SHIELD',40X,
WRITE 600
POINT
CALL TIME
Write
0835 0836
(22)
A,/)
)GFILE
( iX, 'GEOMETRY
OUTPUT
FILE
-
',A
)
0837 To
0838
read
Header
from
.GEM
2-8-91
file
0839 0840
READ
0841
WRITE
0842
DO
( 22,END=630 610
0843
READ
0844
WRITE
0845
610
( i0, J
=
'(
) AA //A
)'
) AA
1,3
(22) ( i0,
CC,BB '(
2A
)'
) CC, BB
CONTINUE
0846 0847
WRITE
( i0, ' ( /A, I4) ')
' Threat
0848
WRITE
( 10, ' ( /A, I4) ')
' Environment
0849
WRITE
( i0, ' ( 5X, A, I5)
')
(I Debris, '
0850
B-13
(i Number
2 Old, of
Meteoroid) 2
New)
Threats
',IType ',IEnv ',NT
and
the
Appendix
0851
C
Write inclin to summary IF ( ITYPE.EQ.I.AND.IEnv.EQ.2
0852 0853 0854
520
FORMAT
file ) WRITE
( lX,'INCLINATION
(DEGREES)
(
READ DO
0858 0859
22
620
)
J-
READ IF
0860
AA
Rotation
1,3
( 22 (IXASC
+ IF
'
)
IXASC(J)
Maximum
(I0,
' (/IX,
A) ')
Axes
(IXASC(J).NE.0.AND.ROTANG(J).NE.0) (J) .EQ. 1)
IX(J)-'X'
0863 0864
IF (IXASC
(J) .EQ. 2)
IX(J)
IF (IXASC
(J) .EQ. 3)
IX(J)ffi'Z'
0865
WRITE
( 10,619
0866
WRITE
( 6,619
0867
619
FORMAT END IF
620
CONTINUE GO TO 650
0871
630
WRITE
0872
640
FORMAT
650
IF
0868 0869 0870
( 2X,
) )
'
( i0,640
)
( /2X,
'---
0875
C
0876
C
0877
C
0878
_'Y
(IBOTHS.EQ.I)
Axisffi',A3,';
Header
(J) Angleffi',F8.2
following
.Gem
RETURN
Close
file
CLOSE
( UNITf22,STATUSf'KEEP'
)
RETURN END C
0884
SUBROUTINE
0885
CHARACTER*80
0886
OPEN(UNITffiI7,FILE='DIRECTORY.LIST',STATUS='OLD',ERRffi30) REWIND 17
0887 5
0889
DIRLIST
0890 20
REWIND
0892
LINE
READ(17,10,ERRffi20) WRITE(6,11) GOTO 5
LINE
LINE
17
i0
CLOSE(UNIT=17,STATUS='KEEP') FORMAT (A80)
0894
ii
FORMAT
0895
30
RETURN
0893
0896
I)
IX(J),ROTANG(J)
C
0882
J.EQ.
'
C
0880
0891
Angles
rotations
C
0874
0888
and
THEN
IX (J) , ROTANG
No
3
(J) .NE. 0 .AND.
IF (IXASC
0883
)INCLIN )
AA
0862
0881
',F6.1
,ROTANG(J)
(J) .NE. 0 .AND.ROTANG
WRITE
0861
0879
( 10,520 -
C
0857
0873
Listing
C
0855 0856
B - A_SURF
( IX, A80)
END
B-14
file
)
---'
)
SD_SURF
Appendix
C.
User's
R_PLOT5
C-i
Manual
Source
Code
R_PLOT5
Listing
Listing from Language Systems FORTRAN (Version 3.0 Tue, Nov 19, 1991) Sat, Sep 12, 1992 11:50 AM Options OFF: A ANSI BKG=O CASE CCD CCX CRAY DYNE EXTENDED F77 I2 LINEFEED MC68040 MC68881 NOIMPLICIT OV R SANE SAVEALL SYM T72 TRACE U VAX W X Z Options ON: C L S 0001 0002 0003 0004 0005 0006 0007 0008 0009 0010 0011 0012 0013 0014 0015 0016 0017 0018 0019 0020 0021 0022 0023 0024 0025 0026 0027 0028 0029 0030 0031 0032 0033 0034 0035 0036 0037 0038 0039 0040 0041 0042 0043 0044 0045 0046 0047 0048
MC68020
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C R,_PLOT5 VER 1.6 8/23/92 C C C C MARTIN MARIEI-FA C C MANNED SPACE SYSTEMS C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C C R_PLOTS VER 1.6 reads the BUMPERII vl.3 - RESPONSE binary output C and converts it to formatted output. The data is put out at C 5 degree increments rather than the 15 degree increments originally C used by RPLOT. The output velocity increment is set to 0.5 km/s for C spreadsheet use. Commas are used to detimit the output to make it C more easily read by a spreadsheet program such as Microsoft EXCEL. C The output may then be utilized by SD_Surf for EXCEL to perform C probability analysis. C C The code requires the output file of the RESPONSE portion of the C BUMPERii code. This file contains the ballistic limit C (minimum diameter to penetrate) as a function of velocity C and obliquity. C C The RESREAD subroutine is taken directly from BUMPERII vl.Za - Shield C version 1.2a except for the COMMONRP.BLK instead of COHMONZ.BLK C BUMPER was developed under the NASA contract 'Integrated Wail Design C Guide and Penetration Control Plan' by M.A.Wright & A.R.Coronado. C C SD_Surf was developed under the NASA contract 'Structural Damage C Prediction and Analysis for Hypervelocity Impacts Study' under the C direction of N. E1fer. C C Version 1.6was modified to read BUMPERII version 1.3 file C C Include module COMMONRP variable list C C C aden = combined areal density of shields C conf = text description of wall configuration C diam = critical diameter , cm C idens - debris density, i- constant density, 2-size function C ienv = environment type, I- 35C 2_I&6_, Z- 7/90 memo C it = current threat case C itype = analysis type , 1- debris, Z-meteoroids C nb = number of angles in the response array C nc = number of wall configurations in the response array
C-I
R_PLOT5
0049 0050 0051 0052 0053 0054 0055 0056 _57 _58 _59 _60 _61 _62 _63 _64 _65 0066 0067 _68 _69 _70 0071 0072 _73 _74 _75 _76 _77 _78 _79 _80 _81 _82 _83 _84 _85 _86 _87 _88 _89 0090 0091 009Z 0093 0094 0095 _96 0097 0098 _99 0100 0_I 0_2
C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C
nv = number of velocities in the response array pid = the property id associated with all elements of the rootfile = file name rsfile = name of R_Plot5 output file for response surface units = english or metric vr = impact (relative) velocity , km/sec vinc = impact (relative) velocity increment , km/sec
adar
= areal density g/cm**2 = array containing the current pid number to process response = array containing the values of the critical diameter a function of impact angle and velocity. (vr,beta,pid) standm = shield stand-off, cm shden = shield density, g/cc shthkm = shield thickness, cm vwden = vessel wall density, g/cc vwthkm = vessel wall thickness, cm pids
Main
Program
Variable
List
Scalers
first = logical if only one PID encountered in requested answer = user input ob = Obliquity for the current threat/element.
FIRST
CHARACTER*80
INCLUDE
ANSWER
'COMHONRP.BLK'
Initialize
the Velocity
increment
and number
of velocities
VINCI0.25 NV168 C C C
ranges
Arrays
LOGICAL
C C C
Listing
Initialize
the Obliquity
increment
BINCI5.0 NB=I9
IBATCON
= 0
C-2
and number
of angles.
ranges
as
R_PLOT5
0103 0104 0105 0106 0107 0108 0109 0110 0111 0112 0113 0114 0115 0116 0117 0118 0119 0120 0121 0122 0123 0124 0125 0126 0127 0128 0129 0130 0131 0132 0133 0134 0135 0136 0137 0138 0139 0140 0141 0142 0143 0144 0145 0146 0147 0148 0149 0150 0151 0152 0153 0154 0155 0156
C C
Write
header
to
screen
and read
Listing
in
orbital
parameters
CALL RPINPUT C C C
SPECIFY
DEBRIS
ANALYSIS
FOR RESREAD
ITYPE = 1 C C C C
Read in the RESPONSE output file. the RESREAD subroutine in BUMPER.
This
is identical
to
CALL RESREAD C C C C C
Determine
the RESPONSE
If number
of cases (NC) is only one then
IF (NC.EQ.1) THEN PIDS(1)=I WRITE ( 6,20 ) 20 FORMAT (/IX,'The C C C C
For multiple Write number
PIDs to process. proceed.
one case in the RESPONSE
file will
be used'
)
PIDs select one or all. If only one, NC is set to 1. of PIDs and first PID in A_SURF to screen.
ELSE 25
WRITE ( 6,25 ) NC FORMAT (/lX,'The Number
30
WRITE ( FORMAT
of
PIDs
in
the
RESPONSE file
is
',I4)
C
2
6,30 ) (//lX,'Enter /1X,'Enter
the PID number to use o specific to use all PIDs.')
C 35
40
45
READ ( 5,35 ) ANSWER FORMAT (A) IF C ANSWERCI:I).EQ.' DO 40 II=I,NC PIDS(I1) - I1 CONTINUE ELSE READ ( ANSWER(1:80),45 FORMAT ( BN,12 ) NC=I ENDIF
' .OR. ANSWER.EQ.'A')
THEN
)PIDS(1)
C ENDIF C C C C C
Print out the flux x area x time surface delimited
format
in a comma
to be read by spreadsheets.
CALL RP5TEXT C C
C-3
PID.',
R_PLOT5
0157 0158 0159 0160 0161 0162 0163 0164 0165 0166 0167 0168 0169
C C C C
Close
summary
Listing
file
CLOSE ( UNIT-Ie,STATUS-'KEEP' ) WRITE( 6,64_03 ) RSFILE FORNATC /' The R_PLOT5 file is complete.'/ ' filename: ',A ) C C C
Finished
Symbol Table The following
END
for: HACll_R_plotS_165main symbols
were
defined
but NOT referenced: Data
Type
Offset
Size
RefCnt
Symbol Name Other Information
Class
ADAR [Cmn ] Dims ADEN [Cmn ] BATCOM
Array
REAL*4
144
192436
0
Scalar
REAL*4
4
192432
0
Scalar
REAL*4
4
192624
0
Scalar
CHARACTER
12
192632
0
Scalar
REAL*4
4
4
0
Scalar Scalar
LOGICAL*4 INTEGER*2
4 2
192628
0 0
Scalar
INTEGER*2
2
192430
0
Scalar
INTEGER*2
2
192428
0
1:36
[_n ] CONF [Cmn ] DIAH [Cmn ] FIRST IBOTHS [Cmn ] IDENS [Cmn ] IENV [Cmn ] IT
Scalar
INTEGER*2
2
8
0
[Cmn ] ITYPEIN
Scalar
INTEGER*2
2
192630
0
ECmn ] NT
Scalar
INTEGER*4
4
16
0
[Cmn ] PID
Scalar
INTEGER*2
2
114
0
[Cmn ] RANGE
Scalar
REAL*4
4
30
0
Array
REAL*4
191520
188
0
Scalar
CHARACTER
40
192584
0
[Cmn ] RESPONSE [Cmn ] Dims ROOTFILE
1:70
1:19
1:36
[Cmn ]
C-4
R_PLOT5
SHDEN [Cmn ] SHTHK [Cmn ] SHTHK_ [Cmn ] STANDt4 [Cmn ] UNITS [Cmn ] VR CCmn ] VWDEN [Cmn ] VW/HK VWTHKH [Cmn ]
Dims -
1:36
Dims
1:36
Dims -
REAL*4
144
191852
0
Scatar
REAL*4
4
192588
0
Array
REAL*4
144
191788
0
Array
REAL*4
144
192284
0
Scatar
CHARACTER
12
192644
0
Scalar
REAL*4
4
2Z
0
Array
REAL*4
144
192140
e
Scalar
REAL*4 REAL*4
4 144
191996
g g
1:36 Array
Dims -
1:36
Alphabetic List: Symbol Name Other Information
ADAR [Cmn ] Dims ADEN [Cmn ] ANSWER BATCOM [Cmn ] BINC [Cmn ] CONF [Cmn ] DIAM [Cmn ] FIRST I1 IBATCOM IBOTHS [Cmn ] IDENS [Cmn ] IENV [Cmn ] ISHLDS 36 IT [Cmn ] ITYPE [Cmn ] ITYPEIN [Cmn ]
Array 1:36
Dims -
-
Listing
Class
Data
Type
Size
Offset
RefCnt
Array
REAL*4
144
192436
0
Scalar
REAL*4
4
192432
0
Scatar Scalar
CHARACTER REAL*4
88 4
38 192624
4 0
Scalar
REAL*4
4
0
i
Scalar
CHARACTER
Scalar
REAL*4
4
4
0
Scalar Scalar Scalar Scalar
LC_ICAL*4 INTEGER*4 INTEGER*2 INTEGER*2
4 4 2 2
12e 118 192628
0 4 I 0
Scalar
INTEGER*2
2
192430
0
Scalar
INTEGER*2
2
192428
0
Parameter
INTEGER*4
Scalar
INTEGER*2
2
8
0
Scalar
INTEGER*2
2
10
I
Scalar
INTEGER*2
2
1:36
C-5
12
192632
0
8=
192630
0
R_.PLOT5
NB
Scalar
[Cmn l NC
Scalar
[O,,n ] NT
]
Scalar
EO,,n]
Scalar
NV PID
Scalar
[Cmn ] PIDS [Cmn ] Dims RANGE [Cmn ] RESPONSE [Cmn ] Dims ROOTFILE FCmn ] RSFILE [Cmn ] SHDEN
-
1:36
Scalar Array
-
1:70 1:19 1:36 Scalar Scalar
[Cmn ] Dims SHTHK
-
1:36
Array Scalar
[Cmn ] SHTHKM [Cmn ] Dims STANDH [Cmn ] Dims UNITS
1:36 1:36
Array Array Scalar
[Cmn ] VINC
Scalar
[Cmn ] VWDEN [Cmn ] Dims VWTHK VWTHKM
-
[Cmn ] Dims _$CMNBASES
-
1:36
1:36
INTEGER*2
2
12
1
INTEGER*2
2
14
4
INTEGER*4
4
16
0
INTEGER*2
2
20
I
INTEGER*2
2
114
0
INTEGER*2
72
116
4
REAL*4
4
30
0
REAL*4
191520
188
0
CHARACTER
40
CHARACTER
80
34
1
REAL*4
144
191852
0
REAL*4
4
192580
0
REAL*4
144
191708
0
REAL*4
144
192284
0
12
192644
0
REAL*4
4
26
1
REAL*4
4
22
0
REAL*4
144
192140
0
REAL*4 REAL*4
4 144
191996
0 0
124
1
CHARACTER
Scalar
[Cmn ] VR
Stack
Array
Usting
Array Scalar Array
192584
Cmn Hndls
0
Frame Information:
Temporaries List: Symbol Name Other Information
_$TEHPI {Work Area} _$TLBO _$TILAsgnO _$TLB1 _$TLB2
Class
Data Type
Scalar
INTEGER*4
Scalar Scalar Scalar Scalar
LOGICAL*I INTEGER*4 LOGICAL*I LOGICAL*I
C-6
Size
Offset
4 24 1 4 1 1
0 4 28 32 36 37
RefCnt
2 4 2 2 2
R_PLOT5
Variable List: Symbol Name Other Information
Class
ANSWER IBATCOM 11 _$CMNBASES
Scalar Scalar Scalar Cmn Hndls
Variables
in
Blank
Common:
Listing
Data
Type
Size
CHARACTER INTEGER*2 INTEGER*4
Offset
80 2 4
RefCnt
38 118 120 124
4 1 4 1
Size:192656 Size
Offset
RefCnt
Symbol Name Other Information
Crass
BINC [Cmn ] DIAM [Cmn ] IT
Scalar
REAL*4
4
0
1
Scalar
REAL*4
4
4
0
Scalar
INTEGER*2
2
8
0
Scatar
INTEGER*2
2
10
1
Scalar
INTEGER*2
2
12
1
Scalar
INTEGER*2
2
14
4
Scalar
INTEGER*4
4
16
0
Scalar
INTEGER*2
2
20
1
Scalar
REAL*4
4
22
0
Scalar
REAL*4
4
26
1
Scalar
REAL*4
4
30
0
Scalar
CHARACTER
80
34
1
Scalar
INTEGER*2
2
114
0
Array
INTEGER*2
72
116
4
Array
REAL*4
191520
188
0
Array
REAL*4
144
191708
0
Array
REAL*4
144
191852
0
Array
REAL*4
144
191996
0
[Cmn ] ITYPE [Cmn ] NB [Cmn ] NC [Cmn ] NT [Cmn ] NV [Cmn ] VR [Cmn ] VINC [Cmn ] RANGE [Cmn ] RSFILE [Cmn ] PID [Cmn ] PIDS [Cmn ] Dims RESPONSE [Cmn ] Dims SHTHKM [Cmn ] Dims SHDEN [Cmn ] Dims VWI'HKM [Cmn ] Dims
-
1:36
-
1:70 1:19
-
1:36
-
1:36
-
1:36
Data
Type
1:36
C-7
R_PLOT5
VWDEN
Array
REAL*4
144
19Z14_
0
Array
REAL*4
144
192Z84
0
Scalar
INTEGER*2
Z
192428
0
Scalar
INTEGER*Z
2
192430
0
Scalar
REAL*4
4
192432
0
Array
REAL*4
144
192436
0
Scalar
REAL*4
4
192580
0
Scalar
CHARACTER
40
192584
0
Scalar
REAL*4
4
19Z624
0
Scalar
INTEGER*2
2
192628
0
Scalar
INTEGER*2
2
192630
0
Scalar
CHARACTER
12
192632
0
Scalar
CHARACTER
12
19Z644
0
1:36
[Cmn ] Dims STANDM [Cmn ] Dims IENV [Cmn ] IDENS [Cmn ] ADEN [Cmn ] ADAR [Cmn ] Dims SHTHK [Cmn ] ROOTFILE [Cmn ] BATC(_4 [Cmn ] IBOTHS [Cmn ] ITYPEIN [Cmn ] CONF [Cmn ] UNITS [Cmn ]
1:36
1:36
Local Stackframe size: Local Symbols: 60
0170 0171 017Z 0173 0174 0175 0176 0177 0178 0179 0180 0181 0182 0183 0184 0185 0186 0187 0188 018g 0190 0191
Usting
IZ8
C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C SUBROUTINE RP5TEXT C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C C
Write
RESPONSE array
C C C C C
Temporary
Variable
RPLINE
for
5 degree
and 0.5 km/sec
increments
list array
= Contains diameters to penetrate at the obliquity increments Identical to Response at a particular velocity and PID.
INCLUDE 'CO_NRP.BLK' REAL*4 RPLINE(19) C C C DO 700 I=I,NC
C-8
R_PLOT5
0192 0193 0194 0195 0196 0197 0198 0199 0200 0201 0202 0203 0204 0205 0206 0207 0208 0209
Listing
WRITE (10,625) PIDS(I) FORMAT (1H1,/lX,'RESPONSE PID=', ' , ',I9) WRITE ( 10,630 ) PIDS(1),(3,J=0,90,5) FORMAT (I6.2,19(',',I12.2) )
625 630 C
DO 690 K=2,NV,2 DO 640 3=1,19 RPLINE(J)=RESPONSE(K,3,PIDS(I)) CONTINUE VR=K*VINC WRITE(10,650) VR,(RPLINE(J),3=l,19) FORMAT (F6.2,19(',',FI2.8) ) CONTINUE CONTINUE
640
650 690 700 C
RETURN C
Symbol Table
END
for: RPSTEXT
The following
symbols
were defined
but NOT referenced: Size
RefCnt
Class
Data Type
Array
REAL*4
144
192436
0
Scalar
REAL*4
4
192432
0
Scalar
REAL*4
4
192624
0
Scalar
REAL*4
4
0
0
Scalar
CHARACTER
12
192632
0
Scalar
REAL*4
4
[Cmn ] IBATCOM IBOTHS
Scalar Scalar
INTEGER*2 INTEGER*2
2 2
192628
0 0
[Cmn ] IDENS
Scalar
INTEGER*2
2
192430
0
[Cmn ] IENV
Scalar
INTEGER*2
2
192428
0
[Cmn ] IT
Scalar
INTEGER*Z
2
8
0
[Cmn ] ITYPE
Scalar
INTEGER*2
2
10
0
[Cmn ] ITYPEIN
Scalar
INTEGER*Z
Z
Symbol Name Other Information
Offset
...........................
......................
ADAR [Cmn ] Dims ADEN [Cmn ] BATCOM [Cmn ] BINC [Cmn ] CONF [Cmn ] DIAM
1:36
[Cmn ]
C-9
4
192630
0
0
R_PLOT5
NB [Cmn 3 NT [Cmn ] PID [¢mn ] RANGE [Cmn 1 ROOTFILE [Cmn ] RSFILE [Cmn ] SHDEN [Cmn ] Dims SHTHK [Cmn ] SHTHKM [Cmn ] Dims STANDM [Cmn ] Dims UNITS [Cmn ] VWDEN [Cmn ] Dims VWTHK VWTHKH [Cmn 3 Dims
-
1:36
-
1:36
[Cmn ] CONF [Cmn ] DIAM [Cmn ] I IBATCOH IBOTHS [Cmn ] IDENS [Cmn ] IENV [Cmn 3
INTEGER*2
2
iZ
0
Scalar
INTEGER*4
4
16
0
Scalar
INTEGER*2
Z
114
0
Scalar
REAL*4
4
30
0
Scalar
CHARACTER
40
19Z584
0
Scalar
CHARACTER
80
34
0
Array
REAL*4
144
191852
0
Scalar
REAL*4
4
192580
0
Array
REAL*4
144
191708
0
Array
REAL*4
144
192284
0
Scalar
CHARACTER
12
19Z644
0
Array
REAL*4
144
19Z140
0
Scalar
REAL*4 REAL*4
4 144
191996
0 0
1:36 Array
-
1:36
Alphabetic List: Symbol Name Other Information
ADAR [Cmn ] Dims ADEN [Cmn 3 BATCON [Cmn ] BINC
Scalar
1:36
-
-
-
Listing
Class
Data
Type
Offset
Size
RefCnt
Array
REAL*4
144
19Z436
0
Scalar
REAL*4
4
19Z43Z
0
Scalar
REAL*4
4
19Z6Z4
0
Scalar
REAL*4
4
Scalar
CHARACTER
Scalar
REAL*4
4
4
0
Scalar Scalar Scalar
INTEGER*4 INTEGER*Z INTEGER*Z
4 Z Z
144 1926Z8
5 0 0
Scalar
INTEGER*Z
2
19Z43g
0
Scalar
INTEGER*Z
Z
19Z428
0
1:36
C-10
12
g 19263Z
0 0
R_PLOT5
ISHLDS 36 IT [Cmn ] ITYPE [Cmn ] ITYPEIN [Cmn ] 3 K NB [Cmn ] NC [Cmn ] NT [Cmn ] NV [Cmn ] PID [Cmn ] PIDS [Cmn ] Dims RANGE [Cmn ] RESPONSE [Cmn ] Dims ROOTFI L E [Cmn ] RPLINE Dims RSFILE [Cmn ] SHDEN [Cmn ] SHTHK [Cmn ] SHTHKH [Cmn ] STANDM [Cmn ] UNITS [Cmn ] VINC [Cmn ] VR [Cmn ] VWDEN [Cmn ] VWTHK VWTHKM
-
-
Parameter
INTEGER*4
Scalar
INTEGER*2
2
8
0
Scalar
INTEGER*2
2
i@
0
Scalar
INTEGER*2
2
19263@
0
Scalar Scalar Scalar
INTEGER*4 INTEGER*4 INTEGER*2
4 4 2
148 152 12
10 4 0
Scalar
INTEGER*2
2
14
1
Scalar
INTEGER*4
4
16
0
Scalar
INTEGER*2
2
28
1
Scalar
INTEGER*2
2
114
0
Array
INTEGER*2
72
116
6
Scalar
REAL*4
4
Array
REAL*4
19152@
188
2
Scalar
CHARACTER
48
192584
0
Array
REAL*4
76
68
4
Scalar
CHARACTER
8@
34
0
Array
REAL*4
144
191852
0
Scalar
REAL*4
4
19258@
0
Array
REAL*4
144
191788
0
Array
REAL*4
144
192284
0
Scalar
CHARACTER
12
192644
0
Scalar
REAL*4
4
26
1
Scalar
REAL*4
4
22
Z
Array
REAL*4
144
192148
0
Scalar Array
REAL*4 REAL*4
4 144
191996
0 0
l_
1:36
1:70 1:19
3@
0
1:36
1:19
Dims
-
1:36
Dims
-
1:36
Dims
-
1:36
Dims
-
[Cmn ] Dims _$CMNBASES Stack
Listing
1:36
1:36 Cmn Hndls
Frame Information:
C-11
156
1
R_PLOT5 Listing
Temporaries List: Symbol Name .Other Information
Class
_$TEMPD4 _$TEMP3 _$TEMPZ _$TEMPD3 _$TEMP1 {Work Area} _$TILO _$TILAsgne _STILl _$TEe Variabte
Data
Type
Scalar Scalar Scalar Scalar Scalar
INTEGER*4 INTEGER*4 INTEGER*4 INTEGER*4 INTEGER*4
Scalar Scalar Scalar Scalar
INTEGER*4 INTEGER*4 INTEGER*4 EXTENDED*12
4 4 4 4 4 Z4 4 4 4 12
RefCnt
0 4 8 1Z 16 ze 44 48 5Z 56
2 Z Z Z Z 8 2 Z 2
List:
Symbol Name Other Information
Crass
RPLINE Dims 1:19 I J K _$CMNBASES
Array
REAL*4
Scalar Scalar Scalar Cmn Hndls
INTEGER*4 INTEGER*4 INTEGER*4
Variables
Offset
Size
in
Blank
Data
Type
Offset
Size
RefCnt
76
68
4
4 4 4
144 148 152 IS6
5 10 4 1
Size:192656
Common:
Data
Offset
RefCnt
Class
BINC [Cmn ] DIAM [Cmn ] IT [Cmn 1 ITYPE [Cmn l NB CCmn ] NC
Scalar
REAL*4
4
0
0
Scalar
REAL*4
4
4
0
Scalar
INTEGER*2
2
8
0
Scatar
INTEGER*Z
Z
10
0
Scatar
INTEGER*Z
2
12
0
Scatar
INTEGER*2
2
14
1
[Cmn ] NT
Scatar
INTEGER*4
4
16
0
[Cmn ] NV
Scatar
INTEGER*Z
2
20
1
[Cmn ] VR [Cmn ]
Scatar
REAL*4
4
2Z
Z
C-12
Type
Size
Symbol Name Other Information
R_PLOT5 Listing VlNC [Cmn ] RANGE [Cmn ] RSFILE [Cmn ] PID [Cmn ] PIDS -
1:36
-
I:70 1:19 1:36
-
1:36
-
1:36
-
1:36
-
1:36
[Cmn ] Dims IENV [Cmn ] IDENS [Cmn ] ADEN [Cmn ] ADAR [Cmn ] Dims SHTHK [Cmn ] ROOTF ILE [Cmn ] BATCOM [Cmn ] IBOTHS [Cmn ] ITYPEIN [Cmn ] CONF [Cmn ] UNITS [Cmn ]
1:36
[Cmn ] Dims RESPONSE [Cmn ] Dims SHTHKM [Cmn ] Dims SHDEN [Cmn ] Dims VWTHKM [Cmn ] Dims VWDEN [Cmn ] Dims STANDM
0210 0211 0212 0213
REAL*4
4
26
i
Scalar
REAL*4
4
3@
0
Scalar
CHARACTER
80
34
0
Scalar
INTEGER*2
2
114
0
Array
INTEGER*2
72
i16
6
Array
REAL*4
191520
188
2
Array
REAL*4
144
1917@8
0
Array
REAL*4
144
191852
0
Array
REAL*4
144
191996
0
Array
REAL*4
144
192140
0
Array
REAL*4
144
192284
0
Scalar
INTEGER*2
2
192428
0
Scalar
INTEGER*2
2
192430
0
Scalar
REAL*4
4
192432
0
Array
REAL*4
144
192436
0
Scalar
REAL*4
4
19258e
0
Scalar
CHARACTER
40
192584
0
Scalar
REAL*4
4
192624
0
Scalar
INTEGER*2
2
192628
0
Scalar
INTEGER*2
2
192630
0
Scalar
CHARACTER
12
192632
0
Scalar
CHARACTER
12
192644
0
1:36
Local Stackframe Local Symbols:
Scalar
size:
160
67
C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C SUBROUTINE RESREAD
C-13
R_PLOT5 Listing 0214 0215 0216 0217 0218 0219 0220 0221 O222 0223 0224 0225 8226 0227 0228 0229 0230 0231 0232 0233 0234 0235 0236 0237 0238 0239 024Q 0241 0242 0243 0244 0245 0246 0247 0248 0249 0250 0251 0252 0253 0254 0255 0256 0257 0258 0259 0260 0261 0262 0263 0264 0265 0266 0267
C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C Resreod reads in the output from the RESPONSE code. This output C consists of the critical diameter data as a function of property C id, impact angle, and impact velocity. C C C note: for variables contained in the common block refer to the main C listing for definition. C C C Variable list C C answer = character string represnting user input C itf 1 analysis type for rfile C rfile = response output filenome C ienvr = environment for response file C CHARACTER LENGTH*2 CHARACTER*80 ANSWER,RFILE,Form CHARACTER*46 A46 CHARACTER*I2 BI2A, BIZB CHARACTER*8 C8A, C8B CHARACTER*2 D2 REAL*4 STND C INTEGER*4 Ill:, ITA, IC, ICT, ICB, IPF C l!! CHANGES TO BE COMPATIBLE WITH BUMPERII verl.3 INTEGER*2 ITF, ITA, IC, ICT, ICB, IPF, IPFUNC3 INTEGER*2 IENVR C INCLUDE 'COMMONRP.BLK' IF (IBOTHS.EQ.2) GOTO 6@ C C Read in the RESPONSE output filename , set default to resp.dat C IF (INDEX(ROOTFILE,' ') .EQ. O)ROOTFILE='STATION.' ANSWER=ROOTFILE(I:INDEX(ROOTFILE,' '))//'RSP' JOT = INDEX( ROOTFILE, ' ' ) WRITE (LENGTH, '(I2)' )30T+3 FORM='(/lX,"RESPONSE OUTPUT FILENAME (=",A'//LENGTH// |
,
If)
>
W|
,$)'
10
WRITE ( 6,FORM )ANSWER READ ( 5,30 ) RFILE 30 FORMAT (A) C C C
!l!!
THIS OPEN FOR THE MAC WILL GIVE THE NORMAL BOX. THE DIRLIST METHOD IS SKIPPED
FINDER
DIALOG
IF (RFILE(I:I).EQ.'?') THEN OPEN ( UNIT=23,FILE=*,STATUS='OLD',FORM='UNFORMA1-FED',ERR=4Q INQUIRE(UNIT=Z3,NAME=RFILE) GOTO 6@ END IF
C-14
)
R_PLOT5
0268 0269 0270 0271 0272 0273 0274 0275 0276 0277 0278 0279 0280 0281 0282 0283 0284 0285 0286 0287 0288 0289 0290 0291 0292 0293 0294 0295 0296 0297 0298 0299 0300 0301 0302 0303 0304 0305 0306 0307 0308 0309 0310 0311 0312 0313 0314 0315 0316 0317 0318 0319 0320 0321
C C C
!!!!
Listing
END OF MAC OPEN IF
(RFILE(I:I).EQ.'?') CALL DIRLIST GOTO 10 END IF IF ( RFILECI:4).EQ.' RFILE-ANSWER ELSE
THEN
' ) THEN
ROOTFILE = RFileCI:INDEX(RFile,' ENDIF
• '))
C IF(IBATCOM.EQ.1) WRITE(13,'(A)') RETURN END IF C C C
THEN RFILE
Open the file OPEN
( UNIT=23,FILE-RFILE,STATUS='OLD',FORM='UNFORMA1-FED',ERR=40
)
C GO TO 60 C C C
Error
controI
on open
40 WRITE ( 6,50 ) 50 FORMAT ( /lX,'UNABLE GO TO 10 C C C C C C C
C C C
Read
in the analysis
TO OPEN
FILE'
)
type and the number
of property
cases•
60 READ (23) ITYPEIN,ITF,IDens,NC !!1 NO ERROR CHECKING ON IENVR 60 READ (23) ITYPEIN,ITF,IENVR,IDens,NC WRITE(6,*)'ITYPEIN,ITF,IDens,NC' WRITE(6,*) ITYPEIN,ITF,IDens,NC IF (ITYPEIN.EQ.3.AND.IBOTHS.EQ.1) ITYPE=I IF (IDens.EQ.1) THEN WRITE (6,63) 63 FORMAT (/5X,' Constant density threat') ELSE IF (IDens.EQ.2) THEN WRITE (6,64) 64 FORMAT (/5X,' Variable density threat') END IF Check
that the response
file is the correct
analysis
type
IF ( ITF.NE.ITYPE ) THEN IF ( ITYPE.EQ.1 ) THEN WRITE ( 6,70 ) 70 FORMAT ( /IX,'DEBRIS ANALYSIS SPECIFIED IN GEOMETRY 1 'BUT RESPONSE FILE IS FOR METEOROIDS ')
C-15
FILE
',
R_PLOT5
0322 0323 0324 0325 0326 0327 0328 0329 0330 0331 033Z 0333 0334 0335 0336 0337 0338 0339 0340 0341 034Z 0343 0344 0345 0346 0347 0348 0349 0350 0351 0352 0353 0354 0355 0356 0357 0358 0359 0360 0361 0362 0363 0364 0365 0366 0367 0368 0369 0370 0371 0372 0373 0374 0375
Usting
ELSE WRITE C 6,80 ) FORMAT (/1X,'METEOROID ANALYSIS SPECIFIED IN GEOMETRY FILE', ' BUT RESPONSE FILE IS FOR DEBRIS' ) END IF
80 1 C
WRITE C 6,9e ) FORMAT ( /1X,'DO YOU WISH TO CONTINUE WITH GEOMETRY OPTION '(-NO) > ',$) READ ( 5,30 ) ANSWER
90 1 C
IF
C ANSWERCI:I).EQ.'Y' GO TO 10 ELSE STOP END IF
.OR.
ANSWERCI:I).EQ.'y'
C END IF C C C C C C C C C C C C C
Read in
the
impact
angle
information
READ CZ3) NB,BINC WRITE(6,*) 'NB,BINC' WRITEC6,*) NB,BINC Read in
the impact
velocity
READ CZ3) NV,VINC WRITEC6,*) 'IMPACT WRITE(6,*) NV,VINC Initiatize
RESPONSE to
information
VELOCITY,
0.0
DO 2_ I-I,NC DO 150 3-I,NB DO 1_ K-I,NV RESPONSE ( K,3,I 100 CONTINUE 150 CONTINUE 200 CONTINUE C C C C C C C C C C C
Read
in-the
critical
VEL INCR.'
) -0.
diameter
data
WRITE(6,*) 'NC,NB,NV' WRITE(6,*) NC,NB,NV WRITE(6,*) 'RESPONSECK,3,I)' Loop thru the property id's DO 400 I=I,NC Loop thru the impact DO 300 3=I,NB
angles
Loop thru the impact DO 250 K=I,NV
veracities
C-16
)
THEN
'
R_PLOT5
0376 0377 0378 0379 0380 0381 0382 0383 0384 0385 0386 0387 0388 0389 0390 0391 0392 0393 0394 0395 0396 0397 0398 0399 0400 0401 0402 0403 0404 0405 0406 0407 0408 0409 0410 0411 0412 0413 0414 0415 0416 0417 0418 0419 0420 0421 0422 0423 0424 0425 0426 0427 0428 0429
C 250 300 400 C C C C
i!!!
411 C
Store the READ (23) WRITE(6,*) CONTINUE CONTINUE CONTINUE
Listing
critical diameter RESPONSE(K,3,I) RESPONSE(K,J,I)
in
response
THE NEXT LINE WAS COMMENTED OUT FOR R_PLOT5 IF (INPUTCD.EQ.2) CALL SETDIAMS READ ( 23,END=44@,ERR=440 ) A46 WRITE(6,*) 'A46' WRITE ( 6,'( //lX,A)' ) A46 WRITE ( 10,'( //IX,A)' ) A46 READ ( 23 ) CSA,ITA,C8B,ICB,UNITS WRITE ( 10,'(A,I4)') ' Threat (i Debris, 2 Meteoroid) WRITE ( 10,'(A,I4)') ' Density (1 Constant, 2 Function) WRITE ( 10,'(A,I4)') ' Number of PID Cases WRITE ( 10,'(2A)') ' Units ',UNITS WRITE(6,') 'C8A,C8B' WRITE(6,*) C8A,C8B WRITE ( 6,'(A,I4)') ' Threat (1 Debris, 2 Meteoroid) WRITE ( 6,'(A,I4)') ' Density (I Constant, 2 Function) WRITE ( 6,'(A,I4)') ' Number of PID Cases WRITE ( 6,'(2A)') ' Units ',UNITS DO 420 I-I,ICB READ ( 23 ) ICT,D2,BIZA,BI2B,IPF,IPFUNC3 WRITE (10,411) I FORMAT( /IX,'PID NUMBER ',14 ) IF (ICT.EQ.2) THEN IF ( IPF.EQ.I ) THEN WRITE (18,485) ELSE IF ( IPF.EQ.2 ) THEN WRITE (10,486) ELSE IF ( IPF.EQ.3 ) THEN WRITE (10,487) ELSE IF ( IPF.EQ.4 ) THEN WRITE (10,488) ELSE IF ( IPF.EQ.5 ) THEN WRITE (10,484) ELSE IF ( IPF.EQ.6 ) THEN WRITE (10,489) ELSE IF ( IPF.EQ.7 ) THEN WRITE (10,490) ELSE IF ( IPF.EQ.8 ) THEN WRITE (10,491) ELSE IF ( IPF.EQ.9 ) THEN WRITE (10,492) ELSE IF ( IPF.EQ.10 ) THEN WRITE (10,493) ELSE IF ( IPF.EQ.11 ) THEN WRITE (10,494)
C-17
',ITA ',IDens ',ICB
',ITA ',IDens ',ICB
R_PLOT5
0430 0431 0432 0433 0434 0435 0436 0437 0438 0439 0440 0441 0442 0443 0444 0445 0446 0447 0448 0449 0450 0451 045Z 0453 0454 0455 0456 0457 0458 0459 0460 0461 0462 0463 0464 0465 0466 0467 0468 0469 0470 0471 0472 0473 0474 0475 0476 0477 0478 0479 0480 0481 048Z 0483
485 486 487 488 484 489 490 491 492 493 494 495 496 497 C C
C C C
C C C C
ELSE IF ( IPF.EQ.12 ) THEN WRITE (10,495) ELSE IF ( IPF.EQ.13 ) THEN WRITE (10,496) ELSE IF ( IPF.EQ.14 ) THEN WRITE (10,497) END IF END IF FORMAT ( /IX,'ORIGINAL PENETRATION FUNCTION') FORMAT ( /IX,'PEN4 PENETRATION FUNCTION') FORMAT ( /IX,'REGRESSION PENETRATION FUNCTION') FORMAT ( /IX,'COUR-PALAIS PENETRATION FUNCTION') FORMAT ( /IX, 'BOEING INTERP PENETRATION FUNCTION') FORMAT ( /IX,'DEVELOPHENTAL6, USER INPUT') FORMAT ( /IX,'DEVELOPHENTALT, USER INPUT') FORMAT ( /IX,'DEVELOPMENTALS, USER INPUT') FORMAT ( /IX,'DEVELOPMENTAL9, USER INPUT') FORMAT ( /IX,'DEVELOPMENTALIO, USER INPUT') FORMAT ( /IX,'DEVELOPMENTALII, USER INPUT') FORMAT ( /IX,'DEVELOPHENTALI2, USER INPUT') FORMAT ( /IX,'DEVELOPMENTALI3, USER INPUT') FORMAT ( /IX,'DEVELOPHENTALI4, USER INPUT') WRITE ( 10,'( /A )' ) ' Configuration Shield Wall' WRITE ( 6,* ) 'ICT,DZ,BIZA,BIZB' WRITE ( 6,* ) ICT,DZ,BI2A,BI2B IF (ICT.EQ.I) CONF = 'Single Plate' IF (ICT.EQ.2) CONF = 'Doubte Plate' IF (ICT.EQ.3) CONF = 'Muttiwatt' WRITE ( 10,'( IX,A,4X,ZA )' ) CONF,BIZA,BIZB WRITE ( 6,'( IX,A,4X,ZA )' ) CONF,BIZA,BI2B READ ( 23 ) ShThk,VWThk,STND,ShDen(1),VWDen(1),ADEN WRITE ( 6,* ) 'ShThk,VWThk, STND, ShDen(1) ,VWDen(1),ADEN,'I' WRITE (6,*) ShThk,VWThk,STND,ShDen(1),VWDen(1),ADEN,I IF (ICT.EQ.3) THEN WRITE ( 10,'(A,A,F8.4)') ' Combined Areal Density', ' of All Shietds = ',ADEN + WRITE ( 10,'(A,F8.4)') ' Total Standoff = ' -.I,STND Combined Areal Density', WRITE ( 6,'(A,A,FS.4)') ' + ' of AII Shields = ',ADEN Total Standoff = ' WRITE ( 6,'(A,F8.4)') ' +
C +
410 C
+
C
Usting
, STND GOTO 410 END IF WRITE ( 6,'(A,F8.4)') ' Shield Thickness = IF (SHTHK.NE.O.0) WRITE ( 10,'(A,F8.4)') ' Shietd Thickness WRITE ( 10,'(A,F8.4)') ' Vesset Wall Thickness = WRITE ( 6,'(A,F8.4)') ' Vessel Wall Thickness = IF (ICT.NE.3) THEN IF CSHTHK. NE.0.0.AND. STND. NE .0.0) WRITE ( 10, '(A,F8.4)') ' Standoff WRITE ( 6,'(A,F8.4)') ' Standoff = END IF
C-18
',ShThk = ',ShThk ',VWThk ',VWThk
= ',STND ',STND
R_PLOT5
0484 0485 0486 9487 0488 0489 0490 0491 0492 0493 8494 0495 0496 0497 0498 0499 0500 0501 0502 0503 0504 0505 0506 0507 0508 0509 0510 0511 0512 0513 0514 0515 0516
Symbol
Listing
IF ( Units .EQ. ' ENGLISH ShThkM(I) = ShThk*2.54 VWThkM(I) = VWThk*2.54 ADAR(I)=ADEN/.0142233 ELSE ShThkM(I) = ShThk VWThkM(I) = VWThk ADAR(I)=ADEN END IF
' ) THEN
With or without 30 MLI READ ( 23 ) A46 WRITE ( 10,'( 4X,A)' ) A46 WRITE ( 6,'( 4X,A)' ) A46 420 CONTINUE GO TO 450 440 WRITE ( 10,42 ) 42 FORMAT ( /2X,' No Header following 450 IF (IBOTHS.EQ.1) C C C
.RSP file
' )
RETURN
Close the file and return CLOSE ( UNIT=Z3,STATUS='KEEP'
C C C
Write
500
Rfile to summary
)
file
WRITE ( 10,500 )RFILE FORMAT(lX,'RESPONSE OUTPUT
FILE - ',A )
C RETURN C END
Table for:
The following
RESREAD
symbols
were
defined
but NOT referenced: Offset
RefCnt
Class
BATCOM [Cmn ] DIAM [Cmn ] IC IENV [Cmn ] IT
Scalar
REAL*4
4
Scalar
REAL*4
4
Scalar Scalar
INTEGER*2 INTEGER*2
Z 2
Scalar
INTEGER*2
2
8
0
Scalar
INTEGER*4
4
16
0
[Cmn ] NT
Data Type
Size
Symbol Name Other Information
[Cmn ]
C-19
192624 4
192428
0 0 0 0
R_PLOT5 Usting PID
Scalar
INTEGER*Z
2
114
g
PIDS
Array
INTEGER*2
72
116
0
Scalar
REAL*4
4
3@
0
Scalar
CHARACTER
80
34
g
Array
REAL*4
144
Scalar
REAL*4
4
[Cmn ] Dims RANGE [Cmn ] RSFILE [Cmn ] STANDM [Cmn ] Dims VR
1:36
192284
0
1:36 ZZ
0
[Cmn ] Alphabetic List: Symbol Name Other Information
Data Type
Class
Size
Offset
RefCnt
46 144
338 19Z436
4 4
4
19243Z
4
80 12 12 4
98 384 396 192624
6 2 2 0
4
0
i
CHARACTER CHARACTER CHARACTER
8 8 12
408 416 192632
i i 4
Scalar Scalar
CHARACTER REAL*4
2 4
4Z4 4
i 0
[Cmn ] FORM I IBATCOM IBOTHS
Scalar Scalar Scalar Scalar
CHARACTER INTEGER*4 INTEGER*Z INTEGER*2
80 4 2 2
Z58 454 444 192628
Z 17 i 3
[Cmn ] IC ICB ICT IDENS
Scalar Scalar Scalar Scalar
INTEGER*Z INTEGER*Z INTEGER*2 INTEGER*2
2 2 Z 2
436 434 192430
0 3 7 4
[Cmn ] IENV
Scalar
INTEGER*Z
2
192428
0
[Cmn ] IENVR IPF IPFUNC3 ISHLDS 36
Scalar Scalar Scalar Parameter
INTEGER*2 INTEGER*2 INTEGER*2 INTEGER*4
A46 ADAR [Cmn ] Dims ADEN [Cmn ] ANSWER BIZA BIZB BATCOH [Cmn ] BINC [Cmn ] CBA CBB CONF [Cmn ] DZ DIAM
Scalar Array
CHARACTER REAL*4
Scalar
REAL*4
Scalar Scalar Scalar Scalar
CHARACTER CHARACTER CHARACTER REAL*4
Scalar
REAL$4
Scalar Scalar Scalar
1:36
C-20
442 438 440
1 15 1 8=
R_PLOT5 Listing IT [Cmn ] ITA ITF ITYPE [Cmn ] ITYPEIN [Cmn ] 3 30T K LENGTH NB [Cmn ] NC [Cmn ] NT [Cmn ] NV [Cmn ] PID [Cmn ] PIDS [Cmn ] Dims RANGE [Cmn ] RESPONSE [Cmn ] Dims RFILE ROOTFILE [Cmn ] RSFILE [Cmn ] SHDEN [Cmn ] Dims SHTHK [Cmn ] SHTHKM [Cmn ] Dims STANDM [Cmn ] Dims STND UNITS [Cmn ] VINC [Cmn ] VR [Cmn ] VWDEN [Cmn ] Dims VWTHK VWTHKM [Cmn ] Dims _$CMNBASES
-
-
-
-
-
Scalar
INTEGER*2
2
8
0
Scalar Scalar Scalar
INTEGER*2 INTEGER*2 INTEGER*2
2 2 2
432 430 10
2 2 3
Scalar
INTEGER*2
2
192630
2
Scalar Scalar Scolor Scolor Scolor
INTEGER*4 INTEGER*4 INTEGER*4 CHARACTER INTEGER*2
4 4 4 2 2
458 450 462 96 12
6 2 6 2 3
Scolor
INTEGER*2
2
14
3
Scolor
INTEGER*4
4
16
0
Scolor
INTEGER*2
2
20
3
Scolor
INTEGER*Z
2
114
0
Arroy
INTEGER*2
72
116
0
Scalar
REAL*4
4
30
0
Array
REAL*4
191520
188
4
Scalar Scalar
CHARACTER CHARACTER
80 40
178 192584
11 6
Scalar
CHARACTER
80
34
0
Array
REAL*4
144
191852
2
Scalar
REAL*4
4
192580
6
Array
REAL*4
144
191708
4
Array
REAL*4
144
192284
0
Scalar Scalar
REAL*4 CHARACTER
4 12
426 192644
4 3
Scotar
REAL*4
4
26
1
Scalar
REAL*4
4
22
0
Array
REAL*4
144
192140
2
Scalar Array
REAL*4 REAL*4
4 144
446 191996
4 4
466
1
1:36
1:70
1:19
1:36
1:36
1:36 1:36
1:36
1:36 Cmn Hndls
C-21
R_PLOT5
Stack
Listing
Frame Information:
Temporaries List: Symbol Name Other Information
_$TEMP10 _$TEMP9 _$TEMP8 _$TEMP7 _$TEMP6 _$TEMP5 _$TEMP4 _$TEMP3 _$TEMP2 _$TEMP1 {Work Area} _$TLBO _$TILO _$TLB1 _$TLBZ _$TILAsgnO _STILl _$TEO
Class
Data
Type
Scatar Scalar Scalar Scalar Scalar Scalar Scalar Scalar Scalar Scalar
INTEGER*4 INTEGER*4 INTEGER*4 INTEGER*4 INTEGER*4 INTEGER*4 INTEGER*4 DYNCHAR DYNCHAR DYNCHAR
Scatar Scalar Scalar Scatar Scatar Scalar Scalar
LOGICAL*I INTEGER*4 LOGICAL*I LOGICAL*I INTEGER*4 INTEGER*4 EXTENDED*12
Variabte List: Symbol Name Other Information
Class
LENGTH ANSWER RFILE FORM A46 B12A B12B C8A C8B D2 STND ITF ITA ICT ICB IPF IPFUNC3 IENVR IBATCOM VWTHK JOT I
Scalar Scalar Scatar Scalar Scalar Scalar Scalar Scatar Scalar Scalar Scalar Scalar Scalar Scalar Scalar Scalar Scalar Scalar Scatar Scalar Scalar Scatar
Data
Type
CHARACTER CHARACTER CHARACTER CHARACTER CHARACTER CHARACTER CHARACTER CHARACTER CHARACTER CHARACTER REAL*4 INTEGER*2 INTEGER*2 INTEGER*2 INTEGER*Z INTEGER*2 INTEGER*2 INTEGER*2 INTEGER*2 REAL*4 INTEGER*4 INTEGER*4
C-22
Size
Offset
4 4 4 4 4 4 4 4 4 4 24 1 4 1 1 4 4 12
Size
RefCnt
_ 4 8 12 16 2e 24 28 32 36 40 64 68 72 73 76 80 84
Offset
2 80 80 80 46 12 12 8 8 2 4 2 2 2 2 2 2 2 2 4 4 4
96 98 178 258 338 384 396 408 416 424 426 430 432 434 436 438 440 442 444 446 458 454
2 2 2 2 2 2 2 3 3 3 72 14 6 6 14 2 6
RefCnt
2 6 11 2 4 2 2 1 1 i 4 2 2 7 3 15 1 1 1 4 2 17
R_PLOT5 3 K
Scalar Scalar Cmn Hndls
-$CMNBASES
Variables Symbol Other
Listing
in Blank
Common:
INTEGER*4 ZNTEGER*4
458 462 466
6 6 1
Size:192656
Name Information
BINC
Class
Data
Type
Size
Offset
RefCnt
---
[Cmn ] DIAM
Scalar
REAL*4
[Cmn ] IT
Scalar
REAL*4
[Cmn ] TTYPE
Scalar
INTEGER*2
[Cmn ] NB
Scalar
INTEGER*2
[Cmn ] NC
Scalar
INTEGER*2
[Cmn ] NT
Scalar
INTEGER*2
[Cmn ] NV
Scalar
[Cmn ] VR
Scalar
INTEGER*2
[Cmn ] VINC
Scalar
REAL*4
[Cmn ] RANGE _Cmn ] RSFILE
Scalar
REAL*4
Scalar
REAL*4
[Cmn ] PID
Scalar
[Cmn ] PIDS
Scalar
-
1:70 1:19 1:36
Array
REAL*4
[Cmn ] Dims SHDEN
-
1:36
Array
REAL*4
[Cmn ] Dims V_FFHKM
-
1:36
Array
REAL*4
[Cmn ] Dims VWDEN
-
1:36
Array
[Cmn ] Dims STANDM
-
1:36
1
4
4
0
2
8
e
2
10
3
2
12
3
2
14
3
4
16
0
2
20
3
4
22
0
4
26
1
4
30
0
80
34
0
2
114
0
72
116
0
188
4
INTEGER*2
[Cmn ] Dims SHTHKM
1:36
Array
e
CHARACTER
-
[Cmn ] Dims IENV
1:36
4
INTEGER*4
[Cmn ] Dims RESPONSE
[Cmn ]
4 4
INTEGER*2
191520 144
191708
4
144
191852
2
144
191996
4
144
192140
2
144
192284
0
2
192428
e
REAL*4
Array
REAL*4
Array
REAL*4
Scalar
INTEGER*2
C-23
R_.PLOT5Listing IDENS
Scalar
INTEGER*2
2
192430
4
Scalar
REAL*4
4
19243Z
4
Array
REAL*4
144
192436
4
Scalar
REAL*4
4
192580
6
Scalar
CHARACTER
4@
192584
6
Scalar
REAL*4
4
192624
@
Scalar
INTEGER*2
Z
192628
3
Scalar
INTEGER*Z
Z
192630
2
Scalar
CHARACTER
12
192632
4
Scalar
CHARACTER
12
192644
3
] ADEN [Cmn ] ADAR [Cmn ] Dims SHTHK [Cmn 3 ROOTFILE [Cmn ] BATCON [Cmn ] IBOTHS COrm ] ITYPEIN [Cmn ] CONF [Cmn ] UNITS [Cmn ]
-
1:36
Local Stackframe size: 470 Local Symbols: 178
0517 0518 0519 0520 9521 0522 8523 0524 0525 o526 0527 0528 0529 0530 0531 0532 0533 0534 0535 0536 0537 0538 053g 0540 0541 0542 0543 0544
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C SUBROUTINE RPINPUT C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C PSINPUT writes the program header to the screen and reads in the C summary ouput filename. It also determines the spacecraft exposure C time and operating altitude. C C C C note: for variables contained in the common block refer to the main C listing for definition C C Variable list C C answer = character string representing user input C Psfile = output filenanme C C C INCLUDE 'COMMONRP.BLK' C CHARACTER*ZO BUM'R74 CHARACTER*80 ANSWER
C-24
R_PLOT5 Listing 0545 0546 0547 0548 0549 0550 0551 0552 0553 0554 0555 0556 0557 0558 0559 0560 0561 0562 0563 0564 0565 0566 0567 8568 0569 0570 0571 0572 0573 0574 0575 0576 0577 0578 0579 0580 0581 0582 0583 0584 0585 0586 0587 0588 0589 0590 0591 0592 0593 0594 0595 0596 0597 0598
C C C C C
Write
header
to screen
WRITE , 6,10 ) 10 FORMAT ***********************************:IX**** 1 'Space Debris SURFace',/lX,9X,'R_PLOT5' 2 /1X,SX,'Ver. 1.6 8/23/92',/1X,SX,'FOR 3
C C C
BUMPERIIvl.3',/IX,
******************************
Read in output
filename,
15 WRITE ( 6,20 ) 20 FORMAT ( /IX,'OUTPUT READ , 5,30 )RSFILE 30 FORMAT (A)
set default
FILENAME
to R_PLOTS.RS
,CR=R_PLOT5.RS)>',$)
C IF ( RSFILE,I:I).EQ.' C C C C C
Open
' ) RSFILE='R_PLOTS.RS'
rsfile
!!!! PUT CREATOR='XCEL' OPEN *
OR
'MSWD' IN OPEN
STATEMENTS
ON MAC
( UNITflO,FILEfRSFILE,STATUS='NEW',IOSTAT=IER, CREATOR='XCEL',ERR=40,RECL-256 )
C GO TO 70 C C C
Error control 40 IF ( IER.EQ.2013 ) THEN WRITE ( 6,50 ) 50 FORMAT ( /IX,'FILE ALREADY READ , 5,30 ) ANSWER
EXISTS OK TO OVERWRITE
,,R-YES,S)>')
C IF ( ANSWER(I:I).EQ.'Y' .OR. ANSWERCI:I).EQ.' ' ) THEN OPEN , UNIT=lO,FILE=RSFILE,STATUS='UNKNOWN',IOSTAT-IER, I ERR=40) REWIND 10 ELSE GO TO 15 END IF ELSE WRITE , 6,60 ) 60 FORMAT (/lX,'UNABLE TO OPEN FILE ' ) GO TO 15 END IF C 70 CONTINUE C C WRITE
( 10,75 )
C-25
R_PLOT5 Listing 0599 0600 0601 0602 0603 0604 0605 0606
SymboI
75 FORMAT ***************************************** 1 ' R_PLOT5', 2 /1X,5X,'Ver. 1.6 8/23/92',/1X,5X,'FOR 3 ********************************
BUMPERIIv1.3',/1X,
RETURN END
Table
for: RPINPUT symbols
The following
were
defined
but NOT referenced:
SymboI Name Other Information
Class
ADAR [Cmn ] Dims ADEN [Cmn ] BATCOM [Cmn ] BINC [Cmn ] BUMTTM CONF [Cmn ] DIAM [Cmn ] IBATCOM IBOTHS [Cmn ] IDENS [Cmn ] IENV [Cmn ] IT [Cmn ] ITYPE [Cmn ] ITYPEIN [Cmn ] NB [Cmn ] NC [Cmn ] NT [Cmn ] NV [Cmn ] PID [Cmn ]
Array
REAL*4
144
192436
0
Scalar
REAL*4
4
192432
0
Scalar
REAL*4
4
192624
0
Scalar
REAL*4
4
0
0
Scalar Scalar
CHARACTER CHARACTER
20 12
192632
0 0
Scalar
REAL*4
4
4
0
Scalar Scalar
INTEGER*2 INTEGER*2
2 2
192628
0 0
Scalar
INTEGER*2
2
192430
0
Scalar
INTEGER*2
2
192428
0
Scalar
INTEGER*2
2
8
0
Scalar
INTEGER*2
2
10
0
Scalar
INTEGER*2
2
192630
0
Scalar
INTEGER*2
2
12
0
Scalar
INTEGER*2
2
14
0
Scalar
INTEGER*4
4
16
0
Scalar
INTEGER*2
2
20
0
Scalar
INTEGER*2
2
114
0
-
Data Type
Size
RefCnt
Offset
1:36
C-26
R__PLOT5Listing PIDS [Cmn ] Dims RANGE [Cmn ] RESPONSE [Cmn ] Dims ROOTFILE [Cmn ] SHDEN [Cmn ] Dims SHTHK [Cmn ] SHTHKM [Cmn ] Dims STANDM [Cmn ] Dims UNITS [Cmn ] VINC [Cmn ] VR [Cmn ] VWDEN [Cmn ] Dims VWTHK VWTHKM
-
[Cmn ] Dims
-
-
-
-
-
[Cmn ] ANSWER BATCOM [Cmn ] BINC [Cmn ] BUMTTM CONF [Cmn ] DIAM [Cmn ] IBATCOM IBOTHS [Cmn ] IDENS [Cmn ] IENV [Cmn ]
-
116
0
4
30
0
188
0
192584
0
144
191852
0
REAL*4
4
192580
0
Array
REAL*4
144
191708
0
Array
REAL*4
144
192284
0
Scalar
CHARACTER
12
192644
0
Scalar
REAL*4
4
26
0
Scalar
REAL*4
4
22
0
Array
REAL*4
144
192140
0
Scalar Array
REAL*4 REAL*4
4 144
191996
0 0
INTEGER*2
Scalar
REAL*4
Array
REAL*4
Scalar
CHARACTER
Array
REAL*4
Scalar
1:36
191520
1:70 I:19 1:36 44)
1:36
1:36 1:36
1:36
1:36
Alphabetic List: Symbol Name Other Information
ADAR [Cmn ] Dims ADEN
72
Array
Class
Data
Type
Offset
Size
RefCnt
192436
0
4
192432
0
80 4
27 192624
3 0
4
0
0
20 12
192632
0 0
Array
REAL*4
144
Scalar
REAL*4
Scalar Scalar
CHARACTER REAL*4
Scalar
REAL*4
Scalar Scalar
CHARACTER CHARACTER
Scalar
REAL*4
4
Scalar Scalar
INTEGER*2 INTEGER*2
2 2
192628
0 0
Scalar
INTEGER*2
2
192430
0
Scalar
INTEGER*2
2
192428
0
1:36
0-27
4
0
R._PLOTS IER ISHLDS 36 IT [Cmn ] ITYPE [Cmn ] ITYPEIN [Cmn ] NB
U_ir_
Scalar Parameter
INTEGER*4 INTEGER*4
4
107
Scalar
INTEGER*2
2
8
%
Scalar
INTEGER*2
2
10
0
Scalar
INTEGER*2
Scalar
INTEGER*2
Scalar
[Cmn ] NV
[Cmn ] NC [Cmn ] NT
2
3 8=
192630
e
2
12
8
INTEGER*2
2
14
0
Scalar
INTEGER*4
4
16
e
[Cmn ] PID
Scalar
INTEGER*2
2
20
[Cmn ] PIDS
Scalar
INTEGER*2
2
114
0
72
116
e
4
3e
0
188
0
[Cmn ] Dims RANGE
-
1:36
[Cmn] RESPONSE [Cmn ] Dims ROOTFILE [Cmn ] RSFILE [Cmn ] SHDEN
-
1:?% 1:19 1:36
Array
INTEGER*2
Scalar
REAL*4
Array
REAL*4
Scalar
[Cmn ] Dims SHTHK
-
1:36
[Cmn ] SHTHKM
CHARACTER
Scalar
CHARACTER
Array
REAL*4
Scalar
REAL*4
[Cmn ] Dims STANDM
-
1:36
Array
REAL*4
[Cmn ] Dims UNITS
-
1:36
Array
REAL*4
Scalar
CHARACTER
Scalar
REAL*4
[Cmn ] VINC [Cmn ] VR [Cmn ] VWDEN [Cmn ] Dims V_FFHK V_FTHKM
-
[Cmn ] Dims ° _$CMNBASES Stack
1:36
1:36
Scalar
REAL*4
Array
REAL*4
Scalar Array
REAL*4 REAL*4
Cmn Hndls
Frame Information:
Temporaries
List:
C-28
191529 4_
192584
0
80
34
5
144
191852
e
4
192580
0
144
191708
144
192284
e
12
192644
0
4
26
0
4
22
%
144
192148
0
191996
0 0
111
1
4 144
R_PLOT5
Class
Symbol Name Other Information
{Work Area} _$TLBO _$TLB1 _$TLBZ
Data Type
Scalar Scalar Scalar
Class
ANSWER IER _$CMNBASES
Scalar Scalar Cmn Hndls
in
Blank
Size
Data Type
0 24 Z5 Z6
Offset
Size
CHARACTER INTEGER*4
RefCnt
Offset
Z4 I I I
LOGICAL*I LOGICAL*I LOGICAL*I
Variable List: Symbol Name Other Information
Variables
Listing
80 4
6 Z Z
RefCnt
Z7 107 iii
3 3 I
Size:lgZ656
Common:
Symbol Name Other Information
Class
BINC [Cmn DIAM [Cmn IT [Cmn ITYPE [Cmn NB [Cmn NC
Scalar
REAL*4
4
0
0
Scalar
REAL*4
4
4
0
Scalar
INTEGER*Z
Z
8
0
Scalar
INTEGER*Z
Z
i@
0
Scalar
INTEGER*Z
Z
IZ
0
Scalar
INTEGER*Z
2
14
0
Scalar
INTEGER*4
4
16
0
Scalar
INTEGER*Z
Z
ZO
0
Scalar
REAL*4
4
ZZ
0
Scalar
REAL*4
4
Z6
0
Scalar
REAL*4
4
30
0
Scalar
CHARACTER
80
34
5
Scalar
INTEGER*Z
Z
114
0
Array
INTEGER*Z
72
116
0
Data Type
Offset
Size
RefCnt
] ] ] 1 ]
[Cmn ] NT [Cmn ] NV [Cmn ] VR [Cmn ] VINC [Cmn ] RANGE [Cmn ] RSFILE [Cmn ] PID [Cmn ] PIDS [Cmn ] Dims
1:36
C-29
R_PLOT5
RESPONSE [Cmn ] Dims SHTHKM [Cmn 3 Dims SHDEN [Cmn ] Dims VWTHKM [Cmn ] Dims VWDEN [Cmn ] Dims STANCH [Cmn 3 Dims IENV [Cmn 3 IDENS [Cmn ] ADEN [Cmn ] ADAR [Cmn ] Dims SHTHK [Cmn ] ROOTF ILE [Cmn ] BATCOM
-
1:70 1:1g 1:36
-
1:36
-
1:36
-
1:36
-
1:36
-
1:36
[Cmn 3 IBOTHS [Cmn 3 ITYPEIN [Cmn ] CONF [Cmn ] UNITS [Cmn ]
Local Stackframe
0607 0608 0609 0610 0611 0612 0613 0614 0615 0616 0617 0618 0619 0620
188
0
144
191708
0
REAL*4
144
191852
0
Array
REAL*4
144
191996
0
Array
REAL*4
144
19Z140
0
Array
REAL*4
144
192284
0
Scalar
INTEGER*2
2
192428
0
Scalar
INTEGER*2
Z
192430
0
Scalar
REAL*4
4
192432
0
Array
REAL*4
144
192436
0
Scalar
REAL*4
4
192580
0
Scalar
CHARACTER
40
192584
0
Scalar
REAL*4
4
192624
0
Scalar
INTEGER*Z
2
19Z628
0
Scalar
INTEGER*2
2
192630
0
Scalar
CHARACTER
12
192632
0
Scalar
CHARACTER
12
192644
0
Array
REAL*4
Array
REAL*4
Array
_1520
1:36
-
Local Symbols:
Listing
size:
116
58
C C
5
20 10 ii
SUBROUTINE DIRLIST CHARACTER*80 LINE OPENCUNIT=IT,FILE='DIRECTORY.LIST',STATUS='OLD',ERR=30) REWIND 17 READCI7,10,ERR=ZO) LINE WRITEC6,11) LINE GOTO 5 REWIND 17 CLOSECUNIT_IT,STATUS='KEEP')
FORMAT (A80) FORMAT
(1X,A80)
C-30
R_PLOT5
0621 0622
30
Symbol
Table
RETURN END
for: DIRLIST
Alphabetic List: Symbol Name Other Information
Class
LINE
Scalar
Stack
Listing
Data
Type
Size
CHARACTER
80
Class
Data
Type
Size
24
Variable List: Symbol Name Other Information
Class
LINE
Scalar
CHARACTER
Class
Result
Stackframe size: Symbols: 8
Symbol Name
2
Data
Type
Size
RefCnt
@
Offset
80
RefCnt
24
2
104
Table
Blank Common DIRLIST
Blank Proc Main Proc Proc Proc
MAClI_R_plotS_16$main RESREAD RPSTEXT RPINPUT 0 serious
Offset
24
{Work Area}
Global Symbol
RefCnt
Frame Information:
Temporaries List: Symbol Name Other Information
Local Local
Offset
errors
Cmn Sub Prgm Sub Sub Sub
detected.
0 warning messages 754 lines compiled.
generated.
C-31
Type Size
192656
Other
Information
SD_SURF
Appendix
D.
BUMPERH
User's
Manual
Modifications
D-i
for
the
Macintosh
Appendix
D
Limitations These
are
the
limits
in ver._ion
3.0 of the
31
significant
characters
255
characters
in each
5100 409 -5500
l,anguage
Systems
in a symbolic source
ccxle
symbols compile
(program,
3200
dimensions combined
in a single array
32
levels
of nested
5O
nesting
depth
32
arguments
512 2O
nested
function
2O
nested
repeat
1500
characters
dimensions DO loops
and IF
nested
function
common
but riot
block
names)
names, statemem l-_bcls, and compiler--generated.
module
implied
DO loops
per CALL reference calls
and
factors
subscript
format
complex
32767
character
constants
21,17483647
maximum statement
record
size
for multiple
21,i7483647
maximum
record
size
for formatted
in a program
constants
references
in a format
in a packed
in any
lines,
definition
32767
fields
72)
statements
real constants
2147.183647
and
in a program
32767
341
77 allows
6)
array
in a statement
arguments
77 allows
all continuation
subprogram
for block
actual
(ANSI
line (ANSI
local symbols (including all symbolic subprogram and function references temporaries) 7
32767
narne
character,'_ in a statement (counting counting comment lines) global in one
FOR'VILCa"q Compiler:
module
a program
module
in a program
module
ifl
items
in an unfomlatted
I/O
I/O
stmct-ure
maximum
size of a STRUCTURE
maximum
iterations
element
for a DO loop
D-1
in an array
of RECORDS
Appendix VAX
EXTENSIONS
THAT
ARE
D
ACCEPTED
BUT
NOT
EXECUTED:
options
CLOSE
'SUBMIT'
DISPOSE/DISP=
'SUBMIT/DELETE' DEFINEFILE DELETE DICTIONARY FIND options
INQUIRE
DEFAULTFILE KEYED ORGANIZATION=
'SEQUENTIAL' 'RELATIVE' 'INDEXED'
'FIXED'
RECORDTYPE=
'VARIABLE' 'STREAM OPEN
CR'
options ACCESS=
'KEYED'
ASSOCIATEVARIABLE BLOCKSIZE BUFFERCOUNT DEFAULTFILE DISPOSE/DISP
'SUBMIT'
=
'SUBMIT/DELETE' EXTENDSIZE INITIALSIZE KEY NOSPANBLOCKS ORGANIZATION=
tSEQUENTIAL' 'RELATIVE' 'INDEXED'
RECORDTYPE=
'FIXED' 'VARIABLE' 'STREAM
CR'
'STREAM
LF'
SHARED USEROPEN READ/WRITE KEYID REWRITE UNLOCK
VAX
VMS
FORTRAN
FEATURES
Octal
constant
Extended
Range
Indexed
DO
"77
BLOCK
loops
EXTERNAL
*v
FUNCTION
DATA
subprogram.s
[,*v]
NAI4E
_V
a [unction
D-2
SUPPORTED: PARAMETER
or "77" (man),
constant.s
External
CAI.,L
notation:
files
Radix-50
TYPE
of the form:
statements
PARAMETER
NOT
olher
p=c,
[p=c]
ty|x_.c;o[ octal notation
arc
avail,qblc)
SD_SURF
Appendix
E.
E-i
User's Manual
SD SURF Macro
Appendix E - SD_SURF..MACRO
A
B
1 2 3._ 4
Summary Title:
Space Debris NAS8-38856
Contract: Version:
C
Information Surfaces
Macro
v1.1
m
7
Programmer:
Norman
Corporation: Creation Date:
Martin !Ver
Notice
This
Elfer,
Ph.D.
Marietta
(504)-257-3162
Manned
1.1 - Feb.
Space
Systems
14_ 1992
8 9 11) 1 1
series
of EXCEL
Macros
were
wdttan
in
support of contract NAS8-38856 from NASA-Marshall Space Flight Center.
12 13 1 4
SUBROUTINES
1 5
NAME
16
Auto_op
PURPOSE en
Calls Opening
17 18
Opens
Dialog
Function
Box.
Macro Sheets.
19
Auto= close A surf 0.5
20
Close., Macro
Close SD Surf from menu bar. Continue
21
Function
22
Open/Save/Set... BL PASTE
Pastes
named
23
RPLOT..Open
Opens
FORTRAN
25
Auto
Auto
26
SD_Surf.name
km/s
Deletes menu. Changes A_surf
output
described
and
arrays
to 0.5
krn/s
sets variable from
Ballistic
text file output
increments EXCEL name. Umit
to
PNP
and pastes to PNP
24 open
open
=GET.DOCUMENT(I) =MESSAGE(1,'SD
27
Pull Down Menu added.
Unhide
macro to
modify.')
28
=CALCULATION(3)
29_
=SHORT.MENUS(FALSE)
30
=ADD.MENU(t
31
=DIALOG.BOX(Intro_Dialog_box)
,SD)
=ALERT('Automatic Options. 32
Recalculation
If you don1 want
Apple-period
to
stop
was
to recalculate
turned when
off.
See
saving
use
re-calc.',2)
33
=MESSAGE(0,)
34
=HIDE()
35
-E R ROR(2,OP EN_SD_FUNCTION
36
=OPEN ?('SD_Function_Mecros',,TRUE)
37
=HIDE()
38
=ERROR(I)
39
=Open_PNP_Template0
40
=RETURN()
m
MACRO)
41 42
auto.close
close
Macro
43
=ACTIVATE(SD_Surf.name)
44
=CLOSE()
45
=RETURN(}
46 47
auto
close
auto
close
48
=DELETE.MENU(1
49
=RETURN()
,'SD')
50 51
command
I Open_SD
Function
Macro
E-1
Open
Function
Macro
Calls
a subrou_ne
Sheets
Appendix E- SD_SURF_MACRO A
B =MESSAGE(1,'Please
52 53 54 55 56 57 58 59 60
=OPEN
It
General open box used
?("SD_Function_M
acros",,TRUE)
=HIDE() =MESSAGE(0,) = Open_P
Calls
N P_Templ ate 0
a subroutine
=RETURN()
command
Open
61
BL Template
TEMPLATE. open_template
open
a Ballistic
Limit Macro
Genera/open user
Template",TRUE)
=MESSAGE(0,) Name
wants
a previously
= IF(open_BLTemplate BL_ Template.
box used if
Default is READ ONLY.')
=OPEN?('BL
64 65
Macro Sheet.
will be hidden.')
=MESSAGE(1,'Please
62 63
C
open the SD Funtion
=FALSE,HALT(),)
to select modified
template.
=GET.DOCUMENT(I)
66
=FULL(TRUE)
67
=FORMUI.A.GOTO0$A$1
68
-RETURNI)
,TRUE)
69 70
command
Save
BL Template
=MESSAGE(1,'Save 71
with new or old name.
Apple-.
to stop
!recalculation.') -SAVE.AS?0
73
=IF(B73
74
-MESSAGE(0,)
75
command
Set
BL
=FALSE,HALT(),) Template
76
=SET.VALUE(BL_Tamplate.Name,GET.DOCUM
77
=RETURN()
ENT(1 ))
78 79
command
Open
PNP
Template
=MESSAGE(1,'Please 8O
open a PNP TEMPLATE.
Default is
81
=OPEN?("PNPIFLUX
82
= IF(B81
user
Template',TRUE)
=FALSE,HALT0,)
PNP_ Template.
Nam
wants
a previously
=MESSAGE(0,) 84
Genera/open
box used if
READ ONLY.') to select modified
template.
=GET.DOCUMENT(I) =FULL(TRUE)
86
=FORMULA.GOTO(I$A$1
87
=RETURN
,TRUE)
0
88 8g
command
Save
PNP
Template
=MESSAGE(1,'Sava 90 91
=SAVE.AS?0
g2
=
g3
=MESSAGE(0,)
94
with new or old name.
Apple-.
to stop
recalculation.')
command
Set
IF(B92
PNP
=FALSE,HALT0,)
Template
95
= SET.VALUE(PN
g6
=RETURN()
P_Template.Name,GET.DOCUM
ENT(1 ))
97 98
command
Open_Area_Maker
99
=MESSAGE(1,'Please
open Area
100
=OPEN?('AREA_MAKE
R_MACRO",,TRUE)
101
= IF(B100
=FALSE,HALT(),)
Maker Macro
sheet.')
Genera/open user a
E-2
wants
previously
box used if to select modified
Appendix E - SD..SURF_MACRO
A
B
102 103
=MESSAGE(0,) =RUN(IAuto_C) pen,FALSE) =RETURN()
!105 106 A surf 107 108 109 110
0.5
A_Alert 111 112 11.._33 114 115 116 117
,=====,i
125 1261 m
1281 m
129__. 13o 131 1 321 m
m
13.._.g3 134 m
13_.._s 136 137 m
138m 139 m
140m 141i 1421 1431 BL.PASTE
144 145 1461 7...._ 1 4 148: 149 150 151 m
m
1 53!
A surf • •
0.5
km
s
Use with ASURF output. Cut from 0.25 km/s and add
to 0.5 km/s multiples • Start on first row (0.25 km/s) =ALERT('This will delete every other 0.25km/s A_Surf entry. You must have selected the 1st row to delete.',1) =IF(A_Alert'=FALSE(),HALT(}_) =ECHO(FALSE)
Alert
Alert.
You can change
your mind here. Speeds up Macro Adds first row to second and deletes first row
=S ELECT('R[+0]C2:R[+0]C20") ,,COPY()
11..._S.S 119 120 12_.j._counter 122 123 124
BLP
km/s
C template.
=SELECT('R[+I]C2") =PASTE.SPECIAL(3,2,FALSE,FALSE) =SELECT('R[-1]') =EDIT.DELETE(2) ,FOR('counter',1,33,1 ) = SELECT('R[+I]') = INSERT(2) = SELECT('RC2:RC20") = FORMULA('0.5") = FILLRIGHT() = SELECT('R[+I]C2:R[+I]C20") - COPY() = SELECT('R[-1]C2:R[-1]C20") , PASTE.SPEClAL(314tFALSEtFALSE ) , COPY() SELECT('R[-1]C2") = PASTE.SPECIAL(3,2,FALSE,FALSE) = SELECT('R[+3]C2") PASTE.SPECIAL(3,2,FALSE,FALSE) ,, SELECT('R[-2]') = EDIT.DELETE(i) = EDIT.DELETE(2) -NEXT() .ECHO(TRUE) -RETURN 0 BL PASTE =ALERT('This will paste the ballistic limit surface on the PNP/FLUX WS Template. The ballistic limit template must be active.',1) =IF(BLP Alert .FALSE(),HALT(),) =ECHO(FALSE) =CALCULATION(3mm ) =SET.NAME('BL.name',GET.DOCUMENT(1)) =ACTIVATE(BL.name) =SELECT('r6c5:rl 6c6") =COPY() =ACTIVATE(P N P_Template.Name) =SELECTf'R16C7"_
E-3
Counter Loop Sets up to divide by2
Divides odd cells by 2
Adds to cells above and below
Alert.
You can change
your mind here. Speed= up Macm
Copy am¢lpaste Header
E. SD._SURF._MACRO A 154
1s_._66 lS._Z 158 159 160
s.__! _ 162 163 164 m
16_.._s _1s....!6 6._.Z : i169 m
17o 171 172 173 174 175 RPLOT_Open
7_...66 1 177'
7._..88 1 179
lS..._o 1J.! _182 m
1s-4 104 los
s1....e6 187 mmmm
e._..ae 1 189
0..._o0 1 19._.!
10._.Z= 9....j3 I 19__4 195 197 m
lo-4 2o-4
2o._.11 20-4 20-4 2O4
B -PASTI_,SPI_CIAL(3.1 .FALSE,FALSE) =ACTIVATE(BL.name) .SELECT('R6C10:R10Cll")
C Copy and paste Sample caic
-COPY() -ACTIVATE (P N P_Tem plate. Name) .SELECT('R17C21 ") -PASTE.SPECIAL(3,1 ,FALSE,FALSE) -SELECT('R15C21 ") .FORMULA("Worksheet calculation')
-SELECT('R1
6C21")
-FORMULA(BLname) Copy and pasle BL data
.ACTIVATE(BLname) .SELECT('R21 C4:R52C22") .COPY() .ACTIVATE(PNP_Template.Name) -SELECT('R44C4") .PASTE.SPECIAL(3.1 .FAI.,_E.FAI._) -CALCULATION(I) -ECHO(TRUE)
-RETVRNn iRPLOT..OpIn -WORKSPACE(,,,,TRUE,..) .MESSAGE(TRUE,'Directions will be listed HERE.') -ALERT('See MESSAGE box at lower left.'_2) -MESSAGE(TRUE,'Open the RPLOT Dm Rle') .OPEN?(,,,2) .IF(B182 -FALSE0,HALT(),) .SET.NAME('RPLOTJlm',GET.DOCUMENT(1)) tx MESSAGE(TRUE, Open the RPLOT SUMMARY File') ix OPEN?(,,,2) x IF(B1?8 -FALSE0,HALT0,) x SET.NAME_RSUM.nm'_GET.DOCUMENT(1)) -MESSAGE(TRUE,'Open the R-PLOT Ballistic Urnit Templet= (BL-RPLOT)') -OPEN?CBL-RPLOT'.TRUE,,) -IF(8190 -FALSE0,HALT0,) .SET.NAME('13Lname',GET.DOCUMENT(1)) =ACTIVATE(RPLOT.name) .SELECT(INPUT('Select first diameter on Response Table. (v .0.25 & obl . 0.)',8,'Response Table',,,)) -IF(B194 .FALSE,HALT(),) -SELECT('RC:R[35]C[20]') -COPY0 -ACTIVATE(BLneme) -SELECT('R17C3") -PASTE.SPECIAL(311tFALSEIFALSE) x ACTIVATE(RSUM.nerne) -ACTIVATE(RPLOT.name) .SELECT(INPUT('Salict up to 12 lines of description to be paste to the Template.',8,'Response Descdptlon','R1C13:RC22",,)) .IF(B203 .FALSE,HALT(),) -COPY0 =ACTIVATE(BL.name)
E.4
Sets up message box.
Open Respenu/R-PIot
x Open Response Summary The summary file is not used due to current output
structure.
Open BL Template
Tanster ResponseData
Transfer Response Summary
Appendix E- SD_SURF_MACRO A 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223' 224i 225 226 227
22_ 229 230 231 232 23:3 234 235 236 23"/
23_ 239 24.__0 1241
24._22 243 244 245 246 247 248 249
2so 251 :252 253
B
C
=SELECT('R2C5") =PASTE.SPECIAL(3rlrFALSErFALSE) = SE LECT('R4C11 ") =FORMULA("RESPONSE OUTPUT FILES') x SELECT('RSC11") x FORMULA(RSUM.name) -SELECT('R6C11")
Record Fi/e Names
=FORM ULA(RPLOT.narne) =ACTIVATE (RPLOT.name) =CLOSE(FALSE) x ACTIVATE(RSUM.name) x CLOSE(FALSE) =ECHO(I'RUE) =MESSAGE('I'RUE,'Enter
CLOSE TEXT FILES
SAVE BALLISTIC LIMIT WS Name
to Save
Ballistic
Limit in EXCEL Format.') -SAVE.AS ?(RPLOT.name,1 ,",FALSE) =MESSAGE(FALSE) =IF(B221 =FALSE0,HALT0,) -SET.NAM E('BL.name'tGET.DOCUM ENT(1 )) =ALERT('OK to do PNP Calculation? This will take a few minutest',l ) .IF(Bi25 =FALSE(),HALT0, ) =MESSAGE(TRUE,'Opan the PNP/FLUX WS TEMPLATE. Cancel to usa current Template.') .OPEN?('PNP/FLUX WS TEMPLATE'.TRUE,,) =IF(B228 -FALS E,ACTIVATE (PNP._Template.Name),) .MESSAGE(TRUE,'Transfering data and calculating flux and PNP.') =SET.NAME('PNP.name'TGET.DOCUM ENT(1 )) =ECHO(FALSE) =CALCULATION(3m.,.) -ACTIVATE(BL.name) -SELECT('r2c5:r14c6") =COPY() =ACTIVATE(PNP.name) =SELECT(" R15C4") =PAST_.SPECIAI_(3,1 ,FAI, SE,FAI,$1_) =ACTIVATE(BL.name) ,SELECT('R4Cl 1 :R10C11 ") =COPY() -ACTIVATE(PNP.name) =SELECT('R15Cl 9") =PASTE.SPECIAL(3,1 ,FALSE,FALSE) =ACTIVATE(BLname) -SELECT('R19C4:R50C22") =COPY() =ACTIVATE(PNP.name) =SELECT('R44C4") =PASTE.SPECIAL(3,1 ,FALSE,FALSE) -CALCULATION(t )
=ECHO(TRUE) =MESSAGE(FALSE) =RETURN()
2s...._4 255
E-5
CONTINUE TO PNP?
Open PNP Template
Speed up macro Copy & paste BL_Header_I
Copy 4 pa=te BL_Header..2
Copy and pare BL data
Perform
Calculations
Appendlx E- SD_SURFMACRO
m
I F I G I"l'lJI
1 2
type
3___ OK Button 4 S___Text Text
7_!_ Text e..!. Text 9___Text Text I__90 Text I...1.1 I..._2Text
le m
19 20
K
L
M
DIALOGS
Text Text Text Text Text Text Text Text
Ix I 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5
init/resull Jwid_highltext 337 374 SD Surf for EXCEL 3.0 330 64 CK 8 Space Debris Surfaces. SD SURF MACRO 31 Ver 1.1 - Feb. 14, 1992 50 69 STRUCTUI=IN. DAMAGE PREDICTION AM:) 92 ANALYSIS FOR HYPERVELOCITY IMPACTS 115 Contract NAS8-38856 134 153 NASA - Marshall Space Flight Center 176 Technical Monitors: 199 Grag Olsen Jennifer Robinson 222 245 Joel WilllamNn 272 Martin Marietta Manned Space Systems 295 Program Manager: 318 Norman Elfar 341 (504).257.3162
JY 242 26 26 26 26 26 26 26 26 26 56 56 56 26 26 56 56
we.
eei.lee
o.eeeel
oe.eeo,
eo.eeoe
g_D
weee*ao
e...ele
lee.eel
leeeloeeo,
E-6
eeewei.
eole
names
Appendix E- SD_SURF_MACRO
ol 1 2
3_
P
MENUS Nan Command SO
4
Open R_Piot Output
O
IRl
m
7_
T
Help RPLOT...Cpon Opens R_PLOT output and Pastes to PNP/Rux TeTlate
Open BL Template Open_BLJemplate Ballistic Limit to PNP..1 BL_Paste
Keeps track of file to use as Ballistic Umit template. Copy Ballistic Limit from Active Template
9 10
Open Area Maker Macro
1__.+ +..._2_2 1__3
Open PNP/Flux Template Open_PNP_Template Save PNPIFlux Template Save_PNP_Templste Set PNP/Flux Template Set_PNPJemplate
14 15
Close SD Surf Macro
m
I
IK Status Bar Text
Macro
s_L 6
s
to PN_Jemplat
Open_Area_Maker Area Maker opens A_SURF output & creates new g_ometdes
m
close_macro
I
Keeps track of which file to use as PNP/FLUX template. Keeps track of which file to use as PNP/FLUX template. Keeps track of which file to use as PNP/FLUX template. Closes Macro and deletes SD menu.
E-7
/
AmendbE- SO_SURF_MACRO W
V m
2 3
VARIABLES Auto_Close
s__ Auto_Open 7
j_
A_Alert A surf_0.5_km
s
m.mm, BLP_Alert
1_3.o BL,_P_qTE
=$BS48:$SS49 =$BS2S'.SBS38 -$A$111 :$B$111 .$B$106:$B$141 .'R_PLOTS.Fm xr' =$._144:$B$144 -$B$144.'SB$173
11 BL_Templete.Name 12 closeMicro 1 2 counter
.;B;_S
14 Intro_Dialog box 1_.5_5 Open_Ares_Maker 1_6.60pen_BL_Templata 1.__7.7 Open_PNP_Templata
,,:_FTp3.'_L_20 zSB$99:$B$104 mSB$61:$BSa :$B$80:$B$87
18 OPEN_SO_FUNCnON.MACRO
=_;S_J31:_;S$$7 z"PNP TEMPLATE" =$B$84
1._9.9PNP.name 20 PNP..Temptnte.Name 21 Print_Area 2...2.2RPLOTJ,,n,
2__.s Fa:tJOT_Omn 2 4 RSUM.rmme 2__5 Seve_BL_Templete 2__6 Save_PNP_Template
2_Z.? so 28 29 30
SD_Sud.nmme Set_BL_Templste Set PNP_Templste
=$S$_.'$854S -34
._;A_I :$C$255.$E$1 .'$M$20.$0
-'R_PLOTSJ_S" =$SS175:SS$253 ='R063125M.sum" =$B$71 :$13577 =$B$O0:$B$g6 =$P$3:$T$15 =$8$26 -$BS76:$B$77
.$8$os:_e_us
E-8
X TYPE 0 2 0 2 0 0 2 0 0 0 0 2 2 2 2 0 0 0 0 2 0 2 2
SD_SURF
Appendix F.
User's Manual
AREA_MAKER
F-i
Macro
Appendix F - AREA MAKER
A
LISTING
B
1 2
Title:
Summary Information Space Debris Surfaces-Area
3
Contract:
NAS8-38856
4
Version:
!v1.1
5...L Programmer: 6 Corporation: Creation
MACRO
Norman
Elfer,
Martin
Date:
Notice
February
This
Ph.D.
Marietta
sedas
support
1"1 12
NASA-Marshall
Macro
Space
Systems
1992
of EXCEL
10
Maker
(504)-257.3162
Manned
14,
C
of contract
Macros
were
NAS8-38856
Space
Right
written
in
from
Center.
13 14
COMMAND
15
NAME
MACROS PURPOSE
1._.._6Auto_open 17
Calls Opening
18
Opens Function Macro Sheets. Removes Menus
Opens
19
Auto
close
20
Set
Template
21
Template,
Open Save
Dialog Box and adds Pull Down Menus.
Area Template
Identifies
active
Worksheet
document
as
the Area
Template
Identifies
it as the
Area
!Opens
the Area Template
Saves
Template
22
Template
23
Clear
24
Rectangle
Clears AreaArray & .Descriptions on Template Creates Area Array and descriptions on Area
25
Disk
Template.
26
Cylinder Disk
for user
Area
Array
and
This is done by opening dialog input and
to Rotate
and
creating
than
facets
which
Template
boxes are
sent
to AreaMatrix.
27
2_.8Cone Sphere 29 30 31
AreaS_to_PNP
32
A_Plot
33 34
Close_AreaS
Copies
Area_array
and
Description_Array
to
PNP Template. Mani_ulatio
FUNCTION
Open
A.Piot
Closes
text file and
macro
which
will
compresses start
for EXCEL
auto_close.
MACROS
Input
/output
3S 36
Area_Matrix
Adds facets
Rotate
Rotates
to area
array on Area
Template
Area, Phi, Theta
37 38
facet
/Total Pro/ected Area Phi, Thata, Pitch,
orientation
39
4_9.o Velocity_
Oist
41
Calculates normalized.
probability
distribution,
f(v).
Needs
to be
Yaw / Phi_ Theta velocity, orbital
inclination
/ r(v/
42 44
Auto
open
AMaker_name
Auto
open
Open
Area
Template
,,GET.DOCUMENT(I)
45
m
4..._6
=HIDE()
47
=ADD.MENU(1
Addrnenu ,AreaS)
4_.g_e
=Reset_initial_values0
49
=CALCULATION(3)
s._9.o
-SHORT.MENUS(FALSE)
Sl
-DI/U...OG.BOX(INTRIO
S2
=S ET.VALU
m
Reset Dialog inital values Turn Calculation Off DIALOG_BOX)
E(A.Templatew'Area
Template') F-1
Appendix F - AREA MAKER MACRO us'rING A 63 m
s_j.4 5._3.s 56 67
B
C Make sure Area Template is open
.ERROR(2,TEM PLATE_OPEN) .ACTIVATE(A_TsmplatQ) -ERROR(I) -FULL(TRUE) .FORMULA.GOTO(ISA$1 ,TRUE) -RETURN()
58 59 command
TEMPLATE_OPEN
Called by Auto_open Error
60
s__!1 62
m
6._3 open_template 64 65 66 67
A_TEMPLATE
s_.._.e
.FULL(TRUE)
69 70 71 72
.FORMULA.GOTO(I$A$1,TRUE) -RETURN(} command
7.__.3 74
75.5 m
.MESSAGE(1,'Please open any AREA TEMPLATE. is READ ONLY.') .ERROR(2) .OPEN?('Area Template'.TRUE) =MESSAGE(0,) =ERROR(I) - IF(open_template-FALSE,Alert1(),) .GET.DOCUMENT(I)
76 77 78 79 80 81 82
Default General open box used ff user wants to select a previously modified template.
;Template Save -MESSAGE(1.'Save with new or old name. Apple-. to stop recalculation.') .SAVE.AS?0 . IF(B74 .FALSE,HALT(),) ,,MESSAGE(0,) -S ET.VALUE(A_Template,G ET.DOCU MENT(1 )) .RETURN()
command
Set.Template =SET.VALUE(A Template,GET.DOCUMENT(1 )) .RETURN()
command
ClaN_macro
83 84 85
.ACTIVATE(A_Maker_name) .CLOSE()
m
Does notsave changes.
e_.s 87 88
.RETURN()
89 command
m
90 91 92 93 94 95 96 97 98
auto_close -DELETE.MENU(1 ,'AreaS') X SAVE.AS(,0) -RETURN()
command
Does not save changes.
Reset initial.values -ACTIVATE(A_Maker_narne) .SELECT(IDlalog .boxes. default.values) F-2
Appendix F - AREA MAKER MACRO us'rING
A
B
99
C
-COPY()
100
=SELECT(!Dialog_boxes_lnitial_values)
101
=PASTE.SPECIAL(3,1,FALSE,FALSE)
102
=RETURN()
i
103 104
command
Clear
Area
Array
Clears
105
=ERROR(2,Alert2)
106
-ACTIVATE(A_Template)
i
'107
-ERROR(I)
108
-ECHO(FALSE)
and
Area
array
descriptions
=SE LECT(!Area_array) 109 -FORMULA.FILL(0) 110 111
=SELECT(IArea_Descriptlons)
112
,CLEAR(3)
113
-ECHO(TRUE)
114
-RETURN()
115 116
Rectangle
11..__7
=ERROR(2,Alert2)
118
= ACTIVATE(A_Template)
119
=ERROR(I)
120
=FORMULA.GOTO(
IAxm,TRU E)
,DIALOG.
BOX(Rectangle_Dialog_box)
122
=IF(B121
-FALSE,HALT0,)
123
-ECHO(FALSE)
i
124
Area.Rec
=Area_Multlplier.rec'Langth.rac'Height.r
Phi.Rec 12.___5 126
=Rotate(90,0,
Theta.Rec
127
=Area
128
=COUNT(OFFSET(IArea_Dascrtptions,0,0,1 =B128+1
i
Matrtx(Area.RectPhi
129
Number.Rec
130 131
Description.Rec No.
=Number.Rec
132
Geom
Rectangle
133
i
L1
=Langth.rec
134
L2
=Height.re¢
135
L3
136
Multiplier
i
137
Pitch
138
Yaw
Calculations
=Rotate(90,0,Pitch.rec,Yaw.rac)
i
i
ac
Pltch.rac,Yaw.rec)
.Rec_Thata.Rec)
Call ))
AreaMatrix
Get current no. of geom.
In area array.
=Area Multlpllar.rec -Pitch.rat =Yaw.rec
i
139 i
tat.
Start
140
I.at.
Rnish
i
141
Incr.
Long. Start 14_...22 14_.._3 Long. Finish 144 Surf Area [m^2] 14.__ss 14_.._6 14__! i148 i
=Area.Rec -ACTIVATE(A_Maker_name)
Description
=SELECT(Description.Rat) =COPY() ,ACTIVATE(A_Template) -SELECT(OFFSET((IAre-..Delmdptiona),Numbar.Rec-
14.._.99
1,0.1.1))
+150
,,PASTE.SPECIAL(3,1,FALSE.TRUE) -ECHO(TRUE)
115= 53
,=RETURN()
F-3
Paste
F-AREA MAKER MACRO USTING A
B
154
Disk
5.._.._s 1
=ERROR(2,Aiert2)
156
=
lS..__Z7
=ERROR(I)
.158
=FORMULA.GOTO(IAxu,TRUE)
C
ACTIVATE(A_Template)
=DIALOG.BOX(Disk_Dialog_box) 160
=IF(B15g
161
-ECHO(FALSE)
162
Area.risk
=FALSE,HALT(),)
=Area_Muftipller.dsk
=Rotate(90,0,Pitch.dsk,Yaw.d
165
-Area
166
=COUNT(OFFSET(IArea_Descriptions,0,0,,1 =B166+1
167 168 __
Description. No. lS_.._Eg
Call
) ))
Area
Matrix
Get current no. of geom. in area
array.
=Number.dsk =Radius.dlk
173
L3 Multiplier Pitch
Surf
Phi.dskrTheta.dsk
Disk
174
182
Matrlx{Area.dskr
dsk
17o Geom 17_ L1 L2 17..._2
Long. 18._._0 Long. IS.__j.1
sk)
=Rotate(90,0,Pltch.dsk,Yaw.dsk)
Number. dsk
Yaw 17..._.6S Lat. 17_..._7 I.at. 17..._._S 179 Incr.
Catcufatione
°Pi() "Radlu s.dsk*2
PhLdsk 16..._33 164 Theta.dsk
=Area_Multiplier.dsk =Pitch.dsk =Yaw.dek
Start Finish
Start Finish Area
[m^21
=Area.dsk
183
=ACTIVATE(A_Maker_name)
lS._.__4
=SELECT(Descrlption.dsk)
185
-COPY()
s._.6s 1
Description
Paste
.ACTIVATE(A_Template) -SELECT(OFFSET((IArea_Descdptions),Number.dsk-
s.._Z 1
1,0,1,1))
188
=PASTE.SPECIAI.(3,1,FALSE,TRUE)
189
-ECHO(TRUE)
lg0
-RETURN(}
il
191 192
Command
Cylinder =ERROR(2,Alert2)
19:; 194
= ACTIVATE(A_Template)
195
=ERROR(I)
196
=FORM U LA.GOTO(=J_xm,TRUE)
197
=DIALOG.BOX(Cylinder_Dialog_Box)
9._.8a 1
=IF(B197
199
-ECHO(FALSE)
2oo
_angle.cyl-etart_angle.cyl)l180 num. of. facets, cyl
=(finilh_angle.cyl-start_angle.cyl)/facat_angle.cyl =IF(OR(B201 =
INT(B201
ALERT('Finish_angle
203
Start
angle,
2o4 2o_._5 +Facet.Area. 2o.__6
-
C)_dw0
207
Calculations
=Area_Multlpller.cyl'Pl()'radius.cyl*Length.cyl'(finish
Area. cyl 201
-FALSE,HALT(),)
and the
),B201 180,360-B260,(IF(B260180),-90,90)
to
+x direction.
=
Rotate(InitlaI.Phi.cone,lnitlaI.Theta.cone,Pitch.cone,Yaw.
2s(;
cone) Theta. r.'tno Rotata(InitlaI.Phi.cone,
lnttlaI.Theta.cone,Pitch.cone,Yaw.
cone)
26._._4 265
.
266
,,NEXT() =COUNT(OFFSET(IArea_Descriptions,0,0,,1
267 1268
Phi.cone,Th
eta.con e) ))
Call
AreaMatrix
Get current no. of geom. in area array.
Description.cone ,Number.cone
No.
Cane
27.._1 Geom 27._2 L1 273
Matrix(Facet.Area.cone,
=B267+1
Number. cone
26._sDescription.cone 270
Area
-Radius_for.cone
L2
,Radius
aft.cone
27._._4L3 275 Multiplier Pitch 27.._..66
.Length.cone -AreaMultiplier.cone
277
Yaw
=Yaw.cone
270
LaL
Start
.Start_angle.cone
279
LaL
Finish
=Finish_angle.cone
=Pitch.cone
28(1 incr.
.facet_angle.cone
Long. Start 28.....!1
28(1 Long. 2831 Surf
Finish Area
=Area.cone
[m^21
2s4__ j__ 2e 28(;
=COPY()
287
•,ACT IVATE(A_Templata)
Description
-ACTIVATE(A_Maker_name) .SELECT(Description.cone)
=SELECT(OFFSET((iArea_Descriptions),Number.cone-
2e(;
1,0,I
289
=PASTE.SPECIAL(3,1,FALSE,TRUE)
29(; 291
,ECHO(TRUE)
,I))
=RETURN()
292 293
Whole,.Sphere =INPUT('Entor
Area. WS
starting
2g.._._4
=tF(B294 Radius. HIS
297
AreaMultiplier.
the area of the sphere
-4"PI0"(Radtus)^2",,)
=FALSE,HALT(),)
=SQRT(Area.WS/P WS
[m^2] or a formula
an equal
sign.',l,'Sphero',"
295 296
with
-INPUT('Entor
10/4) an
area
multlpllor',l,'multiplier',l,,)
29(; 299
=IF(B297
=FALSE,HALT(),)
3O(;
=SET.VALUE(Area.WS,Area.WS'Area_MultlplIer.WS)
=ECHO(FALSE) =FORMULA.ARRAY("-R
3o_.!1 302
one Sphere_areaS'Area.WSI4"
Whole_Sphere_areaS) =ACTIVATE(A_Maker_name)
m
30(1
=S ELECT(IWhole_Sphere_areaS)
304
=COPY()
3O5
=ERROR(2,Alert2) F-6
Paste
Appendix F - AREA MAKER MACRO A
LISTING
B
306
-ACT
30._.._7
=ERROR(I)
3O8
C
IVATE(A_Template)
=SELECT(IArea_array)
309
=PASTE.SPEClAL(3t2wFALSErFALS
310 311
=COU NT(OFFSET(IArea -B310+1
Number. WS
312
Description. N_ 31.__=_ 314
WS
Geom
Whole
L1 31._._S L2
317
L3
T
))
Get current no. of g,.om. in area
array.
Sphere
=Radius.We
Multiplier 31.__s 319
,
'DescripUon. W$ =Number.We
m
316
E)
Descriptlons,O,O,,1
-Area_Multlpller.WS
Pitch
Yaw 32.._.00 32"1 tat. Start !tat. Finish 32..__.22 323 Incr. 32._.._4Long. Start 32._._.55 Long. Finish 326 Surf Area rm^21
=Area.We
327
-ACTIVATE(A_Maker_name)
328
=S ELECT(Description.We)
32..._9 330
=COPY()
Description
Paste
-ACTIVATE(A_Template) =SELECT(OFFSET((!Area_Descdptions),Number.WS.
331'
=1,0,1,1))
332
=PASTE.SPECIAL(3.1,FALSE,TRUE)
333
-ECHO(TRUE)
334
-RETURN(}
335 336
command
Sphere
337
-ERROR(2,Alert2)
338
= ACTIVATE(A_Template)
339
-ERROR(I)
340
=FOR MULA.GOTO(=J_xes, TRUE)
341
=DIALOG.BOX(Sphere_Dialog_Box)
342
-IF(B341
343
=ECHO(FALSE)
!34_ num./at.facets. inure.long.facets. 3j..___5
Sph
346
34"1
Sph
-FALSE,HALT(),) Calculations
=(flnish_lat.sph-start_Lat.sph)/lacet_angle.sph =(finish_Long.sph-start_Long.sph)/facet_angle.sph -IF(OR(start_Lat.sph
90,))
=
must be equal or between
+90
ALERT('tatitudes
=
349
=END.IFO =IF(OR(B345INT(B345),B344INT(B344),B345 180 in Area_Matrix
Macro',2),)
3O
=IF(ABS(F25)>180,HALT(),) Theta.rad.AM
=Theta.AM*Pl()/180 =Phi.AM
3___2 Phi.deg.AM 33
and
Input
24
31
area
=Theta.AM
Theta.deg.AM
34
Parameters
1
VlNC.AM
0.5 [km/s]
36
AINC.AM
5 [deg.]
37
Orb. VeI.AM
8
35
38
=
.39
=FOR('Velocity.AM',0.5,2*Orb.VeI.AM,VINC.AM)
4._9_0 Velo
city. A M
41
Velocity.AM
=
FOR("threat.quad.AM",*l,l,2)
Threat.quad.AM
=
43
Threat.Ang.rad.AM
=
44 Obliquity.rad.
AM
45
ACTIVATE(A_Template)
=
42
Start loop the threat
through velocities
Loop once for Port & once Starboard
threat.quad.AM threat.quad.AM'ACOS(Velo_ty.AMI2/Orb.Vel.AM)
=
IF(ABS(Theta.AM-Threat.Ang.rad.AM*1801PlO)7
1/3)
CriL diam [cm]
=(((t rear wd.cm.WPL'(Material_Strongb%ksi.WPIJ40)A0.5 +t_bumper.c:m.WPL)/(1.248" Den $ity_Pmj.WPL^0.5" COS(Obl quity.deg.WPL'Pl(y180)))^(18/19)'(1.T5-V_normaLWPtJ4)) 3PLATE
filename?
>PLATE
ON ON
EDGE.ASB
EDGE.AST
ONE Area Fraction Table will be created from ALL of the ranges of element IDs selected. INPUT THE STARTING AND ENDING ELEMENT ID FOR EACH ENTER D OR WHEN DONE
RANGE
RANGE 1 IN THE TABLE. STARTING ELEMENT ID : 1 ENDING ELEMENT ID : 2 RANGE STARTING GEOMETRY
2 IN THE ELEMENT OUTPUT
Debris
RANGES-
F_
(-STATION.GEM
) > PLATE. GEM
Analysis
JSC-7/90 Processing
TABLE. ID :
M_-o
Property 1 First
PID-
ID
4 4
EFF.
AREA-
0.35370
The Area Surface file is complete. binary filename: PLATE ON EDGE.ASB text filename: PLATE ON EDGE.AST STOP
I-I
Feb
17,
1992
Output
from
Space
Debris
Ver.
OUTPUT
MacII_P_surf
1.4
Mon,
SURFace 1/31/92
FILENAME
(CR-SDSURF.PS)>PLATE
ON
EDGE.PSURF
ENVIRONMENT ? 1-JSC 20001&6000 2- 7/90 MEMO ANSWER 1 OR 2 > 2 SOLAR FLUX 1-NOMINAL 2-MINIMUM 3-CONSTANT ANSWER 1-3 DATE SPACE
TO
LEVEL
>
BEGIN
EXPOSURE
( 1994-2025
STATION
EXPOSURE
TIME
OPERATING OR ENTERAN
RESPONSE
Binary
OUTPUT
Constant one
case
Output
FILENAME density
in the
equal
file
2 I
3 I
4 I
increments
5 I
6 I
0
..,........,..............,.........,,o....,.,,.,,.,,..,..,..,.,
5
..,.....a.eo......i.o....J..W.O...,..,..J.Ieoie........,..,.,..,
0
,,.,,...,............................,........t.....,.,.....,,.,
5
..............e................i........,.-..-........-,......,,
70 65 60 55 50 45 40 35 30 25 20 15 10
will
:PLATE > ONE
be
RESPONSE FILE: ONE EDGE.ASB Total Flux x Area
IMPACT 1 I
>1
(_I0.0)
(_388.92) IN 5KMILES >
File'
RESPONSE
CONTOURS
Obl Deg
>
ON
EDGE.ASB
RESPONSE.RSP
threat
ON
at
(-1995)
(-STATION.RSP}
RESPONSE PID: 1 A SURF FILE: PLATE PNP(%)= 99.99709 .12345
)
(YEARS)
ALTITUDE(100.-500.km) "E" OR "e" TO ENTER
Area_Surface
The
?
from
VELOCITY 7 I
8 I
9 I
used
RESPONSE.RSP x
Time
(NAT)
0 to max
- 0.29084E-04
NAT
-
12 I
13 I
0.20709E-05
km/s 10 I
11 I
14 I
15 I
16 I
....................................................... 111 ...... ..................................................... 11111 ...... ................................................... 2111 ......... ................................................ 14.11 ........... ........................................... 11.5111 .............. ........................................ 21.31.11 ................ .................................... 14..31 ...................... ................................. 32.12 .......................... .......................... 11.15..11 ............................. ...................... 21..51.11 ................................. .................. 21..31 ........................................ ............. 12...11 ............................................ ..... 11..11..11 ................................................. 1-2
Feb
17,
1992
Output from MacII_P_surf 5
.II..Ii..ii
Mon, Feb 17, 1992
.....................................................
elle,o,e,,,e...e,ee,eoeeo,oe,e,oeeee,,aoe,eeeeoe,e,eeee,ee,e,,e,
The PNP Surface filename: PLATE
file is complete. ON EDGE.PSURF
I-3
PLATE_ON_EDGE.AST.xl GE]4_t
_
FILE -
GZZ_4L"ID_ IW_t'E:RS
l)IAIg.$
:
z,w.t - M._Zt.UNZ Ou¢_t - z,umc._,m Sumaxy o m.SUH ttzr=_
(1 2 Her.mm_
Envlz_uu
2 Hew)
1 2
Nmdo_ of Threats INCI.INkT'ION _ MI4GES- 1 Flrst 1
eZl)..
_/lmD
1
28.5 4 I_;'F. N_k 1
0
5
-
0.353/0
2 )O
15
0.25 O.5
O.0OE_OO O.nnr,,G) 1._E-.O4 1.11E-O4
0.0(3:4.00 O.0CE*O0
O.OOE_O O.OOE_O
O.(X_O0 O.00E40)
0.75
2.9_...05
O.0CE,IO0
O.OOE400
O.OOEH)O O.0CE_O
1 1.2S
O.OOE4QO O.01ZtO0 O.0CE*O0 O.COE_O
O.OOE*CO O.GCE_0 O.(X_*OO 0.00_0
O.OOE4CO O.0OE4OO O.OOE_O O.0CE_O O.0(]E*GO O.0OE_O0
O.COE4_O
1.).1E,.4)4
1.5
I ._
45 -
1.95E-05
25
30
)5
40
O.0OE_O O.OCE4QO O.OCE+O0 O.OCE*O0 O.0OE4OO O.0CE*O0 O.00G*OO O.OOE_O
O.OOl:_O
5)
55
ID
_5
10
_
gO
O.OOE*O0 O.OOE+OO O.OOE*O0 O.00E*O0 0.00E4400 O.OOE_O O.(X)l:4OO O.O0_*(X) O.OOC*O0 O.OCE_O
O.0OE*OO 0.0(O_0
O.00E400
0.CO1_40)
O.00E_)O
O.OOE_O
O.OOE_O0 O.0OE_O0
O.OCE*O0
O.OOEH)O O.0CE44)0 O.(X_OO O.CCE_O0
O.OOEH)O O.COE:*O0 O.OCE_O O.OOE*G) O.OOE_I) O.OCE_O O.OOE_O O.O01E4G) O.GOE*OO O.OOE_00
O.OOE*O0 O.OOE_)O O.OOEH)0 O.GOE400
O.OOEH)O O.0OE_O O.OI_4G0 0.00E44)O O.OOE4QO O.OOE_)O
O.GOl:*O0
O.OOE:_O
O.00_OO
O.OOE+O0 O.(X_E_O
O.OIXH)O
O.00E.mO O.CO_*O0
O.OOE,_O0 O.0OE.oO) O.COE+O0 O.0OE4G)
O.OOE4G)
O.00_*O0 0.00£4.00
0.O0t_O O.0_E_O O.OOE_0 O.COE',O0 O.0CEoQO O.OOE+G)
2.25 2.5
O.OOE*OO O.OOE+G) O.OOE:+O0 2.13_.-09
0.OCE+O0 6.60E-04
O.OOE*O0 O.OOEH)O O.0CEtOO O.00E_G) O.0OEtO0 O.OOE+O0 O.OOE_O O.OOE+O0 O.OOE;+O0 O.OOE+O0 O.OOE+O0 O.OOE+O0 O.OOC+O0 O.OOE_O
2.75
O.OOE:+O0 4.90E:-09
l.g_F,.-.03
O.OQE400
3 3.2'3
0.0_*00 O.COE_O O.OOE+O0 O.O0_tO0
O.(X_*O0 0.001:,00
O.OOC*O0 O.OOE:_O0 O.OOE_O O.OOC+O0 O.OOC+(X) 0.00_*00
0.00_0 0.0(_*00
3.5 3.75
O.(X)_O O.OOE;GO O.OOE*OO O.G)E+O0
1.21£-,O4 1.12E-O)
4.83_-O4 4.47E-O3
O.OOE*O0 O.(X)£:_O O.OOE_O O.CCE_O O.00E*O00.OOE4G)
O.O0£_O
O.OOE4G) O.(X)_O
O.0(_+O0
0.COE_O
O.OCE4OO O.OOE*G)
O.00E,_O O.G]E,_0 O.0OE*O0 O.0Cl_+O0 O.OCE4QO O.O(]E+O0
O.OOE_O
O.00_*00 O.OCE_O O.OOE+O0 O.OOE_O
0.00_)0 O.OOf:+O0 O.OOE_O 0.00_+00 O.OOE400 O.OOiC+O0 O.OOE:+O0 O.OOr.+O0 O.OOE+O0 O.(X_+O0
O.OCE_O
O.OOEtO0
O.OOE,_O O.OOE*O0
O.OCE400
O.OOE*O0
O.OOE4QO O.(X]E_O
O.OOE*OO O.00_*00 O.OOEo00 O.OOE+G)
O.OCE_O O.0CIE+G)
O.OOE400 O.OOE_OO O.OOE+O0 O.OOE+O0 O.OOE_O O.OOE+O0 O.OOE*O0 O.OOE400 O.OCE,_O0 O.OOE+O0 0.0(_+00 O.OOE+O0 O.OOE+O0 O.OOE+O0 O.OOE+O0 O.OOE+O0 O.OOE:+O0 O.OOF.+O0
O.OOE+O0 O.OOE:+O0 O.OOE:+O0 O.OOE+O0 O.OOE:+O0 O.OOE:+O0 O.OOE+O0 O.O01E:',O0 O.OOE+O0 0.001£+(30 0.001_+00
O.OCE*OO O.CCE_O0 O.0OE*G)
6)
O.0OC*O0 O.0OE*O0 O.OOE+(X) O.00_*00
O.OOE_O
O.OOE_GO O.GG_)O O.GOE_00 O.GOE*G)
_
O.(X)E_0 O.COE*O0 O.OOE*O0 0.00_*00
O.00_*O0 4.1_-.O4 1.O4F,-O4 O.OOEK)O O.OOE+O0 O.OC_+O0
O.(X_*OO
O.00£_G)
O.OOE_O
45 O.CO_*O0 O.00_O
1.15 2
4
5.)2E-O4
5n
O.OOE+O0 O.OOE+O0 0.00_+00
O.OOE_O0
O.OCE*O0 O.OOE+O0 O.OOf:_O O.OOE_O0 O.OOE+O0 O.(X)f:*O0 0.0_)0 O.OOE+O0 O.OOC+O0 O.OOE+O0 O.OOE+O0 O.OOF.*O0 O.OOE+O0 O.OOE+O0
O.GO£H)O O.OOE_O
O.I)O£+O0 O.OOEtO0 O.OOE4QO O.OOE_O
O.OOEH)O O.OOE+OO O.GOE*OO O.OOE*O0 O.OOE_O0 O.OOE*O0 O.OCE400 O.OOE+GO O.OOE4OO O.GOE400 O.GOE4{X) O.OOE+O00.OOE+O00.(X)E_O
O.OCE_O
O.0(_400
O.OOE:4QO O.OOE*G)
0.00£400
O.COE*O0 O.0CE_OO
O.(X_*GO
O.00E*O0
O.OOE400
4.2S 4.5
O.OOE:+O0 O.OOE+O0 O.OOE+O0 O.OOE+O0 O.OOE+O0 O.OOE:+O0 O.OOE*O0 O.OOE+O0 O.OOE+O0 O.OOE+O0 O.OOE+O0 O.OOE:,_O0 O.OOE:+O0 O.OOE+O0 O.OOE+O0 0.00_+00 O.OOE:*O0 O.OOE+O0 O.OOE:+O0 O.OOE_O0 O.OOE+O0 O.OOE+O0 O.OOE+O0 O.OOE:+O0 O.OOC+O0 O.OOi_+O0 O.OOE+O0 O.OOF,,+O0 O.OOE+O0 O.OOE*O0 O.OOE*O0 0.00_+00
O.OOE+O0 O.(X)E:+O0 O.O0_tO0 O.OOE:+O0 O.OOE:+O0 O.OOE+O0
4.75
O.(X)C*O0
O.OOE:.lO0 O.OOE:+O0 O.OOE+OQ
O.OOE+O0 O.OOE*O0 4.6_E:-03
6.9gE-03
O.OOE+O0 O.(X)E+O0
0.00(+00
O.OOE:+O0 O.OOE+O0 O.OOE+O0 O.OOE:+O0 O.OOC+O0 O.OOC+O0 O.OOC+O0 0.00_+00
,5 O.0CIE+O0 O.OOE_00 O.00E_00 I.?IE.,04 2.57£.-04 5.25 O.OOC_O0 O.OOC*O0 O.OOE:+O0 O.OOE+O0 0.00(+00
O.OOE_O O.00E*00 O.00E_O0 0.00E+00 O.OOE+O0 O.OOE+O0 O.OOE;+O0 0.00_+00
5.5 5.75
O.OOE:+O0 O.OOE:+O0 O.OOE_O0 O.OOE+O0 0.00£*00 0.00£',00 O.OOE:*O0 O.OOE+O0 O.OOE_O0 1.13(-02
O.OOE:+O0 O.OOE+O0 O.OOE:+O0 O.OOE:+O0 O.OOE+O0 O.OOE+O0 O.OOE:+O0 O.OOE+O0 0.00_+00 7.50£-03 O.OOC+O0 0.00(+00 O.OOE:tO0 O.OOE+O0 O.OOE+O0 O.OOC+O0 0.00£*00 0.00£+00
O.OOE+O0 O.OOE+O0 O.(X)E+O0 O.OOE+O0 O.OOE*O0 O.OOE+O0 O.(X)E+O0 O.OOE+O0 O.OOE+O0 O.OOE:_O0
O.OOc+O0
O.OOE:_O0 O.OOC_O0 9.22(-04
6.15£-04
O.OOC+O0 O.OOC+O0 O.OOC+O0 O.OOE+O0 0.00(+00
0._,_00
0.0(0:+00 0.0(0:+00 O.OOC+O0
6
0.00£+00
O.OOE:+O0 O.OOE+O0 0.00_+00
O.OOE_O0 O.0OE*O0 0.OCIE_O O.OOE+O0 O.00E*OO O.0CIE_O O.OOE*O0 0.00E+O0 O.OCIEK)O O.OOEN]O O.OOE+O0 O.OOE+O0 O.OOE+O0 O.OOE+O0 O.OOE:+O0 O.OOE:+O0 O.OOE:+O0 O.OOE:+O0 O.OOE+O0 O.OOE+O0
O.OOE+O0 0.00£*00
O.OOE:+O0 O.OOc*O0
0.00C+00
6.25
O.OOE:+(X_0.00£+00
6.5 _.';5
O.OOE;,O0 O.OOE+O0 O.OOE;+O0 0.0(0:+00 O.OOE+O0 O.OOC+O0 O.OOE+O0 0.00_+00 O.OOC_O0 0.00_+00
O.OOE+O0 2.31£-0_
0.0(0;+00 0.0(0;+00 O.OOE+O0 0.0(0:+00 O.OOC_O0 O.O01E;+O0O.COF.+O0 O.OOE;,O0 O.OOE_O0 5.T;E-.03 O.OOC+O0 O.OOE*O0 O.OOE+O0 O.OOE+O0 O.OOE:+O0 O.OOC+O0 0.00£+00 0.00£+00
0.00(+00
8.49_-_
2.121[-05
O.OOE:+O0 O.OOE:+O0 0.00_+00
O.OOE+O0 8.0_-09
0.00£_ 7.06(-03
O.OOE+O0 0.00£+00
7
O.OOE_O0
O.OOE+O0 O.OOE+O0 O.OOE+O0 0.00£+00
"/.25 0.0(0;_00 0.0(0;+_ O.O(]F._O0O.OOE+O0 "/.5 O.OOC+O0 O.OOC*O0 0.00(+00 0.00(+00 7._
0.0(0_ 8
0._00
O.O()E+_
O.(]()E_
0.0(0:_ 0.00£+00
0.0(0:+00 0.0(0;+00 3.23E_
2.8]E-.02 0._
0.0(0:+00 0.0(0:+00 O.OOE+O0
O.OOE+O0
0.0(0.+00 0.0(0_+00 O.OOC*O0
0.0(0:+00 O.OOE:_O0 O.OOE_+O0 0.0(0_+00 0.0(0:_ 0.00(+00
O.OOE+O0 O.OOE+O0 0.0(0;_ O.OOE;+O0 O.OOC+O0 O.OOE;+O0
O.OOC+O0 O.OOE:+O0 O.OOE:_O0 O.OOE:*O0 O.OOE+O0 O.OOC+O0 O.OOC*O0 O.OOF.+O0 O.OOE+O0
0.0(0;_ 0.0(0:+00 O.OOE+_ 0._+00 O.OOE_O0 0.0(0_00 O.OOE+O0 O.OOE;+O0 O.O01E_ O.OOE+O0 0.00_+00 O.OOC+O0 O.OOF._O0 O.OOE:+O0 O.OOC*O0 O.OOC+O0 O.(X)E:_O 0.00_+00
0.0(0_+_ O.OOE_ O.OOE:tO0 O.OOE+O0
O.OOE+O0
0.0(0:_
O.O01E_
O.OOE+O0
0.0(0:*00
0.0(0;_00 O.OOE_
0.0(0:+00 0.00£+00
O.OOE*O0
0.00£+00
0.0(0:+00 0.00£+00
O.OOF..tO0O.OOE;*O0 0.0(0:+00 0.0(0:_
0._+00
O.OOE_O0
O.OOE+O0
0.0(0;+00 0.0(0_+00 0.0(0;+00 0.0(0:_00 O.O00tO0
8.25 8.5
0.0(0_+00 O.OOE+O0 O.(]OE_O0 0.0(0;+_
O.OOE+O0 O.OOE+O0
0.00£+00 0.0(0;+_
0.00(+00 0.0(0:+00 O.OOC*O0 O.OOF.+O0 O.OOE;+O0 4.14E-03
O.OOF..+O00.0(0;_00 O.O(O:+O0 O.(X_i;+O00.0(0;+00 O.O(O.tO0 0.0(0:*00 O.OOE+O0 1.66£_ 0.00£+00 0._ 0.00£+00 0.00£+_ 0.0(0:_00 0.00£+00 O.OOE+_
O.OOE+O0 0._+(]0
O.OOC+O0 O.OOE*O0
0.0(0;+00 0.0(0:,'00 O.OOE_ O.OCIE+O0
8{'/5 9
0.00E+00 O.OCE*00
0.00E+00 0.00E+00
0.00£*00 O.00C+00
0.00E+00 0.00£*00
1.3'_-O2 0.00E*00
0.00E_00 0.00E;+00
O.00E*O0 0.00E*O0
0.00E+00 0.00£+00
0.00E*00 0.00£*00
0,OOE+00 0.00C_O0
3.34E-03 O.00£m(X)
0.00£._00 O.OOE+00
O.00C+00 0.00E+00
0.00£+00
O.OOE*O0
0.0(0:+00 0.0(0;+00 O.OOE+O0
0.OOE+00 0.00E,00
0.00E+O0 0.00£+00
0.00E*00 O.00E+O0
O.00E+00 0.00£_00
0.00£_00 0.00£+00
O.OQE+00 O.OOE+G)
PLATE_ON_EDGE.AST.xI 9.25
0.oce:+O0
0.0oE:+_
0.00E:_
0.00E+O0
0.00E+00
0.00E+O0
O.OOE+00
1.07E.-O)
1.61E-.0)
0.00[:+00
0._
0.00£+00
O.0CIE_
0.0CIE:_
O.00E_
0.00E_
0.00£_
0.0CIE_
9.5
O.OOE:+O0
0.0(0:+00
O.OOE+_
0.00E,O0
0.00£H_0
0.00E',00
O.00EH_O
I._--02
1.9(3_:.-02 O.00E+00
0.00E+00
0._+00
0.00E+_
0._+00
O.(X)£H_O
0.00£_
0.00E_
0._
0.00£,_(X)
0.00E+00
0.COg+00
0.00E:+(X)
0.COE+00
0.(X_E+O0
0.0QE+00
0.00(_00
0.00E+O0
0.00E:+O0
0.00£+00
0.(XJ£+O0
0.CX)£+C0
0.00E+(X)
0.00£_
0.00£+00
0._.,O0
0.00EH)0
0._
O.(X]Ii:_4]O
0.(X:]E+00
0.00£+00
0.00E+00
0.00E+00
O.(X3E:+00
0.00F.+O0
0.0(F.+00
0.(X]E+O0
0.00£:.,00
0.00£+00
0.00E+O0
0._
0.00E:+00
0.00i;+00
0._
0.00£+00
0._
0._+00
O.0CIi;_
0.(X_E+(X)
O.00E+O0
0.(X_E:+00
0.00E+00
0.0(F.+O0
0.0(]_:+00
0._+00
0.001_
1.52£-02
1.01£_
0.001;+00
0._+O0
0.00E;+O0
0._
0.00£_
0.0CI£.,00
O.00E:+O0
O.00E:_
0._+00
10.5
O.OOE+O0
0.00E+O0
O.00E,J00
O.00E+O0
0.0CE+O0
0.00E+O0
O.00E+O0
O.00E+O0
4.25£.-0,.1
2.8,_..,03
O.OOE+O0
O.00E+O0
O.00E+O0
O.OOE,K]O
O.00E+O0
O.00E+O0
O.00E+O0
O.00E+00
O.001E+O0
10.75
0.00E+O0
0.00£+00
0.00Et._
0.C0£+00
0.00£_
0.00£+00
0.0(O:+O)
0.00E;+O0
0.00(+00
0.00E+00
0.00E+00
O.0(]F.+00
0.0(0:400
0.00£+00
0._+00
0.0CI£,4]0
0._
0.00E+,00
0._,4_0
11
0.0C£+O0
0.00£+00
O.00E+O0
O.0OE+O0
O.00E+O0
O.00E+00
O.00E+O0
O.00E+O0
O,00E+OO
1,67£-O2
4.18£--03
0.00£+00
O.OOE+O0
O.OOE+O0
O.OOE+O0
O.00E+O0
O.00E+O0
0.00E+O0
0.00¢_0
O.00E+O0
O.00E+O0
O.00E+O0
0.00£+00
O.OOE+O0
0.OOE+O0
0.00E+00
O.00EH00
0.00£+00
7.55£-03
1.89E-0]
0.OOE+O0
0.00£+00
0.00£+00
O.00E,_I0
O.OOE+O0
0.0CIE+00
0.001EH]O
Oo00E+00
11.5
0.00E+_
0.00E+00
O.00E_
0.00E+O0
0.CO[:+O0
0.00[:_
O.CO[:_
0.0Q£_
0.0(]£+00
0.00£:+00
0.00E+00
O.0CE_
0.00E+00
0._
O.001E_
0.00E+00
0.(X)£_
O.00E_
0._
11.75
0.00£+00
0.0(](_,CO
0.00F..+00
0.00E+(X_
0.00E+00
0.(X)E+00
O.(X)E+00
0._
0.C0£_00
1.84£'-O8
2.41E_
0.00E:+00
O.00E;+O0
0.00E+00
0.00£+00
O._H_0
0.(X]EH]O
O.(X_£+00
0.00£_4]0
0.OOE+O0
0.00EHJ0
0.CX]E+00
0.00£+00
0.00£+00
0.OC(+O0
0.00£_,00
0.00E+O0
0.00£+00
4.27£.-09
5.60E_
O.00E+CO
0.00E+00
0.00£+00
0.(X)£+00
O.00E+00
O.0C]E+00
0.00£+00
0._,O0
0.00m_+00
0.(13(+00
0.0(O:+(x)
0.(X_I£+00
0.0CI_H]0
0.0CE_
0.00E+00
0.00E+00
0.0(0_+00
0.(x_£+00
9.46£-.04
3.7g£_
0.0()E+00
0.00£HJO
0.(x31£_
O.(X]EH)0
0.0CE+00
O.0CF.+O)
O.(X)E+O0
0.00E+O0
0.00E_
0.00E+00
0._+00
0.0(][_+00
0.00E+_
O.00E_
0._
0.0()£_
0.00£+00
4.81E_
1.92£-02
O.00E+00
0.00E_
0._
0.0CE+00
0.0CIE_
0.00E_
0.00E_
0._+O0
0.00[:+00
0.00E_
0.00E+00
0.00£+00
0.0CIE_
O.(]Ci£_
0._
O.(X]F..+00 0.00£_
0.00£+00
O.00E+00
0.00E+_
0.0()E;_
0.0(]IE_
0._
0.00E_
0.00E+_
0.00E_
0.(X]£',-00 O.COE+00
0.0C(+00
0.0(](_
0.0QE+00
O.(X]E+O0
0.001_+(X)
0.00E+O0
0.(X_E:+O0
0.00E+00
0.(X]E+(X)
?.6_3
1.14E-O2
O.COE+O0
O.(X_£+O0
O.0CE_
0.00F.+00
0.O0£:+OD
0.00£_q00
O.0CIE+OO
0.00£+O0
O.0OE+QO
O.OOE+O0
O.00£+O0
O.OOE+OO
O.COE+OO
O.O01EHJO
O.OOE+O0
O.OOE+O0
O.O0£+O0
2.4OE--O3
3.61E-O3
O.00£+OO
O.OOE+OO
O.OOE+OO
0.00£+OO
O.OOE+OO
O.OCIE+OO
13.5
O.O0£+O0
O.OOE+OO
O.CO£+O0
O.OOE+O0
O.OCE+O0
O.OCIE+O0
O.O0£+CO
O.00E+O0
O.OOE+O0
O.O0£_00
O.OOE+OO
O.00£+QO
6.51£-O3
4.34E.-O3
O.OOE+O0
O.00£+OO
O.OOEH_O
O.0OE+OO
O.0C£+O0
13.75
0.(X]C+00
0.00F_.+00
O.(X]E+O0
0.(X]I_+O0
0.(X)E+00
0.00£+(X)
0.(X_£+00
0.00E+00
0.(131_+00
0.00E+00
0.00E+00
0.CC£+00
4.I_£_
2.'#7£_3
0.00E+00
0.0011_+00
0.0C_+00
0.(X)E+00
0.00£_,4_0
O.OOE+OO
O.OOE+CO
O.00£+O0
O.OOE+O0
O.COE+O0
O.O0£+O0
O.OOE+O0
O.OOE+OO
O.O0£+00
O.00_+O0
O.0OEH00
O.OOE+OO
O.00£+O0
4.94E-O3
1.23E-O_
O.O0£H00
O.OOE+OO
O.OOE+OO
O.001£+400
O.COE',OO
O.COE+O0
O.OOE+OO
O.OOE+O0
0.00_+00
O.OOE+OO
O.OOE+O0
O.00£+00
0.00£+O0
O.OOE+OO
O.00£+O0
O.OOE+O0
O.OOE+OO
1.87E-O1
5.801[:-OI
O.OOE_O
O.OOE+O0
O.COE+O0
O.OCIE+O0
0.(X3mE+O0
0.(]OE+O0
O.OOE+O0
0.00£+00
O.COE+O0
0.00£+O0
0.O0£+O0
O.OOE+(]O
0.001_+00
O.O01£+400
0.00£+O0
O.OOE+OO
O.OOE+OO
1.331£-O9
8.71£,..O4
O.0CIE+O0
O.O0£+OO
O.OOE+OO
O.O0£_O
O.OOE+CO
O.OOE+O0
O.OOE+O0
O.OOE+CO
O.OOE+OO
O.O0£+O0
O.OOE+O0
O.OOE+O0
O.OOE+O0
O.0OE+O0
O.OOE+O0
O.0CIE+O0
O.00_+O0
O.00£+OO
O.OCIE+00
O.OOE+O0
O.CCE+O0
O.OOE+O0
O.O0£+O0
O.OOE+OO
O.COE+O0
O.0OE+OO
O.OOE+OO
O.001E+O0
O.O0£+OO
O.OOE+O0
O.O0£+00
O.OOE+O0
O.OO1_+O0
0.00£+O0
O.OOE+O0
O.OOE+O0
O.OCIE,,.OO
O.O0£+OO
O.OCIE+OO
O.OOE+00
O.OOE+(]O
O.OOE+O0
15.25
O.OOE+OO
O.OQE+O0
O.0C£+OO
O.0C£+O0
O.00£+OO
O.OOE+OO
O.OOE+O0
0.00£+400
O.OOE+OO
O.OC£+00
O.O0£+O0
O.OOE+O0
O.OOE+OO
O.OOE+00
O.OOE+O0
O.OOE+O0
O.OOE+OO
O.O(F,.+O0
O.0OE+O0
15.5
O.OOE+O0
O.OOE_O0
0.00£+00
O.OOE+O0
O.OOE+O0
O.OOE+O0
O.OOE_O0
O.OCIE+CO
0.00£+00
0.00£+00
O.OOE+O0
O.OOE+O0
O.OOE+O0
O.OOE+O0
O.OOE+O0
O.OOE+O0
O.OOE+O0
O.OOE+O0
O.OOE+O0
0.00£+00
0.00£+00
0.0(0_+00
0.00E:+(X)
0.(X_E+00
O.00E+O0
0.00£+00
O.00E:+00
0.001_+00
0.00£+00
0.(X_¢+O0
0._O0
0.00E:+O0
0.0CE:+00
0._+00
O.0()E:H_0
0.(X]E:+00
0.00£+00
0.00£,,