anchor.diff:        max = max_array(1)
   end subroutine
   
  -!! populates grid -- UNPURE, fix later!
  +!! populates grid
   subroutine populate(pos, dim)
         implicit none
         integer                :: dim
  
positive.code: end subroutine
  
  !! Finds the initial list variable set
  subroutine pop_list(list, dim, pos)
        implicit none
        integer, dimension(:,:) :: list
        integer                 :: dim, i, j, k, swap, tmp
        integer, parameter      :: seed = 1
        real*8                  :: var, theta
        real*8, dimension(:,:)  :: pos
  
        !! Initialization!
        do i = 1, dim
            do j = 1,dim - 1
                list(i,j) = j
            end do 
  
negative.code:       end interface
  
        interface
            subroutine contract_sales(weights, min, dim, beta)
                real*8, intent(inout)  :: weights(:,:)
                real*8, intent(in)     :: beta
                integer, intent(in)    :: dim, min
            end subroutine contract_sales
        end interface
  
        interface
            subroutine expand(pos, min, dim, gamma)
                real*8, intent(inout)  :: pos(:,:)
                integer, intent(in)    :: min, dim
  