No, my initial instinct seems to have been correct...
Starting from a Schwarzschild-like energy Lagrangian
$$\def\Rc{{\cal R}}
\def\Lc{{\cal L}}
\Lc ~:=~ g_{\mu\nu} u^\mu u^\nu
~=~ -\left(1 - \frac{2m}{r} \right) c^2 \dot t^2
~+~ \left(1 - \frac{2m}{r} \right)^{-1} \dot r^2
~+~ \rho r (\dot\theta^2 + \sin^2\!\theta\,\dot\phi^2)~,
~~~~~~~~~ \left[\, m := \frac{GM}{c^2} \,\right] ~,$$$$\mbox{where}~~ \rho ~=~ \rho(r,\Rc) ~,~~~
~~~\mbox{such that}~~ \rho \to r ~~\mbox{for}~ r \ll \Rc \;,
~~~\mbox{and}~~ \rho \to \Rc ~~\mbox{for}~ r \gg \Rc \;.$$The standard method for computing (timelike) geodesic orbits (circular with ##\ddot r = 0##) gives
$$0 ~=~ - c^2 \frac{m}{r^2} ~-~ \left( \frac{2m \rho}{r} + m \rho' \right) \,\dot\phi^2
~+~ \frac12 (\rho + r \rho') \dot\phi^2 ~,$$and there's no Tully-Fisher in there (which would need a relationship between ##\dot\phi^4## and ##M## for large ##r##).
Separately, I also found another issue with Arraut's paper. On p3 where he goes from ##V_1(r)## in eq(13) to ##\nabla V_1(r)## in eq(14) he seems to avoid differentiating ##\gamma## wrt ##r##, even though his ##\gamma## is r-dependent.
OK, enough of this.
