,"" 0F

«*,

1«

CO NBS TECHNICAL NOTE 903

°*EAU of

U.S. DEPARTMENT OF COMMERCE/ National Bureau of Standards

The NIRA Computer

Program Package

otonuclear Data Center)
NATIONAL BUREAU OF STANDARDS

The National Bureau of Standards 1 was established by an act of Congress March 3, 1901. The Bureau's overall goal is to strengthen and advance the Nation's science and technology

and facilitate their effective application for public benefit. To this end, the Bureau conducts

research and provides: (1) a basis for the Nation's physical measurement system, (2) scientific and technological services for industry and government, (3) a technical basis for equity in trade, and (4) technical services to promote public safety. The Bureau consists of the Institute for Basic Standards, the Institute for Materials Research, the Institute for Applied Technology,

the Institute for Computer Sciences and Technology, and the Office for Information Programs.

THE INSTITUTE FOR BASIC STANDARDS provides the central basis within the United

States of a complete and consistent system of physical measurement; coordinates that system

with measurement systems of other nations; and furnishes essential services leading to accurate and uniform physical measurements throughout the Nation's scientific community, industry,

and commerce. The Institute consists of the Office of Measurement Services, the Office of Radiation Measurement and the following Center and divisions:

Applied Mathematics — Electricity — Mechanics — Heat — Optical Physics — Center

for Radiation Research: Nuclear Sciences; Applied Radiation — Laboratory Astrophysics 2 — Cryogenics - — Electromagnetics 2 — Time and Frequency 2

. THE INSTITUTE FOR MATERIALS RESEARCH conducts materials research leading to improved methods of measurement, standards, and data on the properties of well-characterized

materials needed by industry, commerce, educational institutions, and Government; provides

advisory and research services to other Government agencies; and develops, produces, and

distributes standard reference materials. The Institute consists of the Office of Standard

Reference Materials, the Office of Air and Water Measurement, and the following divisions:

Analytical Chemistry — Polymers — Metallurgy — Inorganic Materials — Reactor

Radiation — Physical Chemistry.

THE INSTITUTE FOR APPLIED TECHNOLOGY provides technical services to promote

the use of available technology and to facilitate technological innovation in industry and Government; cooperates with public and private organizations leading to the development of

technological standards (including mandatory safety standards), codes and methods of test; and provides technical advice and services to Government agencies upon request. The Insti- tute consists of the following divisions and Centers: Standards Application and Analysis — Electronic Technology — Center for Consumer

Product Technology: Product Systems Analysis; Product Engineering — Center for Building

Technology: Structures, Materials, and Life Safety; Building Environment; Technical Evaluation and Application — Center for Fire Research: Fire Science; Fire Safety Engineering.

THE INSTITUTE FOR COMPUTER SCIENCES AND TECHNOLOGY conducts research and provides technical services designed to aid Government agencies in improving cost effec- tiveness in the conduct of their programs through the selection, acquisition, and effective utilization of automatic data processing equipment; and serves as the principal focus within

the executive branch for the development of Federal standards for automatic data processing

equipment, techniques, and computer languages. The Institute consists of the following

divisions: Computer Services — Systems and Software — Computer Systems Engineering — Information Technology.

THE OFFICE FOR INFORMATION PROGRAMS promotes optimum dissemination and

accessibility of scientific information generated within NBS and other agencies of the Federal Government; promotes the development of the National Standard Reference Data System and

a system of information analysis centers dealing with the broader aspects of the National Measurement System; provides appropriate services to ensure that the NBS staff has optimum

accessibility to the scientific information of the world. The Office consists of the following

organizational units: Office of Standard Reference Data — Office of Information Activities — Office of Technical

Publications — Library — Office of International Relations — Office of International Standards.

1 Headquarters and Laboratories at Gaithersburg, Maryland, unless otherwise noted; mailing address Washington, D.C. 20234.

2 Located at Boulder, Colorado 80302.
°»STANDARr,, LIBRAE Y

' rfB 1 ? 7976 The NIRA Computer Program Package ^^c

(Photonuclear Data Center) .u*t*3

~

" no,

a/03

J hok- no, SOii ; 976

-4-, -r ^ I

H. J. Vander Molen*

and

Henry M. Gerstenberg

Nuclear Sciences Division

Institute for Basic Standards

National Bureau of Standards

Washington, D.C. 20234

Work supported by the National Science Foundation through the National Academy of Sciences-National Research Council, Committee on Nuclear Science.

*A Nuclear Information Research Associate (NIRA) at the National Bureau of Standards; present address: U.S. Nuclear

Regulatory Commission, Washington, D.C. 20555. *< °F *><

U.S. DEPARTMENT OF COMMERCE, Elliot L Richardson, Secretary Designate

James A. Baker, III, Under Secretary

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

NATIONAL BUREAU OF STANDARDS Ernest Ambler, Acting Director

Issued February 1976
National Bureau of Standards Technical Note 903

Nat. Bur. Stand. (U.S.), Tech. Note 903, 30 pages (Feb. 1976)

CODEN: NBTNAE

U.S. GOVERNMENT PRINTING OFFICE

WASHINGTON: 1976

For sale by the Superintendent of Documents, U.S. Government Printing Office, Washington. D.C. 20402 (Order by SD Catalog No. 03.46:903). Price 85 cents (Add 25 percent additional for other than U.S. mailing).
Preface

The Photo-nuclear Data Center, a part of the Center for Radiation Research, Is one of many data centers within the National Standard Reference Data System receiving support from the National Bureau of Standards' Office of Standard Reference Data, Its function is to collect, abstract, index, and evaluate experimental data published in

the world's technical literature relating to the interaction of electromagnetic radiation with nuclei. It provides information-center services for the technological and

scientific communities concerned with photonuclear reactions. Ever since its creation

in 1963, the Data Center has made considerable use of computer-compatible index,

bibliographic, and data files. These are used with the central computer facilities at the National Bureau of Standards to retrieve specific types of information; to prepare

camera-ready copy for publication of tabular, bibliographical and graphical information;

and, to carry out various manipulations of digital data sets germane to their evaluation.

This Technical Note is the first of several being written to describe all aspects

of the computer-based operations of the Photonuclear Data Center. These notes are meant

both to document the procedures used in the Data Center as well as to serve as 'user's manuals' for individuals having a need to make use of the Data Center's files. The

program library described in this Note has been designed for use with the particular

central computer facilities available at the National Bureau of Standards. No attempt

has been made to write the various programs in a completely general form (facilityindependent). For those desiring to have it, the complete program library (in symbolic

form) is available from the Photonuclear Data Center.

Everett G. Fuller

III

Contents

Page

1. Introduction 1

2. Description of the Library 1

3. Use of the Program Library 2

4. Qualitative Description of the Programs and Their Use „ 2

4.1. Programs Which Create Data Sets 2

4.2. Programs Which "Dress-Up" Data 3 4.3. Programs Which Manipulate Data „ 3 4.4. "Housekeeping" Programs ..... 3 4.5. Programs Which Produce Pinal Copy 3 4.6. Use of Programs 3

5. Some Notes on Efficient Use 4

5.1. Multiple Execution 4

5.2. Temporary File 4

5.3- Breakpoint Pile 5

6. Appendix; NIRA Listings 7 6.1. Programs 7 AVRG 7 CALC 7 CLEAN 7 COMBIN 8 CONCAT 9 DBAL 9 GDR 10

HIST 10

INTERP 10

MONO 11

SMOOTH 11 SORT 11 SQUARE 12

STIFF 12

TRUNC 14 UNHIST 14 UNSTAT 14 YIELD4 15

