Using GR scans with opposite PE to create fieldmaps

No, these files (phasediff and magnitude) come from the measurent made by a different sequence, a gradient-echo sequence, with two echo times, not a gradient-echo EPI image as it is the case for bold.

In the case where you have only GRE-EPI images with reversed PE, you can still try the PEPOLAR method , the correction should work.
Otherwise you can still do SDC with Syn-BOLD DisCo, for wich you don’t need any fieldmap: