c===========================================================
c     Appends a double precision vector to file 'fname'.
c===========================================================

      subroutine dvappend(fname,v,n)
c-----------------------------------------------------------
c     Arguments:
c
c        fname:  (I)    File name
c        v:      (I)    Vector to be written
c        n:      (I)    Length of vector
c-----------------------------------------------------------
         implicit         none

         integer          getu,       indlnb

         character*(*)    fname
         integer          n
         real*8           v(n)

         integer          i,      uto,       rc
         character*1      dummy

         uto = getu()
         open(uto,file=fname(1:indlnb(fname)),
     &        form='formatted',iostat=rc)
         if( rc .ne. 0 ) then
            write(0,*) 'dvto: Error opening ',
     &              fname(1:indlnb(fname))
            return
         end if
c-----------------------------------------------------------
c        Position file at end-of-file
c-----------------------------------------------------------
         do while( .true. ) 
            read(uto,fmt='(a)',end=100) dummy
         end do
 100     continue

         do i = 1 , n
            write(uto,*) v(i)
         end do

         close(uto)

         return

      end