Channel flow models

implemented switch CASE for each type of surface model:
Schumann, Moeng, Moeng non-lin, Wei+Schumann, Wei+Moeng
Channel flow new approach

Using a decomposition of global+deviation for the sfc stress, where the
global part is evaluated from global momentum balance and the deviations
are still using Wei+Schumann form