An integrated modeling workflow capable of finding the steady-state plasma solution with self-consistent core transport, pedestal structure, current profile, and plasma equilibrium physics has been developed and tested against a DIII-D discharge. Key features of the achieved core-pedestal coupled workflow are its ability to account for the transport of impurities in the plasma self-consistently, as well as its use of machine learning accelerated models for the pedestal structure and for the turbulent transport physics. Notably, the coupled workflow is implemented within the One Modeling Framework for Integrated Tasks (OMFIT) framework, and makes use of the ITER integrated modeling and analysis suite data structure for exchanging data among the physics codes that are involved in the simulations. Such technical advance has been facilitated by the development of a new numerical library named ordered multidimensional arrays structure.