Modeling organic iron-binding ligands in a three-dimensional biogeochemical ocean model