-kx is now 2D (ky dependant)
-kperp^2 array for computation reduction
-routines to update the grids