c------------------------------------------------
c out_final
c output the coordinates into Final_conf
	subroutine out_final(n)
	implicit double precision(a-h, o-z)
        parameter(maxa=500)
	integer n
       common/matrix/x(maxa,3),bangle(maxa)
       common/bf_coords/f(maxa,3),v(maxa,3),u(maxa,3)
       common/ischeme/icircular,iele,ihydro,iclose,iknot
     +,neuler

       write(139,5002) ((x(mi,mii),mii=1,3),mi=1,n)
       if ((icircular.eq.1).and.(iclose.eq.1)) then
	  write(139,5002) ((f(mi,mii),mii=1,3),mi=1,n)
	  write(139,5002) ((v(mi,mii),mii=1,3),mi=1,n)
	  write(139,5002) ((u(mi,mii),mii=1,3),mi=1,n)
       endif

5002	format(g13.6,1x,g13.6,1x,g13.6)
	return
	end