Sounds like the simplest approach would be to create a new version of you script which reads in the precalculated dynamical matrix and hamiltonian derivateves, and then computes the IETS:
# Load device configuration
device_configuration = nlread('filename.nc', DeviceConfiguration)[-1]
# Load dynamical matrix
dynamical_matrix = nlread('filename.nc', DynamicalMatrix)[-1]
# Load hamiltonian derivatives
hamiltonian_derivatives = nlread('filename.nc', HamiltonianDerivatives)[-1]
# -------------------------------------------------------------
# Inelastic transmission spectrum
# -------------------------------------------------------------
inelastic_transmission_spectrum = InelasticTransmissionSpectrum(
configuration=device_configuration,
dynamical_matrix=dynamical_matrix,
hamiltonian_derivatives=hamiltonian_derivatives,
...
)