Epochs.plot_psd_topomap() doesn't apply unit scaling

Epochs.plot_psd does scale the measurement units while Epochs.plot_psd_topomap doesn’t, causing confusion

import os.path as op
import mne
from mne.datasets import somato

data_path = somato.data_path()
subject = '01'
task = 'somato'
raw_fname = op.join(data_path, 'sub-{}'.format(subject), 'meg',
                    'sub-{}_task-{}_meg.fif'.format(subject, task))

raw = mne.io.read_raw_fif(raw_fname)
events = mne.find_events(raw, stim_channel='STI 014')

picks = mne.pick_types(raw.info, meg='grad', eeg=False, eog=True, stim=False)
event_id, tmin, tmax = 1, -1., 3.
baseline = (None, 0)
epochs = mne.Epochs(raw, events, event_id, tmin, tmax, picks=picks,
                    baseline=baseline, reject=dict(grad=4000e-13, eog=350e-6),
                    preload=True)

epochs.resample(200., npad='auto')  # resample to reduce computation time

fmin = 0.1
fmax = 40

epochs.plot_psd(fmin=fmin, fmax=fmax, average=True, dB=True, area_mode=None)
epochs.plot_psd_topomap(dB=True)

The displayed data have a scale difference of precisely 13 orders of magnitude.

cc @cbrnr @guiomar

[guiomar]

Thnaks a lot @hoechenberger for the super quick fix!! And @cbrnr for the discussions!!

[cbrnr]

Let’s first discuss if people would like to see this changed. I have several things that I need to finish in the next days, so currently I do not want to take on new tasks right now (but maybe next week or the week after).