MIGHT have fixed the issues with diffusion being done in sigma instead of z. (May have to switch to z to fix) The diffusion result still seemed wrong, added a fudge factor to account for overly strong diffusion.  Also, added plot3d_part(lag_structure,plot every nth particle, plot title)  last arg's are optional