Model code and output for "Decoupled timescales of organic carbon and phosphorus recycling in the global ocean"

Model code, output, and preprocessed input datasets supporting the publication: Sullivan, M.R., Primeau, F.W., Seo, H., Camps-Castellà, J., Inomura, K., & Martiny, A. (2026) Decoupled timescales of organic carbon and phosphorus recycling in the global ocean, PNAS, 123, e2514991123. https://doi.org/10.1073/pnas.2514991123The src directory contains the biogeochemical model that couples the cycling of phosphorus, carbon, and oxygen. The postproc directory contains all MATLAB scripts used to analyze the output and create the figures presented in the publication.The outputs of the model optimization and sensitivity runs are saved in individual .mat files in the output directory. The model output data include: model grid information, model parameters, and 3D fields of dissolved inorganic phosphorus, dissolved organic phosphorus, particulate organic phosphorus, dissolved inorganic carbon, dissolved organic carbon (labile, semi-labile, refractory), particulate organic carbon, alkalinity, particulate inorganic carbon, and oxygen.For each model solution, the sequestration flux partitioned by residence time is stored in the rtime_dist_output directory.