SUBROUTINE setup

IMPLICIT NONE
!=========================================================================================
! This routine obtains initial information need by the propagation and analysis routines.
!   Each routine's purpose is explained below.
!
!  Variables:
!=========================================================================================
!   I N P U T

!=========================================================================================
!   O U T P U T

!=========================================================================================
!    Obtain surface functions for initial wave function and analysis

     CALL basis_functions
!=========================================================================================
!    Calculate and symmetrize initial wave function to be propagated

     CALL initial_wavepacket
!=====================================================================
!  Calculate Bessel functions 

     CALL bessel_coef
!=========================================================================================
!    Calculate Negative Imaginary Potential

     CALL nip 

RETURN
END SUBROUTINE setup