Deadlock when output fields with mpi

[rerdo]

Dear all,

I am unable to output the fields when using Meep with MPI. The relevant section of the python code to output is:

sim.use_output_directory("fields")
sim.run(mp.at_beginning(mp.output_epsilon), 
        mp.at_every(1/fcen/20, mp.output_hfield_z), 
        until=1/fcen)

Here is part of the output when stops:

run 0 finished at t = 400.0 (8000 timesteps)
Meep: using output directory 'fields'
creating output file "fields/test2-eps-000400.00.h5"...

I built HDF5 with MPI support, and here are configuration details:

Features:
---------
                     Parallel HDF5: yes
  Parallel Filtered Dataset Writes: yes
                Large Parallel I/O: yes
                High-level library: yes
Dimension scales w/ new references: no
                  Build HDF5 Tests: yes
                  Build HDF5 Tools: yes
                   Build GIF Tools: no
                      Threadsafety: no
               Default API mapping: v114
    With deprecated public symbols: yes
            I/O filters (external): deflate(zlib)
                  _Float16 support: yes
                     Map (H5M) API: no
                        Direct VFD: no
                        Mirror VFD: no
                     Subfiling VFD: no
                (Read-Only) S3 VFD: no
              (Read-Only) HDFS VFD: no
    Packages w/ extra debug output: none
                       API tracing: no
              Using memory checker: no
                  Use file locking: best-effort
         Strict file format checks: no
      Optimization instrumentation: no

The same script works fine if I don't use MPI.

Any insights or suggestions on how to resolve this issue would be greatly appreciated!

Best regards,
Renato