subroutine out_juxtaed(n,i,id)
	implicit double precision(a-h, o-z)
        parameter(maxa=500,maxn=3*maxa)
	integer n,ntime
       common/matrix/x(maxa,3),bangle(maxa)
       common/param/temp,bendc,r0,beadr,ss1,interval,dtime,gamma,
     +amass,am(maxn),coef,dtime2
      common/bf_coords/f(maxa,3),v(maxa,3),u(maxa,3)

	write(100,*) n,i,id
	
       write(100,*) ((x(mi,mii),mii=1,3),mi=1,n)
       write(100,*) ((f(mi,mii),mii=1,3),mi=1,n)
       write(100,*) ((v(mi,mii),mii=1,3),mi=1,n)
       write(100,*) ((u(mi,mii),mii=1,3),mi=1,n)

	return
	end