Sébastien,
Good news, and I am glad to hear it. It would just be the addition of one new variable in the tubes data file to define a weir coefficient; and a small change to BUSE.F to include that new variable. Hopefully an easy update.
I don't know if you want me to publish my edit to the BUSE.F. I removed much of BUSE.F's true function for a single purpose that suited my needs. But I will describe what I did for others to mimic if they choose.
Instead of applying a standard weir equation, I have 3rd order equations for each spillway defined as:
Q=CE1*S1^3+CE2*S1^2+CS1*S1+CS2
Where CE1,CE2,CS1,CS2 are the variables read in on the existing BUSE.F routine from the existing format for the Tubes Data File.
BUT: In doing this I've lost the ability for the subroutine to check for backflow or oriface flow. This was OK for my work as all my spillways are on dams with no chance of backwater effects; and I have chosen to neglect pressurized flow.