Re: Preliminary advice on finite volume methods

Well, there are many reasons why I dont want to split river links:
- one of these is that the hydraulic module of the simulation program uses these links, and I'd rather use the same links for advection (i did not wrote the hydraulic module, but it is robust, stable, tested...)

besides this, there are some other reasons: - the simulation program uses as input a Digital Elevation Model of a watershed and this is made of very easy_life_making square cells - I dont know anything about mesh refining - the time available to work on the code is almost up - etc...

With my big surprise in 1&1/2 days I managed to set up a simple LaxWendroff+VanLeerlimiter in 1D with constant velocity, and, with my big surprise, it is truly flux-conservative! even with low Courant numbers! Almost magic...
I'll try to insert u(x) in the solution and once I have it working, I'll try to put it into the real code.

Thanks for the help and the suggestions!
A
.