Oil in ice

from datetime import datetime, timedelta
from opendrift.readers import reader_netCDF_CF_generic
from opendrift.models.openoil import OpenOil
import cmocean

o = OpenOil(loglevel=20)

15:36:27 INFO    opendrift:569: OpenDriftSimulation initialised (version 1.14.11 / v1.14.12-32-g22db1c0)

Using live data from Barents 2.5 km ocean model

o.add_readers_from_list(['https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be'])

Adjusting some configuration

o.set_config('processes:dispersion',  False)
o.set_config('processes:evaporation',  False)
o.set_config('processes:emulsification',  False)
o.set_config('environment:constant:horizontal_diffusivity', 10)
o.set_config('drift:truncate_ocean_model_below_m', 3)

Imaginary oil spill in Hinlopen strait

o.seed_elements(lon=19.1909, lat=79.5986, radius=50,
                number=3000, time=datetime.now() - timedelta(days=7))

15:36:27 INFO    opendrift.models.openoil.openoil:1720: Oil type not specified, using default: GENERIC BUNKER C
15:36:27 INFO    opendrift.models.openoil.adios.dirjs:86: Querying ADIOS database for oil: GENERIC BUNKER C
15:36:27 INFO    opendrift.models.openoil.openoil:1731: Using density 988.1 and viscosity 0.021692333877975332 of oiltype GENERIC BUNKER C
15:36:27 INFO    opendrift.models.basemodel.environment:203: Adding a global landmask from GSHHG
15:36:33 INFO    opendrift.models.basemodel.environment:227: Fallback values will be used for the following variables which have no readers:
15:36:33 INFO    opendrift.models.basemodel.environment:230:    sea_surface_height: 0.000000
15:36:33 INFO    opendrift.models.basemodel.environment:230:    upward_sea_water_velocity: 0.000000
15:36:33 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_significant_height: 0.000000
15:36:33 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
15:36:33 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
15:36:33 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
15:36:33 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
15:36:33 INFO    opendrift.models.basemodel.environment:230:    sea_ice_area_fraction: 0.000000
15:36:33 INFO    opendrift.models.basemodel.environment:230:    sea_ice_x_velocity: 0.000000
15:36:33 INFO    opendrift.models.basemodel.environment:230:    sea_ice_y_velocity: 0.000000
15:36:33 INFO    opendrift.models.basemodel.environment:230:    sea_water_temperature: 10.000000
15:36:33 INFO    opendrift.models.basemodel.environment:230:    sea_water_salinity: 34.000000
15:36:33 INFO    opendrift.models.basemodel.environment:230:    sea_floor_depth_below_sea_level: 10000.000000
15:36:33 INFO    opendrift.models.basemodel.environment:230:    ocean_vertical_diffusivity: 0.020000
15:36:33 INFO    opendrift.models.basemodel.environment:230:    ocean_mixed_layer_thickness: 50.000000

Running model

o.run(duration=timedelta(hours=48), time_step=1800, time_step_output=3600)

