Indeed, setting the level is not ideal, unless you extend your model much furthr dowstream.
Also, you can appreciate that 1D river-taylored models such as HEC-RAS are much simpler than 2D generic models such as TELEMAC-2D, which makes the coding of certain boundary types mush easier, given the 1D assumption.
Nonetheless, if you were to give me the theoretical principles of this boundary type found in HEC-RAS, I would be happy to provide you with an example of implementation of the SL subroutine in TELEMAC-2D.
Further, you can also program a free boundary in a number of ways. For instance, one could set the velocity of the node at your outlet boundary to be the velocity one or two row of triangles inside the mesh, assuming flow will always goes one way.
With the source code available, there is always a way ...
Hope this helps.
Sebastien.