finished calcequildensity routine, calc velocities doesn't work properly yet