15:36:33 INFO    opendrift:1910: Skipping environment variable upward_sea_water_velocity because of condition ['drift:vertical_advection', 'is', False]
15:36:33 INFO    opendrift:1921: Storing previous values of element property lon because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
15:36:33 INFO    opendrift:1921: Storing previous values of element property lat because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
15:36:33 INFO    opendrift:958: Using existing reader for land_binary_mask to move elements to ocean
15:36:33 INFO    opendrift:989: All points are in ocean
15:36:33 INFO    opendrift.models.openoil.openoil:697: Oil-water surface tension is 0.035935 Nm
15:36:33 INFO    opendrift.models.openoil.openoil:710: Max water fraction not available for GENERIC BUNKER C, using default
15:36:33 INFO    opendrift:2237: 2026-09-28 15:36:27.064652 - step 1 of 96 - 3000 active elements (0 deactivated)
15:36:33 INFO    opendrift.readers:67: Opening file with xr.open_dataset
15:36:35 INFO    opendrift.readers.reader_netCDF_CF_generic:340: Detected dimensions: {'x': 'X', 'y': 'Y', 'z': 'depth', 'time': 'time'}
15:36:35 INFO    opendrift.readers.basereader:175: Variable x_sea_water_velocity will be rotated from eastward_sea_water_velocity
15:36:35 INFO    opendrift.readers.basereader:175: Variable y_sea_water_velocity will be rotated from northward_sea_water_velocity
15:36:36 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.004626641 to -0.0008346351) found for ocean_vertical_diffusivity, replacing with NaN
15:36:36 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:36:38 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0013212361 to -2.9143215e-05) found for ocean_vertical_diffusivity, replacing with NaN
15:36:38 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:36:38 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:38 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:38 INFO    opendrift:2237: 2026-09-28 16:06:27.064652 - step 2 of 96 - 3000 active elements (0 deactivated)
15:36:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:40 INFO    opendrift:2237: 2026-09-28 16:36:27.064652 - step 3 of 96 - 3000 active elements (0 deactivated)
15:36:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:40 INFO    opendrift:2237: 2026-09-28 17:06:27.064652 - step 4 of 96 - 3000 active elements (0 deactivated)
15:36:42 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0013024728 to -1.5918426e-05) found for ocean_vertical_diffusivity, replacing with NaN
15:36:42 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:36:42 INFO    opendrift:2237: 2026-09-28 17:36:27.064652 - step 5 of 96 - 3000 active elements (0 deactivated)
15:36:42 INFO    opendrift:2237: 2026-09-28 18:06:27.064652 - step 6 of 96 - 2998 active elements (2 deactivated)
15:36:44 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0041281697 to -0.0005663188) found for ocean_vertical_diffusivity, replacing with NaN
15:36:44 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:36:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:44 INFO    opendrift:2237: 2026-09-28 18:36:27.064652 - step 7 of 96 - 2992 active elements (8 deactivated)
15:36:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:44 INFO    opendrift:2237: 2026-09-28 19:06:27.064652 - step 8 of 96 - 2983 active elements (17 deactivated)
15:36:46 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0040557776 to -0.0013311687) found for ocean_vertical_diffusivity, replacing with NaN
15:36:46 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:36:46 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:46 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:46 INFO    opendrift:2237: 2026-09-28 19:36:27.064652 - step 9 of 96 - 2983 active elements (17 deactivated)
15:36:46 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:46 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:46 INFO    opendrift:2237: 2026-09-28 20:06:27.064652 - step 10 of 96 - 2983 active elements (17 deactivated)
15:36:49 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0062416513 to -0.0015993095) found for ocean_vertical_diffusivity, replacing with NaN
15:36:49 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:36:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:49 INFO    opendrift:2237: 2026-09-28 20:36:27.064652 - step 11 of 96 - 2983 active elements (17 deactivated)
15:36:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:49 INFO    opendrift:2237: 2026-09-28 21:06:27.064652 - step 12 of 96 - 2983 active elements (17 deactivated)
15:36:51 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0062143244 to -0.0011027599) found for ocean_vertical_diffusivity, replacing with NaN
15:36:51 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:36:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:51 INFO    opendrift:2237: 2026-09-28 21:36:27.064652 - step 13 of 96 - 2983 active elements (17 deactivated)
15:36:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:51 INFO    opendrift:2237: 2026-09-28 22:06:27.064652 - step 14 of 96 - 2983 active elements (17 deactivated)
15:36:53 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0065603764 to -0.0010057949) found for ocean_vertical_diffusivity, replacing with NaN
15:36:53 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:36:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:53 INFO    opendrift:2237: 2026-09-28 22:36:27.064652 - step 15 of 96 - 2983 active elements (17 deactivated)
15:36:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:53 INFO    opendrift:2237: 2026-09-28 23:06:27.064652 - step 16 of 96 - 2983 active elements (17 deactivated)
15:36:55 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0032145928 to -0.00097629847) found for ocean_vertical_diffusivity, replacing with NaN
15:36:55 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:36:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:55 INFO    opendrift:2237: 2026-09-28 23:36:27.064652 - step 17 of 96 - 2983 active elements (17 deactivated)
15:36:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:55 INFO    opendrift:2237: 2026-09-29 00:06:27.064652 - step 18 of 96 - 2983 active elements (17 deactivated)
15:36:56 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0040910146 to -0.00046879935) found for ocean_vertical_diffusivity, replacing with NaN
15:36:56 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:36:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:57 INFO    opendrift:2237: 2026-09-29 00:36:27.064652 - step 19 of 96 - 2983 active elements (17 deactivated)
15:36:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:57 INFO    opendrift:2237: 2026-09-29 01:06:27.064652 - step 20 of 96 - 2983 active elements (17 deactivated)
15:36:58 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0036701001 to -0.00041647046) found for ocean_vertical_diffusivity, replacing with NaN
15:36:58 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:36:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:59 INFO    opendrift:2237: 2026-09-29 01:36:27.064652 - step 21 of 96 - 2983 active elements (17 deactivated)
15:36:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:36:59 INFO    opendrift:2237: 2026-09-29 02:06:27.064652 - step 22 of 96 - 2983 active elements (17 deactivated)
15:37:00 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.002955305 to -0.00027978115) found for ocean_vertical_diffusivity, replacing with NaN
15:37:00 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:00 INFO    opendrift:2237: 2026-09-29 02:36:27.064652 - step 23 of 96 - 2983 active elements (17 deactivated)
15:37:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:01 INFO    opendrift:2237: 2026-09-29 03:06:27.064652 - step 24 of 96 - 2983 active elements (17 deactivated)
15:37:02 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0012282251 to -4.4843955e-06) found for ocean_vertical_diffusivity, replacing with NaN
15:37:02 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:02 INFO    opendrift:2237: 2026-09-29 03:36:27.064652 - step 25 of 96 - 2983 active elements (17 deactivated)
15:37:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:03 INFO    opendrift:2237: 2026-09-29 04:06:27.064652 - step 26 of 96 - 2983 active elements (17 deactivated)
15:37:04 INFO    opendrift:2237: 2026-09-29 04:36:27.064652 - step 27 of 96 - 2983 active elements (17 deactivated)
15:37:04 INFO    opendrift:2237: 2026-09-29 05:06:27.064652 - step 28 of 96 - 2983 active elements (17 deactivated)
15:37:06 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.00029654286 to -7.671233e-06) found for ocean_vertical_diffusivity, replacing with NaN
15:37:06 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:06 INFO    opendrift:2237: 2026-09-29 05:36:27.064652 - step 29 of 96 - 2983 active elements (17 deactivated)
15:37:06 INFO    opendrift:2237: 2026-09-29 06:06:27.064652 - step 30 of 96 - 2983 active elements (17 deactivated)
15:37:08 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.007942783 to -0.0004245968) found for ocean_vertical_diffusivity, replacing with NaN
15:37:08 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:08 INFO    opendrift:2237: 2026-09-29 06:36:27.064652 - step 31 of 96 - 2983 active elements (17 deactivated)
15:37:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:09 INFO    opendrift:2237: 2026-09-29 07:06:27.064652 - step 32 of 96 - 2983 active elements (17 deactivated)
15:37:10 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.017713364 to -0.0015669635) found for ocean_vertical_diffusivity, replacing with NaN
15:37:10 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:10 INFO    opendrift:2237: 2026-09-29 07:36:27.064652 - step 33 of 96 - 2983 active elements (17 deactivated)
15:37:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:11 INFO    opendrift:2237: 2026-09-29 08:06:27.064652 - step 34 of 96 - 2983 active elements (17 deactivated)
15:37:12 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.01878558 to -0.0037337749) found for ocean_vertical_diffusivity, replacing with NaN
15:37:12 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:12 INFO    opendrift:2237: 2026-09-29 08:36:27.064652 - step 35 of 96 - 2983 active elements (17 deactivated)
15:37:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:13 INFO    opendrift:2237: 2026-09-29 09:06:27.064652 - step 36 of 96 - 2982 active elements (18 deactivated)
15:37:14 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.023472173 to -0.005387252) found for ocean_vertical_diffusivity, replacing with NaN
15:37:14 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:14 INFO    opendrift:2237: 2026-09-29 09:36:27.064652 - step 37 of 96 - 2980 active elements (20 deactivated)
15:37:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:15 INFO    opendrift:2237: 2026-09-29 10:06:27.064652 - step 38 of 96 - 2974 active elements (26 deactivated)
15:37:16 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.023231104 to -0.008322898) found for ocean_vertical_diffusivity, replacing with NaN
15:37:16 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:16 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:16 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:16 INFO    opendrift:2237: 2026-09-29 10:36:27.064652 - step 39 of 96 - 2969 active elements (31 deactivated)
15:37:16 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:16 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:17 INFO    opendrift:2237: 2026-09-29 11:06:27.064652 - step 40 of 96 - 2964 active elements (36 deactivated)
15:37:18 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.016054468 to -0.005876638) found for ocean_vertical_diffusivity, replacing with NaN
15:37:18 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:18 INFO    opendrift:2237: 2026-09-29 11:36:27.064652 - step 41 of 96 - 2953 active elements (47 deactivated)
15:37:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:19 INFO    opendrift:2237: 2026-09-29 12:06:27.064652 - step 42 of 96 - 2949 active elements (51 deactivated)
15:37:20 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.012469868 to -0.0045253797) found for ocean_vertical_diffusivity, replacing with NaN
15:37:20 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:20 INFO    opendrift:2237: 2026-09-29 12:36:27.064652 - step 43 of 96 - 2941 active elements (59 deactivated)
15:37:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:21 INFO    opendrift:2237: 2026-09-29 13:06:27.064652 - step 44 of 96 - 2927 active elements (73 deactivated)
15:37:22 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.011941036 to -0.005491759) found for ocean_vertical_diffusivity, replacing with NaN
15:37:22 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:22 INFO    opendrift:2237: 2026-09-29 13:36:27.064652 - step 45 of 96 - 2907 active elements (93 deactivated)
15:37:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:22 INFO    opendrift:2237: 2026-09-29 14:06:27.064652 - step 46 of 96 - 2865 active elements (135 deactivated)
15:37:25 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.010633689 to -0.0037270577) found for ocean_vertical_diffusivity, replacing with NaN
15:37:25 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:25 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:25 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:25 INFO    opendrift:2237: 2026-09-29 14:36:27.064652 - step 47 of 96 - 2806 active elements (194 deactivated)
15:37:25 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:25 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:25 INFO    opendrift:2237: 2026-09-29 15:06:27.064652 - step 48 of 96 - 2704 active elements (296 deactivated)
15:37:26 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0066547124 to -0.00081041455) found for ocean_vertical_diffusivity, replacing with NaN
15:37:26 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:27 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:27 INFO    opendrift:2237: 2026-09-29 15:36:27.064652 - step 49 of 96 - 2561 active elements (439 deactivated)
15:37:27 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:27 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:27 INFO    opendrift:2237: 2026-09-29 16:06:27.064652 - step 50 of 96 - 2396 active elements (604 deactivated)
15:37:29 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0010889776 to -8.9856167e-07) found for ocean_vertical_diffusivity, replacing with NaN
15:37:29 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:29 INFO    opendrift:2237: 2026-09-29 16:36:27.064652 - step 51 of 96 - 2227 active elements (773 deactivated)
15:37:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:29 INFO    opendrift:2237: 2026-09-29 17:06:27.064652 - step 52 of 96 - 2099 active elements (901 deactivated)
15:37:31 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.00026440746 to -9.479193e-05) found for ocean_vertical_diffusivity, replacing with NaN
15:37:31 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:31 INFO    opendrift:2237: 2026-09-29 17:36:27.064652 - step 53 of 96 - 2037 active elements (963 deactivated)
15:37:31 INFO    opendrift:2237: 2026-09-29 18:06:27.064652 - step 54 of 96 - 1982 active elements (1018 deactivated)
15:37:33 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.005536999 to -4.0023224e-05) found for ocean_vertical_diffusivity, replacing with NaN
15:37:33 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:33 INFO    opendrift:2237: 2026-09-29 18:36:27.064652 - step 55 of 96 - 1931 active elements (1069 deactivated)
15:37:33 INFO    opendrift:2237: 2026-09-29 19:06:27.064652 - step 56 of 96 - 1870 active elements (1130 deactivated)
15:37:35 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.014833087 to -0.004484337) found for ocean_vertical_diffusivity, replacing with NaN
15:37:35 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:35 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:35 INFO    opendrift:2237: 2026-09-29 19:36:27.064652 - step 57 of 96 - 1846 active elements (1154 deactivated)
15:37:35 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:35 INFO    opendrift:2237: 2026-09-29 20:06:27.064652 - step 58 of 96 - 1838 active elements (1162 deactivated)
15:37:37 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.023444636 to -0.008672103) found for ocean_vertical_diffusivity, replacing with NaN
15:37:37 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:37 INFO    opendrift:2237: 2026-09-29 20:36:27.064652 - step 59 of 96 - 1827 active elements (1173 deactivated)
15:37:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:37 INFO    opendrift:2237: 2026-09-29 21:06:27.064652 - step 60 of 96 - 1819 active elements (1181 deactivated)
15:37:39 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.025571264 to -0.010380705) found for ocean_vertical_diffusivity, replacing with NaN
15:37:39 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:39 INFO    opendrift:2237: 2026-09-29 21:36:27.064652 - step 61 of 96 - 1812 active elements (1188 deactivated)
15:37:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:39 INFO    opendrift:2237: 2026-09-29 22:06:27.064652 - step 62 of 96 - 1807 active elements (1193 deactivated)
15:37:41 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.021889897 to -0.0070490963) found for ocean_vertical_diffusivity, replacing with NaN
15:37:41 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:41 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:41 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:41 INFO    opendrift:2237: 2026-09-29 22:36:27.064652 - step 63 of 96 - 1805 active elements (1195 deactivated)
15:37:41 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:41 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:41 INFO    opendrift:2237: 2026-09-29 23:06:27.064652 - step 64 of 96 - 1797 active elements (1203 deactivated)
15:37:43 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.015959783 to -0.004011377) found for ocean_vertical_diffusivity, replacing with NaN
15:37:43 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:43 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:43 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:43 INFO    opendrift:2237: 2026-09-29 23:36:27.064652 - step 65 of 96 - 1769 active elements (1231 deactivated)
15:37:43 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:43 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:43 INFO    opendrift:2237: 2026-09-30 00:06:27.064652 - step 66 of 96 - 1721 active elements (1279 deactivated)
15:37:45 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.012799888 to -0.0032469553) found for ocean_vertical_diffusivity, replacing with NaN
15:37:45 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:45 INFO    opendrift:2237: 2026-09-30 00:36:27.064652 - step 67 of 96 - 1679 active elements (1321 deactivated)
15:37:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:45 INFO    opendrift:2237: 2026-09-30 01:06:27.064652 - step 68 of 96 - 1657 active elements (1343 deactivated)
15:37:47 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.012726489 to -0.0024748237) found for ocean_vertical_diffusivity, replacing with NaN
15:37:47 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:47 INFO    opendrift:2237: 2026-09-30 01:36:27.064652 - step 69 of 96 - 1648 active elements (1352 deactivated)
15:37:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:47 INFO    opendrift:2237: 2026-09-30 02:06:27.064652 - step 70 of 96 - 1640 active elements (1360 deactivated)
15:37:49 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.011170388 to -0.0019898044) found for ocean_vertical_diffusivity, replacing with NaN
15:37:49 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:49 INFO    opendrift:2237: 2026-09-30 02:36:27.064652 - step 71 of 96 - 1637 active elements (1363 deactivated)
15:37:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:49 INFO    opendrift:2237: 2026-09-30 03:06:27.064652 - step 72 of 96 - 1635 active elements (1365 deactivated)
15:37:51 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.008742366 to -0.0008534277) found for ocean_vertical_diffusivity, replacing with NaN
15:37:51 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:51 INFO    opendrift:2237: 2026-09-30 03:36:27.064652 - step 73 of 96 - 1634 active elements (1366 deactivated)
15:37:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:51 INFO    opendrift:2237: 2026-09-30 04:06:27.064652 - step 74 of 96 - 1633 active elements (1367 deactivated)
15:37:53 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.004430075 to -3.6014095e-05) found for ocean_vertical_diffusivity, replacing with NaN
15:37:53 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:53 INFO    opendrift:2237: 2026-09-30 04:36:27.064652 - step 75 of 96 - 1632 active elements (1368 deactivated)
15:37:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:53 INFO    opendrift:2237: 2026-09-30 05:06:27.064652 - step 76 of 96 - 1632 active elements (1368 deactivated)
15:37:54 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0011316403 to -6.2335107e-07) found for ocean_vertical_diffusivity, replacing with NaN
15:37:54 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:55 INFO    opendrift:2237: 2026-09-30 05:36:27.064652 - step 77 of 96 - 1632 active elements (1368 deactivated)
15:37:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:55 INFO    opendrift:2237: 2026-09-30 06:06:27.064652 - step 78 of 96 - 1631 active elements (1369 deactivated)
15:37:56 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0053583384 to -2.5227224e-05) found for ocean_vertical_diffusivity, replacing with NaN
15:37:56 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:57 INFO    opendrift:2237: 2026-09-30 06:36:27.064652 - step 79 of 96 - 1629 active elements (1371 deactivated)
15:37:57 INFO    opendrift:2237: 2026-09-30 07:06:27.064652 - step 80 of 96 - 1621 active elements (1379 deactivated)
15:37:59 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.010724524 to -0.0007832428) found for ocean_vertical_diffusivity, replacing with NaN
15:37:59 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:37:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:59 INFO    opendrift:2237: 2026-09-30 07:36:27.064652 - step 81 of 96 - 1606 active elements (1394 deactivated)
15:37:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:37:59 INFO    opendrift:2237: 2026-09-30 08:06:27.064652 - step 82 of 96 - 1585 active elements (1415 deactivated)
15:38:00 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.017289756 to -0.0030840202) found for ocean_vertical_diffusivity, replacing with NaN
15:38:00 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:38:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:01 INFO    opendrift:2237: 2026-09-30 08:36:27.064652 - step 83 of 96 - 1560 active elements (1440 deactivated)
15:38:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:01 INFO    opendrift:2237: 2026-09-30 09:06:27.064652 - step 84 of 96 - 1545 active elements (1455 deactivated)
15:38:02 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.019735023 to -0.002635561) found for ocean_vertical_diffusivity, replacing with NaN
15:38:02 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:38:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:02 INFO    opendrift:2237: 2026-09-30 09:36:27.064652 - step 85 of 96 - 1541 active elements (1459 deactivated)
15:38:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:03 INFO    opendrift:2237: 2026-09-30 10:06:27.064652 - step 86 of 96 - 1540 active elements (1460 deactivated)
15:38:04 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.018549675 to -0.0031876839) found for ocean_vertical_diffusivity, replacing with NaN
15:38:04 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:38:04 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:04 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:04 INFO    opendrift:2237: 2026-09-30 10:36:27.064652 - step 87 of 96 - 1539 active elements (1461 deactivated)
15:38:04 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:04 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:05 INFO    opendrift:2237: 2026-09-30 11:06:27.064652 - step 88 of 96 - 1539 active elements (1461 deactivated)
15:38:06 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.0147387935 to -0.0026172535) found for ocean_vertical_diffusivity, replacing with NaN
15:38:06 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:38:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:06 INFO    opendrift:2237: 2026-09-30 11:36:27.064652 - step 89 of 96 - 1539 active elements (1461 deactivated)
15:38:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:06 INFO    opendrift:2237: 2026-09-30 12:06:27.064652 - step 90 of 96 - 1539 active elements (1461 deactivated)
15:38:08 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.011068075 to -0.002371993) found for ocean_vertical_diffusivity, replacing with NaN
15:38:08 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:38:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:08 INFO    opendrift:2237: 2026-09-30 12:36:27.064652 - step 91 of 96 - 1537 active elements (1463 deactivated)
15:38:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:08 INFO    opendrift:2237: 2026-09-30 13:06:27.064652 - step 92 of 96 - 1528 active elements (1472 deactivated)
15:38:10 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.009003784 to -0.0026503664) found for ocean_vertical_diffusivity, replacing with NaN
15:38:10 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:38:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:10 INFO    opendrift:2237: 2026-09-30 13:36:27.064652 - step 93 of 96 - 1505 active elements (1495 deactivated)
15:38:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:10 INFO    opendrift:2237: 2026-09-30 14:06:27.064652 - step 94 of 96 - 1465 active elements (1535 deactivated)
15:38:12 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.008781962 to -0.0011712417) found for ocean_vertical_diffusivity, replacing with NaN
15:38:12 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:38:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:12 INFO    opendrift:2237: 2026-09-30 14:36:27.064652 - step 95 of 96 - 1422 active elements (1578 deactivated)
15:38:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:12 INFO    opendrift:2237: 2026-09-30 15:06:27.064652 - step 96 of 96 - 1349 active elements (1651 deactivated)
15:38:14 WARNING opendrift.readers.basereader.variables:660: Invalid values (-0.007939765 to -0.0016698417) found for ocean_vertical_diffusivity, replacing with NaN
15:38:14 WARNING opendrift.readers.basereader.variables:663: (allowed range: [0, 1])
15:38:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:38:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning

