def METHOD_NAME(state: State, variable, vmin=None, vmax=None):
    z = state.data_2d[variable]
    avg = vcm.zonal_average_approximate(state.data_2d.lat, z)
    vmin = None if vmin is None else float(vmin)
    vmax = None if vmax is None else float(vmax)
    avg.plot(x="time", vmin=vmin, vmax=vmax)
    state.tape.save_plot()