anchor.diff:      integer :: node1, node2
       REAL(SP), dimension(ndim) :: distance
       REAL(SP), dimension(2**ndim) :: weights, values
  +    REAL(SP), dimension(2) :: denominator
       integer :: i, j
       integer :: npoint
   
  
positive.code:     integer :: node1, node2, node3, node4, node5
      REAL(SP), dimension(ndim) :: distance
      REAL(SP), dimension(2**ndim) :: weights, values
      integer :: i, j, k, l, m
      integer :: npoint
  
  
negative.code:     integer :: node1, node2,node3,node4,node5,node6,node7
      REAL(SP), dimension(ndim) :: distance
      REAL(SP), dimension(2**ndim) :: weights, values
      REAL(SP), dimension(7) :: denominator
      integer :: i, j, k, l, m, n, o
      integer :: npoint
  
  