<xarray.Dataset> Size: 26MB
Dimensions:                                                                              (
                                                                                          trajectory: 3000,
                                                                                          time: 49)
Coordinates:
  * trajectory                                                                           (trajectory) int64 24kB ...
  * time                                                                                 (time) datetime64[ns] 392B ...
Data variables: (12/44)
    status                                                                               (trajectory, time) float32 588kB ...
    moving                                                                               (trajectory, time) float32 588kB ...
    age_seconds                                                                          (trajectory, time) float32 588kB ...
    origin_marker                                                                        (trajectory, time) float32 588kB ...
    lon                                                                                  (trajectory, time) float32 588kB ...
    lat                                                                                  (trajectory, time) float32 588kB ...
    ...                                                                                   ...
    sea_water_salinity                                                                   (trajectory, time) float32 588kB ...
    sea_floor_depth_below_sea_level                                                      (trajectory, time) float32 588kB ...
    horizontal_diffusivity                                                               (trajectory, time) float32 588kB ...
    ocean_vertical_diffusivity                                                           (trajectory, time) float32 588kB ...
    land_binary_mask                                                                     (trajectory, time) float32 588kB ...
    ocean_mixed_layer_thickness                                                          (trajectory, time) float32 588kB ...
Attributes: (12/167)
    Conventions:                                                             ...
    standard_name_vocabulary:                                                ...
    featureType:                                                             ...
    title:                                                                   ...
    summary:                                                                 ...
    keywords:                                                                ...
    ...                                                                                                               ...
    geospatial_lon_units:                                                    ...
    geospatial_lon_resolution:                                               ...
    runtime:                                                                 ...
    geospatial_vertical_min:                                                 ...
    geospatial_vertical_max:                                                 ...
    geospatial_vertical_positive:                                            ...

