anchor.diff:    i=i+1
    else if (trim(arg)=='-peratom') then
     peratom=.true.
  +  i=i+1
  + else if (trim(arg)=='-peratomx') then
  +    peratom=.true.
  +    jfirst=1
  +    jlast=1
  +  i=i+1
  + else if (trim(arg)=='-peratomy') then
  +    peratom=.true.
  +    jfirst=2
  +    jlast=2
  +  i=i+1
  + else if (trim(arg)=='-peratomz') then
  +    peratom=.true.
  +    jfirst=3
  +    jlast=3
     i=i+1
    else
     write(*,*)
  
positive.code: double precision,dimension(:,:),allocatable :: dip
  double precision,dimension(:,:,:),allocatable :: vel
  double precision,dimension(:),allocatable :: wsave
  
  width=10.d0
  skip=0
  
negative.code:   stop
   end if
  !IR
   do j=jfirst,jlast
    mean=sum(dip(j,:))/dble(nsteps)
    dip(j,:)=dip(j,:)-mean
    v=(0.d0,0.d0)
  