anchor.diff:  * $Id: nwpw_unfold.f 19707 2010-10-29 17:59:36Z d3y133 $
   *
   
  -      subroutine nwpw_unfold(ni,r1,r2)
  +      subroutine nwpw_unfold(ni,unita,r1,r2)
   *     ======================================================= 
   *	This routine unfolds r2
   *     ======================================================= 
         implicit none
         integer ni
  +      real*8  unita(3,3)
         real*8  r1(3,*)
         real*8  r2(3,*)
   
  
positive.code: 	   
     24   IF ((f2(1)-f1(1)) .LE. (-0.5d0)) THEN
  *          WRITE (*,*) 'a1<', I, R2A1, DA(1)/2.0d0
             r2(1,i) = r2(1,i) + lattice_unita(1,1)
             r2(2,i) = r2(2,i) + lattice_unita(2,1)
             r2(3,i) = r2(3,i) + lattice_unita(3,1)
  
              f2(1) =  b(1,1) * r2(1,i)
       >            +  b(2,1) * r2(2,i)
  
negative.code:       real*8  a(3,3),b(3,3),volume
        integer i,j
  
  
  *     ***** Determine the unit lattice vectors and distances ******
        do j=1,3
        do i=1,3
          a(i,j) = unita(i,j)
        end do
        end do
  
  