Starting from the Bethe Ansatz solution of the open integrable spin-1 XXZ quantum spin chain with diagonal boundary terms, we derive a set of nonlinear integral equations (NLIEs), which we propose to describe the boundary supersymmetric sine-Gordon model BSSG+ with Dirichlet boundary conditions on a finite interval. We compute the corresponding boundary S matrix, and find that it coincides with the one proposed by Bajnok, Palla and Takács for the Dirichlet BSSG+ model. We derive a relation between the (UV) parameters in the boundary conditions and the (IR) parameters in the boundary S matrix. By computing the boundary vacuum energy, we determine a previously unknown parameter in the scattering theory. We solve the NLIEs numerically for intermediate values of the interval length, and find agreement with our analytical result for the effective central charge in the UV limit and with boundary conformal perturbation theory.