I think I have identified the bug in grain resonstructions that causes
squared grain boundaries. It has to do with how I am applying periodic
boundary conditions maybe. I am in the process of fixing it but this
commit is broken code.