Dear Sir/Madam,

I am getting the following netCDF4 error when I try to run Gusto:

----------------------

UFL:WARNING Discontinuous Lagrange element requested on interval * interval, creating DQ element.

Traceback (most recent call last):

File "examples/sk_nonlinear.py", line 112, in <module>

stepper.run(t=0, tmax=tmax)

File "/Users/pb412/firedrake/src/gusto/gusto/timeloop.py", line 169, in run

state.dump(t, pickup=False)

File "/Users/pb412/firedrake/src/gusto/gusto/state.py", line 345, in dump

self.pointdata_output.dump(self.fields, t)

File "/Users/pb412/firedrake/src/gusto/gusto/state.py", line 112, in dump

var[idx, :] = vals

File "netCDF4/_netCDF4.pyx", line 1908, in netCDF4._netCDF4.Dataset.__exit__ (netCDF4/_netCDF4.c:14289)

File "netCDF4/_netCDF4.pyx", line 2012, in netCDF4._netCDF4.Dataset.close (netCDF4/_netCDF4.c:16381)

File "netCDF4/_netCDF4.pyx", line 1996, in netCDF4._netCDF4.Dataset._close (netCDF4/_netCDF4.c:16289)

File "netCDF4/_netCDF4.pyx", line 1581, in netCDF4._netCDF4._ensure_nc_success (netCDF4/_netCDF4.c:12601)

RuntimeError: NetCDF: HDF error

--------------------------------------------------------------------------

MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD

with errorcode 1.

----------------------

There error occurs independent of the model I run (e.g. Rayleigh_Taylor, sk_nonlinear, etc) or the branch I am using (e.g. master, boussinesq_2d_lab, etc). So it seems there is something wrong with the way I have setup the netCDF4 software on my machine...?

Please advise what the best way forward is.

Regards,

Paul