Combine active cell and max value loops
Move bed shear stress arrays into roughness module
Refactor VINTERP function
Update vegetation input and module for new Cauchy approach
Fix style and move constant values to module parameters
Remove unecessary modules and add TODO list