The DynACof model was designed to model coffee agroforestry systems and study the trade-offs to e.g. optimize the system facing climate changes. The model simulates net primary productivity (NPP), growth, yield, mortality, energy and water balance of coffee agroforestry systems according to shade tree species and management. Several plot-scale ecosystem services are simulated by the model, such as production, canopy cooling effect, or potential C sequestration. DynACof uses metamodels derived from a detailed 3D process-based model (MAESPA) to account for complex spatial effects, while running fast. It also includes a coffee flower bud and fruit cohort module to better distribute fruit carbon demand over the year, a key feature to obtain a realistic competition between sinks. We compared the model outputs with a highly comprehensive database on a coffee agroforestry farm in Costa Rica. The fluxes simulated by the model were close to the measurements over a 5-year period (RMSE= 1.60 gC m-2 d-1 for gross primary productivity; 0.63 mm d-1 for actual evapo-transpiration, 1.34 MJ m-2 d-1 for sensible heat flux and 1.88 MJ m-2 d-1 for net radiation), and DynACof satisfactorily simulated the yield, NPP, mortality and carbon stock for each coffee organ type over a 35-year rotation.