ADAS Subroutine argam

C
       function argam( l   , a   )                                       
c-----------------------------------------------------------------------
c
c  ****************** fortran77 function: argam   **********************
c
c  purpose:  calculates arggamma(l+1+i*a) where l is an integer not 
c            less than zero
c
c
c  subroutine:
c
c  input : (i*4)  l         = angula momentum >=0 
c  input : (r*8)  a         = real argument 
c                                           
c  output: (r*8)  argam     = arggamma(l+1+i*a)
c                               
c  routines:
c          none
c
c  author:  h. p. summers, university of strathclyde
c           ja7.08
c           tel. 0141-548-4196
c
c  date:    06/06/02
c
c  update:
c
C VERSION: 1.1                          DATE: 04-07-96
C MODIFIED: WILLIAM OSBORN
C               - FIRST VERSION.
C VERSION: 1.2                          DATE: 14-06-2000
C MODIFIED: Martin O'Mullane
C               - changed definition from real argam*8().
C
C VERSION: 1.3                          DATE: 19-12-01
C MODIFIED: Martin O'MULLANE
C               - Removed junk from > column 72.
C
C VERSION: 1.4                          DATE: 18-03-03
C MODIFIED: Hugh Summers
C               - Re-written and updated. 
c
c-----------------------------------------------------------------------
c-----------------------------------------------------------------------
      INTEGER             L
      REAL*8              A