Functional ergodic limits for occupation time processes of site-dependent branching Brownian motions in ℝ