program TEST_RANMAR

c  Author:  Glen Cowan
c  Date:  13 October, 1999    
c  Test program for RANMAR (CERNLIB V113) to generate uniform random numbers.

      implicit        NONE

c  Needed for HBOOK routines

      integer         hsize
      parameter       (hsize = 100000)

      integer         hmemor (hsize)
      common  /pawc/  hmemor

c  Local variables

      character*80    outfile

      integer         i
      integer         icycle
      integer         ierror
      integer         istat
      integer         lun
      integer         NTOTIN          !  0<=NTOTIN<=999 999 999
      integer         NTO2IN          !  0<=NTOT2N<<999 999 999
      integer         num_points
      integer         IJKLIN          !  0<=IJKLIN<=900 OOO OOO

      real            rvec(1)   ! vector of random nos.  (here 1 element only)

c  Initialize HBOOK, open histogram file, book histograms.

      call HLIMIT (hsize)
      lun = 20
      outfile = 'test_ranmar.his'
      call HROPEN (lun, 'histog', outfile, 'N', 1024, istat)
      if (istat .ne. 0) write (*, *) 'HROPEN error, istat = ', istat
      call HBOOK1 (1, 'uniform dist', 100, 0., 1., 0.)

      write (*, *) 'Enter initial seed between 0 and 900 000 000'
      read (*, *) IJKLIN
      NTOTIN = 0
      NTO2IN = 0
      call RMARIN(IJKLIN,NTOTIN,NTO2IN)

C  Generate values and enter into histogram.

      write (*, *) 'enter number of values to generate'
      read (*, *) num_points
      do i = 1, num_points
        call RANMAR (rvec, 1)
        call HF1 (1, rvec(1), 1.)
      end do

c  Store histogram and close.

      call HROUT (0, icycle, ' ')
      call HREND ('histog')

      stop
      END