diff --git a/m_fdm_calc.f90 b/m_fdm_calc.f90
index 5909692..b3ee762 100644
--- a/m_fdm_calc.f90
+++ b/m_fdm_calc.f90
@@ -119,9 +119,11 @@ module m_fdm_calc
     offset = 2*4+6*8 + 2*4+3*8 + 3*(2*4+2*8) + 4 + 2*nx*ny*nz_all*8 + myid*nx*ny*nz*8
     call MPI_FILE_READ_AT_ALL(fh, offset, w_, count, MPI_REAL8, mpi_status, mpi_err)
 
-    count = 2 * nx * ny * nz
-    offset = 2*4+6*8 + 2*4+3*8 + 3*(2*4+2*8) + 4 + 3*nx*ny*nz_all*8 + myid*2*nx*ny*nz*8
-    call MPI_FILE_READ_AT_ALL(fh, offset, y1, count, MPI_REAL8, mpi_status, mpi_err)
+    offset = 2*4+6*8 + 2*4+3*8 + 3*(2*4+2*8) + 4 + 3*nx*ny*nz_all*8 + myid*nx*ny*nz*8
+    call MPI_FILE_READ_AT_ALL(fh, offset, y1(1,1,1,1), count, MPI_REAL8, mpi_status, mpi_err)
+
+    offset = 2*4+6*8 + 2*4+3*8 + 3*(2*4+2*8) + 4 + 4*nx*ny*nz_all*8 + myid*nx*ny*nz*8
+    call MPI_FILE_READ_AT_ALL(fh, offset, y1(1,1,1,2), count, MPI_REAL8, mpi_status, mpi_err)
 
     call MPI_FILE_CLOSE(fh, mpi_err)
     call MPI_INFO_FREE(mpi_info, mpi_err)
@@ -194,9 +196,11 @@ module m_fdm_calc
       offset = 2*4+6*8 + 2*4+3*8 + 3*(2*4+2*8) + 4 + 2*nx*ny*nz_all*8 + myid*nx*ny*nz*8
       call MPI_FILE_WRITE_AT(fh, offset, w_, count, MPI_REAL8, mpi_status, mpi_err)
 
-      count = 2 * nx * ny * nz
-      offset = 2*4+6*8 + 2*4+3*8 + 3*(2*4+2*8) + 4 + 3*nx*ny*nz_all*8 + myid*2*nx*ny*nz*8
-      call MPI_FILE_WRITE_AT(fh, offset, y1, count, MPI_REAL8, mpi_status, mpi_err)
+      offset = 2*4+6*8 + 2*4+3*8 + 3*(2*4+2*8) + 4 + 3*nx*ny*nz_all*8 + myid*nx*ny*nz*8
+      call MPI_FILE_WRITE_AT(fh, offset, y1(1,1,1,1), count, MPI_REAL8, mpi_status, mpi_err)
+
+      offset = 2*4+6*8 + 2*4+3*8 + 3*(2*4+2*8) + 4 + 4*nx*ny*nz_all*8 + myid*nx*ny*nz*8
+      call MPI_FILE_WRITE_AT(fh, offset, y1(1,1,1,2), count, MPI_REAL8, mpi_status, mpi_err)
 
 
       call MPI_FILE_CLOSE(fh, mpi_err)