我有一组 3D 数据v(i,j,k)，我想计算平均 2D 切片，该切片沿其第三维vmean(i,j)求和nz该组的切片。v(i,j,k)我写下了这段 FORTRAN90 代码，但它没有产生正确的结果。任何人都可以建议我的解决方案？提前致谢

gbl_sum = 0.0
do j    = mystarty,myendy
   do i = mystartx,myendy

      lcl_sum = sum(v(i,j,mystartz:myendz))

      call MPI_REDUCE(lcl_sum, gbl_sum(i,j), &
      1, real, mpi_sum, root,mpi_comm_world,err)

      gbl_sum(i,j) = lcl_sum

   enddo
enddo

call mpi_BCAST(gbl_sum,gbl_sumSize,real,root,&
mpi_comm_world,err)

do j    = mystarty,myendy
   do i = mystartx,myendy
      vmean(i,j) = gbl_sum(i,j)/real(nz)
   enddo
enddo