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. 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