M_0 instead of M for channel flow problem
correction calcul poussée