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)

10:34:31 INFO    opendrift:576: OpenDriftSimulation initialised (version 1.14.8 / v1.14.8-2-g680a748)

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('drift: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))

10:34:31 INFO    opendrift.models.openoil.openoil:1708: Oil type not specified, using default: GENERIC BUNKER C
10:34:31 INFO    opendrift.models.openoil.adios.dirjs:86: Querying ADIOS database for oil: GENERIC BUNKER C
10:34:31 INFO    opendrift.models.openoil.openoil:1717: Using density 988.1 and viscosity 0.02169233387797564 of oiltype GENERIC BUNKER C
10:34:31 INFO    opendrift.models.basemodel.environment:203: Adding a global landmask from GSHHG
10:34:35 INFO    opendrift.models.basemodel.environment:227: Fallback values will be used for the following variables which have no readers:
10:34:35 INFO    opendrift.models.basemodel.environment:230:    sea_surface_height: 0.000000
10:34:35 INFO    opendrift.models.basemodel.environment:230:    upward_sea_water_velocity: 0.000000
10:34:35 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_significant_height: 0.000000
10:34:35 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
10:34:35 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
10:34:35 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
10:34:35 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
10:34:35 INFO    opendrift.models.basemodel.environment:230:    sea_ice_area_fraction: 0.000000
10:34:35 INFO    opendrift.models.basemodel.environment:230:    sea_ice_x_velocity: 0.000000
10:34:35 INFO    opendrift.models.basemodel.environment:230:    sea_ice_y_velocity: 0.000000
10:34:35 INFO    opendrift.models.basemodel.environment:230:    sea_water_temperature: 10.000000
10:34:35 INFO    opendrift.models.basemodel.environment:230:    sea_water_salinity: 34.000000
10:34:35 INFO    opendrift.models.basemodel.environment:230:    sea_floor_depth_below_sea_level: 10000.000000
10:34:35 INFO    opendrift.models.basemodel.environment:230:    ocean_vertical_diffusivity: 0.020000
10:34:35 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)

