SPECTER: Radiation Damage Calculations ---------------------------------------- L.R. Greenwood, Battelle Pacific Northwest National Laboratory April 2002 The input is as follows (all free format): l. Title (80 characters) 2. ITYP,ISIG,IGP,IPKA,ACNM,TIME (Free Format) ITYP = 0 IF NO FLUX UNCERTAINTIES, = 2 OTHERWISE ISIG = 1 TO PRINT COLLAPSED CROSS SECTIONS, =0 TO STOP, = 2 For 1-group Thermal Neutron (2200 m/s) IGP = 0 IF GROUP DIFFERENTIAL FLUX (n/cm2-MeV, = output) = 1 IF TRUE DIFFERENTIAL FLUX (n/cm2-MeV) = 2 IF TRUE GROUP FLUX (n/cm2) IPKA = 0 TO CALCULATE PKA DISTRIBUTIONS, = 1 NO PKA'S ACNM = ACTIVITY NORMALIZATION (USUALLY = 1.0 UNLESS FILE IS GENERATED BY OTHER CODES, eg STAYSL) TIME = TIME IN SECONDS TO GET FLUENCE AND DPA, etc. ANY NORM CAN BE USED HERE IF DESIRED 3. NPT (FREE) = NUMBER OF POINTS OR GROUPS 4. ENERGIES (FREE) = LOW TO HIGH IN MEV 5. FLUXES(FREE) = FLUX VALUES NOTE: IF GROUP FLUXES, THEN YOU NEED ONE MORE ENERGY THAN FLUX VALUE TO DEFINE THE UPPER LIMIT OF THE LAST ENERGY BIN. OPTIONAL: 6. COVARIANCE MATRIX (FREE) ONLY HALF OF MATRIX IS ENTERED SINCE MATRIX MUST BE SYMMETRIC. FIRST ENTRY OF EACH LINE IS THE DIAGONAL TERM. Example: -------- HFIR PTP POSITION 100 MW 0 0 0 1 1.0000E+00 9.386E+6 100 1 1.000E-10 1.000E-09 1.000E-08 2.300E-08 5.000E-08 7.600E-08 1.150E-07 1.700E-07 2.550E-07 3.800E-07 5.500E-07 8.400E-07 1.275E-06 1.900E-06 2.800E-06 4.250E-06 6.300E-06 9.200E-06 1.350E-05 2.100E-05 3.000E-05 4.500E-05 6.900E-05 1.000E-04 1.350E-04 1.700E-04 2.200E-04 2.800E-04 3.600E-04 4.500E-04 5.750E-04 7.600E-04 9.600E-04 1.275E-03 1.600E-03 2.000E-03 2.700E-03 3.400E-03 4.500E-03 5.500E-03 7.200E-03 9.200E-03 1.200E-02 1.500E-02 1.900E-02 2.550E-02 3.200E-02 4.000E-02 5.250E-02 6.600E-02 8.800E-02 1.100E-01 1.350E-01 1.600E-01 1.900E-01 2.200E-01 2.550E-01 2.900E-01 3.200E-01 3.600E-01 4.000E-01 4.500E-01 5.000E-01 5.500E-01 6.000E-01 6.600E-01 7.200E-01 7.800E-01 8.400E-01 9.200E-01 1.000E+00 1.200E+00 1.400E+00 1.600E+00 1.800E+00 2.000E+00 2.300E+00 2.600E+00 2.900E+00 3.300E+00 3.700E+00 4.100E+00 4.500E+00 5.000E+00 5.500E+00 6.000E+00 6.700E+00 7.400E+00 8.200E+00 9.000E+00 1.000E+01 1.100E+01 1.200E+01 1.300E+01 1.400E+01 1.500E+01 1.600E+01 1.700E+01 1.800E+01 1.900E+01 2.000E+01 1.248E+21 9.270E+21 1.994E+22 2.193E+22 1.660E+22 9.215E+21 3.570E+21 1.093E+21 3.994E+20 2.551E+20 1.731E+20 1.223E+20 8.586E+19 5.630E+19 3.706E+19 2.475E+19 1.709E+19 1.184E+19 7.959E+18 5.340E+18 3.681E+18 2.441E+18 1.631E+18 1.158E+18 8.840E+17 6.938E+17 5.405E+17 4.232E+17 3.336E+17 2.646E+17 2.043E+17 1.579E+17 1.224E+17 9.455E+16 7.549E+16 5.841E+16 4.460E+16 3.442E+16 2.648E+16 2.058E+16 1.573E+16 1.233E+16 9.753E+15 7.859E+15 6.154E+15 4.816E+15 3.951E+15 3.234E+15 2.646E+15 2.173E+15 1.801E+15 1.565E+15 1.387E+15 1.248E+15 1.137E+15 1.059E+15 9.996E+14 9.440E+14 8.684E+14 7.833E+14 7.069E+14 6.575E+14 6.416E+14 6.323E+14 6.232E+14 6.016E+14 5.616E+14 5.186E+14 4.764E+14 4.353E+14 3.898E+14 3.441E+14 3.085E+14 2.765E+14 2.484E+14 2.203E+14 1.833E+14 1.466E+14 1.170E+14 9.011E+13 6.825E+13 5.260E+13 3.947E+13 2.835E+13 2.050E+13 1.403E+13 9.140E+12 5.911E+12 3.764E+12 2.325E+12 1.065E+12 3.211E+11 1.362E+11 6.138E+10 2.634E+10 1.088E+10 4.755E+09 2.185E+09 1.223E+09 9.596E+08 Logical Unit Assignments ------------------------- The following is used for the VAX computer: $ASSIGN NAME.IN FOR005 input file (as above) $ASSIGN PKAMIN.DAT FOR010 pka file $ASSIGN SIGD.DAT FOR015 dpa file $ASSIGN MACKLIB.DAT FOR017 Kerma data file $ASSIGN GAS157.DAT FOR018 Gas Production File ENDF 533 $ASSIGN COMPOUND.DAT FOR019 Limited compound dpa's $ASSIGN NIGAS.DAT FOR020 NI58 AND NI59 CROSS SECTIONS $ASSIGN FEASTMDPA.DAT FOR021 FE ASTM E693-01 DPA CROSS SECTION $ASSIGN NBGIN.DAT FOR023 (n,g) and beta data file $ASSIGN NAME.OUT FOR006 output file (DAM.OUT on a PC) $RUN SPECTER $PRINT NAME.OUT to run on a VAX type in: @SPECTER NAME(your input file name) OR : SUBMIT SPECTER.COM/PAR=NAME For PC Operation: All the above assignments are embedded in the code. Run by double clicking the executable file (specter.exe) or put a shortcut to this file on your desktop. The code will ask you for the name of your input file. Output automatically is saved to a file DAM.OUT Rename if you want to save it! (file overwritten each time) Notes: ------- l. If PKA's aren't needed, then no PKA file need be specified and dpa's can be computed for all 38 elements by using SIGD in unit 15 This greatly reduces the running time of the code and is strongly recommended! 2. If PKA's are needed, then the file PKAMIN will provide recoils for all 38 cases. At present this file only gives the net recoil distribution. Spectra for each sub reaction requires another version of SPECTER (much larger). 3. Note that in the above command procedure unit assignments are made at execution time. On some computers, it may be necessary to make such assignments in the program via open statements or by JCL or other control language. 4. Compound dpa cross sections can be added for other materials using the SPECOMP computer code. In the present file, displacement threshold energies were the same as for the elements (see manual); however, this may not be a valid assumption. The file compound.txt describes the compound cross sections currently available.