Trace elements and isotopes (TEIs) are important tools to study the changes in the ocean environment both today and in the past. However, the behaviors of TEIs in marine sediments are poorly known, limiting our ability to use them in oceanography. Here we present a modeling framework that can be used to generate and run models of the sedimentary cycling of TEIs assisted with advanced numerical tools in the Julia language, lowering the coding barrier for the general user to study marine TEIs.