[netcdfgroup] netcdf4 file creation in fortran90 with mpi

• To: netcdfgroup@xxxxxxxxxxxxxxxx
• Subject: [netcdfgroup] netcdf4 file creation in fortran90 with mpi
• From: Momme Butenschön <mommebu@xxxxxxxx>
• Date: Tue, 13 Jul 2010 13:43:41 +0000 (GMT)

I'm having problems in the creation of files for parallel IO:
using the standard file creation: nf90_create(FILE_NAME, NF90_NETCDF4, ncid, 
comm = MPI_COMM_WORLD,info = MPI_INFO_NULL)
I receive the following error:
"Can't add HDF5 file metadata"
However, an empty file gets created and when I run the program (the parallel 
test from the netcdf tutorials) a second time, everything works out.

I've found the following information:
http://www.unidata.ucar.edu/support/help/MailArchives/netcdf/msg08863.html

where it is suggested, that instead of using NF90_NETCDF4 one should use 
NC_MPIIO or NC_MPIPOSIX, however, those are part of the c-interface and not of 
the nf90 interface and corresponding constants for the f90 interface don't 
exist?!

Now, if I declare them myself locally as in: 
   integer, parameter :: NF90_MPIIO=8192,NF90_MPIPOSIX=16384 
correspnding to netcdf.h:
   #define NC_MPIIO       0x2000 
   #define NC_MPIPOSIX    0x4000
and use them instead of NF90_NETCDF4 the small test from the tutorial works 
neatly.
So now I'm wondering:
-why don't the constants corresponding to NC_MPIIO and NC_MPIPOSIX exist in the 
F90 interface?
-why do the parallel examples in tutoial and distribution examples as of 
netcdf-4.1.1 use the NF90_NETCDF4 mode, if it doesn't work?

I'm sure I'm missing something here or getting something wrong, so anybody that 
could shed some light on this, please answer...

Thanks,
Momme