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

nu instead of NU
corrected a comment

The perpendicular pressure effects are already added above
so this commented part is useles