Add latest SIR-compatible version of tracer-advection

Work on PSyclone and Dusk/Dawn now means that we can keep all of the compute in a single subroutine. However, local scalars must be replaced by their values and one loop nest that writes to `(ji,jj,jk+1)` must be transformed so as to write to `(ji,jj,jk)`.