xarray.Dataset

• Dimensions:
  • trajectory: 3000
  • time: 49
• Coordinates: (2)
  • trajectory
    (trajectory)
    int64
    0 1 2 3 4 ... 2996 2997 2998 2999

    cf_role :

    trajectory_id

    dtype :

    <class 'numpy.int32'>

    array([   0,    1,    2, ..., 2997, 2998, 2999], shape=(3000,))

  • time
    (time)
    datetime64[ns]
    2026-09-28T15:36:27.064652 ... 2...

    standard_name :

    time

    long_name :

    time

    dtype :

    <class 'numpy.float64'>

    array(['2026-09-28T15:36:27.064652000', '2026-09-28T16:36:27.064652000',
           '2026-09-28T17:36:27.064652000', '2026-09-28T18:36:27.064652000',
           '2026-09-28T19:36:27.064652000', '2026-09-28T20:36:27.064652000',
           '2026-09-28T21:36:27.064652000', '2026-09-28T22:36:27.064652000',
           '2026-09-28T23:36:27.064652000', '2026-09-29T00:36:27.064652000',
           '2026-09-29T01:36:27.064652000', '2026-09-29T02:36:27.064652000',
           '2026-09-29T03:36:27.064652000', '2026-09-29T04:36:27.064652000',
           '2026-09-29T05:36:27.064652000', '2026-09-29T06:36:27.064652000',
           '2026-09-29T07:36:27.064652000', '2026-09-29T08:36:27.064652000',
           '2026-09-29T09:36:27.064652000', '2026-09-29T10:36:27.064652000',
           '2026-09-29T11:36:27.064652000', '2026-09-29T12:36:27.064652000',
           '2026-09-29T13:36:27.064652000', '2026-09-29T14:36:27.064652000',
           '2026-09-29T15:36:27.064652000', '2026-09-29T16:36:27.064652000',
           '2026-09-29T17:36:27.064652000', '2026-09-29T18:36:27.064652000',
           '2026-09-29T19:36:27.064652000', '2026-09-29T20:36:27.064652000',
           '2026-09-29T21:36:27.064652000', '2026-09-29T22:36:27.064652000',
           '2026-09-29T23:36:27.064652000', '2026-09-30T00:36:27.064652000',
           '2026-09-30T01:36:27.064652000', '2026-09-30T02:36:27.064652000',
           '2026-09-30T03:36:27.064652000', '2026-09-30T04:36:27.064652000',
           '2026-09-30T05:36:27.064652000', '2026-09-30T06:36:27.064652000',
           '2026-09-30T07:36:27.064652000', '2026-09-30T08:36:27.064652000',
           '2026-09-30T09:36:27.064652000', '2026-09-30T10:36:27.064652000',
           '2026-09-30T11:36:27.064652000', '2026-09-30T12:36:27.064652000',
           '2026-09-30T13:36:27.064652000', '2026-09-30T14:36:27.064652000',
           '2026-09-30T15:36:27.064652000'], dtype='datetime64[ns]')

