Added the TFIM 2D code to this repo.  It will be
rewritten to so it used the Heisenberg Hamiltonian
for an incomplete *bilayer* system