10:34:35 INFO    opendrift:1810: Skipping environment variable upward_sea_water_velocity because of condition ['drift:vertical_advection', 'is', False]
10:34:35 INFO    opendrift:1821: Storing previous values of element property lon because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
10:34:35 INFO    opendrift:1821: Storing previous values of element property lat because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
10:34:35 INFO    opendrift:952: Using existing reader for land_binary_mask to move elements to ocean
10:34:35 INFO    opendrift:982: All points are in ocean
10:34:35 INFO    opendrift.models.openoil.openoil:692: Oil-water surface tension is 0.035935 Nm
10:34:35 INFO    opendrift.models.openoil.openoil:705: Max water fraction not available for GENERIC BUNKER C, using default
10:34:35 INFO    opendrift:2118: 2026-01-15 10:34:31.485059 - step 1 of 96 - 3000 active elements (0 deactivated)
10:34:35 INFO    opendrift.readers:63: Opening file with xr.open_dataset
10:34:37 INFO    opendrift.readers.reader_netCDF_CF_generic:332: Detected dimensions: {'x': 'X', 'y': 'Y', 'z': 'depth', 'time': 'time'}
10:34:37 INFO    opendrift.readers.basereader:176: Variable x_sea_water_velocity will be rotated from eastward_sea_water_velocity
10:34:37 INFO    opendrift.readers.basereader:176: Variable y_sea_water_velocity will be rotated from northward_sea_water_velocity
10:34:39 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0071644546 to -0.0010733843) found for ocean_vertical_diffusivity, replacing with NaN
10:34:39 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:34:41 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.006821165 to -0.0011060748) found for ocean_vertical_diffusivity, replacing with NaN
10:34:41 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:34:41 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:41 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:42 INFO    opendrift:2118: 2026-01-15 11:04:31.485059 - step 2 of 96 - 3000 active elements (0 deactivated)
10:34:44 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.006721994 to -0.0013423299) found for ocean_vertical_diffusivity, replacing with NaN
10:34:44 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:34:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:44 INFO    opendrift:2118: 2026-01-15 11:34:31.485059 - step 3 of 96 - 3000 active elements (0 deactivated)
10:34:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:44 INFO    opendrift:2118: 2026-01-15 12:04:31.485059 - step 4 of 96 - 3000 active elements (0 deactivated)
10:34:46 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0058689085 to -0.0009474999) found for ocean_vertical_diffusivity, replacing with NaN
10:34:46 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:34:46 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:46 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:46 INFO    opendrift:2118: 2026-01-15 12:34:31.485059 - step 5 of 96 - 3000 active elements (0 deactivated)
10:34:46 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:46 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:47 INFO    opendrift:2118: 2026-01-15 13:04:31.485059 - step 6 of 96 - 3000 active elements (0 deactivated)
10:34:49 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.005230071 to -0.0012837136) found for ocean_vertical_diffusivity, replacing with NaN
10:34:49 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:34:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:49 INFO    opendrift:2118: 2026-01-15 13:34:31.485059 - step 7 of 96 - 3000 active elements (0 deactivated)
10:34:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:49 INFO    opendrift:2118: 2026-01-15 14:04:31.485059 - step 8 of 96 - 3000 active elements (0 deactivated)
10:34:51 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0060158144 to -0.0012293737) found for ocean_vertical_diffusivity, replacing with NaN
10:34:51 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:34:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:51 INFO    opendrift:2118: 2026-01-15 14:34:31.485059 - step 9 of 96 - 3000 active elements (0 deactivated)
10:34:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:52 INFO    opendrift:2118: 2026-01-15 15:04:31.485059 - step 10 of 96 - 3000 active elements (0 deactivated)
10:34:54 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.008228783 to -0.0014546507) found for ocean_vertical_diffusivity, replacing with NaN
10:34:54 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:34:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:54 INFO    opendrift:2118: 2026-01-15 15:34:31.485059 - step 11 of 96 - 3000 active elements (0 deactivated)
10:34:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:55 INFO    opendrift:2118: 2026-01-15 16:04:31.485059 - step 12 of 96 - 3000 active elements (0 deactivated)
10:34:57 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.008884608 to -0.0012681276) found for ocean_vertical_diffusivity, replacing with NaN
10:34:57 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:34:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:57 INFO    opendrift:2118: 2026-01-15 16:34:31.485059 - step 13 of 96 - 3000 active elements (0 deactivated)
10:34:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:57 INFO    opendrift:2118: 2026-01-15 17:04:31.485059 - step 14 of 96 - 3000 active elements (0 deactivated)
10:34:59 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0092358 to -0.0012371199) found for ocean_vertical_diffusivity, replacing with NaN
10:34:59 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:34:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:34:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:00 INFO    opendrift:2118: 2026-01-15 17:34:31.485059 - step 15 of 96 - 3000 active elements (0 deactivated)
10:35:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:00 INFO    opendrift:2118: 2026-01-15 18:04:31.485059 - step 16 of 96 - 3000 active elements (0 deactivated)
10:35:02 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0086290045 to -0.001291035) found for ocean_vertical_diffusivity, replacing with NaN
10:35:02 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:02 INFO    opendrift:2118: 2026-01-15 18:34:31.485059 - step 17 of 96 - 3000 active elements (0 deactivated)
10:35:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:03 INFO    opendrift:2118: 2026-01-15 19:04:31.485059 - step 18 of 96 - 3000 active elements (0 deactivated)
10:35:05 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0074652107 to -0.0013493352) found for ocean_vertical_diffusivity, replacing with NaN
10:35:05 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:05 INFO    opendrift:2118: 2026-01-15 19:34:31.485059 - step 19 of 96 - 3000 active elements (0 deactivated)
10:35:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:05 INFO    opendrift:2118: 2026-01-15 20:04:31.485059 - step 20 of 96 - 3000 active elements (0 deactivated)
10:35:08 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.00733511 to -0.0014242479) found for ocean_vertical_diffusivity, replacing with NaN
10:35:08 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:08 INFO    opendrift:2118: 2026-01-15 20:34:31.485059 - step 21 of 96 - 3000 active elements (0 deactivated)
10:35:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:09 INFO    opendrift:2118: 2026-01-15 21:04:31.485059 - step 22 of 96 - 3000 active elements (0 deactivated)
10:35:11 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0073391083 to -0.0012233868) found for ocean_vertical_diffusivity, replacing with NaN
10:35:11 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:11 INFO    opendrift:2118: 2026-01-15 21:34:31.485059 - step 23 of 96 - 3000 active elements (0 deactivated)
10:35:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:11 INFO    opendrift:2118: 2026-01-15 22:04:31.485059 - step 24 of 96 - 3000 active elements (0 deactivated)
10:35:13 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0070821294 to -0.00095224794) found for ocean_vertical_diffusivity, replacing with NaN
10:35:13 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:13 INFO    opendrift:2118: 2026-01-15 22:34:31.485059 - step 25 of 96 - 3000 active elements (0 deactivated)
10:35:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:14 INFO    opendrift:2118: 2026-01-15 23:04:31.485059 - step 26 of 96 - 3000 active elements (0 deactivated)
10:35:15 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.005810605 to -0.00079130434) found for ocean_vertical_diffusivity, replacing with NaN
10:35:15 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:16 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:16 INFO    opendrift:2118: 2026-01-15 23:34:31.485059 - step 27 of 96 - 3000 active elements (0 deactivated)
10:35:16 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:16 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:16 INFO    opendrift:2118: 2026-01-16 00:04:31.485059 - step 28 of 96 - 3000 active elements (0 deactivated)
10:35:18 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0048449333 to -0.00089195895) found for ocean_vertical_diffusivity, replacing with NaN
10:35:18 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:18 INFO    opendrift:2118: 2026-01-16 00:34:31.485059 - step 29 of 96 - 3000 active elements (0 deactivated)
10:35:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:18 INFO    opendrift:2118: 2026-01-16 01:04:31.485059 - step 30 of 96 - 3000 active elements (0 deactivated)
10:35:21 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0050115865 to -0.0007800259) found for ocean_vertical_diffusivity, replacing with NaN
10:35:21 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:21 INFO    opendrift:2118: 2026-01-16 01:34:31.485059 - step 31 of 96 - 3000 active elements (0 deactivated)
10:35:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:21 INFO    opendrift:2118: 2026-01-16 02:04:31.485059 - step 32 of 96 - 3000 active elements (0 deactivated)
10:35:23 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0065748426 to -0.00060247787) found for ocean_vertical_diffusivity, replacing with NaN
10:35:23 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:24 INFO    opendrift:2118: 2026-01-16 02:34:31.485059 - step 33 of 96 - 3000 active elements (0 deactivated)
10:35:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:24 INFO    opendrift:2118: 2026-01-16 03:04:31.485059 - step 34 of 96 - 3000 active elements (0 deactivated)
10:35:26 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.008386606 to -0.0012679261) found for ocean_vertical_diffusivity, replacing with NaN
10:35:26 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:26 INFO    opendrift:2118: 2026-01-16 03:34:31.485059 - step 35 of 96 - 3000 active elements (0 deactivated)
10:35:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:26 INFO    opendrift:2118: 2026-01-16 04:04:31.485059 - step 36 of 96 - 3000 active elements (0 deactivated)
10:35:28 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.009547254 to -0.0015198219) found for ocean_vertical_diffusivity, replacing with NaN
10:35:28 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:29 INFO    opendrift:2118: 2026-01-16 04:34:31.485059 - step 37 of 96 - 3000 active elements (0 deactivated)
10:35:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:29 INFO    opendrift:2118: 2026-01-16 05:04:31.485059 - step 38 of 96 - 3000 active elements (0 deactivated)
10:35:31 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.009827602 to -0.0014912843) found for ocean_vertical_diffusivity, replacing with NaN
10:35:31 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:32 INFO    opendrift:2118: 2026-01-16 05:34:31.485059 - step 39 of 96 - 3000 active elements (0 deactivated)
10:35:32 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:32 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:32 INFO    opendrift:2118: 2026-01-16 06:04:31.485059 - step 40 of 96 - 3000 active elements (0 deactivated)
10:35:34 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.008239541 to -0.0014198819) found for ocean_vertical_diffusivity, replacing with NaN
10:35:34 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:34 INFO    opendrift:2118: 2026-01-16 06:34:31.485059 - step 41 of 96 - 3000 active elements (0 deactivated)
10:35:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:34 INFO    opendrift:2118: 2026-01-16 07:04:31.485059 - step 42 of 96 - 3000 active elements (0 deactivated)
10:35:37 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.00732667 to -0.0010894116) found for ocean_vertical_diffusivity, replacing with NaN
10:35:37 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:37 INFO    opendrift:2118: 2026-01-16 07:34:31.485059 - step 43 of 96 - 3000 active elements (0 deactivated)
10:35:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:37 INFO    opendrift:2118: 2026-01-16 08:04:31.485059 - step 44 of 96 - 3000 active elements (0 deactivated)
10:35:39 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.006300209 to -0.000713309) found for ocean_vertical_diffusivity, replacing with NaN
10:35:39 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:39 INFO    opendrift:2118: 2026-01-16 08:34:31.485059 - step 45 of 96 - 3000 active elements (0 deactivated)
10:35:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:40 INFO    opendrift:2118: 2026-01-16 09:04:31.485059 - step 46 of 96 - 3000 active elements (0 deactivated)
10:35:42 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.006163765 to -0.0008119836) found for ocean_vertical_diffusivity, replacing with NaN
10:35:42 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:42 INFO    opendrift:2118: 2026-01-16 09:34:31.485059 - step 47 of 96 - 3000 active elements (0 deactivated)
10:35:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:42 INFO    opendrift:2118: 2026-01-16 10:04:31.485059 - step 48 of 96 - 3000 active elements (0 deactivated)
10:35:44 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0056592785 to -0.000644047) found for ocean_vertical_diffusivity, replacing with NaN
10:35:44 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:44 INFO    opendrift:2118: 2026-01-16 10:34:31.485059 - step 49 of 96 - 3000 active elements (0 deactivated)
10:35:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:45 INFO    opendrift:2118: 2026-01-16 11:04:31.485059 - step 50 of 96 - 3000 active elements (0 deactivated)
10:35:46 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.004569642 to -0.0005027029) found for ocean_vertical_diffusivity, replacing with NaN
10:35:46 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:46 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:46 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:47 INFO    opendrift:2118: 2026-01-16 11:34:31.485059 - step 51 of 96 - 3000 active elements (0 deactivated)
10:35:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:47 INFO    opendrift:2118: 2026-01-16 12:04:31.485059 - step 52 of 96 - 3000 active elements (0 deactivated)
10:35:49 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0032275752 to -0.0003251524) found for ocean_vertical_diffusivity, replacing with NaN
10:35:49 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:49 INFO    opendrift:2118: 2026-01-16 12:34:31.485059 - step 53 of 96 - 3000 active elements (0 deactivated)
10:35:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:49 INFO    opendrift:2118: 2026-01-16 13:04:31.485059 - step 54 of 96 - 3000 active elements (0 deactivated)
10:35:51 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0028038435 to -0.0003599533) found for ocean_vertical_diffusivity, replacing with NaN
10:35:51 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:51 INFO    opendrift:2118: 2026-01-16 13:34:31.485059 - step 55 of 96 - 3000 active elements (0 deactivated)
10:35:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:51 INFO    opendrift:2118: 2026-01-16 14:04:31.485059 - step 56 of 96 - 3000 active elements (0 deactivated)
10:35:53 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0037742488 to -0.00026288972) found for ocean_vertical_diffusivity, replacing with NaN
10:35:53 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:54 INFO    opendrift:2118: 2026-01-16 14:34:31.485059 - step 57 of 96 - 3000 active elements (0 deactivated)
10:35:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:54 INFO    opendrift:2118: 2026-01-16 15:04:31.485059 - step 58 of 96 - 3000 active elements (0 deactivated)
10:35:56 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.005928685 to -0.0007636284) found for ocean_vertical_diffusivity, replacing with NaN
10:35:56 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:56 INFO    opendrift:2118: 2026-01-16 15:34:31.485059 - step 59 of 96 - 3000 active elements (0 deactivated)
10:35:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:56 INFO    opendrift:2118: 2026-01-16 16:04:31.485059 - step 60 of 96 - 3000 active elements (0 deactivated)
10:35:58 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.007581459 to -0.0010445701) found for ocean_vertical_diffusivity, replacing with NaN
10:35:58 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:35:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:58 INFO    opendrift:2118: 2026-01-16 16:34:31.485059 - step 61 of 96 - 3000 active elements (0 deactivated)
10:35:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:35:58 INFO    opendrift:2118: 2026-01-16 17:04:31.485059 - step 62 of 96 - 3000 active elements (0 deactivated)
10:36:01 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.00801247 to -0.0010620297) found for ocean_vertical_diffusivity, replacing with NaN
10:36:01 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:01 INFO    opendrift:2118: 2026-01-16 17:34:31.485059 - step 63 of 96 - 3000 active elements (0 deactivated)
10:36:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:02 INFO    opendrift:2118: 2026-01-16 18:04:31.485059 - step 64 of 96 - 3000 active elements (0 deactivated)
10:36:04 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0073418277 to -0.0010712218) found for ocean_vertical_diffusivity, replacing with NaN
10:36:04 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:04 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:04 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:06 INFO    opendrift:2118: 2026-01-16 18:34:31.485059 - step 65 of 96 - 2999 active elements (1 deactivated)
10:36:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:06 INFO    opendrift:2118: 2026-01-16 19:04:31.485059 - step 66 of 96 - 2999 active elements (1 deactivated)
10:36:08 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.006410585 to -0.0011650668) found for ocean_vertical_diffusivity, replacing with NaN
10:36:08 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:08 INFO    opendrift:2118: 2026-01-16 19:34:31.485059 - step 67 of 96 - 2998 active elements (2 deactivated)
10:36:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:09 INFO    opendrift:2118: 2026-01-16 20:04:31.485059 - step 68 of 96 - 2998 active elements (2 deactivated)
10:36:11 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.005330579 to -0.0010146017) found for ocean_vertical_diffusivity, replacing with NaN
10:36:11 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:11 INFO    opendrift:2118: 2026-01-16 20:34:31.485059 - step 69 of 96 - 2997 active elements (3 deactivated)
10:36:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:11 INFO    opendrift:2118: 2026-01-16 21:04:31.485059 - step 70 of 96 - 2997 active elements (3 deactivated)
10:36:13 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0047869203 to -0.00094596273) found for ocean_vertical_diffusivity, replacing with NaN
10:36:13 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:13 INFO    opendrift:2118: 2026-01-16 21:34:31.485059 - step 71 of 96 - 2995 active elements (5 deactivated)
10:36:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:13 INFO    opendrift:2118: 2026-01-16 22:04:31.485059 - step 72 of 96 - 2995 active elements (5 deactivated)
10:36:15 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0044553755 to -0.0008255452) found for ocean_vertical_diffusivity, replacing with NaN
10:36:15 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:15 INFO    opendrift:2118: 2026-01-16 22:34:31.485059 - step 73 of 96 - 2995 active elements (5 deactivated)
10:36:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:15 INFO    opendrift:2118: 2026-01-16 23:04:31.485059 - step 74 of 96 - 2994 active elements (6 deactivated)
10:36:17 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0034397705 to -0.00045716547) found for ocean_vertical_diffusivity, replacing with NaN
10:36:17 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:17 INFO    opendrift:2118: 2026-01-16 23:34:31.485059 - step 75 of 96 - 2994 active elements (6 deactivated)
10:36:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:17 INFO    opendrift:2118: 2026-01-17 00:04:31.485059 - step 76 of 96 - 2993 active elements (7 deactivated)
10:36:19 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0058933976 to -0.0004851404) found for ocean_vertical_diffusivity, replacing with NaN
10:36:19 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:20 INFO    opendrift:2118: 2026-01-17 00:34:31.485059 - step 77 of 96 - 2992 active elements (8 deactivated)
10:36:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:20 INFO    opendrift:2118: 2026-01-17 01:04:31.485059 - step 78 of 96 - 2991 active elements (9 deactivated)
10:36:22 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0052557033 to -0.00035678656) found for ocean_vertical_diffusivity, replacing with NaN
10:36:22 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:22 INFO    opendrift:2118: 2026-01-17 01:34:31.485059 - step 79 of 96 - 2989 active elements (11 deactivated)
10:36:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:22 INFO    opendrift:2118: 2026-01-17 02:04:31.485059 - step 80 of 96 - 2987 active elements (13 deactivated)
10:36:24 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.005823848 to -0.000521282) found for ocean_vertical_diffusivity, replacing with NaN
10:36:24 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:24 INFO    opendrift:2118: 2026-01-17 02:34:31.485059 - step 81 of 96 - 2983 active elements (17 deactivated)
10:36:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:25 INFO    opendrift:2118: 2026-01-17 03:04:31.485059 - step 82 of 96 - 2981 active elements (19 deactivated)
10:36:27 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.008607096 to -0.0005827284) found for ocean_vertical_diffusivity, replacing with NaN
10:36:27 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:27 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:27 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:27 INFO    opendrift:2118: 2026-01-17 03:34:31.485059 - step 83 of 96 - 2977 active elements (23 deactivated)
10:36:27 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:27 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:27 INFO    opendrift:2118: 2026-01-17 04:04:31.485059 - step 84 of 96 - 2974 active elements (26 deactivated)
10:36:29 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.012093728 to -0.0013803246) found for ocean_vertical_diffusivity, replacing with NaN
10:36:29 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:29 INFO    opendrift:2118: 2026-01-17 04:34:31.485059 - step 85 of 96 - 2963 active elements (37 deactivated)
10:36:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:29 INFO    opendrift:2118: 2026-01-17 05:04:31.485059 - step 86 of 96 - 2951 active elements (49 deactivated)
10:36:32 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.014822403 to -0.0012083176) found for ocean_vertical_diffusivity, replacing with NaN
10:36:32 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:32 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:32 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:32 INFO    opendrift:2118: 2026-01-17 05:34:31.485059 - step 87 of 96 - 2928 active elements (72 deactivated)
10:36:32 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:32 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:32 INFO    opendrift:2118: 2026-01-17 06:04:31.485059 - step 88 of 96 - 2908 active elements (92 deactivated)
10:36:34 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0155025255 to -0.001618828) found for ocean_vertical_diffusivity, replacing with NaN
10:36:34 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:34 INFO    opendrift:2118: 2026-01-17 06:34:31.485059 - step 89 of 96 - 2874 active elements (126 deactivated)
10:36:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:34 INFO    opendrift:2118: 2026-01-17 07:04:31.485059 - step 90 of 96 - 2845 active elements (155 deactivated)
10:36:38 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.015079487 to -0.0014143856) found for ocean_vertical_diffusivity, replacing with NaN
10:36:38 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:38 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:38 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:38 INFO    opendrift:2118: 2026-01-17 07:34:31.485059 - step 91 of 96 - 2821 active elements (179 deactivated)
10:36:38 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:38 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:38 INFO    opendrift:2118: 2026-01-17 08:04:31.485059 - step 92 of 96 - 2803 active elements (197 deactivated)
10:36:40 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.013186717 to -0.00084414514) found for ocean_vertical_diffusivity, replacing with NaN
10:36:40 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:40 INFO    opendrift:2118: 2026-01-17 08:34:31.485059 - step 93 of 96 - 2798 active elements (202 deactivated)
10:36:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:41 INFO    opendrift:2118: 2026-01-17 09:04:31.485059 - step 94 of 96 - 2791 active elements (209 deactivated)
10:36:43 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.013383632 to -0.00079648214) found for ocean_vertical_diffusivity, replacing with NaN
10:36:43 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:43 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:43 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:43 INFO    opendrift:2118: 2026-01-17 09:34:31.485059 - step 95 of 96 - 2784 active elements (216 deactivated)
10:36:43 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:43 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:43 INFO    opendrift:2118: 2026-01-17 10:04:31.485059 - step 96 of 96 - 2778 active elements (222 deactivated)
10:36:45 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.012178458 to -0.0012396388) found for ocean_vertical_diffusivity, replacing with NaN
10:36:45 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
10:36:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
10:36:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning

