ENH: Construct resampled data array in F order

Ensures that each volume is contiguous in memory, which should keep it
more efficient for both writing resampled volumes into the array and
writing the array to disk.

fmriprep/interfaces/resampling.py

@@ -372,7 +372,7 @@ async def resample_series_async(
 
     semaphore = asyncio.Semaphore(max_concurrent)
 
-    out_array = np.zeros(coordinates.shape[1:] + data.shape[-1:], dtype=output_dtype)
+    out_array = np.zeros(coordinates.shape[1:] + data.shape[-1:], dtype=output_dtype, order='F')
 
     tasks = [
         asyncio.create_task(