My first thought is that you should run topup on the magnitude data and then use applytopup on the phase images. Does that make sense? (Also pinging @mattcieslak, who might have ideas)
I’m also wondering- what will you use the phase data for? I am familiar with using the phase data with dwidenoise, but I didn’t realize it was useful for DWI data outside of that.