## MOLSCAT v. 14

### Bugfixes

The following problems have been identified since the original
release of MOLSCAT, version 14, in Aug 1994. Affected calculations are
briefly described, along with the action taken to correct the problem. The
"version date" (`PDATE`

) reported by the program has not
been changed unless indicated.

- 1 Nov 94: base.f
`ITYPE=4`

cases incorrectly calculated parity for assignment to parity subblocks.

Corrected and replaced base.f in v14.f and v14.f.Z- 10 Mar 95: gaushp.f
- Gauss-Hermite routine provided incorrect weight if only a single
integration point was requested

Corrected and replaced gaushp.f in v14.f and v14.f.Z - 21 Mar 95: EUNITS
- Allowing the NAMELIST variable EUNITS to accept either an integer
or a character string value, which worked on some earlier machines, rarely
works on
modern computers. Subroutines DRIVER, BASE, ECNV, ECNVX, and STSRCH have
been modified so that EUNITS is restricted to integer values, and a new
CHARACTER*4 NAMELIST variable EUNITC has been added, with a default
value of blanks. A nonblank input value for EUNITC will override any
value input in EUNITS. The version date (
`PDATE`

) has been changed.

Modified routines DRIVER, BASE, ECNV, ECVNX, and STSRCH have been put into v14.f and v14.f.Z and are also available separately here.

Section 2 and Section 3 of the documentation have been modified. - 28 Mar 95: findrm.f
- In trying to estimate a starting distance for the propagation findrm.f
sometimes encounters difficulties. Often these are due to an interaction
potential which is poorly behaved in the classically forbidden region.
In cases where a classically allowed region was encountered and an inner
turning point could not be located or where the origin was reached before
a suitable staring point was found, the original
algorithm chose a starting value of R=0. but this is rarely desirable.
The new routine chooses 0.8 times the input RMIN value for these cases
which generally is reasonable.

Corrected and replaced findrm.f in v14.f and v14.f.Z - 31 Jul 95: thrj.f
- Problems with the S-matrix post-processing program for
differential
scattering cross sections were traced to a faulty initialization of
the factorials in thrj.f.

Function thrj.f has been corrected and replaced in v14.f and v14.f.Z - 3 Aug 95: storag.f
- A bug was found in allocating storage for the WKB integrator in cases
which use numerical derivatives.

Subroutine storag.f has been corrected and replaced in v14.f and v14.f.Z - 30 Aug 95: Driver.f, output.f, restrt.f
- This was not a bug, but these routines were modified so that the program
would give the number of channels and number of open channels as part
of the summary line printed for each set of coupled equations.

The modified routines are available separately and have also been replaced in v14.f and v14.f.Z - 5 Sept 95: couple.f, restrt.f, sig6.f
- Corrected bugs associated with unitialized variables in some error condition
write statements. None of these should have caused erroneous results.

The modified routines have been replaced in v14.f and v14.f.Z - 21 Sept 95: plm.f
- It was pointed out by Dr. Marie-Lise Dubernet <mld@ipp-garching.mpg.de>
that the algorithm used to obtain associated Legendre polynomials
was inaccurate and caused problems for
high values of the indices (higher than typically used by MOLSCAT).
She suggested an alternate routine which has been tested and
incorporated into the code.

**NOTE:**the version posted 28 Jun 95 had bugs. If you don't have a version dated 21 Sept 91, you should update.

New plm.f has been incorporated in v14.f and v14.f.Z - 19 Oct 95: potenl.f
- Some changes were made to get the print out to reflect correctly whether and how the projection mechanism was used. It is believed that actual calculations were not affected by these bugs, but only the print out. New potenl.f has been incorporated in v14.f and v14.f.Z

Return to MOLSCAT Homepage