6.2. Control Decks 16

DIGIT 16

LIST 18

PPLOT 18

A. Program Calling and Post-Processors 19

B. Plot Control Cards 19

C. Parameters and Defaults 20

D. Examples 22

SPCDST 24

The NIRA Computer Program Package

(Photonuclear Data Center)

H. J. Vander Molen and Henry M. Gerstenberg

The Photonuclear Data Center's NIRA library of programs, executable

from mass storage on the National Bureau of Standards' central computer

facility, is described. Detailed instructions are given (with examples),

for the use of the library to analyze, evaluate, synthesize, and produce

for publication, camera-ready tabular and graphical presentations of digital photonuclear reaction cross-section data. NIRA is the acronym for

Nuclear Information Research Associate.

Introduction

There are two aspects of the Photonuclear Data Center's operations which make

use of the NBS computer facilities. The

first is related to the abstracting and indexing of experimental papers published in

the world's technical literature. This

provides the input both for the Center's

computer-searchable index to the available

experimental data for the field of photonuclear physics, as well as the annotated

bibliographies or data indexes which the

Center regularly publishes and updates [l]. 1

The second computer-oriented activity

of the Center relates to its library of

digital cross section data and the systematic evaluations being made of these data.

To do this, it has been necessary to generate procedures for handling, storing,

classifying, evaluating, and presenting

data. The program library described here

is an implementation of these objectives on

the National Bureau of Standards' main computer facility, a UNIVAC 1108 computer operating under the EXEC-8 system. 2

2. Description of the Library

The NIRA program library consists

primarily of pr Dgrams which handle data ,

rather than index or bibliographic information. Most of these programs are quite

simple, both conceptually and in implementation. All can be executed from a demand

terminal, although ordinary batch usage was

the primary intention. All are mutually

compatible; data processed by any program

can be read by any other program without

human manipulation. It is likely that, in

the future, need will arise for even more

programs. However, the basic operations

described here should suffice for most

procedures

.

The programs reside, as addressable

elements, in a file called "NIRA*PROGM.".

(The asterisk and period are essential parts

of the name.) This file is stored as a "read-only" file. It is therefore not possible to alter it by any error in normal use.

(It is possible to delete it, however, so be

cautious when using DELETE statements.) Two or more runs may simultaneously access this

file. As a backup, the entire program

library (in symbolic form), as well as a

complete listing, is kept in the Photonuclear Data Center. At present, this library

deck consists of more than 3,000 punched

cards with control cards in place for re- loading to mass storage.

Most of the individual programs are designed to operate on a deck of punched

cards which define a "data set". Similarly,

most of the programs output either real or

figures in brackets indicate the literature references given on p. 6.

2Use of brand names in this report is for identification purposes only and is not to be

construed as a recommendation of a particular product or company.
virtual punched cards. In all such cases, the program adds a comment card to the deck

which states what the program has done to

the data. This card is listed on the user's

printout for the run In which the program Is called. For reasons of economy (I.e., to

limit the amount of the computer's memory

used) , most of the programs will accept a maximum of 500 data points and 100 comment

cards. This can be changed without great

difficulty, but this quantity of memory has

been more than sufficient thus far.

As defined by the Photonuclear Data

Center, a "data set" consists of the results

of the measurement of an excitation function, angular distribution, or reaction

product spectrum under specified conditions

for a single target. The set consists of a Header Card which contains the acquisition

number, number of data points, reference,

nucleus, and reaction, plus some data source Information; two Label Cards which give the

axis labels, usually "Photon Energy (MeV)"

and "Cross Section (mb)"; an arbitrary num- ber of Comment Cards ; and then Data Cards

