The consideration of marine biogeochemistry is essential for simulating the carbon cycle in an Earth system model. Here we present the implementation and evaluation of a marine biogeochemical model, the Model of Oceanic Pelagic Stoichiometry (MOPS) in the Flexible Ocean and Climate Infrastructure (FOCI) climate model. FOCI-MOPS enables the simulation of marine biological processes, i.e. the marine carbon, nitrogen, and oxygen cycles with prescribed or prognostic atmospheric CO2 concentration. A series of experiments covering the historical period (1850-2014) were performed following the DECK (Diagnostic, Evaluation and Characterization of Klima) and CMIP6 (Coupled Model Intercomparison Project 6) protocols. Overall, modelled biogeochemical tracer distributions and fluxes, transient evolution in surface air temperature, air-sea CO2 fluxes, and changes in ocean carbon and heat contents are in good agreement with observations. Modelled inorganic and organic tracer distributions are quantitatively evaluated by statistically derived metrics. Results of the FOCI-MOPS model, including sea surface temperature, surface pH, oxygen (100-600 m), nitrate (0-100 m), and primary production, are within the range of other CMIP6 model results. Overall, the evaluation of FOCI-MOPS indicates its suitability for Earth climate system simulations.