Fixed CFL condition for Darcy - u_max was squared, so the original term sqrt(u_max) was actually correct - changed so that u_max is squared when it is calculated to avoid confusion