docs: 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])