A numerical routine for the seasonal evolution of open-ocean water column temperature for use in biogeochemical studies