changed output slightly
Works for 1D *if* the unit cell is changed to 2 sites (that means all the embeddings need to be changed accordingly, which is done in order16rectangles_1D.dat).
How will this translate to 2D?  2x2 unit cell???