goreduce appears to work except that I need to get pymol to guess the bonds in the new structure from its coordinates. There are several flags and print statements for debugging that must be taken away when the latter problem has been solved and goreduce is shown to work correctly