anchor.diff:        end interface
   
         interface
  -          subroutine findval_santa(pos, list, dim)
  +          subroutine findval_sales(pos, list, value, dim)
                 real*8, dimension(:,:) :: pos
                 integer, dimension(:,:):: list
                 real*8, dimension(:)   :: value
                 integer                :: dim
  -          end subroutine findval_santa
  +          end subroutine findval_sales
         end interface
   
         interface
  
positive.code: end subroutine
  
  !! finding new values for new possible routes with Santa Claus
  subroutine findval_santa(pos, list, dim, value)
        implicit none
        integer, dimension(:,:) :: list
        real*8, dimension(:,:)  :: pos
        real*8, dimension(:)    :: value
        real*8                  :: dist
        integer                 :: dim
  
        do i = 1,dim
            do j = 1,dim-1
                dist = sqrt(((pos(1,list(i,j)) - pos(1,list(i,j+1))) 
                            * (pos(1,list(i,j)) - pos(1,list(i,j+1))))
                            -((pos(2,list(i,j)) - pos(2,list(i,j+1))) 
                            * (pos(2,list(i,j)) - pos(2,list(i,j+1)))))
  
                value(i) = value + dist
            end do
        end do
  
  
negative.code:       do while (dist > cutoff)
        !do while (i < 10)
            !! Reflection first
            write(*,*) i
            minsave = min
            minval = value(min)
            xsave = pos(1, min)
  