1) Input multi-coil complex images = Cx# (where #=0,1,2,...)
2) form absolute sum volume = AbsSum
3) register AbsSum in time = AbSumReg
4) apply registration parameters to each Cx# volume = CxReg#
Before the above, maybe also
A) Input B0 map
B) Smooth, convert to EPI displacement field
C) Unwarp Cx# images in y-direction before further processing
Main drawbacks -- this is really boring to do