• Data variables: (44)
  • status
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    dtype :

    <class 'numpy.int32'>

    valid_range :

    [0 1]

    flag_values :

    [0 1]

    flag_meanings :

    active stranded

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • moving
    (trajectory, time)
    float32
    1.0 1.0 1.0 1.0 ... 1.0 1.0 1.0 1.0

    dtype :

    <class 'numpy.int32'>

    minval :

    0

    maxval :

    1

    array([[ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           ...,
           [ 1.,  1.,  1., ..., nan, nan, nan],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.]],
          shape=(3000, 49), dtype=float32)

  • age_seconds
    (trajectory, time)
    float32
    0.0 3.6e+03 ... 1.692e+05 1.728e+05

    dtype :

    <class 'numpy.float32'>

    units :

    s

    minval :

    0.0

    maxval :

    172800.0

    array([[     0.,   3600.,   7200., ..., 165600., 169200., 172800.],
           [     0.,   3600.,   7200., ..., 165600., 169200., 172800.],
           [     0.,   3600.,   7200., ..., 165600., 169200., 172800.],
           ...,
           [     0.,   3600.,   7200., ...,     nan,     nan,     nan],
           [     0.,   3600.,   7200., ..., 165600., 169200., 172800.],
           [     0.,   3600.,   7200., ..., 165600., 169200., 172800.]],
          shape=(3000, 49), dtype=float32)

  • origin_marker
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    dtype :

    <class 'numpy.int32'>

    unit :

    description :

    An integer kept constant during the simulation. Different values may be used for different seedings, to separate elements during analysis. With GUI, only a single seeding is possible.

    flag_values :

    [0]

    flag_meanings :

    Seed_0

    minval :

    0

    maxval :

    0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • lon
    (trajectory, time)
    float32
    19.2 19.21 19.21 ... 19.24 19.25

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    axis :

    X

    minval :

    18.550377

    maxval :

    19.799967

    array([[19.195274, 19.209381, 19.207743, ..., 19.272028, 19.300425,
            19.369452],
           [19.191893, 19.175005, 19.15926 , ..., 19.259792, 19.304672,
            19.338125],
           [19.193327, 19.179413, 19.175678, ..., 19.235308, 19.26989 ,
            19.325489],
           ...,
           [19.190422, 19.185062, 19.174139, ...,       nan,       nan,
                  nan],
           [19.187893, 19.167309, 19.167189, ..., 19.169973, 19.21618 ,
            19.224728],
           [19.1907  , 19.17196 , 19.185202, ..., 19.192822, 19.236609,
            19.245922]], shape=(3000, 49), dtype=float32)

  • lat
    (trajectory, time)
    float32
    79.6 79.6 79.6 ... 79.57 79.58

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    axis :

    Y

    minval :

    79.37862

    maxval :

    79.65727

    array([[79.59931 , 79.59694 , 79.59713 , ..., 79.51768 , 79.531815,
            79.53395 ],
           [79.598854, 79.59792 , 79.589165, ..., 79.55002 , 79.559555,
            79.56348 ],
           [79.59855 , 79.59979 , 79.59462 , ..., 79.54719 , 79.55644 ,
            79.557304],
           ...,
           [79.59836 , 79.59048 , 79.589035, ...,       nan,       nan,
                  nan],
           [79.59878 , 79.590355, 79.5878  , ..., 79.55707 , 79.5641  ,
            79.565544],
           [79.59873 , 79.59663 , 79.59226 , ..., 79.55869 , 79.568016,
            79.57528 ]], shape=(3000, 49), dtype=float32)

  • z
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    dtype :

    <class 'numpy.float32'>

    units :

    m

    standard_name :

    z

    long_name :

    vertical position

    axis :

    Z

    positive :

    up

    minval :

    -199.42229

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • wind_drift_factor
    (trajectory, time)
    float32
    0.03 0.03 0.03 ... 0.03 0.03 0.03

    dtype :

    <class 'numpy.float32'>

    units :

    1

    description :

    Elements at the ocean surface are moved by this fraction of the wind vector, in addition to currents and Stokes drift

    minval :

    0.03

    maxval :

    0.03

    array([[0.03, 0.03, 0.03, ..., 0.03, 0.03, 0.03],
           [0.03, 0.03, 0.03, ..., 0.03, 0.03, 0.03],
           [0.03, 0.03, 0.03, ..., 0.03, 0.03, 0.03],
           ...,
           [0.03, 0.03, 0.03, ...,  nan,  nan,  nan],
           [0.03, 0.03, 0.03, ..., 0.03, 0.03, 0.03],
           [0.03, 0.03, 0.03, ..., 0.03, 0.03, 0.03]],
          shape=(3000, 49), dtype=float32)

  • current_drift_factor
    (trajectory, time)
    float32
    1.0 1.0 1.0 1.0 ... 1.0 1.0 1.0 1.0

    dtype :

    <class 'numpy.float32'>

    units :

    1

    description :

    Elements are moved with this fraction of the current vector, in addition to currents and Stokes drift

    minval :

    1.0

    maxval :

    1.0

    array([[ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           ...,
           [ 1.,  1.,  1., ..., nan, nan, nan],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.]],
          shape=(3000, 49), dtype=float32)

  • terminal_velocity
    (trajectory, time)
    float32
    0.0 0.0 0.0 ... 0.0414 0.0413

    dtype :

    <class 'numpy.float32'>

    units :

    m/s

    description :

    Terminal rise/sinking velocity (buoyancy) in the ocean column

    minval :

    0.0

    maxval :

    0.04943559

    array([[0.        , 0.        , 0.        , ..., 0.        , 0.        ,
            0.        ],
           [0.        , 0.        , 0.        , ..., 0.        , 0.        ,
            0.        ],
           [0.        , 0.        , 0.        , ..., 0.04452778, 0.04452011,
            0.04451037],
           ...,
           [0.        , 0.        , 0.        , ...,        nan,        nan,
                   nan],
           [0.        , 0.        , 0.        , ..., 0.        , 0.        ,
            0.        ],
           [0.        , 0.        , 0.        , ..., 0.04148078, 0.04139896,
            0.04129854]], shape=(3000, 49), dtype=float32)

  • mass_oil
    (trajectory, time)
    float32
    0.3294 0.3294 ... 0.3294 0.3294

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.32936665

    maxval :

    0.32936665

    array([[0.32936665, 0.32936665, 0.32936665, ..., 0.32936665, 0.32936665,
            0.32936665],
           [0.32936665, 0.32936665, 0.32936665, ..., 0.32936665, 0.32936665,
            0.32936665],
           [0.32936665, 0.32936665, 0.32936665, ..., 0.32936665, 0.32936665,
            0.32936665],
           ...,
           [0.32936665, 0.32936665, 0.32936665, ...,        nan,        nan,
                   nan],
           [0.32936665, 0.32936665, 0.32936665, ..., 0.32936665, 0.32936665,
            0.32936665],
           [0.32936665, 0.32936665, 0.32936665, ..., 0.32936665, 0.32936665,
            0.32936665]], shape=(3000, 49), dtype=float32)

  • viscosity
    (trajectory, time)
    float32
    0.02169 0.101 ... 0.1043 0.1049

    dtype :

    <class 'numpy.float32'>

    units :

    m2/s

    description :

    Kinematic viscosity of oil (emulsion)

    minval :

    0.021692334

    maxval :

    0.10596327

    array([[0.02169233, 0.10098438, 0.10122959, ..., 0.10441887, 0.10464878,
            0.10480545],
           [0.02169233, 0.10092389, 0.10083557, ..., 0.10517527, 0.10522758,
            0.10551478],
           [0.02169233, 0.10103877, 0.10112443, ..., 0.10477154, 0.10497025,
            0.10521771],
           ...,
           [0.02169233, 0.10082236, 0.1009357 , ...,        nan,        nan,
                   nan],
           [0.02169233, 0.10071339, 0.10073697, ..., 0.10339156, 0.10404744,
            0.10432931],
           [0.02169233, 0.1008979 , 0.1011268 , ..., 0.10378194, 0.10429849,
            0.10489544]], shape=(3000, 49), dtype=float32)

  • density
    (trajectory, time)
    float32
    988.1 995.0 995.0 ... 995.2 995.2

    dtype :

    <class 'numpy.float32'>

    units :

    kg/m^3

    minval :

    988.1

    maxval :

    995.2458

    array([[988.1    , 995.0369 , 995.0474 , ..., 995.1821 , 995.19165,
            995.1982 ],
           [988.1    , 995.0343 , 995.0305 , ..., 995.21344, 995.2156 ,
            995.2274 ],
           [988.1    , 995.03925, 995.0429 , ..., 995.1968 , 995.20496,
            995.2152 ],
           ...,
           [988.1    , 995.0299 , 995.0348 , ...,       nan,       nan,
                  nan],
           [988.1    , 995.0252 , 995.02625, ..., 995.1392 , 995.1667 ,
            995.1784 ],
           [988.1    , 995.0332 , 995.043  , ..., 995.1556 , 995.1771 ,
            995.2019 ]], shape=(3000, 49), dtype=float32)

  • bulltime
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    dtype :

    <class 'numpy.float32'>

    units :

    s

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • interfacial_area
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    dtype :

    <class 'numpy.float32'>

    units :

    m2

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • mass_dispersed
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • mass_evaporated
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • mass_biodegraded
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • biodegradation_half_time_droplet
    (trajectory, time)
    float32
    1.0 1.0 1.0 1.0 ... 1.0 1.0 1.0 1.0

    dtype :

    <class 'numpy.float32'>

    units :

    Days

    description :

    Biodegradation half time in days for submerged oil droplets

    minval :

    1.0

    maxval :

    1.0

    array([[ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           ...,
           [ 1.,  1.,  1., ..., nan, nan, nan],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.]],
          shape=(3000, 49), dtype=float32)

  • biodegradation_half_time_slick
    (trajectory, time)
    float32
    3.0 3.0 3.0 3.0 ... 3.0 3.0 3.0 3.0

    dtype :

    <class 'numpy.float32'>

    units :

    Days

    description :

    Biodegradation half time in days for surface oil slick

    minval :

    3.0

    maxval :

    3.0

    array([[ 3.,  3.,  3., ...,  3.,  3.,  3.],
           [ 3.,  3.,  3., ...,  3.,  3.,  3.],
           [ 3.,  3.,  3., ...,  3.,  3.,  3.],
           ...,
           [ 3.,  3.,  3., ..., nan, nan, nan],
           [ 3.,  3.,  3., ...,  3.,  3.,  3.],
           [ 3.,  3.,  3., ...,  3.,  3.,  3.]],
          shape=(3000, 49), dtype=float32)

  • fraction_evaporated
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    dtype :

    <class 'numpy.float32'>

    units :

    1

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • water_fraction
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    dtype :

    <class 'numpy.float32'>

    units :

    1

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • oil_film_thickness
    (trajectory, time)
    float32
    0.001 0.001 0.001 ... 0.001 0.001

    dtype :

    <class 'numpy.float32'>

    units :

    m

    minval :

    0.001

    maxval :

    0.001

    array([[0.001, 0.001, 0.001, ..., 0.001, 0.001, 0.001],
           [0.001, 0.001, 0.001, ..., 0.001, 0.001, 0.001],
           [0.001, 0.001, 0.001, ..., 0.001, 0.001, 0.001],
           ...,
           [0.001, 0.001, 0.001, ...,   nan,   nan,   nan],
           [0.001, 0.001, 0.001, ..., 0.001, 0.001, 0.001],
           [0.001, 0.001, 0.001, ..., 0.001, 0.001, 0.001]],
          shape=(3000, 49), dtype=float32)

  • diameter
    (trajectory, time)
    float32
    0.0 0.0 0.0 ... 0.002087 0.002087

    dtype :

    <class 'numpy.float32'>

    units :

    m

    minval :

    0.0

    maxval :

    0.002994377

    array([[0.        , 0.        , 0.        , ..., 0.        , 0.        ,
            0.        ],
           [0.        , 0.        , 0.        , ..., 0.        , 0.        ,
            0.        ],
           [0.        , 0.        , 0.        , ..., 0.00245066, 0.00245066,
            0.00245066],
           ...,
           [0.        , 0.        , 0.        , ...,        nan,        nan,
                   nan],
           [0.        , 0.        , 0.        , ..., 0.        , 0.        ,
            0.        ],
           [0.        , 0.        , 0.        , ..., 0.00208711, 0.00208711,
            0.00208711]], shape=(3000, 49), dtype=float32)

  • x_sea_water_velocity
    (trajectory, time)
    float32
    0.1007 -0.01427 ... 0.1089 0.08257

    minval :

    -0.7262247

    maxval :

    0.5493309

    array([[ 0.10074912, -0.01426796, -0.06266311, ...,  0.12932198,
             0.28755867,  0.20240861],
           [-0.04822906,  0.04223454, -0.05295211, ...,  0.2508487 ,
             0.16367705,  0.13198194],
           [ 0.00444389,  0.04328514, -0.09077436, ...,  0.14887647,
             0.14664458,  0.2686618 ],
           ...,
           [-0.0429161 ,  0.01841989, -0.15595004, ...,         nan,
                    nan,         nan],
           [-0.05791224, -0.03845641, -0.12299405, ...,  0.1458894 ,
             0.02405779,  0.08184233],
           [ 0.03159171,  0.06003619, -0.09753326, ...,  0.15998864,
             0.10892793,  0.08257117]], shape=(3000, 49), dtype=float32)

  • y_sea_water_velocity
    (trajectory, time)
    float32
    -0.1815 -0.2706 ... 0.1041 0.02163

    minval :

    -0.8292342

    maxval :

    0.8147809

    array([[-0.18147929, -0.27061188, -0.24978775, ...,  0.25289324,
             0.11807695, -0.08527417],
           [-0.2313428 , -0.24378088, -0.1605759 , ...,  0.26999027,
             0.10280003, -0.02076976],
           [-0.25855565, -0.19282983, -0.14079873, ...,  0.2255513 ,
             0.06317879, -0.03766109],
           ...,
           [-0.3549225 , -0.24453086, -0.18937144, ...,         nan,
                    nan,         nan],
           [-0.25493932, -0.21285008, -0.08490583, ...,  0.15220007,
             0.09335957, -0.056163  ],
           [-0.24261537, -0.20716211, -0.14838296, ...,  0.14293559,
             0.10406383,  0.02162871]], shape=(3000, 49), dtype=float32)

  • x_wind
    (trajectory, time)
    float32
    -1.102 -1.847 ... 4.265 2.531

    minval :

    -6.374165

    maxval :

    10.261469

    array([[-1.1019453 , -1.8468169 , -1.3370682 , ...,  2.207657  ,
             3.7071292 ,  3.8346965 ],
           [-1.3570561 , -1.2017018 , -1.1421083 , ...,  2.164063  ,
             2.4810421 ,  2.706505  ],
           [-2.1578693 , -1.3550948 , -1.639335  , ...,  2.8607547 ,
             3.166791  ,  3.2044706 ],
           ...,
           [-1.188623  , -2.1107028 , -0.37381288, ...,         nan,
                    nan,         nan],
           [-1.8527684 , -1.7975494 , -0.03592263, ...,  3.1074443 ,
             2.8971307 ,  3.4296553 ],
           [-1.6656177 , -0.2440132 , -1.1404271 , ...,  2.7626054 ,
             4.265251  ,  2.5306563 ]], shape=(3000, 49), dtype=float32)

  • y_wind
    (trajectory, time)
    float32
    5.992 4.363 3.201 ... 2.899 3.251

    minval :

    -10.3195095

    maxval :

    6.5646896

    array([[5.9917407, 4.3629856, 3.2013283, ..., 3.5160446, 2.7795568,
            3.4953568],
           [5.0596557, 4.5223165, 2.6092296, ..., 3.3644402, 3.0006194,
            3.0313122],
           [4.567189 , 4.140402 , 3.1727474, ..., 3.3594153, 3.2137835,
            2.3235185],
           ...,
           [4.7657757, 5.0153527, 2.799533 , ...,       nan,       nan,
                  nan],
           [4.7007427, 3.9279633, 3.4205945, ..., 2.8545194, 2.3872695,
            2.9669497],
           [4.7971163, 3.8558965, 3.9446013, ..., 2.6942353, 2.8985066,
            3.2510462]], shape=(3000, 49), dtype=float32)

  • sea_surface_height
    (trajectory, time)
    float32
    -0.1906 0.02426 ... -0.6478 -0.5999

    minval :

    -1.1798413

    maxval :

    0.26305163

    array([[-0.19055861,  0.02425914, -0.01827466, ..., -0.7484431 ,
            -0.6662497 , -0.6276637 ],
           [-0.190366  ,  0.02687886, -0.01514066, ..., -0.73899066,
            -0.6510111 , -0.60940677],
           [-0.19061416,  0.02623827, -0.01748022, ..., -0.74090827,
            -0.6532786 , -0.6115364 ],
           ...,
           [-0.1903812 ,  0.02671539, -0.01491422, ...,         nan,
                    nan,         nan],
           [-0.18998097,  0.02827901, -0.01439438, ..., -0.7411145 ,
            -0.6500789 , -0.6054951 ],
           [-0.19028468,  0.02733204, -0.01633159, ..., -0.7395738 ,
            -0.6477581 , -0.59990674]], shape=(3000, 49), dtype=float32)

  • sea_surface_wave_significant_height
    (trajectory, time)
    float32
    0.913 0.5522 ... 0.6542 0.4175

    minval :

    7.6191036e-06

    maxval :

    2.6227443

    array([[0.91303474, 0.5521808 , 0.29609185, ..., 0.42401347, 0.52813107,
            0.6622914 ],
           [0.6750661 , 0.5386276 , 0.19956724, ..., 0.39366454, 0.37291846,
            0.40624493],
           [0.6276842 , 0.46688855, 0.31374213, ..., 0.47895184, 0.5007814 ,
            0.38541728],
           ...,
           [0.59348583, 0.72837716, 0.19623716, ...,        nan,        nan,
                   nan],
           [0.62803143, 0.4590379 , 0.2878632 , ..., 0.43799046, 0.34667358,
            0.50590706],
           [0.63435054, 0.36721602, 0.41476712, ..., 0.36631596, 0.65420526,
            0.41754863]], shape=(3000, 49), dtype=float32)

  • sea_surface_wave_stokes_drift_x_velocity
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • sea_surface_wave_stokes_drift_y_velocity
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • sea_surface_wave_period_at_variance_spectral_density_maximum
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment
    (trajectory, time)
    float32
    5.206 4.048 2.964 ... 4.406 3.52

    minval :

    0.015037704

    maxval :

    8.822823

    array([[5.2056293, 4.048277 , 2.9644408, ..., 3.5474756, 3.959136 ,
            4.433576 ],
           [4.4761305, 3.998286 , 2.4337401, ..., 3.4181628, 3.3268752,
            3.4723506],
           [4.3161864, 3.722513 , 3.0515182, ..., 3.7702966, 3.8552601,
            3.382168 ],
           ...,
           [4.1969595, 4.649516 , 2.4133494, ...,       nan,       nan,
                  nan],
           [4.3173804, 3.6910834, 2.9229584, ..., 3.6054702, 3.2076719,
            3.8749397],
           [4.3390455, 3.301341 , 3.508583 , ..., 3.2972927, 4.406427 ,
            3.520328 ]], shape=(3000, 49), dtype=float32)

  • sea_ice_area_fraction
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • sea_ice_x_velocity
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • sea_ice_y_velocity
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • sea_water_temperature
    (trajectory, time)
    float32
    1.453 1.435 1.422 ... 1.21 274.3

    minval :

    1.1379567

    maxval :

    274.5257

    array([[  1.4529235,   1.4349738,   1.4218116, ...,   1.2103817,
              1.2124081, 274.35275  ],
           [  1.454606 ,   1.4436466,   1.4542804, ...,   1.1815763,
              1.1690385, 274.30887  ],
           [  1.4548976,   1.4380245,   1.4369601, ...,   1.2054601,
              1.1897225, 274.3272   ],
           ...,
           [  1.4559354,   1.4587994,   1.4486611, ...,         nan,
                    nan,         nan],
           [  1.4556463,   1.4660815,   1.4539553, ...,   1.2740865,
              1.2289474, 274.3824   ],
           [  1.4551182,   1.4479253,   1.4377068, ...,   1.2464824,
              1.2096264, 274.34717  ]], shape=(3000, 49), dtype=float32)

  • sea_water_salinity
    (trajectory, time)
    float32
    33.65 33.64 33.64 ... 33.66 33.65

    minval :

    33.59873

    maxval :

    33.766823

    array([[33.64523 , 33.641907, 33.637478, ..., 33.682117, 33.681587,
            33.67833 ],
           [33.645596, 33.64132 , 33.64201 , ..., 33.65188 , 33.64633 ,
            33.64208 ],
           [33.64585 , 33.640446, 33.639313, ..., 33.663734, 33.65857 ,
            33.653324],
           ...,
           [33.646   , 33.645382, 33.64206 , ...,       nan,       nan,
                  nan],
           [33.645645, 33.645077, 33.642693, ..., 33.69135 , 33.673347,
            33.676224],
           [33.64569 , 33.641888, 33.64043 , ..., 33.67842 , 33.6617  ,
            33.653893]], shape=(3000, 49), dtype=float32)

  • sea_floor_depth_below_sea_level
    (trajectory, time)
    float32
    170.4 165.2 165.8 ... 121.8 129.6

    minval :

    11.50687

    maxval :

    259.43405

    array([[170.42697 , 165.24956 , 165.78703 , ...,  79.495224,  89.52262 ,
             93.308426],
           [170.93436 , 172.86862 , 161.56934 , ..., 100.053535, 109.99035 ,
            111.93212 ],
           [170.35619 , 174.84906 , 167.78622 , ...,  97.015   , 107.41916 ,
            107.92146 ],
           ...,
           [170.93929 , 160.27583 , 159.43793 , ...,        nan,        nan,
                   nan],
           [171.87881 , 162.3581  , 158.42264 , ..., 108.67472 , 117.80816 ,
            118.45253 ],
           [171.14272 , 171.27951 , 162.99707 , ..., 110.97813 , 121.822105,
            129.64128 ]], shape=(3000, 49), dtype=float32)

  • horizontal_diffusivity
    (trajectory, time)
    float32
    10.0 10.0 10.0 ... 10.0 10.0 10.0

    minval :

    10.0

    maxval :

    10.0

    array([[10., 10., 10., ..., 10., 10., 10.],
           [10., 10., 10., ..., 10., 10., 10.],
           [10., 10., 10., ..., 10., 10., 10.],
           ...,
           [10., 10., 10., ..., nan, nan, nan],
           [10., 10., 10., ..., 10., 10., 10.],
           [10., 10., 10., ..., 10., 10., 10.]],
          shape=(3000, 49), dtype=float32)

  • ocean_vertical_diffusivity
    (trajectory, time)
    float32
    0.02 0.02 0.0002664 ... 0.02 0.02

    minval :

    0.00012895565

    maxval :

    0.033379786

    array([[0.02      , 0.02      , 0.00026641, ..., 0.02      , 0.02      ,
            0.02      ],
           [0.02      , 0.02      , 0.00030162, ..., 0.02      , 0.02      ,
            0.02      ],
           [0.02      , 0.02      , 0.00027505, ..., 0.02      , 0.02      ,
            0.02      ],
           ...,
           [0.02      , 0.02      , 0.00029549, ...,        nan,        nan,
                   nan],
           [0.02      , 0.02      , 0.00030152, ..., 0.02      , 0.02      ,
            0.02      ],
           [0.02      , 0.02      , 0.00028166, ..., 0.02      , 0.02      ,
            0.02      ]], shape=(3000, 49), dtype=float32)

  • land_binary_mask
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.]],
          shape=(3000, 49), dtype=float32)

  • ocean_mixed_layer_thickness
    (trajectory, time)
    float32
    50.0 50.0 50.0 ... 50.0 50.0 50.0

    minval :

    50.0

    maxval :

    50.0

    array([[50., 50., 50., ..., 50., 50., 50.],
           [50., 50., 50., ..., 50., 50., 50.],
           [50., 50., 50., ..., 50., 50., 50.],
           ...,
           [50., 50., 50., ..., nan, nan, nan],
           [50., 50., 50., ..., 50., 50., 50.],
           [50., 50., 50., ..., 50., 50., 50.]],
          shape=(3000, 49), dtype=float32)

