rocsy

The rocsy package is designed to accurately and efficiently compute all marine carbonate system variables from dissolved inorganic carbon (DIC), total alkalinity (Alk), temperature, and salinity as well as phosphate and silica concentrations. The pH scale and thermodynamic equilibrium constants used in rocsy are those recommended by Best-Practices Guide (Dickson et al., 2007). Also included are options to use more recent formulations for K1 and K2 from Millero (2010) and for total boron from Lee et al. (2010).

The rocsy package is designed for users with some programming experience, preferably in Fortran or python. Others may well prefer the facility of packages such as the Excel version of CO2sys. But rocsy is computationally much faster, being designed for big data sets and ocean models. For models, rocsy offers improved accuracy by avoiding 3 common approximations. Computed variables generally agree with those from CO2sys within 0.005% (Orr and Epitalon, 2013).