old_contents: """
  Plot MEG inverse solution
  =========================
  
  Data were computed using mne-python (http://martinos.org/mne)
  
  """
  import os
  import numpy as np
  
  from surfer import Brain
  from surfer.io import read_stc
  
  print(__doc__)
  
  """
  define subject, surface and hemisphere(s) to plot
  """
  subject_id, surface = 'fsaverage', 'inflated'
  hemi = 'split'
  
  """
  create Brain object for visualization
  """
  brain = Brain(subject_id, hemi, surface, size=(800, 400))
  
  """
  label for time annotation in milliseconds
  """
  
  
  def time_label(t):
      return 'time=%0.2f ms' % (t * 1e3)
  
  """
  read MNE dSPM inverse solution
  """
  for hemi in ['lh', 'rh']:
      stc_fname = os.path.join('example_data/meg_source_estimate-' +
                               hemi + '.stc')
      stc = read_stc(stc_fname)
  
      """
      data and vertices for which the data is defined
      """
      data = stc['data']
      vertices = stc['vertices']
  
      """
      time points (in seconds)
      """
      time = np.linspace(stc['tmin'], stc['tmin'] + data.shape[1] * stc['tstep'],
                         data.shape[1])
  
      """
      colormap to use
      """
      colormap = 'hot'
  
      brain.add_data(data, colormap=colormap, vertices=vertices,
                     smoothing_steps=10, time=time, time_label=time_label,
                     hemi=hemi)
  
  """
  scale colormap and set time (index) to display
  """
  brain.set_data_time_index(2)
  brain.scale_data_colormap(fmin=13, fmid=18, fmax=22, transparent=True)
  
  """
  uncomment these lines to use the interactive TimeViewer GUI
  """
  # from surfer import TimeViewer
  # viewer = TimeViewer(brain)
  
new_contents: """
  Plot MEG inverse solution
  =========================
  
  Data were computed using mne-python (http://martinos.org/mne)
  
  """
  import os
  import numpy as np
  
  from surfer import Brain
  from surfer.io import read_stc
  
  print(__doc__)
  
  """
  define subject, surface and hemisphere(s) to plot
  """
  subject_id, surface = 'fsaverage', 'inflated'
  hemi = 'split'
  
  """
  create Brain object for visualization
  """
  brain = Brain(subject_id, hemi, surface, size=(800, 400))
  
  """
  label for time annotation in milliseconds
  """
  
  
  def time_label(t):
      return 'time=%0.2f ms' % (t * 1e3)
  
  """
  read MNE dSPM inverse solution
  """
  for hemi in ['lh', 'rh']:
      stc_fname = os.path.join('example_data/meg_source_estimate-' +
                               hemi + '.stc')
      stc = read_stc(stc_fname)
  
      """
      data and vertices for which the data is defined
      """
      data = stc['data']
      vertices = stc['vertices']
  
      """
      time points (in seconds)
      """
      time = np.linspace(stc['tmin'], stc['tmin'] + data.shape[1] * stc['tstep'],
                         data.shape[1])
  
      """
      colormap to use
      """
      colormap = 'hot'
  
      brain.add_data(data, colormap=colormap, vertices=vertices,
                     smoothing_steps=10, time=time, time_label=time_label,
                     hemi=hemi, initial_time=0.1)
  
  """
  scale colormap
  """
  brain.scale_data_colormap(fmin=13, fmid=18, fmax=22, transparent=True)
  
  """
  uncomment these lines to use the interactive TimeViewer GUI
  """
  # from surfer import TimeViewer
  # viewer = TimeViewer(brain)
  
subject: update example