@agramfort I'm wondering whether this should happen inside the decimate_surface function. Would that be consistent with the write surface API?