Directory: optical/om-data/
File:      totela-sys.readme (10 Jul., 2002)
********************************************


                              totela-sys.tgz

               Code for calculation of total, elastic and reaction
            cross sections at intermediate energies using systematics
                 (provided by T. Fukahori on July 10, 2002)
            **********************************************************

                         Tokio FUKAHORI and Koji NIITA
        Nuclear Data Center, Japan Atomic Energy Research Institute (Japan)



Content
-------
Code totela-sys calculates total, elastic scattering and reaction cross
sections, and angular distribution of elastic scattering for neutron- and
proton-induced reactions in the intermediate energy region from 20 MeV to 3GeV.
The approach is based on the systematics of Pearlstein [1], which was modified
to improve agreement with the available experimental data and the LA150
evaluation. The new  was developed to fill gaps in experimental data in the
intermediate energy region and was used for production of the JENDL High Energy
File.

Plots of the calculated results compared with experimental data and with the
LA150 evaluations are provided in the totela-sys.tgz archive.


Files
-----
totela-sys.tgz contains two subdirectories:

   - totela-sys-source: FORTRAN sources Makefile and two sample input  and
                        respective output files (altogether 25 files)
   - totela-sys-plots : set of plots comparing totela-sys predictions with
                        experimental data (7 files)


Format
------
Tarred and gzipped archive. On UNIX/Linux systems use

tar xvzf totela-sys.tgz
or
gunzip totela-sys.tgz
tar xvf totela-sys.tar

to explode the archive. On MS Windows WinZip should be used instead.


Input description
-----------------
Code requires a simple input file containing:

READ(*,*)  NFILE

   NFILE - name of the output file

READ(*,*)  IZ,IA,NINC,IOPT,NE,NANG

   IZ    - Z of the target
   IA    - A of the target
   NINC  - incident particle (1 for neutron, 2 for proton)
   IOPT  - systematics (0 default, 1 Niita, 2 Pearlstein)
   NE    - number of incident energies (default 55)
   NANG  - number of angles (default 180)

READ(*,*) (EIN(I),I=1,NE)

   EIN(I) - incident energies if NE different from 0

The output file is in the ENDF-6 format.


Example of an input file reads:
'sample1_fe056n.dat'
26,56,1,0,0,0


Platform
--------
PC running Red Hat 7.3 Linux, likely to run on other systems with 
FORTRAN compiler.

Testing
-------
Tested by M. Herman on PC running Red Hat 7.3 Linux with g77 compiler.


Reference
---------
[1] S. Pearlstein, Nucl. Sci. Eng. 95, (1987) 116.


Contact
-------
T. Fukahori, Nuclear Data Center, Japan Atomic Energy Research Institute (Japan)
e-mail: fukahori@ndc.tokai.jaeri.go.jp