program main

      implicit none

      include 'const.inc'
      integer i, j, nevts, counter, ranno
      parameter(nevts=1e2)
      real xlam, HRNDM1, vect(5), const, xinit, yinit, plr,xtest
      double precision x0, y0, z0, bx, by, bz
c
c     -------------------------
c     fill the ntuple 'spin.rz' //// ntuple variables ....
      integer max_paw, iquest
      integer istat, icycle
      parameter(max_paw = 10000000) ! PAW common block size
      real quest, p
      common/quest/iquest(100)
      common/pawc/p(max_paw)
c
c     time1=time()
c     ranno=mod(time1,2000000)
c     ranno=mod(time1,1089926518)
      print *,"asking for seed"
      read(*,*) ranno;
c     ranno=seedc()
c     print *, "ranNo is ",ranno, time1
      call hlimit(max_paw)
      iquest(10) = 6500000
c
c-----------------------------------------------------------------------
      call rluxgo(3,ranno,0,0)  ! RANLUX luxury level 3
c-----------------------------------------------------------------------
c
      do 100 i = 1, nevts
c
         call ranlux(vect,3)
         x0=vect(1)*9.5-4.5
         y0=vect(2)*4-2
         z0=vect(3)*19-9.5

         call mag_bx(x0,y0,z0,bx)
         call mag_by(x0,y0,z0,by)
         call mag_bz(x0,y0,z0,bz)
         write(*,*) x0,y0,z0, bx, by, bz

 100  continue

      return
      end