Terrestrial biosphere models that simulate both leaf dynamics and canopy photosynthesis are required to understand vegetation–climate interactions. A time-stepping scheme is proposed to simulate leaf area index, phenology, and gross primary production via climate variables. The method performs well on simulating deciduous broadleaf forests across the eastern United States; it provides a simplified and improved version of the growing production day model for use in land surface modeling.