<xarray.Dataset> Size: 25MB
Dimensions:                                                                              (
                                                                                          trajectory: 3000,
                                                                                          time: 49)
Coordinates:
  * trajectory                                                                           (trajectory) int64 24kB ...
  * time                                                                                 (time) datetime64[ns] 392B ...
Data variables: (12/43)
    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_temperature                                                                (trajectory, time) float32 588kB ...
    sea_water_salinity                                                                   (trajectory, time) float32 588kB ...
    sea_floor_depth_below_sea_level                                                      (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/165)
    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-01-15T10:34:31.485059 ... 2...

    standard_name :

    time

    long_name :

    time

    dtype :

    <class 'numpy.float64'>

    array(['2026-01-15T10:34:31.485059000', '2026-01-15T11:34:31.485059000',
           '2026-01-15T12:34:31.485059000', '2026-01-15T13:34:31.485059000',
           '2026-01-15T14:34:31.485059000', '2026-01-15T15:34:31.485059000',
           '2026-01-15T16:34:31.485059000', '2026-01-15T17:34:31.485059000',
           '2026-01-15T18:34:31.485059000', '2026-01-15T19:34:31.485059000',
           '2026-01-15T20:34:31.485059000', '2026-01-15T21:34:31.485059000',
           '2026-01-15T22:34:31.485059000', '2026-01-15T23:34:31.485059000',
           '2026-01-16T00:34:31.485059000', '2026-01-16T01:34:31.485059000',
           '2026-01-16T02:34:31.485059000', '2026-01-16T03:34:31.485059000',
           '2026-01-16T04:34:31.485059000', '2026-01-16T05:34:31.485059000',
           '2026-01-16T06:34:31.485059000', '2026-01-16T07:34:31.485059000',
           '2026-01-16T08:34:31.485059000', '2026-01-16T09:34:31.485059000',
           '2026-01-16T10:34:31.485059000', '2026-01-16T11:34:31.485059000',
           '2026-01-16T12:34:31.485059000', '2026-01-16T13:34:31.485059000',
           '2026-01-16T14:34:31.485059000', '2026-01-16T15:34:31.485059000',
           '2026-01-16T16:34:31.485059000', '2026-01-16T17:34:31.485059000',
           '2026-01-16T18:34:31.485059000', '2026-01-16T19:34:31.485059000',
           '2026-01-16T20:34:31.485059000', '2026-01-16T21:34:31.485059000',
           '2026-01-16T22:34:31.485059000', '2026-01-16T23:34:31.485059000',
           '2026-01-17T00:34:31.485059000', '2026-01-17T01:34:31.485059000',
           '2026-01-17T02:34:31.485059000', '2026-01-17T03:34:31.485059000',
           '2026-01-17T04:34:31.485059000', '2026-01-17T05:34:31.485059000',
           '2026-01-17T06:34:31.485059000', '2026-01-17T07:34:31.485059000',
           '2026-01-17T08:34:31.485059000', '2026-01-17T09:34:31.485059000',
           '2026-01-17T10:34:31.485059000'], dtype='datetime64[ns]')

