                  V I B I N I (diatomic codes only)

      VIBINI initializes the acquisition of vibrational wavefunctions.
      It is (optionally) called by VIBRMT or TMAT and is not called 
      directly in the main program.  The wave functions
      may either be read in from an external file or calculated from
      the closed form expressions for Morse functions.
      The format for direct input is consistent with the program of
      Leroy (1980) as modified by J Tennyson and LA Morgan (1985)

      If Morse functions are to be calculated the code requires the
      coefficients in the following parameterization of the potential

          V(R) =E0 + DE*(1-EXP(-BETA*(R-RE)))**2
 
      where E0 is the absolute energy in au and all other parameters 
      are given in au
      Alternatively the Morse potential can be specified in terms of
           W = OMEGAe,  A = OMEGAe/OMEGAeXe

      The following PARAMETERs are used to define array dimensions

      MAXTGT = 10  maximum number of target electronic states
      MAXVIB = 20  maximum number of target vibrational levels

      Namelist input /VIBIN/

        The format of the following is 
       (Name, Type, Dimension, Default, Limits, ! description)
        Unless specified otherwise, the default value applies to all elements
        of array variables.

        EBASE     R   1          0.   [:]
           ! Minimum of ground state potential energy curve.  If omitted,
             code uses minimum energy found in fixed nuclei data

        LUVIB     I   1          0    [0:]
           ! Logical unit from which vibrational wavefunctions are to be read
             If zero, then Morse functions are calculated

        LEROY     I   1          1    [0:]
           ! Indicates which version of Leroy code was used to generate
             wavefunctions

        DE        R   1               [:]
           ! Parameter in Morse potential (see above)

        BETA      R   1               [:]
           ! Parameter in Morse potential (see above)

        RE        R   1               [:]
           ! Parameter in Morse potential (see above)

        E0        R   1               [:]
           ! Parameter in Morse potential (see above)

        W         R   1               [:]
           ! Omega_e  (see above)

        A         R   1               [:]
           ! Omega_e/Omega_e X_e (see above)
