Hi Nadi,
You can extract node coordinates at boundaries using BlueKenue. With the mesh selected, do: Tools>T3 Mesh> Extract Edges (Shorelines). This will only give you the coordinates. However,I believe node IDs are assigned sequentially so you can use that as a base to build a code to find the IDs. Note that the boundary conditions file contains the node IDs of your boundary under NBOR. I'm unsure if you're using the 2D mesh or a 3D results file, however I've used the boundary IDs as an index to find the values at those nodes of a 3D results file using something similar to:
IPOIN3 = (IPLAN-1)*NPOIN2 + INDEX.
You might use this to find the XY.
Laballeuer explains this in more detail here:
www.opentelemac.org/index.php/community-...a-for-certain-layers