We introduce a space inhomogeneous generalization of the dynamics on interlacing arrays considered by Borodin and Ferrari in [7]. We show that for a certain class of initial conditions the point process associated to the dynamics has determinantal correlation functions and we calculate explicitly, in the form of a double contour integral, the correlation kernel for one of the most classical initial conditions, the densely packed. En route to proving this we obtain some results of independent interest on non-intersecting general pure-birth chains, that generalize the Charlier process, the discrete analogue of Dyson's Brownian motion. Finally, these dynamics provide a coupling between the inhomogeneous versions of the TAZRP and PushTASEP particle systems which appear as projections on the left and right edges of the array respectively⋆ ⋆ + 1 ⋆ + 2 ⋆ + 3 ⋆ + 4 ⋆ + 5 ⋆ + 6 Figure 1: A configuration of particles in GT 4 . If the clock of the particle labelled X 3,1 4 rings, which happens at rate λ(⋆ + 1), then the move is blocked since interlacing with X 2,1 4 would be violated. On the other hand, if the clock of the particle X 2,2 4 rings, which happens at rate λ(⋆ + 3), then it jumps to the right by one and instantaneously pushes both X 3,3 4 and X 4,4 4 to the right by one as well, for otherwise the interlacing would break.The right edge particle system is called inhomogeneous PushTASEP, see [7], [6] and also [16] for a q-deformation of the homogeneous case. The fact that this particle system has an underlying determinantal structure is new as well, but is also obtained in the independent work of Petrov [46] that uses different methods.