import matplotlib.pyplot as plt
from tenpy.models import lattice
fig, axes = plt.subplots(2, 2, sharex=True, sharey=True, figsize=(6, 8))
lat = lattice.Honeycomb(3, 3, None, bc='periodic')
order_names=['rings', 'Cstyle', 'snake', 'snakeCstyle']
for order_name, ax in zip(order_names, axes.flatten()):
    lat.plot_coupling(ax, linestyle='-', linewidth=1)
    lat.order = lat.ordering(order_name)
    lat.plot_order(ax, linestyle=':', linewidth=2)
    lat.plot_sites(ax)
    lat.plot_basis(ax, origin=(-0.5, -0.5))
    ax.set_title(f"order={order_name!r}")
    ax.set_xlim(-0.5, 3)
    ax.set_ylim(-0.5, 4)
    ax.set_aspect('equal')
axes[0, 0].set_title("order='default'='rings'")
plt.show()