READ me is a work in progress nii_preprocess now has an optional final argument of phase encoding direction, given as a number 1 to 3