tried 3d sphere (super slow and Im[k] goes in wrong direction), added Nz, Mz, Npmlz read from options