READCOUL.ME : comments on the test files. last changed 7/7/95 (ARB) C========================================================================== These are FORTRAN 77 files for COUL90.FOR, SBESJY.FOR and RICBES.FOR Coulomb functions, spherical Bessel functions and Riccati-Bessel functions. (The '90' refers to a year NOT to FORTRAN 90) The subroutine COUL90 calculates Coulomb and spherical/cylindrical Bessels The other two subroutines are special cases. C========================================================================== Simple version: The driver program COUL90_S.FOR calls the test file COULSMPL.IN and creates the output file COULSMPL.OUT which should be identical to COULSMPL.REF The 1st line of input is text. Each line then contains : eta, x, xm, lrange, kfn in format 3F10.3 , 2I5 xm is the minimum L-value desired (usually zero) lrange is the number of additional L-values desired kfn = 0 gives Coulomb functions kfn = 1 gives spherical Bessel functions kfn = 2 gives cylindrical Bessel functions Where appropriate, calls to SBESJY and RICBES are also made. C========================================================================== Full test version: The driver program COUL90_T.FOR calls the test file COULTEST.IN and creates the output file COULTEST.OUT which should be identical to COULTEST.REF The 4-letter key words are COUL SBES BESS and BESI and signify the functions to be calculated, BESS being cylindrical Bessels and BESI being fractional order cylindrical Bessels. NX is the number of x-values for the given orders. NL is the number of orders for the given x-values. The 1st line of input is text. Each 3 [or 4] lines then contains : 4-letter key word NL, NX, eta, xm in format 1X, A4, 7X, 2I3, 2F10.2 NL integer L-values in free format NX real x-values in free format [text to annotate the output] if NX is negative xm is the minimum L-value desired (usually zero) lrange is the number of additional L-values desired kfn = 0 gives Coulomb functions kfn = 1 gives spherical Bessel functions kfn = 2 gives cylindrical Bessel functions Where appropriate, calls to SBESJY and RICBES are also made. The complete description of the programs is to be found in Chapter 9 of 'Computational Atomic Physics' ed. K Bartschat, which is called 'The Calculation of Spherical Bessel and Coulomb Functions' by A R Barnett C==========================================================================