Remove gravity from position initialization functions, gamma_wt = 0.0, work on Neumann BC