in Fortran format (7(E9. V, ') ,E9.4,1X)

.

A few programs require more elaborate

control cards than this MAP-IN-XQT sequence.

These have been combined with the appropriate individual programs into what in the

following discussion are referred to as "job-control decks". Each deck is named,

and only one card need be entered to call a

deck; e.g., PPLOT is called by:

@ADD,P PF. PPLOT

The individual writeups give the exact

cards to use. (Do not forget the ASG-USE

sequence at the beginning of your run, however ! ) Note that the run cards , tape

identification labels, etc., that are used

in a number of the examples, are given for

purposes of illustration only, and that In making use of the library, the user should

be sure to use the appropriate number and

designation associated with his computer

account. Throughout the listings, square

brackets, i.e., [ ], are used to designate

cards containing information germane to

particular runs only. Brackets are not used in entering information unless specifically called for in the program description.

3. Use of the Program Library

To use the library, it is necessary to

assign the library file to a particular

computer name. To do this, use the two

EXEC~B control cards [2]:

@ASG,A NIRA*PROGM.

@USE PF.,NIRA*PROGM.

The "commercial at" sign is represented by

a 7/8 punch on the card punch, and by "@"

on the Teletype.

As indicated in the individual descriptions given in the Appendix, most of the

programs are called by entering three cards

For example: SORT is called by:

@MAP,INX MAPR,S0RT

IN PF.SORT

@XQT SORT

Similarly, SMOOTH is called by:

@MAP,INX MAPR, SMOOTH

IN PF. SMOOTH

@XQT SMOOTH

Qualitative Description of the Programs

and Their Use

A brief description of the purpose of

the various programs is given here. The

individual descriptions should be consulted

before they are used.

4.1. Programs which Create Data Sets

DIGIT is a job-control deck that takes

data from a tape made by the sonic delayline digitizer, converts it to data-library

format, operates on it with CLEAN, provides

punched cards for the data library, and/or

leaves the data in a computer file element

called DIGOUT for use by other programs.

HIST converts a normal data set into a histogram.

UNHIST converts a histogram into a normal data set

.
4.2. Programs Which "Dress-Up" Data

SMOOTH performs a normal vertical

smoothing on raw data.

STIFF is used to "dress-up" digital

data so that it may be used as input to obtain computer-drawn, smooth, continuous

cross section curves.

AVRG smooths raw data over a pre-set

energy interval.

INTERP parabolically interpolates

smoothed data to give points at "nice"

energies

.

SQUARE squares up histograms.

4.3. Programs Which Manipulate Data

COMBIN algebraically adds two cross

sections. It will read a (y,n) and a

(Yj2n) cross section, and punch a (y>xn) cross section, for example.

DBAL reverses a cross section by means of the principle of detailed balance.

CALC renormalizes data sets.

UNSTAT removes the statistical theory's

multiplicity correction from a (Y 5 sn) cross

section.

4.4. "Housekeeping" Programs

CLEAN checks data for format and makes

several computer-oriented (non-manipulative)

corrections. A data set which has passed

through CLEAN should give no problems when

read by the other programs.

MONO forces energies to increase

monotonically

.

SORT sorts data points by energy.

TRUNC picks out a subset of a data set.

CONCAT joins together several data

sets covering different energy ranges.

4.5. Programs Which Produce Final Copy

YTELD4 is used to print a column-wise

table of cross sections. It is also used

to calculate c_, at 10, 15, 20, 25, and

30 MeV.

x

SPCDST is a deck used to print an index table of spectral and angular distribution data.

LIST is used to produce camera-ready

copy of tables made by either YIELD4 or

SPCDST.

PPLOT is a plotting job-control deck.

It produces plots which are suitable for

mounting on reduction paper.

GDR calculates the integral moments a , a_i, and o_2 of the cross section, and estimates the giant resonance parameters E , a , and r. See Ref. [1] for definitions. m m

4.6. Use of Programs

The use of the programs is completely

flexible. Nevertheless, there is a "normal" or "usual" sequence of steps through which

data pass from journal article to final

copy.

(1) The data are digitized and

assigned an acquisition number. The data

set is then run through DIGIT. If the data

already exist in digital form, they are punched and run through CLEAN. UNHIST may

also be used if the data are given as a histogram.

(2) If CLEAN indicates monotonicity

problems with the data, SORT and/or MONO

may be used to correct the trouble.

(3) The data set is manipulated, if

appropriate, to get it into the proper units

and form for comparison. COMBIN, DBAL,

UNSTAT, TRUNC, and CALC are provided for

such manipulation.

(4) The data are plotted, using the

CPLOT mode of PPLOT. Comparisons are made,

and the best data are selected.

(5) The selected data, if a (y.sn),

(y,xn), or (y,y, ), are run through GDR to calculate a„ , a

' , and a_„

.

also gives estimated values for E

n and T.

This program

(6) The selected data are used to create a smoothed cross section curve. This

smooth curve is the source of the evaluated

data set. The curve is usually drawn by

hand. SMOOTH and AVRG are sometimes used

in conjunction with, or even in place of,

the manual drawing. (SQUARE is used to dress up histograms . ) If the manual method

is used, the hand-drawn line is digitized,
run through DIGIT and MONO, and then run

through STIFF to smooth out the hand-induced

irregularities. If the hand-drawn line was digitized in sections, CONCAT is used to join the sections together before STIFF is

used.

(7) The evaluated data are then run

through INTERP to obtain a compact representation with cross-section values given at

specified energy intervals. Occasionally,

MONO is needed between STIFF and INTERP. The output of INTERP is filed as the final

form of the "evaluated" cross section.

(8) The yields are calculated from the

selected (not evaluated) data using the

YIELD aspects of YIELD4.

(9) The evaluated cross sections, as well as the selected data set, are plotted

using PPLOT.

(10) The evaluated cross sections are tabulated using YIELDS and LIST."

(11) The index to supplementary data

is tabulated using SPCDST.

5. Some Notes on Efficient Use

Almost all of the programs described

in the preceding pages read cards for input

and punch cards for output. It is not necessary to execute just one program per

computer run. There are several "tricks"

which enable the user to perform a sizable

amount of work in one computer run.

5.1. Multiple Execution

Most of the programs are called in the

following manner:

@MAP,INX MAPR,GDR

IN PF.GDR

@XQT GDR

The first two cards direct the MAP

processor to transform the relocatable

element PF.GDR into the absolute (executable) element GDR. Even after execution,

the absolute element GDR is still in the

run's temporary program file (TPF$.). It can be re-executed with another XQT card;.

@MAP,INX MAPR,GDR

IN PF.GDR

@XQT GDR

[1st data set] i @XQT GDR

[2nd data set]

@QXT GDR

[3rd data set]

etc.

5.2. Temporary File

A temporary file is useful when the

same data set is needed for input to more than one program. The file is created by:

@ASG,T TF.

Data are then put into the file element by

element . The file elements are fed to the

processing programs by ADD.

Example : Suppose one started out with a

(y 3 n+np) and a (y,2n) data set. . It is

desired to calculate a (y^sn) and a (Y,xn)'

data set. The cards would be:

@ASG,T TF.

<§ELT,IS TF.l

[(y,n+np) data] ., @ELT,IS TF.2

[(y,2n) data]

@MAP,INX MAPR,COMBIN

IN PF.COMBIN

@XQT COMBIN

0,0,1. ,1.,

@ADD,P TF.l

@ADD,P TF.2

@XQT COMBIN

0,0,1. ,2.,

@ADD,P TF.l

@ADD,P TF.2
5.3- Breakpoint File

A breakpoint file is useful when it is

desired to feed the output of one program

into the input of another program. The file

is created by:

@ASG,T BP.

The punched output of a program is diverted

to the breakpoint file by:

@BRKPT PUNCH$/BF (Note no punctuation!)

Then the first program is executed. After

the last data card of the first program, the

punch channel is re-activated by:

SBRKPT PUNCH$

Now the contents of the breakpoint file can

be fed into the 2nd program by means of:

@ADD,P BF.

Example : Suppose it is desired to calculate a (y,xn) total neutron yield. The

starting data are a (Y>sn) data set

measured via a BF3 detector and a statistical multiplicity correction. It is necessary to first call UNSTAT to transform the

(YjSn) data set into a (y,xn) data set,

then to call YIELDM to calculate the yield.

@ASG,T BF. @MAP,INX MAPR, UNSTAT

IN PF. UNSTAT

ffiRKPT PUNCH$/BF

@XQT UNSTAT

[parameter card]

[(YjSn) data]

@BRKPT PUNCH$

@MAP,INX MAPR, YIELDS

IN PF.YIELD4

@XQT YIELD4

(NONE) YIELD

@ADD,P BF.

The breakpoint file is erased whenever

new gBRKPT PUNCH$/BF card is encountered,

iesides adding data with an @ADD statement,

he contents of the breakpoint file can be

isted by entering:

@ADD,P PF.LIST

as can be seen in the description of SPCDST.

The breakpoint file can also be punched by:

@ELT,ID PUNCH

6ADD,D BF. @END

@PCH,S PUNCH

(The SYM processor of EXEC-8 will also punch

BF, but requires cataloging. This had led

to interference between runs .

)

The breakpoint file is an extremely

powerful tool when used in conjunction with

a temporary file. The following example

begins with a (p,Y ) data set. The data set is inverted and plotted. Simultaneously, a smoothed form of the same cross section is

digitized, smoothed, punched, inverted,

plotted, and tabulated.

@RUN,M/NR VAND41,30850-VANDER,NIRA, 1,50/1000

@ASG,A NIRA*PROGM.

@USE PF.,NIRA*PROGM.

@ASG,T BF. @ASG,T TF. @MAP,INX MAPR,DBAL

IN PF.DBAL

gBRKPT PUNCH$/BF

@XQT DBAL

[parameter card]

[(p,YQ

) data]

(3BRKPT PUNCH$

@ELT,IS TF.l

@ADD,P BF. @ELT,IS DIGIN

[digitizer control cards]

@ASG,LJ0 IN.,8C,VAN02

@SETC 1 @ADD,P PF. DIGIT

@MAP,INX MAPR, MONO

IN PF.M3NO

@BRKPT PUNCH$/BF

@XQT MONO

@ADD,P DIGOUT

@BRKPT PUNCH$

@ELT,IS TF.2
@ADD,P BF. @MAP,INX MAPR, STIFF

IN PF. STIFF

gBRKPT PUNCH$/BF

@XQT STIFF

[parameter card]

@ADD,P TF.2

@BRKPT PUNCH$

@ELT,IS TF.2

@ADD,P BF. @BRKPT PUNCH$/BF

@XQT DBAL

[parameter card]

@ADD,P TF.2

@BRKPT PUNCH$

@ELT,IS TF.2

@ADD,P BF. @MAP,INX MAPR,INTERP

IN PF.INTERP

(9BRKPT PUNCH$/BF

@XQT INTERP

@ADD,P TF.2

@BRKPT PUNCH$

@ELT,IS TF.2

@ADD,P BF. @PCH,S TF.2

@MAP,INX MAPR,YIELD4

IN PF.YIELD4

(2BRKPT PUNCH$/BF

@XQT YIELD4

[blank card]

@ADD,P TF.2

gBRKPT PUNCH$

@ADD,P PF.LIST

@ADD,P PF.PPLOT

@ADD,P TF.l

$CURVE

@ADD,P TF.2

$PPLOT

@ADD,P PF.PRINTR

@MSG,W 30850-VANDER,WRITE ENABLE REEL VAN03

@ASG,TMI 9 . ,8C,VAN03W

@CLOSE 9. @XQT NBS*CALCMP.

@FIN

References

[1] E. G. Fuller, H. M. Gerstenberg, H. Vander Molen, and T. C. Dunn, Photonuclear

Reaction Data, 1973 , Nat. Bur. Stand. (U.S.), Spec. Publ. 380, 125 panes

(March 1973).

[2] UNIVAC 1100 Series, EXEC-8 Hardware/Software Summary UP-7824 Rev. 1 (1971-1974) Sperry Rand Corporation, New York, N. Y.
AVRG

6. Appendix; NIRA Listings

6.1. Programs

AVRG is a means of smoothing a data set. The cross section at each energy, E. is

replaced by the average value over a specified energy interval, e. (This is in

contrast to SMOOTH which averages over 3 points regardless of spacing.) The averaged data

set is then punched. The formula used is:

i A+ £/2

a(E.) = i / x a(E)dE . l e rJ E.- e/2

The integration is rectangular. The averaging interval "e" is constant except at the

low energy end. If E.- E < e/2, e = 2(E.- E ).

The averaging interval "e" is punched as a floating point number in free form,

followed by a comma, on a parameter card immediately preceding the data set. If the card

is omitted, e = 1 MeV. Any number of data sets may be entered. The value given to

epsilon remains in force until a new parameter card is encountered.

@MAP,INX MAPR,AVRG

IN PF.AVRG

@XQT AVRG

[parameter card]

[data set]

CALC

CALC manipulates a data set by the formulas: E.= A E.+ A and a.= B a + B„ . The

first card read by CALC must contain A,, A„, B, , B

?

, and the new acquisition number. If

zero is entered for the new acquisition number, the old acquisition number will be used.

The program will read only one data set

.

Example : If it is desired to change a data set scale from millibarns

to microbarns, the parameter card would be: 1. ,0. ,1000. ,0. ,0

@MAP,INX MAPR,CALC

IN PP. CALC

@XQT CALC

[parameter card]

[data set]

CLEAN

CLEAN reads a data set and checks it for format. It will automatically modify the

set as follows:
(1) If the label cards are missing, "PHOTON ENERGY (MEV)" and "CROSS SECTION (MB)" are inserted.

(2) If the number of data points is missing or incorrect, the correct number is

inserted.

(3) If the header card is missing, one is inserted.

(4) If two adjacent data points have the same energy, the second energy is increased

one part per thousand.

(5) If the data are non-monotonlc, except for the case above, a warning is printed.

(6) If a blank card is encountered, it is skipped over.

(7) If two points have identical x- and y-coordlnates (e.g., due to double-sparking

of digitizer), the second point of the doublet is deleted.

This routine can read any "free" format. A warning Is printed if the data are not in the

standard "(7(E9.4,', ' ),E9.4,1X)" Data Center format. The monotonicity and double-sparking

tests are bypassed if the acquisition number is negative (i.e., special test data).

A data set punched by CLEAN can be handled by all the other programs without problems.

@MAP,INX MAPR, CLEAN

IN PF. CLEAN

@XQT CLEAN

COMBIN

COMBIN is used to algebraically add two cross sections and punch a new data set. The

formula is:

c' = aa + 6a, , a b

'

where a and 8 can be any positive or negative floating point number. The two data sets must be preceded by a parameter card containing NA, NB, a, and 8 in free field format. NA

and NB must be fixed-point positive or negative integers, and a and 8 are floating-point

numbers. The maximum energy in the new data set will be the larger of that in the two sets

combined if NA = 0. Setting NA = -1 results in the maximum energy being the lesser of that

in the two sets. In both cases, NB is the acquisition number for the new data set.

Only two data sets may be entered. The set with the best resolution should be the

second data set. If the resolution differs greatly between the two sets, then the programs

AVRG or SMDOTH should be used to obtain comparable resolution before using COMBIN.

The program makes the resolution of output data similar to that of the first data set. For the region where the data sets overlap in energy, the second set is linearly interpolated to match exactly the energies of the first data set. For the regions where there is no overlap and data exist only in the second set, then the resolution of the first set is

averaged over ten data points nearest the region of interest. This average energy

resolution is then used in the linear interpolation of the second data set to obtain output

data points of the second data set having a resolution similar to that of the first data

set.

@MAP , LNX MAPR, COMBIN

IN PF. COMBIN

@XQT COMBIN

[parameter card]

[data set A]

[data set B]
CONCAT

This program reads an arbitrary number of data sets, and punches one data set which

is a simple concatenation (i.e., joining together) of the input data sets. The header and

label cards of the output are, except for the number of points, copied from those of the

first data set read. The total number of points must not exceed 500. Generally, it is

wise to feed the output of CONCAT into SORT.

@MAP , INX MAPR, CONCAT

IN PF. CONCAT

@XQT CONCAT

[data sets]

DBAL

DBAL reverses a cross section via the principle of detailed balance. It reads a parameter card and a data set.

For reactions of the type A+y rf B+b, the equations used are:

a(Y,b

Q

) = f(A,E

e

) o(b,Y ) where,

The parameters are:

(2J +l)(2s+l) m Mg (E -BE)

f(A,E ) = 7

„

T -= t .. , with

e (2JA+1) Mg+m £2

E = =5-2- E,+ BE or E = E - * — e Mg+m b e Y 2(M 2+E ) A y

q = flag to indicate direction of change m = mass of particle b

s = spin of particle b

M. = mass of nucleus A

%

= mass of nucleus B

J

fl = spin of nucleus A

J„ = spin of nucleus B

ts BE = binding energy of b in A

To turn a ground-state capture cross section into the related photodisintegration

cross section, the parameter card is:

1 , m, s, Mg, J

B

, M

A

, J

A

, BE
Similarly, for the reverse transformation, the card should read:

0., m, s, Mg, J

B

, M

A

, JA, BE

Any number of data sets may be entered, but each must have a parameter card.

@MAP,INX MAPR,DBAL

IN PF.DBAL

@XQT DBAL

[parameter card]

[data set]

GDR

GDR reads a data set, calculates an0-1 , a ,, and a „, and estimates E_, a , and the -2' 0' max' width. Only one data set may be entered.

The integration is done by the trapezoidal rule. The upper limit of integration is

the greatest integer which is less than the max energy of the data set. The value of E_

is a calculated mean energy. The estimates of o and width are, in practice, almost

worthless. A better job can be done by eye.

The data set is repunched with the calculated giant resonance parameters given on a comment card.

@MAP,INX MPR,GDR

IN PF.GDR

@XQT GDR

[data set]

HIST

This program converts a normal data set into a histogram, vertical bars are placed

midway between every two points of the input data set. A slight inclination of the

verticals is Introduced to preserve monotonicity. Any number of data sets may be entered.

CAUTION - The output may look strange if the input is not approximately evenly spaced in

energy

.

@MAP,INX MAPR,RTST

IN PF.HIST

@XQT HIST

[data sets]

INTERP

INTERP reads a data set and punches a parabolic interpolation of the input . The

output, except for the first and last points, is given at every tenth, fifth, or half MeV,

depending on local curvature.

The output is suitable for input into PPIOT and YIELDS. Any number of data sets may

be entered.

10
@MAP,INX MAPR,INTERP

IN PF.INTERP

@XQT INTERP

[data sets]

MONO

1VDN0 is a program that corrects monotonicity problems which result from roundoff er- rors in going from the raw digitized data to the Data Center's data card format. The order

of the points is maintained, but energies are adjusted upwards such that they monotonically

increase. (This is not always the appropriate correction - see SORT.)

If more than one data set is entered, each one will be checked, corrected if necessary, and punched.

@MAP,INX MAPR, MONO

IN PF.MONO

@XQT MONO

[data sets]

SMOOTH

SMOOTH performs one pass of a normal vertical smoothing on a data set. Any number of

data sets may be entered. The energies do not need to be equally-spaced. The algorithm

used is:

1 °i

=

2

a a.

a. + _1+1 ~ i-J- (e. - E. ,) + a. , 1 E. ,, - E. , 1 l-l l-l

l+l l-l

In the special case of equally-spaced energies, this reduces to the familiar form:

1 1.1 . 1

o. =7-a. -, + -a. +T-O.J,

i 4 i-I 2 i 4 i+I

@MAP , INK MAPR, SMOOTH

IN PF. SMOOTH

@XQT SMOOTH

[data sets]

SORT

This program sorts the points of a data set by energy. Its primary uses are:

(1) removing monotonicity problems caused by digitizing data in a region of high-point

count density, and (2) fixing a deck of cards someone dropped on the floor. Note that

SORT and MONO have the same goal but accomplish it in different ways; SORT changes the

order of the points, while MONO changes the value of the points.

If more than one data set is entered, the program will sort each data set.

@MAP,INX MAPR, SORT

IN PF.SORT

@XQT SORT

[data sets]

11
SQUARE

This program is used to "dress up" a histogram-type bar graph. (We use "histogram"

to denote a bar graph, not in the strict statistical sense.) If the "corner points" of a histogram are digitized:

~L

T—

T

• •

f—

n '>

then, in general, the output of the digitizer program will give a graph in which the horizontal and vertical lines are not truly parallel to the axes.

SQUARE reads such a "raw" data set, and adjusts the points such that horizontal lines are truly horizontal. Vertical lines are inclined very slightly (one part per thousand)

to avoid monotonicity problems.

The height of each "squared" horizontal is the average of the heights of the two raw

points defining it, and similarly for the verticals. Any number of data sets may be

entered

.

@MAP , INX MAPR, SQUARE

IN PF. SQUARE

@XQT SQUARE

[data sets]

STIFF

STIFF performs one pass of a non-vertical smoothing routine. This means that energies

are also adjusted to smooth the data. Such smoothing is the appropriate treatment for

digitized curves where errors are both vertical and horizontal.

The program is used to smooth digital data which has been generated from "smooth" curves drawn by hand through actual data points. The resulting "smooth" data are used to

generate computer-drawn curves which approximate the actual data presented in graphical

form. They are also used to generate the tables of "smoothed" cross section data.

The smoothing is carried out by averaging over three successive data points, i, j, and

k, when their coordinates are expressed in inches. The program first uses points i and k

to determine a straight line. The coordinate axes are then rotated to establish an x'-axis

which is parallel to this line, and the average y' -coordinate of the point j is determined

from:

y'(j) = Jy'(i) +|y'(J) + Jy'OO.

12
The coordinates of the point j, x'(j), y'(j) are then referred back to the original un- rotated system to obtain the "smoothed" coordinates of the point j.

Since the program must work with the coordinates of each point expressed in actual

linear dimensions , the first card must contain XLNGTH and YLNGTH ( in inches ) , and the three

fiduciary points of the plot from which the data were digitized.

Example:

-i- 80

60

in 40

in

20

i o

20 30

6.0-

40

The parameter card is then:

6.,5.5,10.,0.,40.,0.,10.,*

More than one data set may be entered. If new parameter cards are not entered along

with the new data sets, the old parameters will remain in effect.

The primary use of STIFF is to "dress up" a smooth curve by reducing the small pointto-point irregularities introduced by hand-digitizing. The output is suitable for input

into INTERP.

@MAP,INX MAPR, STIFF

IN PF. STIFF

@XQT STIFF

[parameter card]

[data sets]

13
TRUNC

This program reads a data set and punches a truncated data set consisting of data

points which fall within a selected energy interval. The upper and lower limits are entered on a parameter card as floating point numbers in free format. The energy interval

is inclusive of the end points. Any number of data sets may be entered. If no new energy

limits are entered, the old limits remain in force. If no limits are entered [0.,30.,]

is assumed.

@MAP,INX MAPR, TRUNC

IN PF.TRUNC

@XQT TRUNC

[E, , E ] T_ower' upper

[data sets]

UNHIST

This program reads a histogram (either "squared" or "raw") and converts it into a normal data set.

Output points are placed in the center of the horizontal bars, except when multiple

bins are detected. No problem should occur if (as is almost always the case) the input

bins are of approximately equal width. Any number of data sets may be entered.

@MAP,INX MAPR, UNHIST

IN PF. UNHIST

@XQT UNHIST

[data sets]

UNSTAT

UNSTAT reads a (y 5 sn) cross section, and calculates a (y>xn) cross section via the

statistical model. Its primary use is to strip out the statistical multiplicity correction

used by many experimenters with Halpern detectors. The relations used, taken from

T. K. Deague, E. G. Muirhead, and B. M. Spicer, Nucl. Phys. A139 , 501 (1969), are:

a(sn) = a(y 5 n) + a(y3 np) + a(y 5 2n) +

a(xn) = a(y,n) + a(y,np) + 2a(y,2n) +

a(sn) = a(xn) F(E)

F(E) = 1+ (l-x)Cl- (l-f(E))exp-f(E)]

f(E) = (E-E2) (a/(E-El))

1/2

, where

x = direct de-excitation fraction

a = level density (in MsV-1

) El = the (y,n) threshold (in MeV)

E2 = the (y,2n) threshold (in NfeV)

14
The first card must contain floating-point numbers for x, a. El. and E2 in free format.

Any number of data sets may be entered. If new parameter cards are not entered along

with the new data sets, the old parameters will remain in effect.

@MAP,INX MAPR,UNSTAT

IN PF.UNSTAT

@XQT UNSTAT

[parameter card]

[data sets]

YIELDS

The program YIELDS is designed to operate on a number of data sets and calculate:

(1) a cross-section table with values for each of the reactions entered being given at evenly-spaced energy intervals (usually, every half MeV) , and/or (2) values of a_, for upper-energy integration limits of 10, 15, 20, 25, and 30 MeV for each of the

reactions entered. A linear interpolation is used to calculate values at specified energy

values. This program is a version of YIELD3 which operates on a single data set to produce a printout in which for each energy E. in the set the following quantities are given

E.l

a./E. a./E?

l l

i X>./e. p/E

j

The operations carried out by YIELDS are determined by entries punched on a parameter

card. These entries must be separated by at least one, and no more than five, blank spaces.

There must be no blanks within an entry, and commas must not be used to separate entries.

This parameter card must be provided even if no entries are punched on it. In its normal

state, the program calculates a cross-section table at half-MeV intervals. The following

entries may be made:

DELTA=XXXXXX

(G,XN)

62Ful

PUNCH

(NONE)

If an interval other than 0.5 MeV is desired,

replace the six X*s with a floating-point

number giving the interval. When the normal

DELTA is not used, this must be the first parameter on the card.

To skip a particular reaction in building the

table, enter the reaction with parentheses.

To skip a particular reference in building the

cross-section table, enter the reference number.

The output is listed in the normal Data Center

cross-section library format, i.e., Energy,

Cross Section, Energy, Cross Section etc. Note : "PUNCH" and "YIELD" are mutually exclusive.

No cross-section table is desired.

15
YIELD The yield table is calculated, i.e., c_, values

at 10, 15, 20, 25, and 30 MeV. Note : "PUNCH" and "YIELD" are mutually exclusive.

PRINT Gives a full printout of YIELD3 for each data set entered. See above discussion.

The cross-section table is punched in the proper format to be listed by LIST. See the

description of LIST for more details. Note : The cross-section table can hold a maximum

of 200 energies and 13 reactions.

@MAP,INX MAPR, YIELDS

IN PF. YIELDS

@XQT YIELDS

[parameter card] This card must be here, even if blank.

[data sets]

6.2. Control Decks

The following are decks of UNIVAC job-control statements each identified by a specific

name. These control decks perform various functions primarily related to the formating

of data for various final output modes (plots, tables, etc.). From the user's point of

view, they differ from the programs described above primarily in the way they are called

from mass storage and executed.

DIGIT

DIGIT is a deck of UNIVAC job-control statements, stored as mass storage, which:

(1) executes the digitizer program;

(2) takes the output of the digitizer program and inputs it into CLEAN;

(3) executes CLEAN;

(4) punches the output of CLEAN complete with headers and labels

( Caution : The labels are always "PHOTON ENERGY (MEV)" and

"CROSS SECTION (MB)");

(5) leaves the output in an element called DIGOUT for use in other

programs.

The use of DIGIT is best illustrated by a skeletal example:

@ELT,IS DIGIN

[digitizer control cards]

@ASG,LJ0 IN., 8C, [user's tape number]

@ADD,P PF. DIGIT

The data in the element DIGOUT may be accessed immediately in the same job by

the statement:

@ADD,P DIGOUT

If it is desired to use the data in the same job and not punch the output of CLEAN, precede

the @ADD,P PF. DIGIT card with:

SSETC 1

16
Two digitizer control cards must be given for each curve recorded on the tape. The

order of these pairs of cards must correspond to the order of the data on the tape. The

first card of each pair is a Header Card; the second a Plot Parameter Card. The Header

Card gives the acquisition number of the data set, the reference number, the nuclide, the

reaction, and any pertinent comments in the following format:

Position

1 - 5

6

7 -11

12

13 -18

19 -23

24 -30

31

32 -43

4445

52 -6

Quantity

Acquisition number, right adjust

Blank

Number of data points (to be provided

by piogram)

Blank

Data Center reference number

Blank

Nuclide/Element identification

25,25 Z, right adjust

26,27 Symbol, left adjust

28-30 A, left adjust

Blank

Reaction with ( ) , left adjust

Blank

Source flap;: blank if digitized at Data

Center; or, 1 if digital data was supplied

by original source Available for comments if needed

The Plot Parameter Card must contain nine numbers separated by commas. The first six,

in floating-point notation, give the coordinates of the three fiducial points used to establish the axis scales for the digitizer programs. The last three numbers in fixed-point

notation give the acquisition number for the plot, an x-axis flag, and a y-axis flag. For

a linear scale the flag is 0; setting the flag equal to 1 indicates that the corresponding

axis scale is logarithmic

.

Example:

h-O

LdCO

CO

COo

cro

200«

>.

, r i i 1 Vf 15

Plot Parameter Card entries:

35

ENERGY, MeV

15. ,0. ,35. ,0. ,15. ,200. ,932,0,0,

17
LIST

LIST is a combination of job-control statements (in PP. LIST) and a simple program

(in PP.LISTR) which lists data in the breakpoint file (BF.) on the slow printer. It is

designed to produce camera-ready copy. Each line is printed twice to ensure dark type.

Whenever a card with "$PAGE" (starting in column one) is encountered, the program skips

to the top of the next page in the printer. The $PAGE card is, of course, not printed.

Lines longer than 80 characters may be printed with this program. The first 80 characters of the long line are entered on an 80-column card in the normal fashion. The rest

of the line is entered in columns 7-46 of a continuation card immediately following the

primary card. The continuation card must have "$IVDRE" punched in columns 1-5. A maximum

of 120 characters may be printed on one line. The primary use of LIST is to build tables

in formats suitable for mounting on H0% reduction paper. The long linelength required by

such forms precludes the use of simple 80-80 listing.

There are two uses for LIST:

(1) to list a deck of cards:

@ASG,T BF. @DATA,I BF.,,LISTR

[card deck to be listed]

@END LISTR

@ADD,P PF.LIST

(2) to list an already-existing breakpoint file called BF:

@ADD,P PF.LIST

Example : LIST may be used to immediately print the output of YTFLD4:

@ASG,T BF.

@MAP,INX MAPR, YIELDS

IN PF. YIELDS

@BRKPT PUNCH$/BF

@XQT YIELD4

[parameter card]

[data sets]

6BRKPT PUNCH$

@ADD,P PF.LIST

PPLOT

PPLOT is a production plotting job-control deck. It will plot points, smooth curves, or zigzag lines. As written, the program produces point plots in a standard format.

Normally, the nuclide symbol is drawn above the plot. Also, whenever a data set is plotted

in POINT mode, the reaction and reference are drawn in the upper right-hand corner of the

plot. Two control cards are provided for altering the legend. Provisions are made,

however, to modify this format both through the use of plot control cards as well as by

entering the desired changes for specific plot parameters which control trie format used for a given plot. An arbitrary number of data sets may be plotted. Plotting may be done on

the printer, the SC 4020, or the Calcomp. It is recommended that a printer plot be made on

the slow printer.

18
A. Program Calling and Post-Processors

The plot program is called by the statement:

@ADD,P PF.PPLOT

A printer plot may be created on the slow printer by also calling

@ADD,P PF.PRINTR

The normal sequence of cards is:

@ADD,P PP.PPLOT

[control cards]

$$PARMS

$ Parameter Specifications ^ (only if needed)

$$END

[data sets]

@ADD,P PF.PRINTR (if printer plot is desired)

@MSG,W [user's account number], WRITE ENABLE REEL [user's reel number]

@ASG,TMJ 9.,8C,[user's reel number]W

gCLOSE 9

@XQT NBS*CALCMP.

B. Plot Control Cards

PPLOT is controlled by a set of control cards. On all of these cards, a dollar sign

must appear in column one. The control cards include:

General Controls

$CURVE

$LINES

$POINT

$CPLOT

The program normally plots data points. This

card switches to a curve-tracing mode.

This card switches to a zigzag line tracing mode.

(This is also useful for drawing histogram bar graphs.)

This card returns the program to point plotting mode.

The program will stay in CPLOT mode until a $FRAME

card is encountered. While in CPLOT mode:

-plots are 12" x 9.5"

-the nuclide ID is suppressed

-each data set is plotted on a new frame

-all these frames are drawn to the same scale

-the data sets are plotted as points connected

by zigzag lines

19
$SMALL A special control card to plot small frames. This

card causes a 4 x 2.5 in. frame to be plotted.

Two such small frames will fit side-by-side on 40%

reduction paper. Or, a small frame may be pasted

onto a regular frame as an insert. The effect of

$SMALL persists until a $FRAME card is encountered.

$FRAME This card marks the beginning of a new frame. The"

program returns to POINT mode.

Legend Controls

$LABEL Anything punched in card columns 7 through 24 of

this card is written as a single line into the plot

legend.

$PPLOT This card causes the current frame to be repeated

(re-drawn) without any annotation. The primary use of this card is to plot graphs to be sent to a drafts- man for lettering. If entered as the first card of

a frame, $PPLOT will cause a "stripped" frame to be

plotted. If entered as the last control card of a

frame, two plots will be drawn, one with annotation

and one without.

Data Modification Controls

The energies and/or cross section values of data sets may be renormalized by the

control cards:

$MUL E[value] Multiplies the energies of the next data set by the

given floating-point value.

$MUL S[value] Multiplies the cross sections of the next data set by the given value.

Similarly, offsets can be added to the data by:

$ADD E[value] Adds the value to the energies of the next data set.

$ADD S[value] Adds the value to the cross sections -of the next- data

set

.

All four of the above manipulate only the next data set. If used in combination, the

multiplication is done first. All values are given as floating-point numbers which must

start immediately after the E or S (i.e., no space).

C. Parameters and Defaults

PPLOT does not require any plot parameters to be specified. The normal plot constants

are:

XDIST =1.25

YDIST =1.25

XLNGTH =9.25

YLNGTH =6.00

FONTB = .09375 -FONT is plot symbol size. .09375 = 3/32 in.

20
t

1 1 1 r t

1 r

CM

i 2

o* lyvisAa:<

\-

(/> x

CMCM O

CM

J l_L

<o

X

I-o

z

_JX

tf)

tf> CM 00 (0 CM CM

I

H19N1A -Kisiqah

E3 o

M t/1 (- 0)

CD 3

C/> rH

<u «

Si >

ex

O o 3

if) o,+j o

<D >-H • X >—I XlrH o g mo • Sg 1/1 M-i o •H ^ II (J 0) hJ q <u <j£ ex

•H T3 C +-> U

XI <L> 4-1 CL> Q

CD X! > ^ u +-» 6a rt

CD +J —1 CD 6 iw3h -

oj U-Q •

Si (/) fn-H i—

I

rt 0) oj !h II uh cd > in «3>h CD .—I U Q

U CX-H +J >H S , f °

rf 4-1 S-l,—

I

« *J OUft in fi, <ni •H ^ h N

ex

<D ex n T3

21
YDOWN SYMBOL

YUP XSTART

YDIV YSTART

YDECPL

FONTH = .09375

CHARB = 0.1875 -CHAR is legend and numeric symbol size.

CHARH = 0.1875 0.1875 = 3/16 in.

SCALB = 0.15625 -SCAL is axis label symbol size.

SCALH = 0.15625 0.15625 = 5/32 in.

XLBL = PHOTON ENERGY (MEV) -If only one data set is plotted, XLBL and

YLBL = CROSS SECTION (MB)

YLBL are coPled f

*rom the label cards of the

data set.

The plot variables are normally calculated from the data. They include

XDOWN

XUP

XDIV

XDECPL

The various plot parameters and variables are defined in the figure given on p. 21.

The plot control switches, and defaults, are:

RULE = 1 Set to zero to inhibit upper right legend.

NTOP = 1 Set to zero to inhibit nuclide ID. MOVE = 1 Set to zero to inhibit frame advance.

RECALC = Set to 1 to re-calculate XDIV, YDIV, XDECPL,

YDECPL, XSTART, and YSTART.

Any and all of these parameters (in fact, the entire SPECS array) can be modified by means of the $$PARMS cards. The constants and variables are meant to be changed. However, it is

recommended that the casual user leave the plot switches alone.

Parameters are entered in FORTRAN "NAMELIST" format. That is, you punch the parameter

name, an equals sign, the numeric value, and a comma.

Example : $$PARMS

$ XDOWN=10 . ,XUP=40 . ,XDTV=4 . ,XDECPL=0 .

,

$$END

Don't forget the dollar signs in column one.' All parameters return to their default values

when "$FRAME" is encountered.

D. Examples :

Actual use of PPLOT is best illustrated by a few examples.

1. To plot a (y,sn), (y,2n), and (y,p) graph in a form ready for the draftsman,

use:

@ADD,P PF. PPLOT

[(y,sn) data]

C( Y ,2n) data]

[(Y,p) data]

22
$CURVE

[(Y,sn) curve]

[(Yj2n) curve]

[(y,p) curve]

$PPLOT

@ADD,P PF.PRINTR

@MSG,W [user's account number], WRITE ENABLE REEL [user's reel number]

@ASG,TMJ 9 .,8C, [user's reel number ]W

@CLOSE 9

.

@XQT NBS*CALCMP.

This will give both a printer plot and a tape for the Calcomp. To do the same thing, but

with a nine-inch, rather than a six-inch plot height, use:

@ADD,P PF.PPLOT

$$PARMS

$ YLNGTH=9 .

,

$$END

[(YjSn), (v,2n), and (y»p) data sets]

$CURVE

[(y,sn), (y,2n), and (y,p) curves]

$PPLOT

@ADD,P PF.PRINTR

etc.

2. To plot a (Y»xn) cross section and a (y 5 p) cross section on separate frames, use:

@ADD,P PF.PPLOT

[(y,xn) points]

$CURVE

[(y,xn) curve]

$FRAME

[(y,p) points]

$CURVE

[(y,p) curve]

@ADD,P PF.PRINTR

etc.

3. To plot a (y»p) cross section and a (Y>Mt

) cross section on separate frames, with

the intention of pasting the (y,u.) graph in the upper right-hand corner of the (y>p) graph

as an insert, and also assuming tne (yjM^) data consist of bar graph corners, use:

@ADD,P PF.PPLOT

$$PARMS

$ XSTART=1.5, (this moves the legend to the upper left comer)

$$END

23
[(y,p) points]

$CURVE

[(y,p) curve]

$PPLOT

$FRAME

$SMALL

$LABEL (G,MU-T) 64ER2

$LINES

[(Y,ut

) data]

$PPLOT

@ADD,P PF.PRINTR

etc.

4. To obtain a set of comparison plots, use:

@ADD 3 P PF.PPLOT

$CPLOT

[all (y,n), (y,sn), (y,xn), (y,v

t

) data]

$FRAME

$CPLOT

[all (y>p) and (y 5 xp) data]

$FRAME

$CPLOT

[all (Y,yt

) data]

$FRAME

$CPLOT

[all (Yjc) data]

etc.

SPCDST

This deck creates an index table of spectrum and angular distribution papers. Use is

made of the Photonuclear Data Index which must be in mass storage in a file called

"DC*INDEX.". The output table is placed in a file called "BF.", where it can be listed by

LIST. The user enters:

1) the nuclide

2) the references and reactions of SPC entries

3) DST

4) the references and reactions of DST entries

5) @EOF

