In this study, we investigate the stationary states of the Glauber-Fock oscillator waveguide array. We begin by transforming the associated Hamiltonian into the form of a quantum harmonic oscillator Hamiltonian, allowing the implementation of a supersymmetric (SUSY) approach. By considering the simplest case for the intertwining operator, the optical ladder operators are straightforwardly constructed and shown to map eigensolutions into eigensolutions of the corresponding Hamiltonian operator, in pretty much the same manner as it is done for the quantum harmonic oscillator case. The ladder of the corresponding (eigen) supermodes is then easily established.
