On May 20, 2011, at 6:23 AM, Jeff Squyres wrote: > Shouldn't ijlena and ijdisp be 1D arrays, not 2D arrays?
Ok, if I convert ijlena and ijdisp to 1D arrays, I don't get the compile error (even though they're allocatable -- so allocate was a red herring, sorry). That's all that "use mpi" is complaining about -- that the function signatures didn't match. use mpi is your friend -- even if you don't use F90 constructs much. Compile-time checking is Very Good Thing (you were effectively "getting lucky" by passing in the 2D arrays, I think). Attached is my final version. And with this version, I see the hang when running it with the "T" parameter. That being said, I'm not an expert on the MPI IO stuff -- your code *looks* right to me, but I could be missing something subtle in the interpretation of MPI_FILE_SET_VIEW. I tried running your code with MPICH 1.3.2p1 and it also hung. Rob (ROMIO guy) -- can you comment this code? Is it correct? -- Jeff Squyres jsquy...@cisco.com For corporate legal information go to: http://www.cisco.com/web/about/doing_business/legal/cri/
x.f90
Description: Binary data