• Data variables: (43)
  • 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.22 19.21 ... 18.85 18.85

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    axis :

    X

    minval :

    17.817108

    maxval :

    19.296974

    array([[19.195274, 19.215483, 19.20522 , ..., 18.45055 , 18.413813,
            18.363558],
           [19.191893, 19.173586, 19.16449 , ..., 18.683363, 18.679735,
            18.670153],
           [19.193327, 19.19992 , 19.187004, ..., 18.600574, 18.597836,
            18.579905],
           ...,
           [19.190422, 19.213375, 19.207834, ...,       nan,       nan,
                  nan],
           [19.187893, 19.165318, 19.151188, ..., 18.104942, 18.05368 ,
            17.994984],
           [19.1907  , 19.200506, 19.222773, ..., 18.865133, 18.850815,
            18.854134]], shape=(3000, 49), dtype=float32)

  • lat
    (trajectory, time)
    float32
    79.6 79.6 79.6 ... 79.68 79.69

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    axis :

    Y

    minval :

    79.58218

    maxval :

    79.875984

    array([[79.59931 , 79.60179 , 79.59995 , ..., 79.75329 , 79.75846 ,
            79.76617 ],
           [79.598854, 79.59973 , 79.60118 , ..., 79.72176 , 79.7263  ,
            79.73102 ],
           [79.59855 , 79.599724, 79.6014  , ..., 79.7198  , 79.7252  ,
            79.72965 ],
           ...,
           [79.59836 , 79.59816 , 79.5977  , ...,       nan,       nan,
                  nan],
           [79.59878 , 79.5978  , 79.60296 , ..., 79.76043 , 79.77339 ,
            79.784676],
           [79.59873 , 79.598976, 79.6026  , ..., 79.68001 , 79.68156 ,
            79.68513 ]], 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 :

    -36.195084

    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 :

    %

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

    dtype :

    <class 'numpy.float32'>

    units :

    m/s

    description :

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

    minval :

    0.0

    maxval :

    0.047993936

    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_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.1643 ... 0.1652 0.1656

    dtype :

    <class 'numpy.float32'>

    units :

    m2/s

    description :

    Kinematic viscosity of oil (emulsion)

    minval :

    0.021692334

    maxval :

    0.16865978

    array([[0.02169233, 0.16429181, 0.16412897, ..., 0.1629603 , 0.16299309,
            0.1630899 ],
           [0.02169233, 0.16410461, 0.16397496, ..., 0.16426033, 0.16428944,
            0.16447215],
           [0.02169233, 0.16423514, 0.1640881 , ..., 0.16335255, 0.16329779,
            0.16329779],
           ...,
           [0.02169233, 0.16422807, 0.16422099, ...,        nan,        nan,
                   nan],
           [0.02169233, 0.16410303, 0.1639671 , ..., 0.15984851, 0.1571463 ,
            0.15432593],
           [0.02169233, 0.16422257, 0.16420525, ..., 0.16526048, 0.16523911,
            0.16558312]], shape=(3000, 49), dtype=float32)

  • density
    (trajectory, time)
    float32
    988.1 997.1 997.1 ... 997.2 997.2

    dtype :

    <class 'numpy.float32'>

    units :

    kg/m^3

    minval :

    988.1

    maxval :

    997.2391

    array([[988.1    , 997.1277 , 997.1235 , ..., 997.09314, 997.094  ,
            997.0965 ],
           [988.1    , 997.12286, 997.1195 , ..., 997.1269 , 997.1276 ,
            997.1323 ],
           [988.1    , 997.1262 , 997.12244, ..., 997.10333, 997.1019 ,
            997.1019 ],
           ...,
           [988.1    , 997.12604, 997.12585, ...,       nan,       nan,
                  nan],
           [988.1    , 997.1228 , 997.11926, ..., 997.01117, 996.93866,
            996.8615 ],
           [988.1    , 997.1259 , 997.1254 , ..., 997.15265, 997.1521 ,
            997.16095]], 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 :

    %

    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 :

    %

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

    dtype :

    <class 'numpy.float32'>

    units :

    m

    minval :

    0.0

    maxval :

    0.0029900854

    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)

  • x_sea_water_velocity
    (trajectory, time)
    float32
    0.1842 0.09822 ... 0.006773 -0.1454

    minval :

    -0.30299678

    maxval :

    0.35651928

    array([[ 0.18417373,  0.09822447,  0.05832989, ..., -0.01310237,
            -0.11932994, -0.01332879],
           [ 0.03641567,  0.05884644, -0.00563928, ...,  0.00587553,
            -0.04166141,  0.06716818],
           [ 0.08822528, -0.03482154, -0.02732509, ..., -0.01114503,
            -0.01293597,  0.00331432],
           ...,
           [ 0.04205337,  0.02579064,  0.0211725 , ...,         nan,
                    nan,         nan],
           [ 0.02848651,  0.08037608,  0.05509579, ..., -0.12797241,
            -0.04618533, -0.08709075],
           [ 0.11669202,  0.16743945,  0.04705049, ..., -0.05152597,
             0.00677292, -0.14541174]], shape=(3000, 49), dtype=float32)

  • y_sea_water_velocity
    (trajectory, time)
    float32
    0.01112 -0.1608 ... -0.02521

    minval :

    -0.31306228

    maxval :

    0.46632648

    array([[ 0.01112498, -0.16084616, -0.11601207, ...,  0.00396374,
             0.01710097,  0.05635867],
           [-0.03927519, -0.05471944, -0.05408569, ...,  0.07021736,
             0.01645594,  0.00635464],
           [-0.06559813, -0.06612746, -0.10998987, ...,  0.11091574,
            -0.05712   ,  0.06244832],
           ...,
           [-0.1626923 , -0.11050987, -0.18946414, ...,         nan,
                    nan,         nan],
           [-0.06412833, -0.0491679 , -0.14968336, ...,  0.28378925,
             0.30354145,  0.14009605],
           [-0.05078913, -0.04332309, -0.08730384, ...,  0.01642409,
             0.0283372 , -0.0252054 ]], shape=(3000, 49), dtype=float32)

  • x_wind
    (trajectory, time)
    float32
    -1.346 -3.488 ... -4.261 -5.147

    minval :

    -8.065187

    maxval :

    0.47352323

    array([[-1.3460633, -3.4883156, -2.482778 , ..., -5.358876 , -6.8437576,
            -6.634155 ],
           [-1.5886765, -3.3568559, -2.1822171, ..., -4.31921  , -4.056411 ,
            -3.9749339],
           [-2.3908858, -1.965581 , -2.852019 , ..., -5.6090674, -5.1697483,
            -5.466607 ],
           ...,
           [-1.4125258, -3.0546725, -1.9329821, ...,        nan,        nan,
                   nan],
           [-2.0729935, -2.1778252, -3.0710924, ..., -5.5486484, -5.653829 ,
            -6.08087  ],
           [-1.8931137, -3.9610903, -2.743093 , ..., -4.511143 , -4.260729 ,
            -5.147438 ]], shape=(3000, 49), dtype=float32)

  • y_wind
    (trajectory, time)
    float32
    4.734 5.302 4.269 ... 3.483 3.456

    minval :

    -1.7923073

    maxval :

    9.228616

    array([[4.733804 , 5.3019047, 4.2688684, ..., 2.8143551, 3.3126333,
            4.00304  ],
           [3.7792525, 5.296665 , 4.4049487, ..., 2.6777213, 2.5615273,
            2.3969426],
           [3.2968335, 4.3909607, 4.4596567, ..., 3.95062  , 4.0414543,
            2.9437196],
           ...,
           [3.4759586, 5.308793 , 5.9775257, ...,       nan,       nan,
                  nan],
           [3.3931162, 5.043725 , 4.826462 , ..., 5.115598 , 5.2001667,
            5.0795913],
           [3.5087416, 4.6184206, 4.9497423, ..., 3.0590727, 3.4825113,
            3.4557524]], shape=(3000, 49), dtype=float32)

  • sea_surface_height
    (trajectory, time)
    float32
    -0.5842 -0.5245 ... -0.6586 -0.6409

    minval :

    -0.91321945

    maxval :

    -0.25864506

    array([[-0.5842121 , -0.52452654, -0.4602727 , ..., -0.7128321 ,
            -0.64891696, -0.63235563],
           [-0.58420265, -0.52391386, -0.4587706 , ..., -0.71872437,
            -0.6526317 , -0.6359518 ],
           [-0.5842197 , -0.52434075, -0.45955685, ..., -0.71862227,
            -0.6521623 , -0.6352624 ],
           ...,
           [-0.58420634, -0.5246022 , -0.46047613, ...,         nan,
                    nan,         nan],
           [-0.5841788 , -0.5237719 , -0.45826536, ..., -0.7206361 ,
            -0.65635514, -0.6382364 ],
           [-0.5841981 , -0.5243726 , -0.4608162 , ..., -0.7292706 ,
            -0.6586011 , -0.6408914 ]], shape=(3000, 49), dtype=float32)

  • sea_surface_wave_significant_height
    (trajectory, time)
    float32
    0.5958 0.9909 ... 0.7449 0.9456

    minval :

    0.00044649135

    maxval :

    2.9287808

    array([[0.5958314 , 0.990852  , 0.59993064, ..., 0.90129846, 1.4221396 ,
            1.4768939 ],
           [0.41344336, 0.9673493 , 0.5944748 , ..., 0.63531387, 0.56619084,
            0.5300177 ],
           [0.40800193, 0.5693435 , 0.6893547 , ..., 1.1578982 , 1.0592673 ,
            0.9483122 ],
           ...,
           [0.34630692, 0.922852  , 0.97089386, ...,        nan,        nan,
                   nan],
           [0.3889393 , 0.74247915, 0.80506796, ..., 1.4011382 , 1.4515848 ,
            1.5443691 ],
           [0.39102063, 0.9106931 , 0.7878028 , ..., 0.7308251 , 0.7449297 ,
            0.9455832 ]], 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
    4.205 5.423 4.22 ... 4.702 5.298

    minval :

    0.115116104

    maxval :

    9.323372

    array([[4.2052445, 5.42293  , 4.2196856, ..., 5.172064 , 6.4968176,
            6.620705 ],
           [3.5029795, 5.358228 , 4.2004547, ..., 4.342339 , 4.0993123,
            3.9662013],
           [3.4798512, 4.110709 , 4.523254 , ..., 5.8622556, 5.6070232,
            5.305243 ],
           ...,
           [3.2059753, 5.233541 , 5.3680367, ...,       nan,       nan,
                  nan],
           [3.397586 , 4.6943097, 4.888165 , ..., 6.448669 , 6.563731 ,
            6.770256 ],
           [3.4066648, 5.1989493, 4.8354664, ..., 4.657323 , 4.70205  ,
            5.2976036]], shape=(3000, 49), dtype=float32)

  • sea_ice_area_fraction
    (trajectory, time)
    float32
    0.0 0.0 0.0 ... 5.184e-06 0.0

    minval :

    0.0

    maxval :

    0.1115536

    array([[0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ..., 1.6243074e-03,
            2.7442509e-03, 2.7136356e-03],
           [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ..., 3.4316283e-05,
            4.4012417e-05, 4.6166075e-05],
           [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ..., 2.4445323e-04,
            1.2303040e-04, 9.7358192e-05],
           ...,
           [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ...,           nan,
                      nan,           nan],
           [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ..., 6.3851185e-02,
            7.7982828e-02, 7.2946236e-02],
           [0.0000000e+00, 0.0000000e+00, 0.0000000e+00, ..., 3.2211469e-06,
            5.1841325e-06, 0.0000000e+00]], 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.2208789

    maxval :

    0.014384421

    array([[ 0.        ,  0.        ,  0.        , ..., -0.09367246,
            -0.10267102, -0.10480528],
           [ 0.        ,  0.        ,  0.        , ...,  0.        ,
             0.        ,  0.        ],
           [ 0.        ,  0.        ,  0.        , ..., -0.00077722,
            -0.0009054 , -0.00061051],
           ...,
           [ 0.        ,  0.        ,  0.        , ...,         nan,
                    nan,         nan],
           [ 0.        ,  0.        ,  0.        , ..., -0.13532062,
            -0.15044807, -0.16206488],
           [ 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.018850487

    maxval :

    0.36290377

    array([[0.        , 0.        , 0.        , ..., 0.10054162, 0.10863838,
            0.11004911],
           [0.        , 0.        , 0.        , ..., 0.        , 0.        ,
            0.        ],
           [0.        , 0.        , 0.        , ..., 0.0015798 , 0.00144787,
            0.00088584],
           ...,
           [0.        , 0.        , 0.        , ...,        nan,        nan,
                   nan],
           [0.        , 0.        , 0.        , ..., 0.2919153 , 0.2777135 ,
            0.2712284 ],
           [0.        , 0.        , 0.        , ..., 0.        , 0.        ,
            0.        ]], shape=(3000, 49), dtype=float32)

  • sea_water_temperature
    (trajectory, time)
    float32
    -1.69 -1.689 ... -1.732 271.4

    minval :

    -1.8546681

    maxval :

    272.40015

    array([[ -1.6904825,  -1.6890671,  -1.6855277, ...,  -1.6419146,
             -1.6444199, 271.50525  ],
           [ -1.6900493,  -1.6818881,  -1.6785452, ...,  -1.6889262,
             -1.6938149, 271.45148  ],
           [ -1.6905925,  -1.6867524,  -1.6815501, ...,  -1.6529294,
             -1.6544192, 271.49713  ],
           ...,
           [ -1.6900767,  -1.6908206,  -1.687326 , ...,         nan,
                    nan,         nan],
           [ -1.689203 ,  -1.6812527,  -1.6764303, ...,  -1.4602901,
             -1.3440121, 271.8577   ],
           [ -1.689869 ,  -1.6873441,  -1.6878756, ...,  -1.7251881,
             -1.7324876, 271.4086   ]], shape=(3000, 49), dtype=float32)

  • sea_water_salinity
    (trajectory, time)
    float32
    33.94 33.94 33.94 ... 33.92 33.92

    minval :

    33.907787

    maxval :

    34.102524

    array([[33.936398, 33.935234, 33.937973, ..., 33.947918, 33.94826 ,
            33.948742],
           [33.936905, 33.94126 , 33.943115, ..., 33.931156, 33.929817,
            33.928684],
           [33.93661 , 33.937557, 33.940754, ..., 33.93984 , 33.939243,
            33.939674],
           ...,
           [33.937077, 33.934948, 33.937107, ...,       nan,       nan,
                  nan],
           [33.93757 , 33.94238 , 33.944386, ..., 33.973663, 33.99073 ,
            33.998726],
           [33.937088, 33.93733 , 33.935555, ..., 33.92342 , 33.922478,
            33.92023 ]], shape=(3000, 49), dtype=float32)

  • sea_floor_depth_below_sea_level
    (trajectory, time)
    float32
    170.4 167.3 168.5 ... 180.8 173.7

    minval :

    51.28601

    maxval :

    298.1138

    array([[170.4251 , 167.30254, 168.46524, ..., 201.8335 , 215.18501,
            224.06992],
           [170.93388, 175.85362, 179.35652, ..., 128.9069 , 121.77753,
            117.02251],
           [170.35571, 169.5801 , 174.01733, ..., 177.11356, 168.67076,
            168.38046],
           ...,
           [170.93738, 165.22348, 166.1994 , ...,       nan,       nan,
                  nan],
           [171.87785, 174.06766, 183.96413, ..., 263.32608, 252.23602,
            253.23837],
           [171.14177, 168.90584, 166.16011, ..., 181.30623, 180.77452,
            173.70924]], shape=(3000, 49), dtype=float32)

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

    minval :

    0.006965164

    maxval :

    0.02

    array([[0.02, 0.02, 0.02, ..., 0.02, 0.02, 0.02],
           [0.02, 0.02, 0.02, ..., 0.02, 0.02, 0.02],
           [0.02, 0.02, 0.02, ..., 0.02, 0.02, 0.02],
           ...,
           [0.02, 0.02, 0.02, ...,  nan,  nan,  nan],
           [0.02, 0.02, 0.02, ..., 0.02, 0.02, 0.02],
           [0.02, 0.02, 0.02, ..., 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)

• Indexes: (2)
  • trajectory
    PandasIndex

    PandasIndex(Index([   0,    1,    2,    3,    4,    5,    6,    7,    8,    9,
           ...
           2990, 2991, 2992, 2993, 2994, 2995, 2996, 2997, 2998, 2999],
          dtype='int64', name='trajectory', length=3000))

  • time
    PandasIndex

    PandasIndex(DatetimeIndex(['2026-01-15 10:34:31.485059', '2026-01-15 11:34:31.485059',
                   '2026-01-15 12:34:31.485059', '2026-01-15 13:34:31.485059',
                   '2026-01-15 14:34:31.485059', '2026-01-15 15:34:31.485059',
                   '2026-01-15 16:34:31.485059', '2026-01-15 17:34:31.485059',
                   '2026-01-15 18:34:31.485059', '2026-01-15 19:34:31.485059',
                   '2026-01-15 20:34:31.485059', '2026-01-15 21:34:31.485059',
                   '2026-01-15 22:34:31.485059', '2026-01-15 23:34:31.485059',
                   '2026-01-16 00:34:31.485059', '2026-01-16 01:34:31.485059',
                   '2026-01-16 02:34:31.485059', '2026-01-16 03:34:31.485059',
                   '2026-01-16 04:34:31.485059', '2026-01-16 05:34:31.485059',
                   '2026-01-16 06:34:31.485059', '2026-01-16 07:34:31.485059',
                   '2026-01-16 08:34:31.485059', '2026-01-16 09:34:31.485059',
                   '2026-01-16 10:34:31.485059', '2026-01-16 11:34:31.485059',
                   '2026-01-16 12:34:31.485059', '2026-01-16 13:34:31.485059',
                   '2026-01-16 14:34:31.485059', '2026-01-16 15:34:31.485059',
                   '2026-01-16 16:34:31.485059', '2026-01-16 17:34:31.485059',
                   '2026-01-16 18:34:31.485059', '2026-01-16 19:34:31.485059',
                   '2026-01-16 20:34:31.485059', '2026-01-16 21:34:31.485059',
                   '2026-01-16 22:34:31.485059', '2026-01-16 23:34:31.485059',
                   '2026-01-17 00:34:31.485059', '2026-01-17 01:34:31.485059',
                   '2026-01-17 02:34:31.485059', '2026-01-17 03:34:31.485059',
                   '2026-01-17 04:34:31.485059', '2026-01-17 05:34:31.485059',
                   '2026-01-17 06:34:31.485059', '2026-01-17 07:34:31.485059',
                   '2026-01-17 08:34:31.485059', '2026-01-17 09:34:31.485059',
                   '2026-01-17 10:34:31.485059'],
                  dtype='datetime64[ns]', name='time', freq='h'))

• Attributes: (165)

  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-01-22 10:34:35.765053

  date_created :

  2026-01-22T10:34:35.765073

  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', 'global_landmask'])

  time_coverage_start :

  2026-01-15 10:34:31.485059

  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: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:horizontal_diffusivity :

  10.0

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

  seed_lon :

  19.1909

  seed_lat :

  79.5986

  seed_radius :

  50.0

  seed_time :

  2026-01-15 10:34:31.485059

  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-01-22 10:34:35.778801

  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_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-01-17 10:34:31.485059

  time_coverage_duration :

  P2DT0H0M0S

  time_coverage_resolution :

  P0DT0H30M0S

  performance :

  -------------------- Reader performance: -------------------- global_landmask 0:00:01.5 total 0:00:00.0 preparing 0:00:01.5 reading 0:00:00.0 masking -------------------- https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be 0:01:49.1 total 0:00:00.0 preparing 0:01:47.8 reading 0:00:01.9 interpolation 0:00:00.1 interpolation_time 0:00:01.2 rotating vectors 0:00:00.0 masking -------------------- Performance: 2:14.1 total time 4.2 configuration 0.0 preparing main loop 0.0 moving elements to ocean 2:09.7 main loop 10.6 updating elements 0.0 oil weathering 0.0 updating viscosities 0.0 updating densities 9.9 vertical mixing 0.0 cleaning up --------------------

  geospatial_bounds_crs :

  EPSG:4326

  geospatial_bounds_vertical_crs :

  EPSG:5831

  geospatial_lat_min :

  79.58218

  geospatial_lat_max :

  79.875984

  geospatial_lat_units :

  degrees_north

  geospatial_lat_resolution :

  point

  geospatial_lon_min :

  17.817108

  geospatial_lon_max :

  19.296974

  geospatial_lon_units :

  degrees_east

  geospatial_lon_resolution :

  point

  runtime :

  0:02:14.107818

  geospatial_vertical_min :

  -36.195084

  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:
--------------------
global_landmask
 0:00:01.5  total
 0:00:00.0  preparing
 0:00:01.5  reading
 0:00:00.0  masking
--------------------
https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
 0:01:49.1  total
 0:00:00.0  preparing
 0:01:47.8  reading
 0:00:01.9  interpolation
 0:00:00.1  interpolation_time
 0:00:01.2  rotating vectors
 0:00:00.0  masking
--------------------
Performance:
 2:14.1 total time
    4.2 configuration
    0.0 preparing main loop
      0.0 moving elements to ocean
 2:09.7 main loop
     10.6 updating elements
        0.0 oil weathering
          0.0 updating viscosities
          0.0 updating densities
        9.9 vertical mixing
    0.0 cleaning up
--------------------
===========================
Model:  OpenOil     (OpenDrift version 1.14.8)
        2765 active Oil particles  (235 deactivated, 0 scheduled)
-------------------
Environment variables:
  -----
  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-01-15 10:34:31.485059 UTC
        Present: 2026-01-17 10:34:31.485059 UTC
        Calculation steps: 96 * 0:30:00 - total time: 2 days, 0:00:00
        Output steps: 49 * 1:00:00
===========================

10:36:56 INFO    opendrift:4667: Saving animation to /root/project/docs/source/gallery/animations/example_oil_ice_0.gif...
10:37:36 INFO    opendrift:3101: Time to make animation: 0:00:51.205247

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-01-15 10:34 to 2026-01-17 10:34 UTC (49 steps)'}>, <Figure size 629.473x1100 with 2 Axes>)