anchor.diff:  double precision,dimension(:,:),allocatable :: dip
   double precision,dimension(:,:,:),allocatable :: vel
   double precision,dimension(:),allocatable :: wsave
  +character*6 vdos_tag
  +integer jfirst,jlast
   
   width=10.d0
   skip=0
  
positive.code:   stop
   end if
  !IR
   do j=1,3
    mean=sum(dip(j,:))/dble(nsteps)
    dip(j,:)=dip(j,:)-mean
    v=(0.d0,0.d0)
  
negative.code:  end do
  !VDOS
   do i=1,nat
  !  do j=1,3
    do j=jfirst,jlast
     v=(0.d0,0.d0)
     v(1:nsteps)=vel(j,i,:)
     call zfftf(pn,v,wsave)
  