Directory: optical/om-get/ File: om-retrieve.readme (July 15, 2002) ********************************************* om-retrieve.tgz Code for retrieving optical model potentials from the RIPL-2 library (provided by P. Young on July 15, 2002) ************************************************************** Content ------- The om-retrieve.for code makes inputs from any potential in the OMP library. The user selects whether SCAT2000 or ECIS96 inputs are to be made, although some potentials are incompatible with one or the other of the codes. The SCAT2000 code only deals with spherical potentials, whereas inputs for coupled-channels rotational, vibrational, and DWBA calculations can be made for ECIS96. On the other hand, om-retrieve.for will only make a SCAT2000 input when a Gaussian form factor is used for Wd. Note that the code puts vibrational or DWBA ECIS96 inputs out on separate files. The user must input reference numbers from the omp library and energies at which the inputs are made. There is a default set of energies that can be used, but these are extensive (0-200 MeV) and are principally useful for making tables of potential values when the ECIS96 (non-DWBA) option is chosen. The code also contains a simple option for making inputs from every potential in the library, which is useful for testing. Additionally, the user can direct the code to zero all negative values of the Wv and Wd potentials. Code should be compiled using: f77 om-retrieve.for -C -O3 -o omretrieve That is, -r8 option should NOT be turned on in order for the dispersive corrections to be computed accurately. Files ----- om-retrieve.tgz contains: - om-retrieve.for : FORTRAN source - om-retrieve.cmb : common block - gs-mass-sp.dat : ground-state masses, spins and parities Input/Output files ------------------ Input omparameter.dat : RIPL-2 Optical Model Potential (OMP) file (om-parameter-u.dat or om-parameter-a.dat) ominput.inp : input file (see input description below) Output sc2.inp : input file for the scat2000 code (modtyp=1) ec96.inp : standard input file for the ecis96 code (modtyp=2) ec96tc.inp : alternate input file for the ecis96 code. (modtyp=2). Useful for generating transmission transmission coefficients. ec96vib : input file for the ecis96 code with vibrational model activated (modtyp=2). ec96dw.inp : input file for the ecis96 code with DWBA model activated (modtyp=3). massinfo.out : descriptive remarks about the ground-state mass data used in the gs-mass-sp.dat file. omp-table : a table of numerical values of each potential made at each incident energy when modtyp=2 (ECIS input). Format ------ Tarred and gzipped archive. On UNIX/Linux systems use tar xvzf om-retrieve.tgz or gunzip om-retrieve.tgz tar xvf om-retrieve.tar to explode the archive. On MS Windows WinZip should be used instead. Input description (ominput.inp file) ------------------------------------ read(5,*) ne ne = number of incident energies = 0 to use a built-in array of incident energies if(ne.gt.0) read(5,*) (en(n),n=1,ne) en(n) = incident energies in MeV in the laboratory system read(5,*,end=990) iztar, iatar, irefget, modtyp abs(iztar) = Z of target nucleus (set iztar negative to provide integer projectile mass in input decks) abs(iatar) = A of target nucleus (set iatar negative to provide integer target mass in input decks) irefget = reference number of optical model potential to be retrieved from the RIPL-2 library. Set irefget = negative izaproj (projectile) to retrieve all spherical potentials in the RIPL-2 library for this izaproj and the inputted izatar (target). abs(modtyp)= 1 to generate SCAT2 input file (sc2.inp) 2 to generate ECIS96 input files (ec96.inp, ec96tc.inp) 3 to generate ECIS96 DWBA input files (ec96dw.inp), using structure information from an external file (deform.dat), which may be the om-deformations.dat*file or a user-provided deform.dat file. if(abs(modtyp).eq.3) read(5,*) kdef Set modtyp negative to force negative imaginary volume and surface potentials to equal zero. If modtyp = -1, then all SCAT2000 inputs are single energy, i.e., the compact energy representation of SCAT2000 is bypassed. kdef = 1 use JENDL-3.2 betas from om-deformations.dat file = 2 use ENSDF(Q) = 3 use ENSDF(BE2) and ENSDF(BE3) = 4 use Raman and Spear for BE2,BE3 = 5 use deform.dat file provided by user See om-deformations.readme file for more information on the kdef options and for the format of deform.dat for using kdef=5. For kdef=1-4, the om-deformations.dat file must be copied to deform.dat. Note that where multiple entries occur for a single state in om-deformations.dat, the code makes an ECIS96 input for each entry. Platform -------- SUN Microsystem: Sunblade 100 computer, Solaris 8 operating system. Testing ------- Code was tested by P. Young while checking the phenomenological OMP of the RIPL-2 library. R. Capote performed testing of the dispersive potentials. Contact ------- P.G.Young, Group T-2, Los Alamos National Laboratory Mail Stop B243, e-mail: pgy@lanl.gov Coding for dispersive potentials and gamma factor by R. Capote Noy, Centro de Estudios Aplicados al Desarrollo Nuclear, Cuba, e-mail: rcapotenoy@yahoo.com