docs: def _plot_rprof_list(sdat, lovs, rprofs, metas, stepstr, rads=None):
      """Plot requested profiles"""
      if rads is None:
          rads = {}
      for vfig in lovs:
          fig, axes = plt.subplots(ncols=len(vfig), sharey=True)
          axes = [axes] if len(vfig) == 1 else axes
          fname = 'rprof_'
          for iplt, vplt in enumerate(vfig):
              xlabel = None
              for ivar, rvar in enumerate(vplt):
                  fname += rvar + '_'
                  rad = rads[rvar] if rvar in rads else rprofs['r']
                  if conf.rprof.depth:
                      rad = rprofs['bounds'][1] - rad
                  axes[iplt].plot(rprofs[rvar], rad,
                                  conf.rprof.style,
                                  label=metas[rvar].description)
                  if conf.rprof.depth:
                      axes[iplt].invert_yaxis()
                  if xlabel is None:
                      xlabel = metas[rvar].kind
                  elif xlabel != metas[rvar].kind:
                      xlabel = ''
              if ivar == 0:
                  xlabel = metas[rvar].description
              if xlabel:
                  _, unit = sdat.scale(1, metas[rvar].dim)
                  if unit:
                      xlabel += ' ({})'.format(unit)
                  axes[iplt].set_xlabel(xlabel)
              if vplt[0][:3] == 'eta':  # list of log variables
                  axes[iplt].set_xscale('log')
              axes[iplt].set_xlim(left=conf.plot.vmin, right=conf.plot.vmax)
              if ivar:
                  axes[iplt].legend()
          ylabel = 'Depth' if conf.rprof.depth else 'Radius'
          _, unit = sdat.scale(1, 'm')
          if unit:
              ylabel += ' ({})'.format(unit)
          axes[0].set_ylabel(ylabel)
          misc.saveplot(fig, fname + stepstr)