def plot_q_flows(network):
    """Plot maximal reactive line load. 
    
    Parameters
    ----------
    network: PyPSA network container
        Holds topology of grid including results from powerflow analysis

    """
    cmap_line = plt.cm.jet
    
    q_flows_max = abs(network.lines_t.q0.abs().max()/(network.lines.s_nom))
    
    ll = network.plot(line_colors = q_flows_max, line_cmap = cmap_line)
    boundaries = [min(q_flows_max), max(q_flows_max)]
    v = np.linspace(boundaries[0], boundaries[1], 101)

    cb = plt.colorbar(ll[1], boundaries=v,
                      ticks=v[0:101:10])

    cb.set_clim(vmin=boundaries[0], vmax=boundaries[1])