Function padded_dims 

fn padded_dims(dims: [usize; 3], pad: [usize; 3]) -> [usize; 3]

Padded grid: pad[a] voxels on both sides of axis a, plus one at the high end where that is odd.