acceleration# acceleration(gradH, rho0, hbaroverm) ACCEL Summary of this function goes here Detailed explanation goes here