• Attributes: (167)

  Conventions :

  CF-1.11, ACDD-1.3

  standard_name_vocabulary :

  CF Standard Name Table v85

  featureType :

  trajectory

  title :

  OpenDrift trajectory simulation

  summary :

  Output from simulation with OpenDrift framework

  keywords :

  trajectory, drift, lagrangian, simulation

  history :

  Created 2026-10-05 15:36:33.304432

  date_created :

  2026-10-05T15:36:33.304443

  source :

  Output from simulation with OpenDrift

  model_url :

  https://github.com/OpenDrift/opendrift

  opendrift_class :

  OpenOil

  opendrift_module :

  opendrift.models.openoil.openoil

  readers :

  odict_keys(['LazyReader: https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be', 'constant_reader', 'global_landmask'])

  time_coverage_start :

  2026-09-28 15:36:27.064652

  time_step_calculation :

  0:30:00

  time_step_output :

  1:00:00

  config_environment:constant:x_sea_water_velocity :

  None

  config_environment:fallback:x_sea_water_velocity :

  None

  config_environment:constant:y_sea_water_velocity :

  None

  config_environment:fallback:y_sea_water_velocity :

  None

  config_environment:constant:x_wind :

  None

  config_environment:fallback:x_wind :

  None

  config_environment:constant:y_wind :

  None

  config_environment:fallback:y_wind :

  None

  config_environment:constant:sea_surface_height :

  None

  config_environment:fallback:sea_surface_height :

  0

  config_environment:constant:upward_sea_water_velocity :

  None

  config_environment:fallback:upward_sea_water_velocity :

  0

  config_environment:constant:sea_surface_wave_significant_height :

  None

  config_environment:fallback:sea_surface_wave_significant_height :

  0

  config_environment:constant:sea_surface_wave_stokes_drift_x_velocity :

  None

  config_environment:fallback:sea_surface_wave_stokes_drift_x_velocity :

  0

  config_environment:constant:sea_surface_wave_stokes_drift_y_velocity :

  None

  config_environment:fallback:sea_surface_wave_stokes_drift_y_velocity :

  0

  config_environment:constant:sea_surface_wave_period_at_variance_spectral_density_maximum :

  None

  config_environment:fallback:sea_surface_wave_period_at_variance_spectral_density_maximum :

  0

  config_environment:constant:sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment :

  None

  config_environment:fallback:sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment :

  0

  config_environment:constant:sea_ice_area_fraction :

  None

  config_environment:fallback:sea_ice_area_fraction :

  0

  config_environment:constant:sea_ice_x_velocity :

  None

  config_environment:fallback:sea_ice_x_velocity :

  0

  config_environment:constant:sea_ice_y_velocity :

  None

  config_environment:fallback:sea_ice_y_velocity :

  0

  config_environment:constant:sea_water_temperature :

  None

  config_environment:fallback:sea_water_temperature :

  10

  config_environment:constant:sea_water_salinity :

  None

  config_environment:fallback:sea_water_salinity :

  34

  config_environment:constant:sea_floor_depth_below_sea_level :

  None

  config_environment:fallback:sea_floor_depth_below_sea_level :

  10000

  config_environment:constant:horizontal_diffusivity :

  10.0

  config_environment:fallback:horizontal_diffusivity :

  0

  config_environment:constant:ocean_vertical_diffusivity :

  None

  config_environment:fallback:ocean_vertical_diffusivity :

  0.02

  config_environment:constant:land_binary_mask :

  None

  config_environment:fallback:land_binary_mask :

  None

  config_environment:constant:ocean_mixed_layer_thickness :

  None

  config_environment:fallback:ocean_mixed_layer_thickness :

  50

  config_general:use_auto_landmask :

  True

  config_drift:current_uncertainty :

  0.05

  config_drift:current_uncertainty_uniform :

  0

  config_drift:max_speed :

  1.3

  config_readers:max_number_of_fails :

  1

  config_general:simulation_name :

  config_general:coastline_action :

  stranding

  config_general:coastline_approximation_precision :

  0.001

  config_general:time_step_minutes :

  60

  config_general:time_step_output_minutes :

  None

  config_seed:ocean_only :

  True

  config_seed:number :

  1

  config_drift:max_age_seconds :

  None

  config_drift:advection_scheme :

  euler

  config_drift:profiles_depth :

  50

  config_drift:wind_uncertainty :

  0.5

  config_drift:relative_wind :

  False

  config_drift:deactivate_north_of :

  None

  config_drift:deactivate_south_of :

  None

  config_drift:deactivate_east_of :

  None

  config_drift:deactivate_west_of :

  None

  config_seed:origin_marker :

  0

  config_seed:z :

  0

  config_seed:wind_drift_factor :

  0.03

  config_seed:current_drift_factor :

  1

  config_seed:terminal_velocity :

  0.0

  config_seed:oil_film_thickness :

  0.001

  config_drift:vertical_advection :

  False

  config_drift:vertical_advection_at_surface :

  False

  config_drift:water_column_stretching :

  False

  config_drift:vertical_advection_correction :

  False

  config_drift:vertical_mixing :

  True

  config_drift:vertical_mixing_at_surface :

  False

  config_vertical_mixing:timestep :

  60

  config_vertical_mixing:diffusivitymodel :

  environment

  config_vertical_mixing:background_diffusivity :

  1.2e-05

  config_vertical_mixing:TSprofiles :

  False

  config_drift:wind_drift_depth :

  0.1

  config_drift:stokes_drift :

  True

  config_drift:stokes_drift_profile :

  Phillips

  config_drift:use_tabularised_stokes_drift :

  False

  config_drift:tabularised_stokes_drift_fetch :

  25000

  config_general:seafloor_action :

  lift_to_seafloor

  config_drift:truncate_ocean_model_below_m :

  3.0

  config_seed:seafloor :

  False

  config_seed:m3_per_hour :

  1

  config_seed:droplet_size_distribution :

  uniform

  config_seed:droplet_diameter_mu :

  0.001

  config_seed:droplet_diameter_sigma :

  0.0005

  config_seed:droplet_diameter_min_subsea :

  0.0005

  config_seed:droplet_diameter_max_subsea :

  0.005

  config_processes:dispersion :

  False

  config_processes:evaporation :

  False

  config_processes:emulsification :

  False

  config_processes:biodegradation :

  False

  config_biodegradation:method :

  Adcroft

  config_processes:update_oilfilm_thickness :

  False

  config_wave_entrainment:droplet_size_distribution :

  Johansen et al. (2015)

  config_wave_entrainment:entrainment_rate :

  Li et al. (2017)

  config_seed:oil_type :

  GENERIC BUNKER C

  opendrift_version :

  1.14.11

  seed_lon :

  19.1909

  seed_lat :

  79.5986

  seed_radius :

  50.0

  seed_time :

  2026-09-28 15:36:27.064652

  seed_number :

  3000.0

  seed_z :

  0.0

  seed_m3_per_hour :

  1.0

  seed_oiltype :

  GENERIC BUNKER C

  seed_geojson :

  {"features": [], "type": "FeatureCollection"}

  simulation_time :

  2026-10-05 15:36:33.424116

  reader_x_sea_water_velocity :

  ['https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be']

  reader_y_sea_water_velocity :

  ['https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be']

  reader_x_wind :

  ['https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be']

  reader_y_wind :

  ['https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be']

  reader_sea_surface_height :

  ['https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be']

  reader_sea_surface_wave_significant_height :

  0

  reader_sea_surface_wave_stokes_drift_x_velocity :

  0

  reader_sea_surface_wave_stokes_drift_y_velocity :

  0

  reader_sea_surface_wave_period_at_variance_spectral_density_maximum :

  0

  reader_sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment :

  0

  reader_sea_ice_area_fraction :

  ['https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be']

  reader_sea_ice_x_velocity :

  ['https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be']

  reader_sea_ice_y_velocity :

  ['https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be']

  reader_sea_water_temperature :

  ['https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be']

  reader_sea_water_salinity :

  ['https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be']

  reader_sea_floor_depth_below_sea_level :

  ['https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be']

  reader_horizontal_diffusivity :

  10.0

  reader_ocean_vertical_diffusivity :

  ['https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be']

  reader_land_binary_mask :

  ['global_landmask']

  reader_ocean_mixed_layer_thickness :

  50

  time_coverage_end :

  2026-09-30 15:36:27.064652

  time_coverage_duration :

  P2DT0H0M0S

  time_coverage_resolution :

  P0DT0H30M0S

  performance :

  -------------------- Reader performance: -------------------- constant_reader 0:00:00.0 total 0:00:00.0 preparing 0:00:00.0 reading 0:00:00.0 masking -------------------- global_landmask 0:00:00.0 total 0:00:00.0 preparing 0:00:00.0 reading 0:00:00.0 masking -------------------- https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be 0:01:28.7 total 0:00:00.0 preparing 0:01:27.8 reading 0:00:01.3 interpolation 0:00:00.1 interpolation_time 0:00:00.9 rotating vectors 0:00:00.0 masking -------------------- Performance: 1:47.5 total time 6.2 configuration 0.1 preparing main loop 0.0 moving elements to ocean 1:41.1 main loop 5.0 updating elements 0.0 oil weathering 0.0 updating viscosities 0.0 updating densities 4.6 vertical mixing 0.0 cleaning up --------------------

  geospatial_bounds_crs :

  EPSG:4326

  geospatial_bounds_vertical_crs :

  EPSG:5831

  geospatial_lat_min :

  79.37862

  geospatial_lat_max :

  79.65727

  geospatial_lat_units :

  degrees_north

  geospatial_lat_resolution :

  point

  geospatial_lon_min :

  18.550377

  geospatial_lon_max :

  19.799967

  geospatial_lon_units :

  degrees_east

  geospatial_lon_resolution :

  point

  runtime :

  0:01:47.517314

  geospatial_vertical_min :

  -199.42229

  geospatial_vertical_max :

  0.0

  geospatial_vertical_positive :

  up