The references must start in Column 1. The user should leave at least one space after the

reference as a separator. The reactions do not need parentheses.

21
Calling sequences and example:

@ASG,T BF. @ADD,P PF.SPCDST

5B 10

68SH6 G,P

DST

68SH6 G,P

63FU3 P,G

65PA1 HE,G

68LI1 HE,G

@EOF

@ADD,P PF.LIST

25 ftU.S. GOVERNMENT PRINTING OFFICE: 1976 210-801/165 1-3
NBS-1 14A (REV. 7-731

U.S. DEPT. OF COMM. BIBLIOGRAPHIC DATA

SHEET

1. PUBLICATION OR REPORT NO. NBS TN-9032. Gov't Accession

No.

3. Recipient's Accession No.

4. TITLE AND SUBTITLE

The NIRA Computer Program Package

(Photonuclear Data Center)

5. Publication Date February 1976

6. Performing Organization Code

7. AUTHOR(S)

H. J. Vander Molen and Henry M. Gerstenberj

8. Performing Organ. Report No.

9. PERFORMING ORGANIZATION NAME AND ADDRESS

NATIONAL BUREAU OF STANDARDS

DEPARTMENT OF COMMERCE

WASHINGTON, D.C. 20234

10. Project/Task/Work Unit No. 2410232

11. Contract/Grant No.

12. Sponsoring Organization Name and Complete Address (Street, City, State, ZIP)

National Science Foundation through the National

Academy of Sciences-National Research Council,

Committee on Nuclear Science

13. Type of Report & Period CoveredFinal

14. Sponsoring Agency Code

15. SUPPLEMENTARY NOTES

16. ABSTRACT (A 200-word or less tactual summary of most significant information. If document includes a significant bibliography or literature survey, mention it here.)

The Photonuclear Data Center's NIRA library of programs, executable from mass storage on the National Bureau of Standard's central computer facility, is

described. Detailed instructions are given (with examples) , for the use of the

library to analyze, evaluate, synthesize, and produce for publication, camera-ready

tabular and graphical presentations of digital photonuclear reaction cross-section

data. NIRA is the acronym for Nuclear Information Research Associate.

17. KEY WORDS (six to twelve entries; alphabetical order; capitalize only the first letter of the first key word unless a proper name; separated by semicolons)

Central computer; computer programs; data manipulation; digital data; mass storage:

NIRA; Photonuclear Data Center.

18. AVAILABILITY [x} Unlimited

_ For Official Distribution. Do Not Release to NTIS

|y [ Order From Sup. of Doc, U.S. Government Printing Office Washington, D.C. 20402, SD Cat. No. CI 3. 46: 903

_J Order From National Technical Information Service (NTIS)

Springfield, Virginia 22I51

19. SECURITY CLASS

(THIS REPORT)

UNCLASSIFIED

20. SECURITY CLASS

(THIS PAGE)

UNCLASSIFIED

21. NO. OF PAGES

30

22. Pr

USCOMM-DC 29042-P74
NBS TECHNICAL PUBLICATIONS

PERIODICALS JOURNAL OF RESEARCH reports National Bureau

of Standards research and development in physics,

mathematics, and chemistry. It is published in two sections, available separately:

• Physics and Chemistry (Section A)

Papers of interest primarily to scientists working in these fields. This section covers a broad range of physical and chemical research, with major emphasis on

standards of physical measurement, fundamental constants, and properties of matter. Issued six times a

year. Annual subscription: Domestic, $17.00; Foreign,

$21.25.

• Mathematical Sciences (Section B)

Studies and compilations designed mainly for the mathematician and theoretical physicist. Topics in mathematical statistics, theory of experiment design, numerical analysis, theoretical physics and chemistry, logical

design and programming of computers and computer

systems. Short numerical tables. Issued quarterly. Annual subscription: Domestic, $9.00; Foreign, $11.25. DIMENSIONS/NBS (formerly Technical News Bulletin) —This monthly magazine is published to inform

scientists, engineers, businessmen, industry, teachers,

students, and consumers of the latest advances in science and technology, with primary emphasis on the work at NBS. The magazine highlights and reviews such

issues as energy research, fire protection, building tech- nology, metric conversion, pollution abatement, health and safety, and consumer product performance. In addition, it reports the results of Bureau programs in measurement standards and techniques, properties of matter and materials, engineering standards and services, instrumentation, and automatic data processing.

Annual subscription: Domestic, $9.45; Foreign, $11.85.

NONPERIODICALS Monographs—Major contributions to the technical liter- ature on various subjects related to the Bureau's scientific and technical activities. Handbooks—Recommended codes of engineering and

industrial practice (including safety codes) developed

in cooperation with interested industries, professional

organizations, and regulatory bodies.

Special Publications—Include proceedings of confer- ences sponsored by NBS, NBS annual reports, and other

special publications appropriate to this grouping such

as wall charts, pocket cards, and bibliographies.

Applied Mathematics Series—Mathematical tables, manuals, and studies of special interest to physicists,

engineers, chemists, biologists, mathematicians, computer programmers, and others engaged in scientific and technical work.

National Standard Reference Data Series—Provides

quantitative data on the physical and chemical properties of materials, compiled from the world's literature and critically evaluated. Developed under a world-wide

program coordinated by NBS. Program under authority

of National Standard Data Act (Public Law 90-396). NOTE: At present the principal publication outlet for these data is the Journal of Physical and Chemical

Reference Data (JPCRD) published quarterly for NBS

by the American Chemical Society (ACS) and the American Institute of Physics (AIP). Subscriptions, reprints, and supplements available from ACS, 1155 Sixteenth

St. N. W., Wash. D. C. 20056. Building Science Series—Disseminates technical infor- mation developed at the Bureau on building materials,

components, systems, and whole structures. The series presents research results, test methods, and performance criteria related to the structural and environmental functions and the durability and safety characteristics of building elements and systems.

Technical Notes—Studies or reports which are complete

in themselves but restrictive in their treatment of a

subject. Analogous to monographs but not so comprehensive in scope or definitive in treatment of the subject area. Often serve as a vehicle for final reports of work performed at NBS under the sponsorship of other government agencies.

Voluntary Product Standards—Developed under pro- cedures published by the Department of Commerce in Part 10, Title 15, of the Code of Federal Regulations.

The purpose of the standards is to establish nationally

recognized requirements for products, and to provide

all concerned interests with a basis for common under- standing of the characteristics of the products. NBS

administers this program as a supplement to the activities of the private sector standardizing organizations. Federal Information Processing Standards Publications (FIPS PUBS)—Publications in this series collectively constitute the Federal Information Processing Standards Register. Register serves as the official source of information in the Federal Government regarding stand- ards issued by NBS pursuant to the Federal Property and Administrative Services Act of 1949 as amended,

Public Law 89-306 (79 Stat. 1127), and as implemented

by Executive Order 11717 (38 FR 12315, dated May 11, 1973) and Part 6 of Title 15 CFR (Code of Federal

Regulations). Consumer Information Series—Practical information, based on NBS research and experience, covering areas

of interest to the consumer. Easily understandable language and illustrations provide useful background

knowledge for shopping in today's technological

marketplace.

NBS Interagency Reports (NBSIR)—A special series of interim or final reports on work performed by NBS for outside sponsors (both government and non-government). In general, initial distribution is handled by the sponsor; public distribution is by the National Technical Information Service (Springfield, Va. 22161) in paper

copy or microfiche form.

Order NBS publications (except NBSIR's and Biblio- graphic Subscription Services) from: Superintendent of Documents, Government Printing Office, Washington,

D.C. 20402. BIBLIOGRAPHIC SUBSCRIPTION SERVICES

The following current-awareness and literature-survey

bibliographies are issued periodically by the Bureau:

Cryogenic Data Center Current Awareness Service A literature survey issued biweekly. Annual subscription: Domestic, $20.00; foreign, $25.00.

Liquefied Natural Gas. A literature survey issued quarterly. Annual subscription: $20.00. Superconducting Devices and Materials. A literature

survey issued quarterly. Annual subscription : $20.00.

Send subscription orders and remittances for the preceding bibliographic services to National Bureau of Standards, Cryogenic Data Center (275.02)

Boulder, Colorado 80302. Electromagnetic Metrology Current Awareness Service Issued monthly. Annual subscription: $24.00. Send

subscription order and remittance to Electromagnetics

Division, National Bureau of Standards, Boulder,

Colo. 80302.
U.S. DEPARTMENT OF COMMERCE

National Bureau of Standards

Washington. DC. 20234

OFFICIAL BUSINESS

Penalty for Private Use. $300

POSTAGE AND FEES PAID

U.S. DEPARTMENT OF COMMERCE COM-215

SPECIAL FOURTH-CLASS RATE BOOK

V/U) YEARS

/^6-l9'fe 1901-197B