At Poznan University of Technology, we have homemade Matlab scripts for extracting data both from devices and BNDF, and then to specify boundary conditions in Abaqus input file or by Fortran subroutine.
Using devices instead of BNDF is useful, when size of the structural element is not comparable to the grid size. This I use for steel structures, see
https://www.sciencedirect.com/science/article/pii/S0263823121005759?via%3Dihub.
Similarly, when you consider surfaces that do not coincide with rectangular grid, devices allow to capture heat exposure from any direction.
Usually we adjust the scripts for particular case studies.