[Errno -101] NetCDF: HDF error: on output netcdf file

[dbonaldo]

Hi,
I am running a Python script which should run a series of OceanDrift simulations and, after each of them, access the netcdf output file and edit some attributes. The problem is that when my script tries to actually edit the file (with the command nc.Dataset(outputfile, 'r+') ), I get the following error message: "[Errno -101] NetCDF: HDF error:". It seems as though the oceandrift instance keeps the file open even if the simulation is finished. Does anybody have any suggestion to solve this problem?
Any help is most appreciated, thank you very much!
Below is the piece of code that gives rise to my problem

` o = OceanDrift(loglevel=20) #
run_reader=reader_ROMS_native.Reader(runs[run][0]+str(year)+'*_zip.nc')

        o.add_reader(run_reader)
        o.set_config("general:use_auto_landmask",False)
        o.set_config('general:coastline_action', 'previous')
        # Seed elements at defined positions, depth and time
        o.seed_elements(lon=lonseedarr2, lat=latseedarr2, radius=0, z=zseedarr2, time=run_reader.start_time+datetime.timedelta(days=starttimedelta))
        o.run(time_step=3600,duration=datetime.timedelta(days=durdays),outfile=cartella_output+'/'+outfilename)
        
        if runs[run][1]=='365_day': 
            ds=nc.Dataset(cartella_output+'/'+outfilename,'r+')# <- here is where I get stuck
            ds.variables['time'].setncattr('calendar','365_day')
            ds.close()`

[gauteh]

You can use ds = o.run(), or ds = o.result. you can also close it manually using: o.result.close()

[dbonaldo]

Brilliant, it works! Thank you so much!