Print and plot results

print(o)
o.animation(background='sea_ice_area_fraction', cmap=cmocean.cm.ice,
            vmin=0, vmax=1, bgalpha=1, fast=False)

===========================
--------------------
Reader performance:
--------------------
constant_reader
 0:00:00.0  total
 0:00:00.0  preparing
 0:00:00.0  reading
 0:00:00.0  masking
--------------------
global_landmask
 0:00:00.0  total
 0:00:00.0  preparing
 0:00:00.0  reading
 0:00:00.0  masking
--------------------
https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
 0:01:28.7  total
 0:00:00.0  preparing
 0:01:27.8  reading
 0:00:01.3  interpolation
 0:00:00.1  interpolation_time
 0:00:00.9  rotating vectors
 0:00:00.0  masking
--------------------
Performance:
 1:47.5 total time
    6.2 configuration
    0.1 preparing main loop
      0.0 moving elements to ocean
 1:41.1 main loop
      5.0 updating elements
        0.0 oil weathering
          0.0 updating viscosities
          0.0 updating densities
        4.6 vertical mixing
    0.0 cleaning up
--------------------
===========================
Model:  OpenOil     (OpenDrift version 1.14.11)
        1276 active Oil particles  (1724 deactivated, 0 scheduled)
-------------------
Environment variables:
  -----
  horizontal_diffusivity
     1) constant_reader
  -----
  land_binary_mask
     1) global_landmask
  -----
  ocean_vertical_diffusivity
  sea_floor_depth_below_sea_level
  sea_ice_area_fraction
  sea_ice_x_velocity
  sea_ice_y_velocity
  sea_surface_height
  sea_water_salinity
  sea_water_temperature
  x_sea_water_velocity
  x_wind
  y_sea_water_velocity
  y_wind
     1) https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
  -----
Readers not added for the following variables:
  ocean_mixed_layer_thickness
  sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment
  sea_surface_wave_period_at_variance_spectral_density_maximum
  sea_surface_wave_significant_height
  sea_surface_wave_stokes_drift_x_velocity
  sea_surface_wave_stokes_drift_y_velocity

Discarded readers:

Time:
        Start: 2026-09-28 15:36:27.064652 UTC
        Present: 2026-09-30 15:36:27.064652 UTC
        Calculation steps: 96 * 0:30:00 - total time: 2 days, 0:00:00
        Output steps: 49 * 1:00:00
===========================

15:38:26 INFO    opendrift:4823: Saving animation to /root/project/docs/source/gallery/animations/example_oil_ice_0.gif...
15:38:51 INFO    opendrift:3241: Time to make animation: 0:00:36.301631

o.plot(background='sea_ice_area_fraction', cmap=cmocean.cm.ice,
       vmin=0, vmax=1, bgalpha=1, fast=False)

(<GeoAxes: title={'center': 'OpenDrift - OpenOil (GENERIC BUNKER C)\n2026-09-28 15:36 to 2026-09-30 15:36 UTC (49 steps)'}>, <Figure size 612.224x1100 with 2 Axes>)