- 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?

for adaptive time stepping
and cleaning

the fields have been adapted to Hermite-Laguerre basis and kr,kz fourier space and the non useful lines and scripts have been deleted.