Remove unecessary module parameter
Use correct finite differencing for second order turbulence terms
Add mixing length model for eddy viscosity
Rename module and consolidate subroutines into one