La principal característica del riego por gravedad es la forma en la que se distribuye el agua en el perfil del suelo y por lo tanto, son necesarios estudiar y desarrollar metodologías para el estudio del proceso. En este trabajo se presenta una solución numérica del modelo hidrodinámico de las ecuaciones de Barré de Saint-Venant para el flujo superficial, acopladas a la ecuación de Green y Ampt para modelar el flujo subsuperficial. Mediante la optimización de los parámetros de conductividad hidráulica a saturación y la presión de succión en el frente de humedad, parámetros característicos de la ecuación de Green y Ampt, se obtiene el gasto de riego óptimo necesario para un valor máximo de uniformidad. La solución de las ecuaciones de Barré de Saint-Venant se aproximan utilizando un esquema Lagrangiano en diferencias finitas