still haven't worked out how to vary rho in load2 or do the mass in derivs