-dirty (lots of uncommented lines)
-unoptimized (could remove aux var)
-ExB factor still not working

- 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