Problem in implementation of Shallow water equation
 

I am trying to implement shallow water equations as described in Finite-Volume Model for Shallow_water Flooding of Arbitrary Topography(Bradford and sanders). I got matlab code from the site of Sanders( http://sanders.eng.uci.edu/matlabcodes.html). For 2D it didn't take into account for bed elevation and wet/dry boundary. I extended the implementation to take account for bed elevation and wet/dry boundary. But the system gets unstable after sometime(Courant no. increases more than 1).

For Wet/Dry Boundary,
instead of F = 0.5(FL + FR - R*A*dW), equation for equating fluxes at the interface
i m doing , F = 0.5*(FL + FR), for extrapolation, is it correct?
and when height goes below 0.0, i m setting F = 0;

Is there anything else i need to do? Can anyone guide me how to apply conditions for wet/dry boundary? or provide me link to the code in which case of wet/dry boundary is implemented.

Thanks in advance

Cheers
Sumit

solved?
 

hi DeliriousPsycho,

did you manage to solve the problem with the wetdrying front? Currrently, I'm working on that problem, too. So I would be happy, if you could share your progress!

Best regards,

Andy