def plot_miri_filters(ax, wl, filters, names, ylim=(0.0, 1.0), leg=True):
    """
    Plot JWST/MIRI photometric filter bands.

    :param :py:obj:`Axis` ax: Matplotlib Axis
    :param array_like wl: Wavelengths
    :param array_like filters: Filter throughput curves
    :param list names: Names of filters
    :param tuple ylim: y-axis limits. Default `(0.0, 1.0)`
    :param bool leg: Add a legend to the plot. Default :py:obj:`True`

    """
    ax1 = ax.twinx()
    for i in range(len(names)):
        ax1.fill_between(wl, filters[:, (i)], label=names[i], alpha=0.3)
        ax1.set_ylabel('Filter Throughput', rotation=270, labelpad=25)
    if leg:
        leg = ax1.legend(loc=2, fontsize=16, ncol=1)
        leg.get_frame().set_alpha(0.0)
    ax1.set_ylim(ylim)
