If you do something like:
step = pm.NUTS()
trace = pm.sample(step=step, tune=0)
It basically turn off all the tuning and just use the bare bone implementation.
However, I think over the years the NUTS sampler improves in many small details and the implementation might not be exactly the same as Algorithm 6 on the paper any more.