Add estimated width to mainstreams vector in m.swim.routing

[Michel Wortmann]

• move or copy the width estimation from m.swim.substats
• also make sure each subbasin only has one mainstream feature