An improved triangular form-based multiple flow direction algorithm for determining the nonuniform flow domain over grid networks