Se presenta la estructura de un algoritmo de solución para las ecuaciones de aguas pandas escritas en forma bidimensional y aplicadas a la propagación de ondas en flujos a superficie libre. El modelo propuesto considera todos los términos en las ecuaciones de aguas pandas y en particular incluye las aceleraciones de carácter convectivo en la ecuación de conservación del momentum lineal. La solución numérica de este modelo utiliza el método de las diferencias finitas.