<div dir="ltr">Hello,<br><br>I have a 4D variable inside a NetCDF file that I wish to distribute over a number of MPI tasks.  The variable will be decomposed over the first 3 dimensions but not the fouth (i.e. the fourth dimension is kept global for all MPI tasks). In other words:<br>
<br>              GLOBAL_FIELD[nx,ny,nz,nv]  ==&gt; LOCAL_FIELD[nx_local,ny_local,nz_local,nv]<br><br>I am trying to achieve via a nfmpi_get_vara_double_all call but the data keeps getting corrupted.  I am sure that my offsets and local domain sizes are correct.  If I modify my code to read only a single 3D slice (i.e. along 1 point in the fourth dimension), the code and input data are correct.<br>
<br>Can parallel-netcdf handle a local dimension being equal to a global dimension?  Or should I be using another call?<br><br>Thanks,<br>Mark<br clear="all"><br>-- <br>Mark Patrick Cheeseman<br><br>Computational Scientist<br>
KSL (KAUST Supercomputing Laboratory)<br>Building 1, Office #126<br>King Abdullah University of Science &amp; Technology<br>Thuwal 23955-6900<br>Kingdom of Saudi Arabia<br>         <br>EMAIL   : <a href="mailto:mark.cheeseman@kaust.edu.sa">mark.cheeseman@kaust.edu.sa</a><br>
PHONE : +966   (2) 808 0221 (office)<br>               +966 (54) 470 1082 (mobile)<br>SKYPE : mark.patrick.cheeseman<br>
</div>