anchor.diff:  
     end function next_position
     
  -end module
  +end module
positive.code:             nd2 = surf(j)%nodes(jj)
              if( nd1==nd2 ) then
                surf(i)%n_neighbor = surf(i)%n_neighbor+1
                call insert_int2array( j, surf(i)%neighbor )
              endif
            enddo
          enddo		
  
negative.code: !                                                                      !
  !      Module Name : lib                                               !
  !                                                                      !
  !                    Written by Xi YUAN                                !
  !                                                                      !
  !      Contact address :  IIS,The University of Tokyo, CISS            !
  !                                                                      !
  