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