c*******************************************************
c init_rand(n)
c
      subroutine init_rand(n,api)
      implicit double precision(a-h, o-z)
      parameter(pi=3.1415926535d0,maxa=500,maxn=maxa*3)
	integer n,n3
      common/param/temp,bendc,r0,beadr,ss1,interval,dtime,gamma,
     +amass,am(maxn),coef,dtime2
      common/myrandom/amean,var(maxn),fvar(maxa)
       common/diffu_tensor/diffu(maxn,maxn),w(maxn,maxn),varr(maxn) 
      common/tor_param/dr,cg,phi0,dlk,phii,subseg
      
         n3=n*3

	 amean = 0.d0
         do i = 1, n3
            var(i) = dtime*2/api
         enddo

          vartmp = 2.d0*dtime
         do i = 1, n3
          varr(i) = vartmp
         enddo

          dr=1.5/api/beadr/beadr
          vartmp = 2.d0*dtime*dr
         do i = 1, n
            fvar(i) = vartmp
         enddo

      return 
      end