Understanding the nuances of the effects of nutrient limitation on oceanic primary production has been the focus of many bioassay experiments by oceanographers. A theme of these investigations is that they identify the currently limiting nutrient at a given location, or in other words they identify the proximate limiting nutrient (PLN). However, the ultimate limiting nutrient (ULN; the nutrient whose supply controls system productivity over extensive timescales) can be different from the PLN. Our motivation is to investigate the identity of the ULN. To facilitate this the <strong>c</strong>arbon-centric <strong>G</strong>rid <strong>En</strong>abled <strong>I</strong>ntegrated <strong>E</strong>arth system model (cGEnIE) <strong>nut</strong>rient cycles have been extended to create NutGEnIE. NutGEnIE incorporates the nutrients nitrogen, phosphorus, and iron. The impacts of diazotrophs, capable of fixing nitrogen, are represented alongside those of other phytoplankton. NutGEnIE is capable of extended model simulations necessary to investigate the ULN while, at the same time, including iron as a potentially limiting nutrient. NutGEnIE will be described, with particular focus on the biogeochemical cycles of iron, nitrogen and phosphorus. Model results are compared to ocean observational data to assess the degree of realism. Model-data comparisons include physical properties, nutrient concentrations, and process rates (e.g., export and nitrogen fixation). These comparisons support the conclusion that NutGEnIE is appropriate for the investigation of the ULN.