I am not completely sure, but I believe you should use the -blip_forward_dset and -blip_reverse_dset options using your fmriAP and fmriPA datasets in addition to those options for spatial warping to the macaque template described in the MACAQUE_DEMO scripts. For the fieldmap distortion correction instead, you could use epi_b0_correct.py as a separate procedure.