1. Bug fixed in compute_force_mesh() near Nyquist frequency. 2. Changed initial velocity step a_v=da from 0.5*da. This gives significantly better matter power spcetrum, reporoducing the original serical code.