anchor.diff:        list_node%id_lagrange(1) = id_lagrange
         countNon0_lagrange = countNon0_lagrange + 1
       else
  +      allocate(id_lagrange_save(num_lagrange))
         id_lagrange_save(1:num_lagrange) = list_node%id_lagrange(1:num_lagrange)
         location = find_locationINarray(id_lagrange,num_lagrange,list_node%id_lagrange)
         if(location /= 0)then
  
positive.code:         if(location /= num_lagrange) list_node%id_lagrange(location+1:num_lagrange) = id_lagrange_save(location:num_lagrange-1)
          countNon0_lagrange = countNon0_lagrange + 1
        endif
      endif
  
    end subroutine insert_lagrange
  
negative.code:     integer(kind=kint) :: countNon0_lagrange !< counter of node-based number of non-zero items
      !< in Lagrange multiplier-related matrix
      integer(kind=kint) :: ierr, num_lagrange, location
      integer(kind=kint), allocatable :: id_lagrange_save(:)
  
      character(len=1)   :: answer
  
  