import matplotlib.pyplot as plt
from tenpy.models import lattice
fig, axes = plt.subplots(2, 2, sharex=True, sharey=True, figsize=(8, 6))
orders = ['Cstyle', 'snakeCstyle', 'Fstyle', 'snakeFstyle']
lat = lattice.Square(5, 3, None, bc='periodic')
for order, ax in zip(orders, axes.flatten()):
    lat.order = lat.ordering(order)
    lat.plot_order(ax, linestyle=':', linewidth=2)
    lat.plot_sites(ax)
    lat.plot_basis(ax, origin=-0.25*(lat.basis[0] + lat.basis[1]))
    ax.set_title(repr(order))
    ax.set_aspect('equal')
    ax.set_xlim(-1)
    ax.set_ylim(-1)
plt.show()