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

iuins
input unit containing blms (output from symgen).
iukmat
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.
iuplt
output unit for differential cross section in form suitable for 2-D plotting program.
iprint
integer, print flag.
title
character string (LEN = 70), job title.
lmax
integer, maximum value of l to include from the K matrices.
lbig
integer, maximum value of l in the Legendre polynomial expansion of the calculated A_L coefficient (and DCS).
nang
integer, number of theta angles where the DCS will be calculated.
lbigl
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.
neng
integer, number of energies to do in this run.
enefnd(i)
real, values of the kinetic energies (in eV) to do in this run.