Orthogonal reduce field that applies reduction function `fun` over slices along orthogonal dimensions of a field iterator.
r = fun(r, fields[i][index]); reduction by i
non empty slice
Orthogonal reduce field that applies reduction function `fun` over slices along orthogonal dimensions of a field iterator.