diff --git a/gather_4.f90 b/gather_4.f90
index b53fecd..d6165bc 100644
--- a/gather_4.f90
+++ b/gather_4.f90
@@ -45,18 +45,15 @@ subroutine gather_4
   ! velocities
   call xFFT3d(-1,1)
   fname = 'u.'//file_ext
-  tmp4(1:nx,1:ny,1:nz) = wrk(1:nx,1:ny,1:nz,1)
-  call gather_tmp4(u_)
+  u_(1:nx,1:ny,1:nz) = wrk(1:nx,1:ny,1:nz,1)
 
   call xFFT3d(-1,2)
   fname = 'v.'//file_ext
-  tmp4(1:nx,1:ny,1:nz) = wrk(1:nx,1:ny,1:nz,2)
-  call gather_tmp4(v_)
+  v_(1:nx,1:ny,1:nz) = wrk(1:nx,1:ny,1:nz,2)
 
   call xFFT3d(-1,3)
   fname = 'w.'//file_ext
-  tmp4(1:nx,1:ny,1:nz) = wrk(1:nx,1:ny,1:nz,3)
-  call gather_tmp4(w_)
+  w_(1:nx,1:ny,1:nz) = wrk(1:nx,1:ny,1:nz,3)
 
   ! scalars
   if (int_scalars) then
diff --git a/main.f90 b/main.f90
index 3db4169..a6abb28 100644
--- a/main.f90
+++ b/main.f90
@@ -130,7 +130,7 @@ program x_code
         end if
 
         if (fdm_sw.ne.0) then  ! Reacting
-            !call gather_4
+            call gather_4
                call fdm_exe
         endif