Using Shallow Water Equations
 

I am trying to use the shallowWaterFoam Solver to extend it to a bigger area (more than 25 meters) to calculate height fields with varying boundary velocity.

However whatever deltaT I set the program produces a floating point error after some time. When I reduce the deltaT, the the execution increases by a little number of loops and still crashes.

This is because the velocity increases more than the amount I assume when calculating the Courant Number.

WHat can I do??
Can someone please help?