Adjunto imagen con resultado.
El procedimiento realizado fue, al raster del factor les resto a cada valor -1. Luego realizo el flujo de dirección, para los valores MDT, dados.
Y luego multiplico cada valor dado, entre cuadrículas. En las esquinas al dar valor 0, directamente en el flujo de dirección no lo tengo en cuenta.