Routines to model ocean carbonate system thermodynamics