The code Figure13.py needs following input files:
phasediagram_landau_OBC_maxo=4_hermTrue_D0=2.1000_B0=12.0000_rtrunc=-1.0000.txt
phasediagram_landau_OBC_maxo=9_hermTrue_D0=2.1000_B0=12.0000_rtrunc=-1.0000.txt
These files contain map of the LLL-like structure validity (1 - valid, 0 - not valid).

They can be created by running
python phasediagram.py
The maxo can be changed by changing max_orbital in the code (see comments inside)

Then, Figure 13 is created by running
python Figure13.py

