anchor.diff:      use m_utilities
       type(tSurfElement), intent(inout) :: surf(:)
       integer :: i, j, ii,jj, nd1, nd2, nsurf
  +    integer :: k, oldsize
  +    integer, pointer :: dumarray(:)
  +
       nsurf = size(surf)
       do i=1,nsurf
         do ii=1, size(surf(i)%nodes)
  
positive.code: !                                                                      !
  !      Module Name : lib                                               !
  !                                                                      !
  !                    Written by X. YUAN                                !
  !                                                                      !
  !      Contact address :  IIS,The University of Tokyo, CISS            !
  !                                                                      !
  
negative.code: 
    end function next_position
    
  end module