Use smarter truncation for Fourier Modes

**Problem**
For the moment, the first $n_{i}$  fourier modes are truncated. 
While this works well for most diffusion like problems, this fairs poorly for advection dominated ones.

**Solution**
Perform an rfftshift and truncate the modes about the center of the array.