- not working yet

-the routines have been tested with the CBC
-looks good, the ExB shear is not tested yet

-FFTW parallelized in ky and x
-local_nkx has been replaced by local_nx

-1D routines development
-NL factor to correct moving grid method
-still an issue with the 1D r2c FFT
-ExB not working yet.

Better default input

- kx was not 2D but it is now since ExB shear implementation

-removed the parameters only modules
-simplified the names of the input namelists

-it is made now in the increasing or decreasing kx order
-the order is decided with the sign of the jump (= sign of the ExBrate)

also improved the way to init from image

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

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