pub fn idx3d(i: usize, j: usize, k: usize, nx: usize, ny: usize) -> usize
Index into a 3D array stored in Fortran order (column-major) index = x + ynx + znx*ny