c******************************************************* c out_juxtapair(ntime) c c np600=1 output the distance vs. time directly c subroutine out_juxtapair(ntime) implicit double precision(a-h, o-z) parameter(maxpair=100,maxa=500,maxjudge=10) common /juxtapair/ npair,npairs(maxpair,2),nflag(maxpair,maxjudge) +,njudge, sjuxta_judge(maxjudge),mp(maxpair) common/matrix/x(maxa,3),bangle(maxa) do i=1,npair rr=0.d0 do j=1,3 rr=rr+(x(npairs(i,1),j)-x(npairs(i,2),j))*(x(npairs(i,1),j) + -x(npairs(i,2),j)) enddo write(mp(i),700) npairs(i,1),npairs(i,2),ntime,dsqrt(rr) enddo 700 format(i5,1x,i5,1x,i8,1x,g11.5) return end