Here, we demonstrate and evaluate the Simulating Optimal FUNctioning (rsofun) R package, which provides a computationally efficient implementation of the P-model for site-scale simulations of ecosystem photosynthesis and the acclimation of photosynthetic traits, complemented with functionalities for Bayesian model-data integration and the estimation of model parameters and prediction uncertainty.