I finally got it to work (with yours and a colleague’s help). Given that the sounddevice library you recommended was able to detect instances of the device with 4-channel outputs, I used that same library to specify and play the 4-channel device as follows:
import sounddevice as sd
import soundfile as sf
data, fs = sf.read(filename, dtype='float32')
sd.play(data, fs, device=10)
status = sd.wait()
The PsychoPy sound library also calls the same sounddevice library so I assume this can also be done using the former?
Thanks again for your help!
Mick