Support `OrdinalRange`s in `ndrange` argument

I think it'd be nice if one could specify e.g.
```julia
kernel(u, v, ndrange=(2:2:size(u, 1), 2:2:size(u,2)))
```
to get kernel indices that don't necessarily start at 1 or iterate in unit steps. This'd also fit better with the name `ndrange`