MIME-Version: 1.0 Content-Location: file:///C:/F2CAA98E/Virgin.htm Content-Transfer-Encoding: quoted-printable Content-Type: text/html; charset="us-ascii"
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
PROGRAM
VIRGIN =
&nb=
sp; =
&nb=
sp;
Virgin
VERSION 76=
-1
(NOVEMBER 1976) =
&nb=
sp; =
Virgin
VERSION 84=
-1
(JUNE 1984) *DOUBLE PRE=
CISION
ENERGY =
Virgin
VERSION 86=
-1
(JANUARY 1986)*FORTRAN-77/H VERSION =
Virgin
VERSION 88=
-1 (JULY
1988) *OPTION...INTERNA=
LLY
DEFINE ALL I/O
Virgin
=
&nb=
sp;
FILE NAMES (SEE, SUBROUTINE FILEIO Virgin
=
&nb=
sp;
FOR DETAILS). =
&nb=
sp;
Virgin
=
&nb=
sp;
*IMPROVED =
BASED
ON USER COMMENTS.
Virgin
VERSION 89=
-1
(JANUARY 1989)*PSYCHOANALYZED BY PROGRAM FREUD TO Virgin
=
&nb=
sp;
INSURE PROGRAM WILL NOT DO ANYTHING Virgin
=
&nb=
sp;
CRAZY. =
&nb=
sp;
Virgin
=
&nb=
sp;
*UPDATED TO USE NEW PROGRAM CONVERT Virgin
=
&nb=
sp;
KEYWORDS. =
&nb=
sp;
Virgin
=
&nb=
sp;
*ADDED
=
&nb=
sp;
CONVENTIONS. =
&nb=
sp;
Virgin
VERSION 92=
-1
(JANUARY 1992)*COMPLETE RE-WRITE =
&nb=
sp;
Virgin
=
&nb=
sp;
*OUTPUT IN PLOTTAB FORMAT =
Virgin
=
&nb=
sp;
*UP TO 2000 THICKNESSES =
Virgin
=
&nb=
sp;
*INCREASED INCORE PAGE SIZE TO 6000 Virgin
=
&nb=
sp;
CROSS SECTION POINTS =
Virgin
=
&nb=
sp;
*ADDED PHOTON CALCULATIONS =
Virgin
=
&nb=
sp;
*ADDED BLACKBODY SPECTRUM =
Virgin
=
&nb=
sp;
*ADDED MULTIPLE LAYERS =
Virgin
=
&nb=
sp;
*ADDED SPATIALLY DEPENDENT DENSITY Virgin
=
&nb=
sp;
*ADDED FORTRAN SAVE OPTION =
Virgin
=
&nb=
sp;
*COMPLETELY CONSISTENT I/O ROUTINES - Virgin
=
&nb=
sp; TO
MINIMIZE COMPUTER DEPENDENCE.
Virgin
VERSION 92=
-2
(MAY 1992) *CORRE=
CTED
TO HANDLE MULTIGROUP CROSS
Virgin
=
&nb=
sp;
SECTIONS AS INPUT IN ENDF/B FORMAT. Virgin
VERSION 96=
-1
(JANUARY 1996) *COMPLETE RE-WRITE =
Virgin
=
&nb=
sp;
*IMPROVED COMPUTER
=
&nb=
sp;
*ALL DOUBLE PRECISION =
Virgin
=
&nb=
sp;
*ON S=
CREEN
OUTPUT =
&nb=
sp;
Virgin
=
&nb=
sp;
*UNIFORM TREATMENT OF ENDF/B I/O
Virgin
=
&nb=
sp;
*IMPROVED OUTPUT PRECISION =
Virgin
=
&nb=
sp;
*DEFINED SCRATCH FILE NAMES =
Virgin
VERSION 99=
-1
(MARCH 1999) *CORRECTED
CHARACTER TO FLOATING Virg=
in
=
&nb=
sp; =
POINT READ FOR MORE DIGITS =
Virgin
=
&nb=
sp;
*UPDATED TEST FOR ENDF/B FORMAT
Virgin
=
&nb=
sp; =
VERSION BASED ON RECENT FORMAT CHANGE Virgin
=
&nb=
sp;
*GENERAL IMPROVEMENTS BASED ON
Virgin
=
&nb=
sp; =
USER FEEDBACK =
=
Virgin
VERS. 2000=
-1
(FEBRUARY 2000)*GENERAL IMPROVEMENTS BASED ON
Virgin
=
&nb=
sp; =
USER FEEDBACK =
&nb=
sp;
Virgin
VERS. 2002=
-1
(MAY 2002)
*OPTIONAL INPUT PARAMETERS Virgin
VERS. 2004=
-1
(MARCH 2004) *ADDED INC=
LUDE
FOR COMMON =
Virgin
=
&nb=
sp;
*UP TO 2000 THICKNESSES =
Virgin
=
&nb=
sp;
*INCREASED INCORE PAGE SIZE TO 60,000 Virgin
VERS. 2007=
-1
(JAN. 2007) *CHECKED AGAINST ALL ENDF/B-=
VII.
Virgin
=
&nb=
sp;
*INCREASED INCORE PAGE SIZE TO
Virgin
=
&nb=
sp; =
240,000 FROM 60,000. =
Virgin
=
&nb=
sp; =
&nb=
sp; =
Virgin
OWNED,
MAINTAINED AND DISTRIBUTED BY =
&nb=
sp;
Virgin
------------------------------------ =
&nb=
sp;
Virgin
THE NUCLEA=
R DATA
SECTION =
&nb=
sp; =
Virgin
INTERNATIO=
NAL
ATOMIC ENERGY AGENCY =
&nb=
sp;
Virgin
P.O. BOX
100 =
&nb=
sp; =
&nb=
sp;
Virgin
A-1400, VI=
ENNA,
AUSTRIA =
&nb=
sp; =
Virgin
EUROPE =
&nb=
sp; =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
ORIGINALLY
WRITTEN BY
=
&nb=
sp; =
Virgin
------------------------------------ =
&nb=
sp;
Virgin
DERMOTT E.
CULLEN =
&nb=
sp; =
&nb=
sp;
Virgin
UNIVERSITY=
OF
CALIFORNIA =
=
&nb=
sp; Virgin
LAWRENCE
LIVERMORE NATIONAL LABORATORY =
&nb=
sp;
Virgin
L-159 =
&nb=
sp; =
&nb=
sp; =
Virgin
P.O. BOX
808 =
&nb=
sp; =
=
Virgin
LIVERMORE,=
CA
94550 =
&nb=
sp; =
&nb=
sp;
Virgin
U.S.A. =
&nb=
sp; =
&nb=
sp; =
Virgin
TELEPHONE<=
span
style=3D'mso-spacerun:yes'> 925-423-7359 =
&nb=
sp; =
Virgin
E. MAIL
WEBSITE
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
PURPOSE
-------
THIS PROGR=
AM IS
DESIGNED TO CALCULATE UNCOLLIDED (I.E. VIRGIN) Virgin
FLUX AND REACTIONS DUE TO
TRANSMISSION OF A MONODIRECTIONAL
Virgin
BEAM OF NE=
UTRONS
THROUGH ANY THICKNESS OF MATERIAL. IN ORDER
Virgin
TO SIMULAT=
E AN
EXPERIMENTAL MEASUREMENT THE RESULTS ARE GIVEN Virgin
AS INTEGRA=
LS
OVER ENERGY TALLY GROUPS (AS OPPOSED TO POINTWISE Virgin
IN ENERGY)=
. BY
TAKING THE RATIO OF REACTIONS TO FLUX IN EACH
Virgin
GROUP AN
EQUIVALENT SPATIALLY DEPENDENT GROUP AVERAGED CROSS
Virgin
SECTION IS
CALCULATED BY THE PROGRAM. =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
EVALUATED
DATA =
&nb=
sp; =
&nb=
sp;
Virgin
-------------- =
=
&nb=
sp; =
Virgin
THE EVALUA=
TED
DATA MUST BE IN THE ENDF/B FORMAT. HOWEVER IT
Virgin
MUST BE
LINEAR-LINEAR INTERPOLABLE IN ENERGY-CROSS SECTION
Virgin
BETWEEN
TABULATED POINTS. SINCE ONLY CROSS SECTIONS (FILE 3 OR 23) Virgin
ARE USED, =
THIS
PROGRAM WILL WORK ON ANY VERSION OF ENDF/B =
Virgin
(I.E. ENDF=
/B-I,
II, III, IV, V OR VI). =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp;
=
Vi=
rgin
RELATED CO=
MPUTER
CODES =
&nb=
sp; =
Virgin
---------------------- =
&nb=
sp; =
Virgin
IN ORDER TO
CONVERT ENDF/B DATA TO THE FORM REQUIRED BY THIS CODE Virgin
THE FOLLOW=
ING
COMPUTER CODES MAY BE USED, =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
LINEAR - CONVERT FROM GENERAL ENDF/B
INTERPOLATION TO LINEAR- =
Virgin
=
LINEAR INTERPOLATION. =
&nb=
sp; =
Virgin
RECENT - ADD THE RESONANCE CONTRIBUTION TO
TABULATED BACKGROUND Vi=
rgin
=
CROSS SECTIONS TO OBTAIN LINEAR-LINEAR INTERPOLABLE Virgin
=
RESULTS. =
&nb=
sp; =
&nb=
sp;
Virgin
SIGMA1 - DOPPLER BROADEN CROSS SECTION TO
OBTAIN LINEAR-LINEAR
Virgin
=
INTERPOLABLE RESULTS. =
&nb=
sp; =
Virgin
MIXER - MIX INDIVIDUAL MATERIALS TO=
GETHER
TO DEFINE COMPOSITE
Virgin
=
MIXTURES, E.G., COMBINE MATERIALS TO DEFINE STAINLESS Virgin
=
STELL. =
&nb=
sp; =
&nb=
sp;
Virgin
=
&nb=
sp;
=
&nb=
sp; =
&nb=
sp; Virgin
IN ORDER T=
O PLOT
THE OUTPUT RESULTS OF THIS CODE USE PROGRAM
Virgin
PLOTTAB. =
&nb=
sp; =
&nb=
sp;
Virgin
=
&nb=
sp;
=
&nb=
sp; =
Virgin
COPIES OF =
ANY OR
ALL OF THESE CODES MAY BE OBTAINED FROM D.E. Virgin
CULLEN AT =
THE
ABOVE ADDRESS. =
&nb=
sp; =
Virgin
=
&nb=
sp; =
=
&nb=
sp; Virgi=
n
OUTPUT
FORMAT =
&nb=
sp; =
&nb=
sp;
Virgin
------------- =
&nb=
sp; =
&nb=
sp;
Virgin
FOR ALL VE=
RSIONS
OF THIS PROGRAM PRIOR TO VERSION 92-1 OUTPUT WAS Virgin
IN TABULAR
FORM. =
&nb=
sp; =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
FOR VERSIO=
N 92-1
AND LATER VERSIONS OF THIS CODE ALL OUTPUT IS IN Virgin
THE PROGRAM
PLOTTAB FORMAT TO ALLOW RESULTS TO BE EASILY PLOTTED. Virgin
FOR A COPY=
OF
PROGRAM PLOTTAB CONTACT D.E. CULLEN AT THE ABOVE Virgin
ADDRESS. =
&nb=
sp; =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
TALLY
GROUPS =
&nb=
sp; =
&nb=
sp;
Virgin
------------ =
&nb=
sp; =
&nb=
sp;
Virgin
THE TALLY =
GROUP
STRUCTURE MAY BE ANY SET OF MONONTONICALLY
Virgin
INCREASING
ENERGY BOUNDARIES. THERE MAY BE UP TO 2000 TALLY
Virgin
GROUPS. BY=
USING
THE INPUT PARAMETERS THE USER MAY SPECIFY ANY Virgin
ARBITRARY =
TALLY
GROUP STRUCTURE OR SELECT ONE OF THE FOLLOWING Virgin
BUILT-IN G=
ROUP
STRUCTURES. =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
(1) TART 1=
75 GROUPS
(INPUT 0) =
&nb=
sp; =
Virgin
(2) ORNL 50 GROUPS (INPUT -1) =
&nb=
sp; =
Virgin
(3) ORNL 1=
26
GROUPS (INPUT -2) =
&nb=
sp; =
Virgin
(4) ORNL 1=
71
GROUPS (INPUT -3) =
&nb=
sp; =
Virgin
(5) SAND-I=
I 620
GROUPS...UP TO 18 MEV (INPUT -4) =
Virgin
(6) SAND-I=
I 640
GROUPS...UP TO 20 MEV (INPUT -5) =
Virgin
(7) WIMS 69
GROUPS (INPUT -6)
=
&nb=
sp; =
Virgin
(8) GAM-I =
68
GROUPS (INPUT -7) =
&nb=
sp; =
Virgin
(9) GAM-II=
99
GROUPS (INPUT -8) =
&nb=
sp; =
Virgin
(10) MUFT 54 GRO=
UPS
(INPUT -9) =
=
&nb=
sp; Virgin
(11) ABBN 28 GRO=
UPS
(INPUT -10) =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
INCIDENT
SPECTRUM =
&nb=
sp; =
=
Virgin
----------------- =
&nb=
sp; =
&nb=
sp;
Virgin
THE INCIDE=
NT
SPECTRUM MAY BE ANY TABULATED FUNCTION THAT IS
Virgin
GIVEN BY A=
SET
OF POINTS THAT IS MONOTONICALLY INCREASING IN Virgin
ENERGY AND
LINEAR-LINEAR INTERPOLABLE IN ENERGY-SPECTRUM =
Virgin
BETWEEN
TABULATED POINTS. THERE IS NO LIMIT TO THE NUMBER OF
Virgin
POINTS USE=
D TO
DESCRIBE THE SPECTRUM. THERE ARE FIVE BUILT-IN Virgin
OPTIONS FO=
R THE
SPECTRUM. =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
(1)
CONSTANT...ENERGY INDEPENDENT (INPUT 0) =
&nb=
sp;
Virgin
(2) 1/E (I=
NPUT
1) =
&nb=
sp; =
&nb=
sp;
Virgin
(3) BLACKB=
ODY -
PHOTON SPECTRUM =
&nb=
sp; =
Virgin
(4) BLACKB=
ODY -
ENERGY SPECTRUM (E TIMES THE PHOTON SPECTRUM) Virgin
(5) TRANSM=
ITTED
SPECTRUM FROM PREVIOUS CASE =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
NORMALIZAT=
ION OF
SPECTRUM =
&nb=
sp; =
Virgin
------------------------- =
&nb=
sp; =
Virgin
ANY INCIDE=
NT
SPECTRUM, EITHER READ AS INPUT OR ONE OF THE =
Virgin
BUILT-IN
SPECTRA, WILL BE NORMALIZED TO UNITY WHEN INTEGRATED Virgin
OVER THEIR
ENTIRE ENERGY RANGE. =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
TRANSMITTED
SPECTRA WILL NOT BE RE-NORMALIZED, SINCE IT ALREADY Virgin
INCLUDES T=
HE
NORMALIZATION OF THE INCIDENT SPECTRUM. =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
NOTE, INCI=
DENT
SPECTRA IS NORMALIZED TO UNITY OVER THEIR ENTIRE Virgin
ENERGY RAN=
GE -
NOT OVER THE ENERGY RANGE OF THE GROUPS. IF THE Virgin
ENERGY RAN=
GE OF
THE GROUPS IS LESS THAN THAT OF THE SPECTRUM
Virgin
ONLY THAT
PORTION OF THE SPECTRUM WILL BE USED AND THIS WILL
Virgin
NOT BE
RE-NORMALIZED TO UNITY. =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
COMPOSITIO=
N OF A
LAYER =
&nb=
sp; =
Virgin
---------------------- =
&nb=
sp; =
Virgin
YOU MAY RUN
PROBLEMS INVOLVING =
&nb=
sp; =
Virgin
1) A LAYER=
OF
UNIFORM DENSITY - DENSITY FOR ATTENUATION IS THAT Virgin
OF THE TOTAL. DENSITY FOR REACTIONS IS THAT OF THE REACTION. Virgin
2) A LAYER=
OF
UNIFORM DENSITY - DENSITY IS THE SUM OF THE TOTAL Virgin
AND REACTION DENSITIES - THE SUM OF THE CROSS SECTIONS IS Virg=
in
USED FOR ATTENUATION AND REACTIONS. =
&nb=
sp;
Virgin
3) A LAYER=
OF
VARYING DENSITY BASED ON A UNIFORM TOTAL DENSITY Virgin
PLUS A VARIATION BETWEEN 0 AND A MAXIMUM BASED ON THE =
Virgin
REACTION DENSITY - 0 AT 0 THICKNESS AND MAXIMUM AT MAXIMUM Virgin
THICKNESS. IN THIS CASE THE AVERAGE REACTION DENSITY IS EQUAL Virgin
TO THE INPUT REACTION DENSITY. THE VARIATION IN REACTION
Virgin
DENSITY CAN BE LINEAR, SQUARE OR CUBIC. =
&nb=
sp;
Virgin
4) A LAYER=
OF
VARYING DENSITY BASED ON A TOTAL DENSITY WHICH
Virgin
VARYING FROM MAXIMUM AT 0 THICKNESS TO 0 AT MAXIMUM THICKNESS Virgin
PLUS A REACTION DENSITY WHICH VARIES FROM 0 AT 0 THICKNESS Virgin
TO MAXIMUM AT MAXIMUM THICKNESS. IN THIS CASE THE AVERAGE
Virgin
DENSITY OF THE TOTAL AND REACTION WILL BOTH BE EQUAL TO THE Virgin
INPUT TOTAL AND REACTION DENSITIES. THE VARIATION IN TOTAL Virgin
AND REACTION DENSITY CAN BE LINEAR, SQUARE OR CUBIC. =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
IN THE FIR=
ST
CASE THE TWO REQUESTED CROSS SECTIONS ARE CONSIDERED Virgin
TO BE
INDEPENDENT - THE TOTAL CROSS SECTION IS USED TO CALCULATE Virgin
ATTENUATIO=
N AND
THE REACTION CROSS SECTION IS USED TO CALCULATE Virgin
REACTIONS,=
E.G.,
TRANSMISSION THROUGH NATURAL URANIUM (THE TOTAL Virgin
CROSS SECT=
ION
SHOULD BE THAT OF NATURAL URANIUM) AND REACTIONS Virgin
IN A U-235
DETECTOR (THE REACTION CROSS SECTION SHOULD BE THAT OF Virgin
U-235).
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
IN THE OTHER THREE CASES THE =
TWO
REQUESTED CROSS SECTIONS ARE Virgin
TREATED AS=
TWO
CONSTITUENTS OF A MIXTURE OF TWO MATERIALS AND Virgin
THE TWO CR=
OSS
SECTIONS ARE USED BOTH TO DEFINE A TOTAL CROSS
Virgin
SECTION FOR
ATTENUATION AND A REACTION CROSS SECTION TO DEFINE Virgin
REACTIONS.=
IN
THESE CASES THE MIXTURE WILL VARY CONTINUOUSLY, Virgin
E.G., IN C=
ASE 4)
HALF WAY THROUGH THE LAYER THE COMPOSITION WILL Virgin
BE 1/2 THE
MATERIAL DEFINED BY THE TOTAL AND 1/2 THE MATERIAL Virgin
BASED ON T=
HE
REACTION. IN THESE CASES RATHER THAN THINKING OF Virgin
THE TWO CR=
OSS
SECTIONS AS A TOTAL AND REACTION CROSS SECTION, Virgin
IT IS BETT=
ER TO
THINK OF THEM AS THE TOTAL CROSS SECTIONS FOR Virgin
MATERIALS =
A AND
B AND THE CALCULATED REACTIONS WILL BE BASED
Virgin
ON THESE T=
WO
TOTAL CROSS SECTIONS. =
&nb=
sp;
Virgin
=
&nb=
sp; =
=
&nb=
sp; Virgin
MULTIPLE
LAYERS =
&nb=
sp; =
&nb=
sp;
Virgin
--------------- =
&nb=
sp; =
&nb=
sp;
Virgin
THIS CODE =
MAY BE
USED TO RUN EITHER A NUMBER OF INDEPENDENT
Virgin
PROBLEMS, =
EACH
INVOLVING TRANSMISSION THROUGH A SINGLE LAYER OF Virgin
MATERIAL, =
OR
TRANSMISSION THROUGH A NUMBER OF LAYERS ONE AFTER Virgin
THE OTHER.=
=
&nb=
sp; =
&nb=
sp;
=
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
IN THE CAS=
E OF
MULTIPLE LAYERS, ONE LAYER AFTER ANOTHER, THE
Virgin
TRANSMITTED
ENERGY DEPENDENT SPECTRUM IS USED AS THE INCIDENT Virgin
SPECTRUM F=
OR THE
NEXT LAYER. THERE IS NO LIMIT TO THE NUMBER
Virgin
OF LAYERS =
WHICH
MAY BE USED - EACH LAYER IS TREATED AS A =
Virgin
COMPLETELY
INDEPENDENT PROBLEM WITH A DEFINED INCIDENT SOURCE, Virgin
AND AS SUCH THE =
CYCLE
OF TRANSMISSION THROUGH EACH LAYER AND
Virgin
USING THE
TRANSMITTED SPECTRUM AS THE INCIDENT SPECTRUM FOR THE Virgin
NEXT LAYER=
MAY
BE REPEATED ANY NUMBER OF TIMES. =
Virgin
=
=
&nb=
sp; =
&nb=
sp; =
Virgin
REMEMBER -=
THE
INCIDENT SPECTRUM IS ASSUMED TO BE LINEARLY
Virgin
INTERPOLAB=
LE IN
ENERGY AND SPECTRUM BETWEEN THE ENERGIES AT
Virgin
WHICH IT I=
S TABULATED.
THE TRANSMITTED SPECTRUM WILL BE TABULATED=
Virgin
AT THE UNI=
ON OF
ALL ENERGIES OF THE INCIDENT SPECTRUM AND CROSS Virgin
SECTIONS (=
TOTAL
AND REACTION). IN ORDER TO INSURE THE ACCURACY Virgin
OF THE RES=
ULT
WHEN PERFORMING MULTIPLE LAYER CALCULATION BE SURE Virgin
TO SPECIFY=
THE
INCIDENT SPECTRUM ON THE FIRST LAYER TO SUFFICIENT Virgin
DETAIL (EN=
OUGH
ENERGY POINTS CLOSELY SPACED TOGETHER) IN ORDER TO Virgin
ALLOW THE
TRANSMITTED SPECTRUM TO BE ACCURATELY REPRESENTED BY Virgin
LINEAR
INTERPOLATION BETWEEN SUCCESSIVE ENERGY POINTS - THERE IS Virgin
NO LIMIT T=
O THE
NUMBER OF POINTS ALLOWED IN THE INCIDENT SPECTRUM, Virgin
SO IF YOU =
ARE IN
DOUBT, SIMPLY USE MORE ENERGY POINTS TO SPECIFY Virgin
THE INCIDE=
NT
SPECTRUM. =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
RESULT OUT=
PUT
UNITS =
&nb=
sp; =
=
Virgin
------------------- =
&nb=
sp; =
&nb=
sp;
Virgin
FLUX =3D EXACTL=
Y AS
CALCULATED =
&nb=
sp;
Virgin
REACTIONS =
=3D 1/CM
OR 1/GRAM =
&nb=
sp; =
Virgin
AVERAGE
CROSS =
&nb=
sp; =
&nb=
sp; =
Virgin
SECTION
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
THICKNESS =
AND
DENSITY =
&nb=
sp; =
Virgin
--------------------- =
&nb=
sp; =
Virgin
THE UNCOLL=
IDED
CALCULATION ONLY DEPENDS ON THE PRODUCT OF =
Virgin
THICKNESS =
AND
DENSITY (I.E. GRAMS PER CM SQUARED). THIS FACT
Virgin
MAY BE USE=
D TO
SIMPLIFY INPUT BY ALLOWING THE THICKNESS AND
Virgin
DENSITY TO=
BE
GIVEN EITHER AS CM AND GRAMS/CC RESPECTIVELY
Virgin
OR ELSE TO=
GIVE
THICKNESS IN GRAMS/(CM*CM) AND INPUT A =
Virgin
DENSITY OF=
1.0 -
OR IN ANY OTHER CONVENIENT UNITS AS LONG AS
Virgin
THE PRODUC=
T OF
THICKNESS AND DENSITY IS IN THE CORRECT GRAMS
Virgin
PER CENTIM=
ETER
SQUARED. =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
GRAMS/(CM*=
CM)
ARE RELATED TO ATOMS/BARN THROUGH THE RELATIONSHIP Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
GRAMS/(CM*CM)=3D(ATOMS/BARN)*(GRAMS/MOLE)*(MOLE/ATOM) =
Virgin
=
&nb=
sp; =
&nb=
sp;
=
&nb=
sp; Virgin
OR... =
&nb=
sp; =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
GRAMS/(CM*CM)=3D(ATOMS/BARN)*(ATOMIC WEIGHT)/0.602 =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
CROSS SECT=
IONS
AT A SPACE POINT AND OPTICAL THICKNESS =
Virgin
----------------------------------------------------- =
Virgin
THIS PROGR=
AM
ALLOWS LAYERS OF EITHER UNIFORM DENSITY OR =
Virgin
CONTINUOUS=
LY
VARYING DENSITY. THE DENSITY CAN BE ONE OF THE
Virgin
FOLLOWING
FORMS, =
&nb=
sp; =
&nb=
sp;
Virgin
1) C =
=3D UNIFORM DENSITY =
&nb=
sp;
Virgin
2)
C*2*(X/T) =
=3D LINEAR VARIATION FROM 0 TO C =
Virgin
3)
C*(2-2*(X/T)) =3D =
LINEAR
VARIATION FROM C TO 0 =
Virgin
4)
C*3*(X/T)**2 =3D
SQUARE VARIATION FROM 0 TO C =
Virgin
5)
C*(3-3*(X/T)**2)/2 =3D SQUARE
VARIATION FROM C TO 0 =
Virgin
6)
C*4*(X/T)**3 =3D
CUBIC VARIATION FROM 0 TO C =
Virgin
7) C*(4-4*=
(X/T)**3)/3 =3D CUBIC VARIATION FROM C TO 0
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
IN ORDER TO
CALCULATE REACTIONS AT A POINT THE MICROSCOPIC
Virgin
REACTION C=
ROSS
SECTION NEED MERELY BE SCALED BY THESE DENSITIES. Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
IN ORDER TO
CALCULATE TRANSMISSION WE MUST DEFINE THE OPTICAL Virgin
PATH LENGTH
WHICH MAY BE DEFINED BY INTEGRATING EACH OF THE
Virgin
ABOVE DENS=
ITY
FORMS TO FIND, =
&nb=
sp; =
Virgin
1) C*X =
&nb=
sp; =
&nb=
sp; =
Virgin
2)
C*X*(X/T) =
&nb=
sp;
=
&nb=
sp; Virgin
3)
C*X*(2-(X/T)) =
&nb=
sp; =
&nb=
sp;
Virgin
4)
C*X*(X/T)**2 =
&nb=
sp; =
&nb=
sp;
Virgin
5)
C*X*(3-(X/T)**2)/2 =
&nb=
sp;
=
Virgin
6)
C*X*(X/T)**3 =
&nb=
sp; =
&nb=
sp;
Virgin
7)
C*X*(4-(X/T)**3))/3 =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
Virgi=
n
IN ORDER TO
CALCULATE TRANSMISSION TO A POINT THE MICROSCOPIC Virgin
TOTAL CROSS
SECTION NEED MERELY BE SCALED BY THESE DENSITIES
Virgin
TO DEFINE =
THE
OPTICAL PATH LENGTH. =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
THE VARIAT=
ION OF
THE DENSITY THROUGH THE LAYER MAY BE DEFINED Virgin
BY SETTING=
X =3D 0
OR X =3D T TO FIND, =
&nb=
sp;
Virgin
X =3D 0 X =3D T
-----
----- =
&nb=
sp; =
&nb=
sp;
Virgin
1) C =
C =
&nb=
sp; =
&nb=
sp;
Virgin
2) 0 2*C =
&nb=
sp; =
&nb=
sp;
Virgin
3) 2*C
0 =
&nb=
sp; =
&nb=
sp;
Virgin
4) 0 =
3*C =
&nb=
sp; =
&nb=
sp;
Virgin
5) 3*C/2 0 =
&nb=
sp; =
&nb=
sp; Virgin
6) 0 =
4*C =
&nb=
sp; =
&nb=
sp;
Virgin
7) 4*C/3 0 =
&nb=
sp; =
&nb=
sp;
Virgin
=
&nb=
sp; =
=
&nb=
sp; =
Virgin
THE OPTICA=
L PATH
THROUGH A LAYER OF THICKNESS T MAY BE DEFINED Virgin
FROM THE A=
BOVE
EXPRESSIONS BY SETTING X=3DT TO FIND THAT IN ALL Virgin
CASES THE =
ANSWER
WILL BY C*T. THE CONSTANTS IN THE ABOVE =
Virgin
EXPRESSION=
S HAVE
BEEN INTRODUCED IN ORDER TO FORCE THIS RESULT. Virgin
WITH THESE
FACTORS THE OPTICAL PATH LENGTH THROUGH THE LAYER
Virgin
WILL EXACT=
LY
CORRESPOND TO AN AVERAGE DENSITY CORRESPONDING TO Virgin
THAT INPUT=
FOR
THE TOTAL AND/OR REACTION, I.E., C CORRESPONDS Virgin
TO THE INP=
UT
DENSITY. =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
=
Virgin
NOTE - FOR=
THE
SAME OPTICAL PATH LENGTHS THROUGH THE LAYER THE Virgin
TRANSMISSI=
ON
WILL BE EXACTLY THE SAME. HOWEVER, VARYING THE
Virgin
DENSITY WI=
LL
ALLOW YOU TO MODIFY THE REACTION RATES AT SPECIFIC Virgin
DEPTHS INT=
O THE
LAYER. =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
COMPUTATIO=
N OF
INTEGRALS =
&nb=
sp; =
Virgin
-----------------------=
- =
&nb=
sp; =
Virgin
STARTING F=
ROM
TOTAL CROSS SECTIONS, REACTION CROSS SECTIONS AND Virgin
A SOURCE
SPECTRUM ALL OF WHICH ARE GIVEN IN TABULAR FORM WITH Virgin
LINEAR INT=
ERPOLATION
BETWEEN TABULATED POINTS ALL REQUIRED
Virgin
INTEGRALS =
CAN BE
DEFINED BY ANALYTICAL EXPRESSIONS INVOLVING
Virgin
NOTHING MO=
RE
COMPLICATED THAN EXPONENTIALS. THE INTEGRALS THAT Virgin
MUST BE
EVALUATED ARE OF THE FORM... =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
FLUX =
&nb=
sp; =
&nb=
sp; =
Virgin
---- =
&nb=
sp;
=
&nb=
sp; =
Virgin
(INTEGRAL =
EK TO
EK+1) (S(E)* EXP(-XCT(E)*Z)*DE) =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
REACTIONS<=
span
style=3D'mso-spacerun:yes'> =
&nb=
sp;
=
&nb=
sp; Virgin
---------<=
span
style=3D'mso-spacerun:yes'> =
&nb=
sp; =
&nb=
sp;
Virgin
(INTEGRAL =
EK TO
EK+1) (S(E)*XCR(E)*EXP(-XCT(E)*Z)*DE) =
Virgin
=
&nb=
sp; =
&nb=
sp;
=
Virgin
WHERE..
EK TO EK+1=
=3D
LONGEST ENERGY INTERVAL OVER WHICH S(E), XCT(E) AND Virgin
=
XCR(E) ARE ALL LINEARLY INTERPOLABLE. =
Virgin
S(E) =3D =
ENERGY
DEPENDENT WEIGHTING SPECTRUM =
Virgin
XCR(E) =3D REACTION CRO=
SS
SECTION =
&nb=
sp;
Virgin
XCT(E) =3D OPTICAL PATH=
LENGTH
(BASED ON TOTAL CROSS SECTION) =
Virgin
Z =
=3D MATERIAL THICKNESS =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
S(E), XCR(=
E) AND
XCT(E) ARE ALL ASSUMED TO BE GIVEN IN TABULAR Virgin
FORM WITH LINEAR
INTERPOLATION USED BETWEEN TABULATED POINTS. Virgin
IN OTHER W=
ORDS
BETWEEN TABULATED POINTS EACH OF THESE THREE IS Virgin
DEFINED BY=
A
FUNCTION OF THE FORM... =
&nb=
sp;
Virgin
=
=
&nb=
sp; =
&nb=
sp; =
Virgin
F(E)=3D((E=
-
EK)*FK+1 + (EK+1 - E)*FK)/(EK+1 - EK) =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
EACH OF TH=
ESE THREE
CAN BE CONVERTED TO NORMAL FORM BY THE
Virgin
CHANGE OF
VARIABLES.... =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
X=3D(E - 0=
.5*(EK+1
+ EK))/(EK+1 - EK) =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
IN WHICH C=
ASE X
WILL VARY FROM -1 (AT EK) TO +1 (AT EK+1) AND Virgin
EACH FUNCT=
ION
REDUCES TO THE NORMAL FORM... =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
F(X)=3D0.5=
*(FK*(1
- X) + FK+1*(1 + X)) =
&nb=
sp;
Virgin
=
=3D0.5*(FK+1 + FK) + 0.5*(FK+1 - FK)*X =
&nb=
sp; Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
BY DEFININ=
G THE
AVERAGE VALUE AND 1/2 THE CHANGE ACROSS THE
Virgin
INTERVAL.<=
span
style=3D'mso-spacerun:yes'> =
&nb=
sp; =
=
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
AVF=3D0.5*=
(FK+1 +
FK) =
&nb=
sp; =
&nb=
sp;
Virgin
DF=3D 0.5*=
(FK+1 -
FK) =
&nb=
sp; =
Virgin
DE=3D 0.5*=
(EK+1 -
EK) =
&nb=
sp; =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
EACH OF THE
THREE FUNCTIONS REDUCES TO THE SIMPLE FORM... =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
F(X)=3DAVF+DF*X =
&nb=
sp; =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
AND THE TWO
REQUIRED INTEGRALS REDUCE TO... =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
FLUX =
&nb=
sp; =
&nb=
sp; =
Virgin
----
=
&nb=
sp; =
&nb=
sp; Virgin
DE*EXP(-AV=
XCT*Z)
* (INTEGRAL -1 TO +1) =
&nb=
sp;
Virgin
((AVS+DS*X)*EXP(-DXCT*Z*X)*DX) =
&nb=
sp; =
Virgin
=
&nb=
sp;
=
&nb=
sp; =
Vi=
rgin
REACTION =
&nb=
sp; =
&nb=
sp;
Virgin
-------- =
&nb=
sp; =
&nb=
sp;
Virgin
DE*EXP(-AV=
XCT*Z)
* (INTEGRAL -1 TO +1) =
&nb=
sp;
Virgin
((AVS*AVXCR+(AVS*DXCR+AVXCR*DS)*X+DS*DXCR*X*X)*EXP(-DXCT*Z*X)*DX)
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
WHERE =
&nb=
sp; =
=
&nb=
sp; Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
AVXCT =3D AVERAGE VALUE OF THE TOT=
AL CROSS
SECTION =
Virgin
AVXCR =3D AVERAGE VALUE OF THE REA=
CTION
CROSS SECTION =
Virgin
AVS =3D AVERAGE VALU=
E OF THE
SOURCE =
&nb=
sp;
Virgin
DXCT =3D 1/2 THE CHANGE IN =
THE
TOTAL CROSS SECTION =
Virgin
DXCR =3D 1/2 THE CHANGE IN =
THE
REACTION CROSS SECTION
Virgin
DS =3D 1/2 TH=
E CHANGE
IN THE SOURCE =
&nb=
sp;
Virgin
DE =3D 1/2 TH=
E CHANGE
IN THE ENERGY =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
NOTE THAT =
IN
THIS FORM THE ENERGY ONLY APPEARS IN FRONT OF THE Virgin
INTEGRALS =
AND
THE INTEGRALS ARE EXPRESSED ONLY IN TERMS OF THE Virgin
TABULATED =
VALUES
OF S(E), XCT(E) AND XCR(E). IN PARTICULAR NO Virgin
DERIVATIVE=
S ARE
USED, SO THAT THERE ARE NO NUMERICAL INSTABILITY Virgin
PROBLEMS I=
N THE
VACINITY OF DISCONTINUITIES IN S(E), XCT(E) OR Virgin
XCR(E). IN=
DEED,
SINCE (EK+1 - EK) APPEARS IN FRONT OF THE INTEGRAL Virgin
POINTS OF =
DISCONTINUITY
AUTOMATICALLY MAKE ZERO CONTRIBUTION TO Virgin
THE
INTEGRALS. =
&nb=
sp; =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
THE REQUIR=
ED
INTEGRALS CAN BE EXPRESSED IN TERMS OF THE THREE Virgin
INTEGRALS =
IN
NORMAL FORM.... =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
F(A,N) =3D
(INTEGRAL -1 TO 1) (X**N*EXP(-A*X)*DX), N=3D0,1 AND 2. Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
THESE THREE
INTEGRALS CAN BE EVALUATED TO FIND... =
Virgin
=
&nb=
sp; =
&nb=
sp;
=
&nb=
sp; Virgin
N=3D0 =
&nb=
sp; =
&nb=
sp; =
Virgin
--- =
&nb=
sp; =
&nb=
sp; =
Virgin
F(A,0) =3D
(EXP(A)-EXP(-A))/A =
&nb=
sp;
=
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
N=3D1 =
&nb=
sp; =
&nb=
sp; =
Virgin
--- =
&nb=
sp; =
&nb=
sp;
Virgi=
n
F(A,1) =3D
((1-A)*EXP(A)-(1+A)*EXP(-A))/(A*A) =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
N=3D2 =
&nb=
sp; =
&nb=
sp; =
Virgin
--- =
&nb=
sp; =
&nb=
sp; =
Virgin
F(A,2) =3D
((2-2*A+A*A)*EXP(A)-(2+2*A+A*A)*EXP(-A))/(A*A*A) =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
HOWEVER THESE EXPRESSIONS ARE NUMER=
ICALLY
UNSTABLE FOR SMALL
Virgin
VALUES OF =
A.
THEREFORE FOR SMALL A THE EXPONENTIAL IN THE =
Virgin
INTEGRALS =
ARE
EXPANDED IN A POWER SERIES... =
&nb=
sp;
Virgin
=
=
&nb=
sp; =
&nb=
sp; =
Virgin
EXP(-AX)=3D1.0-(AX)+(AX)**2/2-(AX)**3/6+(AX)**4/24-........ =
Virgin
=
=3D(SUM K=3D0 TO INFINITY) (-AX)**K/(K FACTORIAL) =
Virgin
=
&nb=
sp;
=
&nb=
sp; =
&nb=
sp; Virgin
AND THE IN=
TEGRAL
REDUCES TO THE FORM.... =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
(SUM K=3D0=
TO
INFINITY) ((-A)**K/(K FACTORIAL)) * =
Virgin
(INTEGRAL =
-1 TO
1) (X**(N+K))*DX =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
WHICH CAN =
BE
ANALYTICALLY EVAULATED TO FIND.... =
Virgin
(K(N) =3D K
FACTORIAL) =
&nb=
sp; =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
N=3D0 =
&nb=
sp; =
=
Virgin
--- =
&nb=
sp; =
&nb=
sp; =
Virgin
F(A,0) =3D
2*(1+(A**2)/K(3)+(A**4)/K(5)+(A**6)/K(7)+.... =
Virgin
=
&nb=
sp; =
&nb=
sp; =
=
Virgin
N=3D1 =
&nb=
sp; =
&nb=
sp; =
Virgin
--- =
&nb=
sp; =
&nb=
sp; =
Virgin
F(A,1) =3D
-2*A*(2/K(3)+4*(A**2)/K(5)+6*(A**4)/K(7)+8*(A**6)/K(9)+.. Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
N=3D2 =
&nb=
sp; =
&nb=
sp; =
Virgin
--- =
&nb=
sp; =
&nb=
sp; =
Virgin
F(A,2) =3D
2*(2/K(3)+3*4*(A**2)/K(5)+5*6*(A**4)/K(7)+ =
Virgin
=
7*8*(A**6)/K(9)+.... =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
THESE EXPA=
NSIONS
ARE USED WHEN THE ABSOLUTE VALUE OF A IS LESS Virgin
THAN 0.1. =
BY
TRUNCATING THE ABOVE SERIES BEFORE A**8 THE ERROR Virgin
RELATIVE T=
O THE
LEADING TERM OF THE SERIES WILL BE 10**(-10), Virgin
YIELDING 10
DIGIT ACCURACY. =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
AFTER EVAL=
UATING
THE ABOVE FUNCTIONS, EITHER DIRECTLY OR BY USING Virgin
THE EXPANS=
ION
THE TWO REQUIRED INTEGRALS CAN BE WRITTEN AS... Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
FLUX =
&nb=
sp; =
&nb=
sp; =
Virgin
---- =
&nb=
sp;
=
&nb=
sp; Virgin
DE*EXP(-AVXCT*Z)*(AVS*F(A,0) + DS*F(A,1)) =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
REACTIONS<=
span
style=3D'mso-spacerun:yes'> =
&nb=
sp; =
=
&nb=
sp;Virgin
---------<=
span
style=3D'mso-spacerun:yes'> =
&nb=
sp; =
&nb=
sp;
Virgin
DE*EXP(-AVXCT*Z)* =
&nb=
sp; =
&nb=
sp;
Virgin
(AVS*AVXCR*F(A,0) + (AVS*DXCR+AVXCR*DS)*F(A,1) + DS*DXCR*F(A,2)) Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
INPUT FILE=
S =
&nb=
sp; =
&nb=
sp;
Virgin
----------=
- =
&nb=
sp; =
&nb=
sp;
Virgin
FILENAME UNIT DESCRIPTION =
&nb=
sp; =
Virgin
-------- ---- ----------- =
&nb=
sp; =
Virgin
INPUT 2
ENDFIN 10 EVALUATED DATA IN ENDF/B FOR=
MAT =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
OUTPUT
FILES =
&nb=
sp; =
&nb=
sp;
Virgin
----------=
-- =
&nb=
sp; =
&nb=
sp;
Virgin
FILENAME UNIT DESCRIPTION =
&nb=
sp; =
Virgin
-------- ---- ----------- =
&nb=
sp; =
Virgin
OUTPUT 3 OUTPUT LISTING =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
SCRATCH
FILES =
&nb=
sp; =
&nb=
sp;
Virgin
------------- =
=
&nb=
sp; =
Virgin
FILENAME UNIT DESCRIPTION =
&nb=
sp; =
Virgin
-------- ---- ----------- =
&nb=
sp; =
Virgin
SCR1 12 REACTION, FLUX AND CROSS SEC=
TION
RESULTS (BCD)
Virgin
=
&nb=
sp;
(SORTED AT END OF RUN AND OUTPUT SEPARATELY)
Virgin
SCR2 13 TALLY GROUP ENERGY BOUNDARIES
(BINARY) =
Virgin
SCR3 14 SOURCE SPECTRUM (BINARY) =
&nb=
sp; Virgin
SCR4 15 TOTAL CROSS SECTION (BINARY)=
=
&nb=
sp;
Virgin
SCR5 16 REACTION CROSS SECTION (BINA=
RY) =
Virgin
=
&nb=
sp; =
&nb=
sp; =
=
Virgin
OPTIONAL
STANDARD FILE NAMES (SEE SUBROUTINE FILIO1 AND FILEIO2) Virgin
---------------------------------------------------------------- Virgin
UNIT FILE NAME FORMAT =
&nb=
sp; =
Virgin
---- ---------- ------ =
&nb=
sp; =
Virgin
2 VIRGIN.INP BCD =
&nb=
sp; =
Virgin
3 VIRGIN.LST BCD =
&nb=
sp; =
Virgin
10 ENDFB.IN BCD =
&nb=
sp; =
Virgin
11-15 (SCRATCH) BINARY =
&nb=
sp; =
Virgin
16 PLOTTAB.CUR PLOTTAB OUTPUT F=
ORMAT
DATA =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
Virgin
INPUT LINE=
S =
&nb=
sp; =
&nb=
sp;
Virgin
----------=
- =
&nb=
sp; =
&nb=
sp;
Virgin
ANY NUMBER=
OF
CASES MAY BE RUN ONE AFTER THE OTHER. AFTER THE Virgin
FIRST CASE=
HAS
BEEN RUN THE FOLLOWING CASES MAY USE THE SAME
Virgin
THICKNESSE=
S,
GROUP STRUCTURE AND SPECTRUM AS THE PRECEDING CASE. Virgin
IN ADDITIO=
N THE
TRANSMITTED SPECTRUM FROM ONE CASE MAY BE USED Virgin
AS THE INC=
IDENT
SPECTRUM IN THE NEXT CASE, TO ALLOW MULTIPLE
Virgin
LAYERS OF
DIFFERENT MATERIALS. =
&nb=
sp; =
Virgin
=
&nb=
sp; =
=
&nb=
sp; =
Virgin
LINE COLS. FORMAT DESCRIPTION =
&nb=
sp; =
Virgin
---- ----- ------ ---------- =
&nb=
sp; =
Virgin
1 1-60 ENDF/B INPUT DATA FILENAME =
&nb=
sp; Virgin
=
&nb=
sp;
(STANDARD OPTION =3D ENDFB.IN) =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
LEAVE THE
DEFINITION OF THE FILENAMES BLANK - THE PROGRAM WILL Virgin
THEN USE
STANDARD FILENAMES. =
&nb=
sp; =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
2-3<=
span
style=3D'mso-spacerun:yes'> 1-72 18A4 TWO LINE TITLE DESCRIBING
PROBLEM Virgin
4 1- 6 I6 ZA (1000*Z+A) OF TARGE=
T FOR
TOTAL =
Virgin
=
7-11 I5
=
12-22 E11.4 DENSITY FOR TOTAL =
&nb=
sp;
Virgin
=
23-28 I6 ZA (1000*Z+A) OF TARGE=
T FOR REACTION =
Virgin
=
29-33 I5 MT OF REACTION =
&nb=
sp;
Virgin
=
&nb=
sp;
=3D 0 - NO REACTION CALCULATION (ONLY FLUX). Virgin
=
&nb=
sp;=3D
GREATER THAN 0 - CALCULATE REACTIONS.
Virgin
=
34-44 E11.4 DENSITY FOR REACTION =
&nb=
sp;
Virgin
=
45-50 I6 NUMBER OF TARGET THICK=
NESSES =
Virgin
=
=
=3D
GREATER THAN 0 =3D READ FROM INPUT =
Virgin
=
&nb=
sp;
(1 TO 2000 ALLOWED) =
&nb=
sp;
Virgin
=
&nb=
sp;
=3D 0 =3D SAME AS LAST=
CASE =
&nb=
sp;
Virgin
=
51-55 I5 NUMBER OF TALLY GROUPS =
&nb=
sp;
Virgin
=
&nb=
sp;
(REMEMBER NUMBER OF GROUP BOUNDARIES =
Virgin
=
&nb=
sp;
IS ONE MORE THAN THE NUMBER OF GROUPS)
Virgin
=
&nb=
sp;
UP TO 2000 GROUPS ARE ALLOWED =
Virgin
=
&nb=
sp;
BUILT-IN GROUP STRUCTURES. =
Virgin
=
&nb=
sp;
=3D GREATER THAN 0 =3D READ FROM INPUT =
Virgin
=
&nb=
sp;
=3D 0 SAME AS LAST CAS=
E =
&nb=
sp;
Virgin
=
&nb=
sp;
=3D -1 TART 175 GROUPS &=
nbsp; &nbs=
p;
Virgin
=
&nb=
sp;
=3D -2 ORNL 50 GROUPS<=
span
style=3D'mso-spacerun:yes'> =
&nb=
sp;
Virgin
=
&nb=
sp;
=3D -3 ORNL 126 GROUPS &=
nbsp;
=
Virgin
=
&nb=
sp;
=3D -4 ORNL 171 GROUPS &=
nbsp; &nbs=
p;
Virgin
=
&nb=
sp;
=3D -5 SAND-II 620 GROUPS..UP TO 18 MEV.
Virgin
=
&nb=
sp;
=3D -6 SAND-II 640 GROUPS..UP TO 20 MEV. Virgin
=
&nb=
sp;
=3D -7 WIMS 69 GROUPS &n=
bsp;  =
;
Virgin
=
&nb=
sp;
=3D -8 GAM-I 68 GROUPS &=
nbsp; &nbs=
p;
Virgin
=
&nb=
sp;
=3D -9 GAM-II 99 GROUPS =
&nb=
sp;
Virgin
=
&nb=
sp;
=3D-10 MUFT 54 GROUPS &n=
bsp;  =
;
Virgin
=
&nb=
sp;
=3D-11 ABBN 28 GROUPS &n=
bsp;  =
;
Virgin
=
56-60 I5 NUMBER OF POINTS IN SO=
URCE SPECTRUM =
Virgin
=
&nb=
sp; (MUST
BE AT LEAST TWO POINTS) &=
nbsp; &nbs=
p;
Virgin
=
&nb=
sp;
=3D GREATER THAN 1 =3D READ FROM INPUT =
Virgin
=
&nb=
sp;
=3D 0 =3D SAME AS LAST=
CASE =
&nb=
sp;
Virgin
=
=
=3D
-1 =3D CONSTANT (ENERGY INDEPENDENT) =
Virgin
=
&nb=
sp;
=3D -2 =3D 1/E &nb=
sp; =
&nb=
sp;
Virgin
=
&nb=
sp;
=3D -3 =3D BLACKBODY - PHOTON SPECTRUM =
Virgin
=
&nb=
sp;
=3D -4 =3D=
BLACKBODY
- ENERGY SPECTRUM =
Virgin
=
&nb=
sp;
=3D -5 =3D TRANSMITTED SPECTRUM FROM LAST CASE Virgin
=
&nb=
sp;
NOTE, ALL SPECTRA, EXCEPT THE TRANSMITTED Virgin
=
&nb=
sp;
SPECTRUM FROM THE LAST CASE, WILL BE =
Virgin
=
&nb=
sp;
NORMALIZED SUCH THAT ITS INTEGRAL OVER
Virgin
=
&nb=
sp;
ENERGY WILL BE UNITY. &n=
bsp;  =
;
Virgin
=
61-64 1X,3I1 SPATIALLY DEPENDENT OUTOUT =
Virgin
=
&nb=
sp;
=3D 0 =3D NO  =
; &n=
bsp;  =
;
Virgin
=
&nb=
sp;
=3D 1 =3D YES &nbs=
p; &=
nbsp; &nbs=
p;
Virgin
=
&nb=
sp;
FOR THE 3 QUANTITIES &nb=
sp;
=
Vi=
rgin
=
&nb=
sp;
COLUMN 67 FLUX &nb=
sp; =
&nb=
sp;
Virgin
=
&nb=
sp;
68 REACTIONS =
&nb=
sp;
Virgin
=
&nb=
sp;
69 AVERAGE CROSS SECTION =
Virgin
=
65-65 I1 ENERGY DEPENDENT OUTOU=
T =
&nb=
sp;
Virgin
=
&nb=
sp;
=3D 0 =3D NONE &nb=
sp; =
&nb=
sp;
Virgin
=
&nb=
sp;
=3D 1 =3D INCIDENT SPECTRUM =
&nb=
sp;
Virgin
=
&nb=
sp;
=3D 2 =3D TRANSMITTED SPECTRUM =
Virgin
=
&nb=
sp;
=3D 3 =3D INCIDENT REACTIONS =
&nb=
sp;
Virgin
=
&nb=
sp;
=3D 4 =3D TRANSMIITED REACTIONS =
Virgin
=
&nb=
sp; =3D
5 =3D TOTAL CROSS SECTION  =
; &n=
bsp;
Virgin
=
&nb=
sp;
=3D 6 =3D REACTION CROSS SECTION =
Virgin
5 1-11 E11.4 BLACKBODY TEMPERATURE IN EV =
Virgin
=
12-22 E11.4 FLUX NORMALIZATION =
&nb=
sp;
Virgin
=
23-33 E11.4 REACTION NORMALIZATION =
&nb=
sp;
Virgin
=
&nb=
sp;
CALCULATIONS WILL BE BASED ON THE SPECTRUM Virgin
=
&nb=
sp;
AND
CROSS SECTIONS AS READ. AT OUTPUT THE Virgin
=
&nb=
sp;
RESULTS WILL BE MULTIPLIED BY THESE =
Virgin
=
&nb=
sp;
NORMALIZATION FACTORS. &=
nbsp; &nbs=
p;
Virgin
=
34-44 I11<=
span
style=3D'mso-spacerun:yes'> DENSITY PROFILE =
&nb=
sp;
Virgin
=
&nb=
sp;
=3D 0 - UNIFORM - BASED ON TOTAL DENSITY
Virgin
=
&nb=
sp;
=3D 1 - UNIFORM - TOTAL + REACTION DENSITY Virgin
=
&nb=
sp;
=3D 2 - TOTAL + LINEAR REACTION =
Virgin
=
&nb=
sp;
=3D 3 - LINEAR (TOTAL + REACTION) =
Virgin
=
&nb=
sp;
=3D 4 - TOTAL + SQUARE REACTION =
Virgin
=
&nb=
sp;
=3D 5 - SQUARE (TOTAL + REACTION) =
Virgin
=
&nb=
sp;
=3D 6 - TOTAL + CUBIC REACTION =
Virgin
=
&nb=
sp;
=3D 7 - CUBIC (TOTAL + REACTION) =
Virgin
6-N<=
span
style=3D'mso-spacerun:yes'> 1-66 6E11.4 TARGET THICKNESSES IN CM =
=
Virgin
=
&nb=
sp;
IF SAME AS LAST CASE THIS SECTION IS NOT Virgin
=
&nb=
sp;
INCLUDED IN THE INPUT. &=
nbsp; &nbs=
p;
Virgin
VARY 1-66 6E11.4 TALLY GROUP ENERGY BOUNDARIES =
Virgin
=
&nb=
sp;
(NUMBER OF BOUNDARIES IS ONE MORE THAN
Virgin
=
&nb=
sp;
THE NUMBER OF TALLY GROUPS) =
Virgin
=
&nb=
sp;
IF THE STANDARD OPTION (-11 TO 0) IS =
Virgin
=
&nb=
sp;
SELECTED THIS SECTION IS NOT INCLUDED
Virgin
=
&nb=
sp;
IN THE INPUT  =
; &n=
bsp;  =
;
Virgin
VARY 1-66 6E11.4 SOURCE SPECTRUM IN ENERGY (EV)-SOU=
RCE
PAIRS Virgin
=
=
(MUST
BE AT LEAST TWO POINTS) =
Virgin
=
&nb=
sp;
IF STANDARD OPTION (-5 TO 0) IS SELECTED THIS Virgin
=
&nb=
sp;
SECTION IS NOT INCLUDED IN THE INPUT =
Virgin
=
=
&nb=
sp; =
&nb=
sp; Virgin
ANY NUMBER=
OF
CASES MAY BE RUN ONE AFTER ANOTHER. =
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
EXAMPLE IN=
PUT
NO. 1 =
&nb=
sp; =
Virgin
------------------- =
&nb=
sp; =
&nb=
sp;
Virgin
CALCULATE =
THE
UNCOLLIDED FLUX AND CAPTURE (MT=3D102) THROUGH
Virgin
30 CM OF I=
RON
(DENSITY 7.87 G/CC). TALLY THE RESULTS USING
Virgin
THE TART 1=
75
GROUP STRUCTURE. THE SOURCE WILL BE CONSTANT =
Virgin
FROM 1 KEV=
TO 20
MEV. USE THE STANDARD ENDF/B INPUT DATA =
Virgin
FILENAME.<=
span
style=3D'mso-spacerun:yes'> =
&nb=
sp; =
=
&nb=
sp; Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
ENDFB.IN =
&nb=
sp; =
&nb=
sp;
Virgin
IRON 0 TO =
30 CM
THICK. =
&nb=
sp;
=
Virgin
CONSTANT S=
OURCE
FROM 1 KEV TO 20 MEV. =
&nb=
sp;
Virgin
2600=
0 1 7.87000+ 0 26000 102 7.87000+ 0 2 0 2 1100 Virgin
0.00=
000+ 0
1.00000+ 0 1.00000+ 0 =
0 0.00000+00  =
;
Virgin
0.00=
000+00
3.00000+01 =
&nb=
sp; =
Virgin
1.00=
00E+03
1.0000E+00 2.0000E+07 1.0000E+00 =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
&nb=
sp;
Virgin
EXAMPLE IN=
PUT
NO. 2 =
&nb=
sp; =
&nb=
sp;
Virgin
------------------- =
&nb=
sp; =
&nb=
sp;
Virgin
CALCULATE =
THE
UNCOLLIDED PHOTON FLUX THROUGH A MIXTURE OF SILICON Virgin
AND IRON F=
OR 100
MEV PHOTONS INCIDENT. THE TRANSMISSION WILL BE Virgin
CALCULATED=
FOR
21 THICKNESSES VARYING BETWEEN 0 AND 1 CM. THERE Virgin
WILL BE ON=
LY 1
TALLY GROUP SPANNING A VERY NARROW ENERGY RANGE Virgin
NEAR 100 M=
EV,
AND THE SOURCE SPECTRUM WILL BE CONSTANT OVER THE Virgin
SAME ENERGY
RANGE. USE THE STANDARD ENDF/B INPUT DATA FILENAME Virgin
BY LEAVING=
THE
FIRST INPUT LINE BLANK. =
&nb=
sp;
Virgin
=
&nb=
sp;
=
&nb=
sp; =
Virgin
(THIS IS A=
BLANK
LINE TO USE THE STANDARD INPUT FILENAME) =
Virgin
100 MEV
PHOTONS =
&nb=
sp; =
&nb=
sp;
Virgin
SILICON + =
5 %
IRON =
=
&nb=
sp; =
Virgin
1400=
0 521 2.30000+ 0 26000 521 1.15000- 1 21 1 2 1000 Virgin
0.00=
000+ 0
1.00000+ 0 1.00000+ 0 =
1 0.00000+00  =
;
Virgin
0.00=
000+00
5.00000-01 1.00000+00 1.50000+00 2.00000+00 2.50000+00 Virgin
3.00=
000+00
3.50000+00 4.00000+00 4.50000+00 5.00000+00 5.50000+00 Virgin
6.00=
000+00
6.50000+00 7.00000+00 7.50000+00 8.00000+00 8.50000+00 Virgin
9.00=
000+00
9.50000+00 1.00000+01 =
=
Virgin
9.99=
000+ 7
1.00100+ 8 =
&nb=
sp; =
Virgin
9.99=
000+ 7
1.00000+ 4 1.00100+ 8 1.00000+ 4 =
&nb=
sp;
Virgin
=
&nb=
sp; =
&nb=
sp; =
Virgin
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D
Virgin