In this work, we consider four f (R) gravity models -the Hu-Sawicki, Starobinsky, Exponential and Tsujikawa models -and use a range of cosmological data, together with Markov Chain Monte Carlo sampling techniques, to constrain the associated model parameters. Our main aim is to compare the results we get when Ω k,0 is treated as a free parameter with their counterparts in a spatially flat scenario. The bounds we obtain for Ω k,0 in the former case are compatible with a flat geometry. It appears, however, that a higher value of the Hubble constant H0 allows for more curvature. Indeed, upon including in our analysis a Gaussian likelihood constructed from the local measurement of H0, we find that the results favor an open universe at a little over 1σ. This is perhaps not statistically significant, but it underlines the important implications of the Hubble tension for the assumptions commonly made about spatial curvature. We note that the late-time deviation of the Hubble parameter from its ΛCDM equivalent is comparable across all four models, especially in the non-flat case. When Ω k,0 = 0, the Hu-Sawicki model admits a smaller mean value for Ω cdm,0 h 2 , which increases the said deviation at redshifts higher than unity. We also study the effect of a change in scale by evaluating the growth rate at two different wavenumbers k † . Any changes are, on the whole, negligible, although a smaller k † does result in a slightly larger average value for the deviation parameter b.