Refactor parameter syntax
Combine active cell and max value loops
Move bed shear stress arrays into roughness module