I think that should be a script, I don't know if it would be useful to have in shore.py. It is useful for theoretical comparison but not for the reconstruction itself. What do you think @Garyfallidis?