
                E I G E N P

      EIGENP calculates eigen-phases from K-matrices.

      Input files:
        Channel and target data          (on unit LUCHAN)
        K-matrices                       (on unit LUKMT )

      Output files:
        Eigenphases                      (on unit LUPHAS)

      The following PARAMETERs are used to define array dimensions

      MAXERN = 10  maximum number of energy subranges (defined in RSOLVE)
                   
      Namelist input /PHASIN/

        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.
        

        EMIN     R   1         0.    [0.:MAXE]
           ! minimum required scattering energy (in units as specified by 
             IEUNIT).  MAXE is the maximum energy on the T-matrix file

        EMAX     R   1         MAXE  [0.:MAXE]
           ! maximum required scattering energy (MAXE is the maximum energy
             on the T-matrix file)

        ICFORM   C   1         'U'   ['U'|'F']
           ! format flag for channel dataset, 'F' for formatted, 'U' for
             unformatted

        ICON     I   1          0    [0,1,2]
           ! eigenphase smoothing switch
                 0 = no smoothing  
                 1 = add multiples of pi to minimize the difference between
                     successive energy points (only within one threshold)
                 2 = use second-order finite differences to ensure
                     continuity (not implemented in present version of codes)

        IEUNIT   I   1          1    [1,2]  
           ! units in which input scattering energies are input
                 1= Ryd, 2= eV

        IKFORM   C   1         'U'   ['U'|'F']
           ! Formatted/unformatted switch for unit LUKMT

        IPRNT    I   6          0    [0|1]
           ! Debug print switches :
             In each case -1 gives less than the default output and +1
             more than the default
                 (1) =1 Print all input data
                 (2) =1 Print eigenvectors corresponding to each eigenphase
                 (3)  not used
                 (4)  not used
                 (5)  not used
                 (6) =1 Print all output data

        IWRITE   I   1          6    [1:]        
           ! logical unit for printed output

        LUCHAN   I   1         10    [1:]        
           ! logical unit for input channel and target datasets

        LUKMT    I   1         11    [1:]
           ! logical unit holding input K-matrices

        LUPHSI   I   1          0    [1:]
           ! logical unit holding eigenphase sums (input, see below)

        LUPHSO   I   1          0    [1:]
           ! logical unit holding eigenphase sums in a form suitable
             for input to a plotting program.  If LUPHSI is non zero
             then the eigenphases already on it will be copied to
             LUPHSO and the current eigenphase sums appended as an
             additional column

        NAME     C   60                          
           ! title for any output

        NCHSET   I   1          1    [1:]        
           ! set number of required channel dataset. 

        NEPRNT   I   1          6    [0:]
           ! Number of eigenphases to be printed at each energy
             NEPRNT=0 give eigenphase sum only.

        NKSET     I   MAXPTS          1    [1:]
           ! Set numbers for input K-matrices 

        SNAME     C   8 
           ! Short name used as column heading on LUPHSO
