Consider 2D-Pencil Decomposition - Improve scalability for extreme scale simulations

- Currently, the code is setup to decompose the domain using y-z slabs which can severely limit the scalability when `N_cpu > Nx`
- Consider porting the MPI version of the code to 2D pencil decomposition for scalability