docs: def plot_data():
      '''Plot sample data up with the fancy colormaps.
  
      '''
  
      var = ['temp', 'oxygen', 'salinity', 'fluorescence-ECO', 'density', 'PAR', 'turbidity', 'fluorescence-CDOM']
      # colorbar limits for each property
      lims = np.array([[26, 33], [0, 10], [0, 36], [0, 6], [1005, 1025], [0, 0.6], [0, 2], [0, 9]])  # reasonable values
      # lims = np.array([[20,36], [26,33], [1.5,5.6], [0,4], [0,9], [0,1.5]]) # values to show colormaps
  
      for fname in fnames:
          fig, axes = plt.subplots(nrows=4, ncols=2)
          fig.set_size_inches(20, 10)
          fig.subplots_adjust(top=0.95, bottom=0.01, left=0.2, right=0.99, wspace=0.0, hspace=0.07)
          i = 0
          for ax, Var, cmap in zip(axes.flat, var, cmaps):  # loop through data to plot up
  
              # get variable data
              lat, lon, z, data = test.read(Var, fname)
  
              map1 = ax.scatter(lat, -z, c=data, cmap=cmap, s=10, linewidths=0., vmin=lims[i, 0], vmax=lims[i, 1])
              # no stupid offset
              y_formatter = mpl.ticker.ScalarFormatter(useOffset=False)
              ax.xaxis.set_major_formatter(y_formatter)
              if i == 6:
                  ax.set_xlabel('Latitude [degrees]')
                  ax.set_ylabel('Depth [m]')
              else:
                  ax.set_xticklabels([])
                  ax.set_yticklabels([])
              ax.set_ylim(-z.max(), 0)
              ax.set_xlim(lat.min(), lat.max())
              cb = plt.colorbar(map1, ax=ax, pad=0.02)
              cb.set_label(cmap.name + ' [' + '$' + cmap.units + '$]')
              i += 1
  
          fig.savefig('figures/' + fname.split('.')[0] + '.png', bbox_inches='tight')