- introduce variables ngX_o2 = ngX/2 where X = p,j,z
- replace all ngX/2 call by a call to ngX_o2
- compiler should optimize it anyway but can we trust it?

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

- not working 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.

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