• Microscopic cross-sections are pre-generated for each isotope and reaction type. • Local concentrations of over 1000 nuclides are calculated using micro depletion. • Time-dependent decay heat power is calculated by summing the isotopic decay heat contribution. • The method is implemented in nodal code DYN3D and verified against Serpent 2.