integer i, j, jj, k, n, repno
      parameter(n=1500)
      double precision x0, y0, z0
      double precision x(n), y(n), z(n), b(n), tmp
      double precision dr(8),rat(8),bp(8), rattot, bmrat, bave
      integer xrep, yrep, zrep, flag
      parameter(xrep=3)
      parameter(yrep=1)