do ie = 1, nelem
       do i = 1, ngp
          lim_inf_x = fe_x(1,lnods(1,ie))
          lim_inf_y = fe_x(2,lnods(1,ie))
          lim_sup_x = fe_x(1,lnods(3,ie))
          lim_sup_y = fe_x(2,lnods(3,ie))

          if( (fd_x(1,i) >= lim_inf_x) .and. (fd_x(2,i) >= lim_inf_y) .and. &
              (fd_x(1,i) <= lim_sup_x) .and. (fd_x(2,i) <= lim_sup_y) )     &
          then
             fd_fe(i) = ie
             !write(*,*) i, fd_fe(i), fd_x(1,i), fd_x(2,i)
          end if
      end do
    end do