如何在Fortran中读取使用MPI-IO编写的二进制数据？

Question

我想从Fortran代码中读出一些输出数据到数据文件以进行后处理。给我的代码使用MPI写入这些数据文件。对于最基本的情况，当处理器的数量应该是1而且Nx, Ny, Nz are all set to 1 and l is 3时，写入数据文件的输出应该都是1.但是，写入数据文件的输出是非某些非 - 可读格式，如下所示：

^@^@^@^@^@^@^@^@^@^@^@^@

下面附有写入数据文件的代码的相关部分。

 do rst=1,1
   fname='phi'
   fname = trim('ufs')//":" // trim(fname)
   write(buffer,"(i3.3)") 2*rst-1
   fname = trim(adjustl(fname))//'.'//trim(adjustl(buffer))
   print *, 'The output file to which all the data is written is ', fname
   call MPI_FILE_OPEN(MPI_COMM_WORLD,fname,MPI_MODE_WRONLY+MPI_MODE_CREATE,MPI_INFO_NULL,ifile,ierr)
   do l=1,numvar
      disp = Nx_MOK*Ny_MOK*Nz_MOK*WP_MOK*(l-1)
      call MPI_FILE_SET_VIEW(ifile,disp,MPI_REAL_SP,view,"native",MPI_INFO_NULL,ierr)
      call MPI_FILE_WRITE_ALL(ifile,phi_xyz(1:Nxp,1:Nyp,1:Nzp,l,2*rst 1),Nxp*Nzp*Nyp, MPI_REAL_SP,status,ierr)
   end do
   call MPI_FILE_CLOSE(ifile,ierr)
   call MPI_BARRIER(MPI_COMM_WORLD, ierr)
end do