Hi Daniel, thank you very much for the useful suggestions. I came up with averaging after it gave me some improvements when I was still using fsl, but now I will pay more attention to subject-specific movement, thanks.
Also, I didn't know about epi_b0_correct.py, I'm glad it has been finally implemented
"I'm surprised you would get that much distortion between EPI and anat though; I would expect only relatively small amounts left, especially after the B0 correction"
Yes, I understand. Since when I tried -maxlev 0 I had very good results for the cortex and poor results only for medulla, I guess that's because we have a problem of fieldmaps quality in the most inferior regions. That's why distortion correction can't do the whole work in there.
I will give a look at the other options you mentioned, anyway, let me say that you're all making a really great job!