Program edcs

General program to calculate differential cross sections (DCS) from K matrices generated by scatstab.

This is a slightly modified version of a code provided by N. Sanna, Jul 1994. This is a reduced version of the EROTVIB code by A. Jain and D. J. Thompson (Comp. Phys. Com. 32, 367 (1984)) which calculated rotational transitions for the Td symmetry group only.

Input data records

  1. iuins, iukmat, iuplt
  2. iprint
  3. title
  4. lmax
  5. lbig
  6. nang
  7. If lbig < 0 then read lbigl, lbig
  8. neng, (enefnd(i), i = 1, neng)

Definition of the input variables

input unit containing blms (output from symgen).
if iukmat > 0 then this is the input unit containing K-matrices. The program searches this unit to find all K-matrices with the energy enefnd(i). If iukmat < 0 then the program reads from ABS(iukmat) previously computed A_L coefficients and generates the corresponding DCS.
output unit for differential cross section in form suitable for 2-D plotting program.
integer, print flag.
character string (LEN = 70), job title.
integer, maximum value of l to include from the K matrices.
integer, maximum value of l in the Legendre polynomial expansion of the calculated A_L coefficient (and DCS).
integer, number of theta angles where the DCS will be calculated.
integer, minimum value of l in the Legendre polynomial expansion to compute in this run. By using different values of lbigl and lbig it is possible to do the calculation in differnent sections and to get the full distribution by combining the results from the different sections by running the program with iukmat < 0.
integer, number of energies to do in this run.
real, values of the kinetic energies (in eV) to do in this run.