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)

08:24:09 INFO    opendrift:568: OpenDriftSimulation initialised (version 1.14.9 / v1.14.9-9-gcde0dcf)

Using live data from Barents 2.5 km ocean model

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

Adjusting some configuration

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

Imaginary oil spill in Hinlopen strait

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

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

08:24:13 INFO    opendrift:1828: Skipping environment variable upward_sea_water_velocity because of condition ['drift:vertical_advection', 'is', False]
08:24:13 INFO    opendrift:1839: Storing previous values of element property lon because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
08:24:13 INFO    opendrift:1839: Storing previous values of element property lat because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
08:24:13 INFO    opendrift:947: Using existing reader for land_binary_mask to move elements to ocean
08:24:13 INFO    opendrift:978: All points are in ocean
08:24:13 INFO    opendrift.models.openoil.openoil:696: Oil-water surface tension is 0.035935 Nm
08:24:13 INFO    opendrift.models.openoil.openoil:709: Max water fraction not available for GENERIC BUNKER C, using default
08:24:13 INFO    opendrift:2136: 2026-03-18 08:24:09.382179 - step 1 of 96 - 3000 active elements (0 deactivated)
08:24:13 INFO    opendrift.readers:64: Opening file with xr.open_dataset
08:24:15 INFO    opendrift.readers.reader_netCDF_CF_generic:340: Detected dimensions: {'x': 'X', 'y': 'Y', 'z': 'depth', 'time': 'time'}
08:24:15 INFO    opendrift.readers.basereader:178: Variable x_sea_water_velocity will be rotated from eastward_sea_water_velocity
08:24:15 INFO    opendrift.readers.basereader:178: Variable y_sea_water_velocity will be rotated from northward_sea_water_velocity
08:24:18 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.020515554 to -0.0035420728) found for ocean_vertical_diffusivity, replacing with NaN
08:24:18 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:20 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.016343934 to -0.003075265) found for ocean_vertical_diffusivity, replacing with NaN
08:24:20 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:21 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:21 INFO    opendrift:2136: 2026-03-18 08:54:09.382179 - step 2 of 96 - 3000 active elements (0 deactivated)
08:24:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:21 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:21 INFO    opendrift:2136: 2026-03-18 09:24:09.382179 - step 3 of 96 - 3000 active elements (0 deactivated)
08:24:23 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.014146863 to -0.0022660033) found for ocean_vertical_diffusivity, replacing with NaN
08:24:23 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:23 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:23 INFO    opendrift:2136: 2026-03-18 09:54:09.382179 - step 4 of 96 - 3000 active elements (0 deactivated)
08:24:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:23 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:23 INFO    opendrift:2136: 2026-03-18 10:24:09.382179 - step 5 of 96 - 3000 active elements (0 deactivated)
08:24:26 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.015302924 to -0.0022114657) found for ocean_vertical_diffusivity, replacing with NaN
08:24:26 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:26 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:26 INFO    opendrift:2136: 2026-03-18 10:54:09.382179 - step 6 of 96 - 3000 active elements (0 deactivated)
08:24:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:26 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:26 INFO    opendrift:2136: 2026-03-18 11:24:09.382179 - step 7 of 96 - 3000 active elements (0 deactivated)
08:24:28 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.01591355 to -0.002327347) found for ocean_vertical_diffusivity, replacing with NaN
08:24:28 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:28 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:28 INFO    opendrift:2136: 2026-03-18 11:54:09.382179 - step 8 of 96 - 3000 active elements (0 deactivated)
08:24:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:29 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:29 INFO    opendrift:2136: 2026-03-18 12:24:09.382179 - step 9 of 96 - 3000 active elements (0 deactivated)
08:24:31 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.019155629 to -0.0018842785) found for ocean_vertical_diffusivity, replacing with NaN
08:24:31 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:31 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:31 INFO    opendrift:2136: 2026-03-18 12:54:09.382179 - step 10 of 96 - 3000 active elements (0 deactivated)
08:24:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:31 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:31 INFO    opendrift:2136: 2026-03-18 13:24:09.382179 - step 11 of 96 - 3000 active elements (0 deactivated)
08:24:34 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.02196795 to -0.0033564668) found for ocean_vertical_diffusivity, replacing with NaN
08:24:34 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:34 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:34 INFO    opendrift:2136: 2026-03-18 13:54:09.382179 - step 12 of 96 - 3000 active elements (0 deactivated)
08:24:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:34 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:34 INFO    opendrift:2136: 2026-03-18 14:24:09.382179 - step 13 of 96 - 3000 active elements (0 deactivated)
08:24:36 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0045618583 to -3.6048554e-05) found for ocean_vertical_diffusivity, replacing with NaN
08:24:36 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:36 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:36 INFO    opendrift:2136: 2026-03-18 14:54:09.382179 - step 14 of 96 - 3000 active elements (0 deactivated)
08:24:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:37 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:37 INFO    opendrift:2136: 2026-03-18 15:24:09.382179 - step 15 of 96 - 3000 active elements (0 deactivated)
08:24:39 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0011200011 to -2.3114642e-05) found for ocean_vertical_diffusivity, replacing with NaN
08:24:39 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:39 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:39 INFO    opendrift:2136: 2026-03-18 15:54:09.382179 - step 16 of 96 - 3000 active elements (0 deactivated)
08:24:39 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:39 INFO    opendrift:2136: 2026-03-18 16:24:09.382179 - step 17 of 96 - 3000 active elements (0 deactivated)
08:24:42 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.008685812 to -0.001130904) found for ocean_vertical_diffusivity, replacing with NaN
08:24:42 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:42 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:42 INFO    opendrift:2136: 2026-03-18 16:54:09.382179 - step 18 of 96 - 3000 active elements (0 deactivated)
08:24:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:42 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:42 INFO    opendrift:2136: 2026-03-18 17:24:09.382179 - step 19 of 96 - 3000 active elements (0 deactivated)
08:24:44 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.03248814 to -0.003870071) found for ocean_vertical_diffusivity, replacing with NaN
08:24:44 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:44 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:44 INFO    opendrift:2136: 2026-03-18 17:54:09.382179 - step 20 of 96 - 3000 active elements (0 deactivated)
08:24:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:45 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:45 INFO    opendrift:2136: 2026-03-18 18:24:09.382179 - step 21 of 96 - 3000 active elements (0 deactivated)
08:24:47 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.048029788 to -0.0052189655) found for ocean_vertical_diffusivity, replacing with NaN
08:24:47 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:47 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:47 INFO    opendrift:2136: 2026-03-18 18:54:09.382179 - step 22 of 96 - 3000 active elements (0 deactivated)
08:24:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:47 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:47 INFO    opendrift:2136: 2026-03-18 19:24:09.382179 - step 23 of 96 - 3000 active elements (0 deactivated)
08:24:49 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.028574307 to -0.003476681) found for ocean_vertical_diffusivity, replacing with NaN
08:24:49 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:50 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:50 INFO    opendrift:2136: 2026-03-18 19:54:09.382179 - step 24 of 96 - 3000 active elements (0 deactivated)
08:24:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:50 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:50 INFO    opendrift:2136: 2026-03-18 20:24:09.382179 - step 25 of 96 - 3000 active elements (0 deactivated)
08:24:52 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.025184317 to -0.0016013231) found for ocean_vertical_diffusivity, replacing with NaN
08:24:52 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:52 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:52 INFO    opendrift:2136: 2026-03-18 20:54:09.382179 - step 26 of 96 - 3000 active elements (0 deactivated)
08:24:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:52 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:52 INFO    opendrift:2136: 2026-03-18 21:24:09.382179 - step 27 of 96 - 3000 active elements (0 deactivated)
08:24:55 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.019246317 to -0.0020095168) found for ocean_vertical_diffusivity, replacing with NaN
08:24:55 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:55 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:55 INFO    opendrift:2136: 2026-03-18 21:54:09.382179 - step 28 of 96 - 3000 active elements (0 deactivated)
08:24:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:55 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:55 INFO    opendrift:2136: 2026-03-18 22:24:09.382179 - step 29 of 96 - 3000 active elements (0 deactivated)
08:24:57 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0159128 to -0.0013982987) found for ocean_vertical_diffusivity, replacing with NaN
08:24:57 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:24:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:58 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:58 INFO    opendrift:2136: 2026-03-18 22:54:09.382179 - step 30 of 96 - 3000 active elements (0 deactivated)
08:24:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:24:58 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:24:58 INFO    opendrift:2136: 2026-03-18 23:24:09.382179 - step 31 of 96 - 3000 active elements (0 deactivated)
08:25:01 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.014308319 to -0.0012376952) found for ocean_vertical_diffusivity, replacing with NaN
08:25:01 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:01 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:01 INFO    opendrift:2136: 2026-03-18 23:54:09.382179 - step 32 of 96 - 3000 active elements (0 deactivated)
08:25:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:01 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:01 INFO    opendrift:2136: 2026-03-19 00:24:09.382179 - step 33 of 96 - 3000 active elements (0 deactivated)
08:25:03 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.014180235 to -0.0023934671) found for ocean_vertical_diffusivity, replacing with NaN
08:25:03 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:03 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:03 INFO    opendrift:2136: 2026-03-19 00:54:09.382179 - step 34 of 96 - 3000 active elements (0 deactivated)
08:25:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:03 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:03 INFO    opendrift:2136: 2026-03-19 01:24:09.382179 - step 35 of 96 - 3000 active elements (0 deactivated)
08:25:06 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.018003292 to -0.002411025) found for ocean_vertical_diffusivity, replacing with NaN
08:25:06 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:06 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:06 INFO    opendrift:2136: 2026-03-19 01:54:09.382179 - step 36 of 96 - 3000 active elements (0 deactivated)
08:25:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:06 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:06 INFO    opendrift:2136: 2026-03-19 02:24:09.382179 - step 37 of 96 - 3000 active elements (0 deactivated)
08:25:08 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.00334509 to -3.2893913e-06) found for ocean_vertical_diffusivity, replacing with NaN
08:25:08 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:08 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:08 INFO    opendrift:2136: 2026-03-19 02:54:09.382179 - step 38 of 96 - 3000 active elements (0 deactivated)
08:25:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:08 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:08 INFO    opendrift:2136: 2026-03-19 03:24:09.382179 - step 39 of 96 - 3000 active elements (0 deactivated)
08:25:12 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:12 INFO    opendrift:2136: 2026-03-19 03:54:09.382179 - step 40 of 96 - 3000 active elements (0 deactivated)
08:25:14 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:14 INFO    opendrift:2136: 2026-03-19 04:24:09.382179 - step 41 of 96 - 2999 active elements (1 deactivated)
08:25:16 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.005410226 to -5.701807e-06) found for ocean_vertical_diffusivity, replacing with NaN
08:25:16 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:17 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:17 INFO    opendrift:2136: 2026-03-19 04:54:09.382179 - step 42 of 96 - 2999 active elements (1 deactivated)
08:25:17 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:17 INFO    opendrift:2136: 2026-03-19 05:24:09.382179 - step 43 of 96 - 2999 active elements (1 deactivated)
08:25:19 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0133799715 to -0.0025298162) found for ocean_vertical_diffusivity, replacing with NaN
08:25:19 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:19 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:19 INFO    opendrift:2136: 2026-03-19 05:54:09.382179 - step 44 of 96 - 2999 active elements (1 deactivated)
08:25:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:19 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:19 INFO    opendrift:2136: 2026-03-19 06:24:09.382179 - step 45 of 96 - 2999 active elements (1 deactivated)
08:25:21 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.035921395 to -0.003031918) found for ocean_vertical_diffusivity, replacing with NaN
08:25:21 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:22 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:22 INFO    opendrift:2136: 2026-03-19 06:54:09.382179 - step 46 of 96 - 2999 active elements (1 deactivated)
08:25:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:22 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:22 INFO    opendrift:2136: 2026-03-19 07:24:09.382179 - step 47 of 96 - 2999 active elements (1 deactivated)
08:25:24 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.025441261 to -0.0027314208) found for ocean_vertical_diffusivity, replacing with NaN
08:25:24 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:24 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:24 INFO    opendrift:2136: 2026-03-19 07:54:09.382179 - step 48 of 96 - 2999 active elements (1 deactivated)
08:25:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:24 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:24 INFO    opendrift:2136: 2026-03-19 08:24:09.382179 - step 49 of 96 - 2999 active elements (1 deactivated)
08:25:27 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.019147234 to -0.0020653924) found for ocean_vertical_diffusivity, replacing with NaN
08:25:27 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:27 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:27 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:27 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:27 INFO    opendrift:2136: 2026-03-19 08:54:09.382179 - step 50 of 96 - 2999 active elements (1 deactivated)
08:25:27 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:27 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:27 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:27 INFO    opendrift:2136: 2026-03-19 09:24:09.382179 - step 51 of 96 - 2999 active elements (1 deactivated)
08:25:29 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.017256832 to -0.002684246) found for ocean_vertical_diffusivity, replacing with NaN
08:25:29 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:29 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:29 INFO    opendrift:2136: 2026-03-19 09:54:09.382179 - step 52 of 96 - 2999 active elements (1 deactivated)
08:25:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:29 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:29 INFO    opendrift:2136: 2026-03-19 10:24:09.382179 - step 53 of 96 - 2999 active elements (1 deactivated)
08:25:32 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.010985414 to -0.00173109) found for ocean_vertical_diffusivity, replacing with NaN
08:25:32 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:32 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:32 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:32 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:32 INFO    opendrift:2136: 2026-03-19 10:54:09.382179 - step 54 of 96 - 2999 active elements (1 deactivated)
08:25:32 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:32 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:32 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:32 INFO    opendrift:2136: 2026-03-19 11:24:09.382179 - step 55 of 96 - 2999 active elements (1 deactivated)
08:25:34 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.007496473 to -0.0014612336) found for ocean_vertical_diffusivity, replacing with NaN
08:25:34 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:34 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:34 INFO    opendrift:2136: 2026-03-19 11:54:09.382179 - step 56 of 96 - 2999 active elements (1 deactivated)
08:25:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:34 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:34 INFO    opendrift:2136: 2026-03-19 12:24:09.382179 - step 57 of 96 - 2999 active elements (1 deactivated)
08:25:37 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.008019248 to -0.0018854325) found for ocean_vertical_diffusivity, replacing with NaN
08:25:37 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:37 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:37 INFO    opendrift:2136: 2026-03-19 12:54:09.382179 - step 58 of 96 - 2999 active elements (1 deactivated)
08:25:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:37 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:37 INFO    opendrift:2136: 2026-03-19 13:24:09.382179 - step 59 of 96 - 2999 active elements (1 deactivated)
08:25:39 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.011362858 to -0.00097478763) found for ocean_vertical_diffusivity, replacing with NaN
08:25:39 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:39 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:39 INFO    opendrift:2136: 2026-03-19 13:54:09.382179 - step 60 of 96 - 2999 active elements (1 deactivated)
08:25:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:39 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:39 INFO    opendrift:2136: 2026-03-19 14:24:09.382179 - step 61 of 96 - 2999 active elements (1 deactivated)
08:25:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:42 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:42 INFO    opendrift:2136: 2026-03-19 14:54:09.382179 - step 62 of 96 - 2999 active elements (1 deactivated)
08:25:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:42 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:42 INFO    opendrift:2136: 2026-03-19 15:24:09.382179 - step 63 of 96 - 2999 active elements (1 deactivated)
08:25:45 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:45 INFO    opendrift:2136: 2026-03-19 15:54:09.382179 - step 64 of 96 - 2999 active elements (1 deactivated)
08:25:45 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:45 INFO    opendrift:2136: 2026-03-19 16:24:09.382179 - step 65 of 96 - 2999 active elements (1 deactivated)
08:25:47 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0017964853 to -8.260468e-05) found for ocean_vertical_diffusivity, replacing with NaN
08:25:47 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:47 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:47 INFO    opendrift:2136: 2026-03-19 16:54:09.382179 - step 66 of 96 - 2999 active elements (1 deactivated)
08:25:47 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:47 INFO    opendrift:2136: 2026-03-19 17:24:09.382179 - step 67 of 96 - 2999 active elements (1 deactivated)
08:25:50 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.00986558 to -0.00091831374) found for ocean_vertical_diffusivity, replacing with NaN
08:25:50 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:50 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:50 INFO    opendrift:2136: 2026-03-19 17:54:09.382179 - step 68 of 96 - 2999 active elements (1 deactivated)
08:25:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:50 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:50 INFO    opendrift:2136: 2026-03-19 18:24:09.382179 - step 69 of 96 - 2999 active elements (1 deactivated)
08:25:52 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.03650338 to -0.0020284287) found for ocean_vertical_diffusivity, replacing with NaN
08:25:52 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:52 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:52 INFO    opendrift:2136: 2026-03-19 18:54:09.382179 - step 70 of 96 - 2999 active elements (1 deactivated)
08:25:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:53 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:53 INFO    opendrift:2136: 2026-03-19 19:24:09.382179 - step 71 of 96 - 2999 active elements (1 deactivated)
08:25:55 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.044042066 to -0.0017460802) found for ocean_vertical_diffusivity, replacing with NaN
08:25:55 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:55 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:55 INFO    opendrift:2136: 2026-03-19 19:54:09.382179 - step 72 of 96 - 2999 active elements (1 deactivated)
08:25:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:55 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:55 INFO    opendrift:2136: 2026-03-19 20:24:09.382179 - step 73 of 96 - 2999 active elements (1 deactivated)
08:25:57 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.020520544 to -0.0022486984) found for ocean_vertical_diffusivity, replacing with NaN
08:25:57 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:25:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:58 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:58 INFO    opendrift:2136: 2026-03-19 20:54:09.382179 - step 74 of 96 - 2999 active elements (1 deactivated)
08:25:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:25:58 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:25:58 INFO    opendrift:2136: 2026-03-19 21:24:09.382179 - step 75 of 96 - 2999 active elements (1 deactivated)
08:26:00 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.011393199 to -0.0016087126) found for ocean_vertical_diffusivity, replacing with NaN
08:26:00 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:26:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:00 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:00 INFO    opendrift:2136: 2026-03-19 21:54:09.382179 - step 76 of 96 - 2999 active elements (1 deactivated)
08:26:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:00 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:00 INFO    opendrift:2136: 2026-03-19 22:24:09.382179 - step 77 of 96 - 2999 active elements (1 deactivated)
08:26:03 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.009209966 to -0.0011883661) found for ocean_vertical_diffusivity, replacing with NaN
08:26:03 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:26:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:03 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:03 INFO    opendrift:2136: 2026-03-19 22:54:09.382179 - step 78 of 96 - 2999 active elements (1 deactivated)
08:26:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:03 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:03 INFO    opendrift:2136: 2026-03-19 23:24:09.382179 - step 79 of 96 - 2999 active elements (1 deactivated)
08:26:06 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0078893835 to -0.0011795226) found for ocean_vertical_diffusivity, replacing with NaN
08:26:06 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:26:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:06 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:06 INFO    opendrift:2136: 2026-03-19 23:54:09.382179 - step 80 of 96 - 2999 active elements (1 deactivated)
08:26:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:06 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:06 INFO    opendrift:2136: 2026-03-20 00:24:09.382179 - step 81 of 96 - 2999 active elements (1 deactivated)
08:26:09 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.008151861 to -0.0011513788) found for ocean_vertical_diffusivity, replacing with NaN
08:26:09 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:26:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:09 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:09 INFO    opendrift:2136: 2026-03-20 00:54:09.382179 - step 82 of 96 - 2999 active elements (1 deactivated)
08:26:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:09 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:09 INFO    opendrift:2136: 2026-03-20 01:24:09.382179 - step 83 of 96 - 2999 active elements (1 deactivated)
08:26:12 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.015570448 to -0.0007537814) found for ocean_vertical_diffusivity, replacing with NaN
08:26:12 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:26:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:12 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:12 INFO    opendrift:2136: 2026-03-20 01:54:09.382179 - step 84 of 96 - 2999 active elements (1 deactivated)
08:26:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:12 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:12 INFO    opendrift:2136: 2026-03-20 02:24:09.382179 - step 85 of 96 - 2999 active elements (1 deactivated)
08:26:15 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.013726724 to -0.00036403805) found for ocean_vertical_diffusivity, replacing with NaN
08:26:15 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:26:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:15 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:15 INFO    opendrift:2136: 2026-03-20 02:54:09.382179 - step 86 of 96 - 2999 active elements (1 deactivated)
08:26:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:15 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:15 INFO    opendrift:2136: 2026-03-20 03:24:09.382179 - step 87 of 96 - 2999 active elements (1 deactivated)
08:26:18 WARNING opendrift.readers.basereader.variables:656: Invalid values (-4.1307896e-05 to -1.6080891e-06) found for ocean_vertical_diffusivity, replacing with NaN
08:26:18 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:26:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:18 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:18 INFO    opendrift:2136: 2026-03-20 03:54:09.382179 - step 88 of 96 - 2999 active elements (1 deactivated)
08:26:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:18 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:18 INFO    opendrift:2136: 2026-03-20 04:24:09.382179 - step 89 of 96 - 2999 active elements (1 deactivated)
08:26:20 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0041307495 to -1.9958284e-06) found for ocean_vertical_diffusivity, replacing with NaN
08:26:20 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:26:21 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:21 INFO    opendrift:2136: 2026-03-20 04:54:09.382179 - step 90 of 96 - 2999 active elements (1 deactivated)
08:26:21 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:21 INFO    opendrift:2136: 2026-03-20 05:24:09.382179 - step 91 of 96 - 2999 active elements (1 deactivated)
08:26:24 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.010589432 to -0.00044187967) found for ocean_vertical_diffusivity, replacing with NaN
08:26:24 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:26:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:24 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:24 INFO    opendrift:2136: 2026-03-20 05:54:09.382179 - step 92 of 96 - 2999 active elements (1 deactivated)
08:26:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:24 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:24 INFO    opendrift:2136: 2026-03-20 06:24:09.382179 - step 93 of 96 - 2999 active elements (1 deactivated)
08:26:28 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.032017127 to -0.0017200754) found for ocean_vertical_diffusivity, replacing with NaN
08:26:28 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:26:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:28 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:28 INFO    opendrift:2136: 2026-03-20 06:54:09.382179 - step 94 of 96 - 2999 active elements (1 deactivated)
08:26:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:28 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:28 INFO    opendrift:2136: 2026-03-20 07:24:09.382179 - step 95 of 96 - 2999 active elements (1 deactivated)
08:26:31 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.04875413 to -0.0024149227) found for ocean_vertical_diffusivity, replacing with NaN
08:26:31 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
08:26:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:31 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice
08:26:31 INFO    opendrift:2136: 2026-03-20 07:54:09.382179 - step 96 of 96 - 2999 active elements (1 deactivated)
08:26:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
08:26:31 INFO    opendrift.models.openoil.openoil:1193: Ice concentration above 30%, using Nordam scheme for advection in ice

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

xarray.Dataset

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

    cf_role :

    trajectory_id

    dtype :

    <class 'numpy.int32'>

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

  • time
    (time)
    datetime64[ns]
    2026-03-18T08:24:09.382179 ... 2...

    standard_name :

    time

    long_name :

    time

    dtype :

    <class 'numpy.float64'>

    array(['2026-03-18T08:24:09.382179000', '2026-03-18T09:24:09.382179000',
           '2026-03-18T10:24:09.382179000', '2026-03-18T11:24:09.382179000',
           '2026-03-18T12:24:09.382179000', '2026-03-18T13:24:09.382179000',
           '2026-03-18T14:24:09.382179000', '2026-03-18T15:24:09.382179000',
           '2026-03-18T16:24:09.382179000', '2026-03-18T17:24:09.382179000',
           '2026-03-18T18:24:09.382179000', '2026-03-18T19:24:09.382179000',
           '2026-03-18T20:24:09.382179000', '2026-03-18T21:24:09.382179000',
           '2026-03-18T22:24:09.382179000', '2026-03-18T23:24:09.382179000',
           '2026-03-19T00:24:09.382179000', '2026-03-19T01:24:09.382179000',
           '2026-03-19T02:24:09.382179000', '2026-03-19T03:24:09.382179000',
           '2026-03-19T04:24:09.382179000', '2026-03-19T05:24:09.382179000',
           '2026-03-19T06:24:09.382179000', '2026-03-19T07:24:09.382179000',
           '2026-03-19T08:24:09.382179000', '2026-03-19T09:24:09.382179000',
           '2026-03-19T10:24:09.382179000', '2026-03-19T11:24:09.382179000',
           '2026-03-19T12:24:09.382179000', '2026-03-19T13:24:09.382179000',
           '2026-03-19T14:24:09.382179000', '2026-03-19T15:24:09.382179000',
           '2026-03-19T16:24:09.382179000', '2026-03-19T17:24:09.382179000',
           '2026-03-19T18:24:09.382179000', '2026-03-19T19:24:09.382179000',
           '2026-03-19T20:24:09.382179000', '2026-03-19T21:24:09.382179000',
           '2026-03-19T22:24:09.382179000', '2026-03-19T23:24:09.382179000',
           '2026-03-20T00:24:09.382179000', '2026-03-20T01:24:09.382179000',
           '2026-03-20T02:24:09.382179000', '2026-03-20T03:24:09.382179000',
           '2026-03-20T04:24:09.382179000', '2026-03-20T05:24:09.382179000',
           '2026-03-20T06:24:09.382179000', '2026-03-20T07:24:09.382179000',
           '2026-03-20T08:24:09.382179000'], dtype='datetime64[ns]')

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

    dtype :

    <class 'numpy.int32'>

    valid_range :

    [0 1]

    flag_values :

    [0 1]

    flag_meanings :

    active stranded

    array([[0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.],
           ...,
           [0., 0., 0., ..., 0., 0., 0.],
           [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., ..., 1., 1., 1.],
           [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., ..., 165600., 169200., 172800.],
           [     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., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.]], shape=(3000, 49), dtype=float32)

  • lon
    (trajectory, time)
    float32
    19.2 19.17 19.17 ... 18.44 18.43

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    axis :

    X

    minval :

    17.904596

    maxval :

    19.234175

    array([[19.195274, 19.170328, 19.169367, ..., 18.30459 , 18.309666,
            18.300161],
           [19.191893, 19.183203, 19.159924, ..., 18.45804 , 18.457468,
            18.43934 ],
           [19.193327, 19.189678, 19.169493, ..., 18.54503 , 18.54474 ,
            18.548805],
           ...,
           [19.190422, 19.18949 , 19.17141 , ..., 18.667616, 18.68844 ,
            18.69379 ],
           [19.187893, 19.187216, 19.163433, ..., 18.36745 , 18.342323,
            18.348803],
           [19.1907  , 19.178617, 19.159544, ..., 18.426527, 18.4401  ,
            18.434914]], shape=(3000, 49), dtype=float32)

  • lat
    (trajectory, time)
    float32
    79.6 79.6 79.6 ... 79.73 79.73

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    axis :

    Y

    minval :

    79.58677

    maxval :

    79.86448

    array([[79.59931 , 79.599815, 79.600426, ..., 79.7387  , 79.7514  ,
            79.75943 ],
           [79.598854, 79.600136, 79.59915 , ..., 79.70362 , 79.70116 ,
            79.70712 ],
           [79.59855 , 79.59786 , 79.60654 , ..., 79.67899 , 79.67671 ,
            79.68229 ],
           ...,
           [79.59836 , 79.60119 , 79.60689 , ..., 79.672585, 79.67504 ,
            79.67937 ],
           [79.59878 , 79.59928 , 79.606415, ..., 79.7045  , 79.70657 ,
            79.71412 ],
           [79.59873 , 79.60263 , 79.60374 , ..., 79.72581 , 79.72777 ,
            79.73107 ]], 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 :

    -277.43335

    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., ..., 0., 0., 0.],
           [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, ..., 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]],
          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., ..., 1., 1., 1.],
           [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.006213 0.006197

    dtype :

    <class 'numpy.float32'>

    units :

    m/s

    description :

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

    minval :

    0.0

    maxval :

    0.048643064

    array([[0.        , 0.        , 0.        , ..., 0.        , 0.        ,
            0.        ],
           [0.        , 0.        , 0.        , ..., 0.        , 0.        ,
            0.        ],
           [0.        , 0.        , 0.        , ..., 0.03908038, 0.03910714,
            0.03910575],
           ...,
           [0.        , 0.        , 0.        , ..., 0.        , 0.        ,
            0.        ],
           [0.        , 0.        , 0.        , ..., 0.03061241, 0.03061863,
            0.0305756 ],
           [0.        , 0.        , 0.        , ..., 0.00621425, 0.0062135 ,
            0.00619665]], 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, ..., 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]], shape=(3000, 49), dtype=float32)

  • viscosity
    (trajectory, time)
    float32
    0.02169 0.1679 ... 0.1625 0.1635

    dtype :

    <class 'numpy.float32'>

    units :

    m2/s

    description :

    Kinematic viscosity of oil (emulsion)

    minval :

    0.021692334

    maxval :

    0.16871561

    array([[0.02169233, 0.16785607, 0.16767345, ..., 0.15286705, 0.15084207,
            0.15116991],
           [0.02169233, 0.16782467, 0.16779973, ..., 0.16512434, 0.16542529,
            0.16563946],
           [0.02169233, 0.16784076, 0.16780697, ..., 0.16668276, 0.16636667,
            0.16635391],
           ...,
           [0.02169233, 0.16784237, 0.16786009, ..., 0.16744366, 0.16741477,
            0.1675762 ],
           [0.02169233, 0.1678456 , 0.1678609 , ..., 0.16244379, 0.16236292,
            0.16290335],
           [0.02169233, 0.16796558, 0.16783835, ..., 0.162457  , 0.16250058,
            0.16353807]], shape=(3000, 49), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    kg/m^3

    minval :

    988.1

    maxval :

    997.2405

    array([[988.1    , 997.2188 , 997.2142 , ..., 996.82104, 996.7642 ,
            996.7735 ],
           [988.1    , 997.218  , 997.21735, ..., 997.1492 , 997.15686,
            997.16235],
           [988.1    , 997.2184 , 997.2175 , ..., 997.189  , 997.18097,
            997.18066],
           ...,
           [988.1    , 997.21844, 997.2189 , ..., 997.2084 , 997.20764,
            997.21173],
           [988.1    , 997.2185 , 997.21893, ..., 997.07965, 997.0775 ,
            997.0916 ],
           [988.1    , 997.22156, 997.2183 , ..., 997.07996, 997.0811 ,
            997.10815]], 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., ..., 0., 0., 0.],
           [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., ..., 0., 0., 0.],
           [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., ..., 0., 0., 0.],
           [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., ..., 0., 0., 0.],
           [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., ..., 0., 0., 0.],
           [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., ..., 1., 1., 1.],
           [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., ..., 3., 3., 3.],
           [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., ..., 0., 0., 0.],
           [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., ..., 0., 0., 0.],
           [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, ..., 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]],
          shape=(3000, 49), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    m

    minval :

    0.0

    maxval :

    0.0029961192

    array([[0.        , 0.        , 0.        , ..., 0.        , 0.        ,
            0.        ],
           [0.        , 0.        , 0.        , ..., 0.        , 0.        ,
            0.        ],
           [0.        , 0.        , 0.        , ..., 0.00208017, 0.00208017,
            0.00208017],
           ...,
           [0.        , 0.        , 0.        , ..., 0.        , 0.        ,
            0.        ],
           [0.        , 0.        , 0.        , ..., 0.00183094, 0.00183094,
            0.00183094],
           [0.        , 0.        , 0.        , ..., 0.00082493, 0.00082493,
            0.00082493]], shape=(3000, 49), dtype=float32)

  • x_sea_water_velocity
    (trajectory, time)
    float32
    0.03202 -0.04709 ... -0.1148

    minval :

    -0.34112322

    maxval :

    0.40525693

    array([[ 0.03202463, -0.04708779, -0.06891973, ...,  0.02359335,
            -0.06042445, -0.08665338],
           [-0.11692698, -0.02192352, -0.05102093, ...,  0.01390223,
            -0.10303815, -0.06052083],
           [-0.06559531, -0.05346442, -0.09168199, ..., -0.02694835,
            -0.10240863, -0.20052585],
           ...,
           [-0.11239264, -0.03267949, -0.09758457, ..., -0.09215479,
            -0.07039911, -0.03695942],
           [-0.12534332,  0.02293452, -0.0751506 , ...,  0.06690645,
             0.05301801, -0.04398245],
           [-0.03699256, -0.00652943, -0.05112452, ...,  0.18013118,
            -0.0348544 , -0.1148179 ]], shape=(3000, 49), dtype=float32)

  • y_sea_water_velocity
    (trajectory, time)
    float32
    0.1532 0.08778 ... 0.2587 0.2378

    minval :

    -0.30323792

    maxval :

    0.4792071

    array([[ 0.15322474,  0.0877753 ,  0.03957657, ...,  0.34001514,
             0.35668433,  0.34956458],
           [ 0.10414992,  0.0965358 ,  0.04653655, ...,  0.27452093,
             0.26696724,  0.30741873],
           [ 0.07881223,  0.13039048,  0.06176583, ...,  0.14149354,
             0.1571109 ,  0.21514678],
           ...,
           [-0.01777082, -0.01224638, -0.03451151, ...,  0.10413116,
             0.17365052,  0.17066222],
           [ 0.07942259, -0.021503  ,  0.0169006 , ...,  0.27106076,
             0.2916237 ,  0.24240382],
           [ 0.09298141,  0.00189674,  0.03220924, ...,  0.27660018,
             0.25870365,  0.2377526 ]], shape=(3000, 49), dtype=float32)

  • x_wind
    (trajectory, time)
    float32
    -2.432 -7.382 ... -3.889 -4.156

    minval :

    -11.146086

    maxval :

    1.2351601

    array([[-2.432051 , -7.381656 , -9.389392 , ..., -2.9949014, -3.4006174,
            -4.7115564],
           [-2.6575665, -7.967294 , -7.519096 , ..., -3.4247038, -3.4035242,
            -3.1127384],
           [-3.4848685, -6.68146  , -8.387659 , ..., -3.2088683, -2.1797013,
            -2.2543406],
           ...,
           [-2.4846785, -7.286793 , -7.6419435, ..., -2.8715546, -3.1530032,
            -2.4219742],
           [-3.1059206, -7.286898 , -8.688975 , ..., -3.8910515, -4.5926456,
            -3.9396358],
           [-2.9547296, -7.5062985, -8.307768 , ..., -3.1357307, -3.8886154,
            -4.155704 ]], shape=(3000, 49), dtype=float32)

  • y_wind
    (trajectory, time)
    float32
    5.262 7.853 7.814 ... 3.459 3.361

    minval :

    -3.2070353

    maxval :

    9.93701

    array([[5.262234  , 7.853217  , 7.813651  , ..., 3.076181  , 4.3414865 ,
            5.0067706 ],
           [4.2625804 , 7.362263  , 8.433201  , ..., 2.08552   , 2.5912364 ,
            2.0618165 ],
           [3.8081596 , 7.3587155 , 7.6211567 , ..., 0.8213629 , 1.7835411 ,
            0.58243436],
           ...,
           [3.945094  , 8.600196  , 7.620029  , ..., 0.823211  , 2.0665977 ,
            1.0787276 ],
           [3.8151362 , 7.5035057 , 7.4798675 , ..., 0.89313996, 2.3121026 ,
            2.3187866 ],
           [3.9755597 , 6.9676995 , 7.1976376 , ..., 3.144321  , 3.4591188 ,
            3.3610237 ]], shape=(3000, 49), dtype=float32)

  • sea_surface_height
    (trajectory, time)
    float32
    -1.009 -0.7824 ... -0.8849 -1.009

    minval :

    -1.0509715

    maxval :

    0.22749454

    array([[-1.0088246 , -0.7824441 , -0.6687845 , ..., -0.54569393,
            -0.9015695 , -1.0175142 ],
           [-1.0089235 , -0.7824282 , -0.6689595 , ..., -0.5107356 ,
            -0.88367903, -1.0149461 ],
           [-1.0090377 , -0.7827779 , -0.66863966, ..., -0.49158853,
            -0.87699527, -1.0154582 ],
           ...,
           [-1.0090597 , -0.7823101 , -0.66860944, ..., -0.47528905,
            -0.86430013, -1.0074153 ],
           [-1.0088935 , -0.78255594, -0.6687182 , ..., -0.5217195 ,
            -0.8964407 , -1.0205055 ],
           [-1.008946  , -0.7821026 , -0.6688165 , ..., -0.52403367,
            -0.884851  , -1.0086664 ]], shape=(3000, 49), dtype=float32)

  • sea_surface_wave_significant_height
    (trajectory, time)
    float32
    0.8267 2.858 ... 0.6663 0.7027

    minval :

    0.0006431289

    maxval :

    4.429563

    array([[0.8267072 , 2.8575819 , 3.6706598 , ..., 0.45343515, 0.7481525 ,
            1.1627563 ],
           [0.62071335, 2.894945  , 3.1403294 , ..., 0.39551857, 0.4501427 ,
            0.34292957],
           [0.6555011 , 2.4302979 , 3.159497  , ..., 0.26989824, 0.19513005,
            0.13336353],
           ...,
           [0.5347399 , 3.1256936 , 2.865018  , ..., 0.21951814, 0.3496215 ,
            0.17292847],
           [0.5953694 , 2.6912766 , 3.2335892 , ..., 0.3920743 , 0.65037996,
            0.51407844],
           [0.6035733 , 2.5803764 , 2.9722948 , ..., 0.48510122, 0.66633606,
            0.7027323 ]], 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., ..., 0., 0., 0.],
           [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., ..., 0., 0., 0.],
           [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., ..., 0., 0., 0.],
           [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.953 9.209 10.44 ... 4.447 4.567

    minval :

    0.13815881

    maxval :

    11.465949

    array([[ 4.9534235,  9.209349 , 10.437624 , ...,  3.6684885,  4.71221  ,
             5.874541 ],
           [ 4.2921524,  9.269359 ,  9.65422  , ...,  3.4262023,  3.6551456,
             3.190304 ],
           [ 4.410789 ,  8.492965 ,  9.683638 , ...,  2.8302813,  2.406532 ,
             1.9895191],
           ...,
           [ 3.9838307,  9.631696 ,  9.221323 , ...,  2.5524943,  3.2212813,
             2.2654936],
           [ 4.203614 ,  8.937349 ,  9.796524 , ...,  3.4112515,  4.3935256,
             3.9061084],
           [ 4.2324767,  8.751271 ,  9.392377 , ...,  3.794423 ,  4.4470935,
             4.566932 ]], shape=(3000, 49), dtype=float32)

  • sea_ice_area_fraction
    (trajectory, time)
    float32
    0.8499 0.8561 ... 0.7343 0.7266

    minval :

    0.23116264

    maxval :

    0.9590042

    array([[0.8498519 , 0.85607433, 0.84903455, ..., 0.6253505 , 0.64532375,
            0.6347474 ],
           [0.8500652 , 0.8550038 , 0.84724116, ..., 0.74189895, 0.7289704 ,
            0.71791613],
           [0.8472686 , 0.8420654 , 0.8660408 , ..., 0.75764847, 0.741105  ,
            0.73392344],
           ...,
           [0.8485416 , 0.8570273 , 0.8666023 , ..., 0.75038826, 0.7315725 ,
            0.72652626],
           [0.8527366 , 0.8509469 , 0.866957  , ..., 0.69510674, 0.6419698 ,
            0.65860486],
           [0.8503361 , 0.86280715, 0.8603245 , ..., 0.73821884, 0.7342602 ,
            0.72662514]], shape=(3000, 49), dtype=float32)

  • sea_ice_x_velocity
    (trajectory, time)
    float32
    -0.06683 -0.07159 ... 0.007738

    minval :

    -0.19806804

    maxval :

    0.123460524

    array([[-0.06682706, -0.0715929 , -0.08298226, ...,  0.01777646,
             0.00396253, -0.00633034],
           [-0.06734639, -0.07081167, -0.08440272, ...,  0.02259549,
             0.00940765,  0.00055943],
           [-0.06779749, -0.07307111, -0.08434071, ...,  0.0387949 ,
             0.0218038 ,  0.01466168],
           ...,
           [-0.06797165, -0.07028689, -0.08441018, ...,  0.0821842 ,
             0.07844236,  0.06885772],
           [-0.0673214 , -0.07127047, -0.08440293, ...,  0.00321473,
            -0.00406215, -0.00534542],
           [-0.06747466, -0.07009612, -0.08368967, ...,  0.02257076,
             0.01205419,  0.00773795]], shape=(3000, 49), dtype=float32)

  • sea_ice_y_velocity
    (trajectory, time)
    float32
    0.0429 0.04467 ... 0.06714 0.08036

    minval :

    -0.05935711

    maxval :

    0.25169706

    array([[0.04289569, 0.0446682 , 0.05524966, ..., 0.14120242, 0.14056651,
            0.14929657],
           [0.04353057, 0.04524523, 0.05544188, ..., 0.04135025, 0.06310836,
            0.07978033],
           [0.04456671, 0.05104783, 0.04913784, ..., 0.03424614, 0.05222256,
            0.05817362],
           ...,
           [0.0445842 , 0.0439512 , 0.04889475, ..., 0.00888555, 0.01255217,
            0.01534255],
           [0.04307368, 0.0470291 , 0.0489747 , ..., 0.10136338, 0.1319844 ,
            0.12727302],
           [0.04365329, 0.04152827, 0.05140836, ..., 0.05637969, 0.06714275,
            0.0803597 ]], shape=(3000, 49), dtype=float32)

  • sea_water_temperature
    (trajectory, time)
    float32
    -1.828 -1.827 ... -1.657 271.5

    minval :

    -1.8604981

    maxval :

    272.60452

    array([[ -1.8283312,  -1.8271378,  -1.8245183, ...,  -1.2113247,
             -1.1850721, 271.98978  ],
           [ -1.8283664,  -1.8262813,  -1.8254111, ...,  -1.7191588,
             -1.7422321, 271.40643  ],
           [ -1.8275379,  -1.8221794,  -1.8292209, ...,  -1.7794179,
             -1.7717047, 271.37903  ],
           ...,
           [ -1.827896 ,  -1.8266244,  -1.8294306, ...,  -1.8109095,
             -1.8157201, 271.33243  ],
           [ -1.8291545,  -1.825124 ,  -1.8294607, ...,  -1.6168998,
             -1.5920829, 271.51254  ],
           [ -1.8284405,  -1.8282821,  -1.8272544, ...,  -1.6093132,
             -1.6571704, 271.48776  ]], shape=(3000, 49), dtype=float32)

  • sea_water_salinity
    (trajectory, time)
    float32
    34.15 34.15 34.15 ... 34.22 34.22

    minval :

    34.136696

    maxval :

    34.42497

    array([[34.14569 , 34.14795 , 34.1475  , ..., 34.28625 , 34.289642,
            34.29449 ],
           [34.14597 , 34.14614 , 34.14876 , ..., 34.21255 , 34.20745 ,
            34.20638 ],
           [34.145603, 34.14421 , 34.150284, ..., 34.201508, 34.198414,
            34.19729 ],
           ...,
           [34.145947, 34.14573 , 34.150208, ..., 34.200966, 34.20072 ,
            34.20082 ],
           [34.146526, 34.145203, 34.151028, ..., 34.22427 , 34.22871 ,
            34.221848],
           [34.146088, 34.147896, 34.150166, ..., 34.227627, 34.221016,
            34.220318]], shape=(3000, 49), dtype=float32)

  • sea_floor_depth_below_sea_level
    (trajectory, time)
    float32
    170.4 176.4 177.5 ... 248.8 248.2

    minval :

    136.18945

    maxval :

    298.257

    array([[170.42508, 176.43799, 177.50842, ..., 286.5962 , 274.48215,
            275.61252],
           [170.93387, 174.09714, 176.9357 , ..., 266.21533, 268.69095,
            271.27972],
           [170.35571, 170.76001, 181.99438, ..., 265.35223, 268.09534,
            266.95206],
           ...,
           [170.93738, 173.21982, 181.66443, ..., 243.12506, 235.52292,
            232.17055],
           [171.87784, 172.40976, 183.60275, ..., 294.08783, 295.29187,
            289.4547 ],
           [171.14177, 177.08955, 183.08224, ..., 259.5948 , 248.81027,
            248.22571]], shape=(3000, 49), dtype=float32)

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

    minval :

    10.0

    maxval :

    10.0

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

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

    minval :

    0.00022164051

    maxval :

    0.072839335

    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, ..., 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]],
          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., ..., 0., 0., 0.],
           [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., ..., 50., 50., 50.],
           [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-03-18 08:24:09.382179', '2026-03-18 09:24:09.382179',
                   '2026-03-18 10:24:09.382179', '2026-03-18 11:24:09.382179',
                   '2026-03-18 12:24:09.382179', '2026-03-18 13:24:09.382179',
                   '2026-03-18 14:24:09.382179', '2026-03-18 15:24:09.382179',
                   '2026-03-18 16:24:09.382179', '2026-03-18 17:24:09.382179',
                   '2026-03-18 18:24:09.382179', '2026-03-18 19:24:09.382179',
                   '2026-03-18 20:24:09.382179', '2026-03-18 21:24:09.382179',
                   '2026-03-18 22:24:09.382179', '2026-03-18 23:24:09.382179',
                   '2026-03-19 00:24:09.382179', '2026-03-19 01:24:09.382179',
                   '2026-03-19 02:24:09.382179', '2026-03-19 03:24:09.382179',
                   '2026-03-19 04:24:09.382179', '2026-03-19 05:24:09.382179',
                   '2026-03-19 06:24:09.382179', '2026-03-19 07:24:09.382179',
                   '2026-03-19 08:24:09.382179', '2026-03-19 09:24:09.382179',
                   '2026-03-19 10:24:09.382179', '2026-03-19 11:24:09.382179',
                   '2026-03-19 12:24:09.382179', '2026-03-19 13:24:09.382179',
                   '2026-03-19 14:24:09.382179', '2026-03-19 15:24:09.382179',
                   '2026-03-19 16:24:09.382179', '2026-03-19 17:24:09.382179',
                   '2026-03-19 18:24:09.382179', '2026-03-19 19:24:09.382179',
                   '2026-03-19 20:24:09.382179', '2026-03-19 21:24:09.382179',
                   '2026-03-19 22:24:09.382179', '2026-03-19 23:24:09.382179',
                   '2026-03-20 00:24:09.382179', '2026-03-20 01:24:09.382179',
                   '2026-03-20 02:24:09.382179', '2026-03-20 03:24:09.382179',
                   '2026-03-20 04:24:09.382179', '2026-03-20 05:24:09.382179',
                   '2026-03-20 06:24:09.382179', '2026-03-20 07:24:09.382179',
                   '2026-03-20 08:24:09.382179'],
                  dtype='datetime64[ns]', name='time', freq='h'))

• Attributes: (167)

  Conventions :

  CF-1.11, ACDD-1.3

  standard_name_vocabulary :

  CF Standard Name Table v85

  featureType :

  trajectory

  title :

  OpenDrift trajectory simulation

  summary :

  Output from simulation with OpenDrift framework

  keywords :

  trajectory, drift, lagrangian, simulation

  history :

  Created 2026-03-25 08:24:13.660116

  date_created :

  2026-03-25T08:24:13.660133

  source :

  Output from simulation with OpenDrift

  model_url :

  https://github.com/OpenDrift/opendrift

  opendrift_class :

  OpenOil

  opendrift_module :

  opendrift.models.openoil.openoil

  readers :

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

  time_coverage_start :

  2026-03-18 08:24:09.382179

  time_step_calculation :

  0:30:00

  time_step_output :

  1:00:00

  config_environment:constant:x_sea_water_velocity :

  None

  config_environment:fallback:x_sea_water_velocity :

  None

  config_environment:constant:y_sea_water_velocity :

  None

  config_environment:fallback:y_sea_water_velocity :

  None

  config_environment:constant:x_wind :

  None

  config_environment:fallback:x_wind :

  None

  config_environment:constant:y_wind :

  None

  config_environment:fallback:y_wind :

  None

  config_environment:constant:sea_surface_height :

  None

  config_environment:fallback:sea_surface_height :

  0

  config_environment:constant:upward_sea_water_velocity :

  None

  config_environment:fallback:upward_sea_water_velocity :

  0

  config_environment:constant:sea_surface_wave_significant_height :

  None

  config_environment:fallback:sea_surface_wave_significant_height :

  0

  config_environment:constant:sea_surface_wave_stokes_drift_x_velocity :

  None

  config_environment:fallback:sea_surface_wave_stokes_drift_x_velocity :

  0

  config_environment:constant:sea_surface_wave_stokes_drift_y_velocity :

  None

  config_environment:fallback:sea_surface_wave_stokes_drift_y_velocity :

  0

  config_environment:constant:sea_surface_wave_period_at_variance_spectral_density_maximum :

  None

  config_environment:fallback:sea_surface_wave_period_at_variance_spectral_density_maximum :

  0

  config_environment:constant:sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment :

  None

  config_environment:fallback:sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment :

  0

  config_environment:constant:sea_ice_area_fraction :

  None

  config_environment:fallback:sea_ice_area_fraction :

  0

  config_environment:constant:sea_ice_x_velocity :

  None

  config_environment:fallback:sea_ice_x_velocity :

  0

  config_environment:constant:sea_ice_y_velocity :

  None

  config_environment:fallback:sea_ice_y_velocity :

  0

  config_environment:constant:sea_water_temperature :

  None

  config_environment:fallback:sea_water_temperature :

  10

  config_environment:constant:sea_water_salinity :

  None

  config_environment:fallback:sea_water_salinity :

  34

  config_environment:constant:sea_floor_depth_below_sea_level :

  None

  config_environment:fallback:sea_floor_depth_below_sea_level :

  10000

  config_environment:constant:horizontal_diffusivity :

  10.0

  config_environment:fallback:horizontal_diffusivity :

  0

  config_environment:constant:ocean_vertical_diffusivity :

  None

  config_environment:fallback:ocean_vertical_diffusivity :

  0.02

  config_environment:constant:land_binary_mask :

  None

  config_environment:fallback:land_binary_mask :

  None

  config_environment:constant:ocean_mixed_layer_thickness :

  None

  config_environment:fallback:ocean_mixed_layer_thickness :

  50

  config_general:use_auto_landmask :

  True

  config_drift:current_uncertainty :

  0.05

  config_drift:current_uncertainty_uniform :

  0

  config_drift:max_speed :

  1.3

  config_readers:max_number_of_fails :

  1

  config_general:simulation_name :

  config_general:coastline_action :

  stranding

  config_general:coastline_approximation_precision :

  0.001

  config_general:time_step_minutes :

  60

  config_general:time_step_output_minutes :

  None

  config_seed:ocean_only :

  True

  config_seed:number :

  1

  config_drift:max_age_seconds :

  None

  config_drift:advection_scheme :

  euler

  config_drift:profiles_depth :

  50

  config_drift:wind_uncertainty :

  0.5

  config_drift:relative_wind :

  False

  config_drift:deactivate_north_of :

  None

  config_drift:deactivate_south_of :

  None

  config_drift:deactivate_east_of :

  None

  config_drift:deactivate_west_of :

  None

  config_seed:origin_marker :

  0

  config_seed:z :

  0

  config_seed:wind_drift_factor :

  0.03

  config_seed:current_drift_factor :

  1

  config_seed:terminal_velocity :

  0.0

  config_seed:oil_film_thickness :

  0.001

  config_drift:vertical_advection :

  False

  config_drift:vertical_advection_at_surface :

  False

  config_drift:water_column_stretching :

  False

  config_drift:vertical_advection_correction :

  False

  config_drift:vertical_mixing :

  True

  config_drift:vertical_mixing_at_surface :

  False

  config_vertical_mixing:timestep :

  60

  config_vertical_mixing:diffusivitymodel :

  environment

  config_vertical_mixing:background_diffusivity :

  1.2e-05

  config_vertical_mixing:TSprofiles :

  False

  config_drift:wind_drift_depth :

  0.1

  config_drift:stokes_drift :

  True

  config_drift:stokes_drift_profile :

  Phillips

  config_drift:use_tabularised_stokes_drift :

  False

  config_drift:tabularised_stokes_drift_fetch :

  25000

  config_general:seafloor_action :

  lift_to_seafloor

  config_drift:truncate_ocean_model_below_m :

  3.0

  config_seed:seafloor :

  False

  config_seed:m3_per_hour :

  1

  config_seed:droplet_size_distribution :

  uniform

  config_seed:droplet_diameter_mu :

  0.001

  config_seed:droplet_diameter_sigma :

  0.0005

  config_seed:droplet_diameter_min_subsea :

  0.0005

  config_seed:droplet_diameter_max_subsea :

  0.005

  config_processes:dispersion :

  False

  config_processes:evaporation :

  False

  config_processes:emulsification :

  False

  config_processes:biodegradation :

  False

  config_biodegradation:method :

  Adcroft

  config_processes:update_oilfilm_thickness :

  False

  config_wave_entrainment:droplet_size_distribution :

  Johansen et al. (2015)

  config_wave_entrainment:entrainment_rate :

  Li et al. (2017)

  config_seed:oil_type :

  GENERIC BUNKER C

  opendrift_version :

  1.14.9

  seed_lon :

  19.1909

  seed_lat :

  79.5986

  seed_radius :

  50.0

  seed_time :

  2026-03-18 08:24:09.382179

  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-03-25 08:24:13.670882

  reader_x_sea_water_velocity :

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

  reader_y_sea_water_velocity :

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

  reader_x_wind :

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

  reader_y_wind :

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

  reader_sea_surface_height :

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

  reader_sea_surface_wave_significant_height :

  0

  reader_sea_surface_wave_stokes_drift_x_velocity :

  0

  reader_sea_surface_wave_stokes_drift_y_velocity :

  0

  reader_sea_surface_wave_period_at_variance_spectral_density_maximum :

  0

  reader_sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment :

  0

  reader_sea_ice_area_fraction :

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

  reader_sea_ice_x_velocity :

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

  reader_sea_ice_y_velocity :

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

  reader_sea_water_temperature :

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

  reader_sea_water_salinity :

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

  reader_sea_floor_depth_below_sea_level :

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

  reader_horizontal_diffusivity :

  10.0

  reader_ocean_vertical_diffusivity :

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

  reader_land_binary_mask :

  ['global_landmask']

  reader_ocean_mixed_layer_thickness :

  50

  time_coverage_end :

  2026-03-20 08:24:09.382179

  time_coverage_duration :

  P2DT0H0M0S

  time_coverage_resolution :

  P0DT0H30M0S

  performance :

  -------------------- Reader performance: -------------------- constant_reader 0:00:00.0 total 0:00:00.0 preparing 0:00:00.0 reading 0:00:00.0 masking -------------------- global_landmask 0:00: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:02:02.6 total 0:00:00.0 preparing 0:02:01.5 reading 0:00:01.6 interpolation 0:00:00.1 interpolation_time 0:00:01.0 rotating vectors 0:00:00.0 masking -------------------- Performance: 2:22.0 total time 4.3 configuration 0.0 preparing main loop 0.0 moving elements to ocean 2:17.7 main loop 5.9 updating elements 0.0 oil weathering 0.0 updating viscosities 0.0 updating densities 5.4 vertical mixing 0.0 cleaning up --------------------

  geospatial_bounds_crs :

  EPSG:4326

  geospatial_bounds_vertical_crs :

  EPSG:5831

  geospatial_lat_min :

  79.58677

  geospatial_lat_max :

  79.86448

  geospatial_lat_units :

  degrees_north

  geospatial_lat_resolution :

  point

  geospatial_lon_min :

  17.904596

  geospatial_lon_max :

  19.234175

  geospatial_lon_units :

  degrees_east

  geospatial_lon_resolution :

  point

  runtime :

  0:02:22.090100

  geospatial_vertical_min :

  -277.43335

  geospatial_vertical_max :

  0.0

  geospatial_vertical_positive :

  up

Print and plot results

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

===========================
--------------------
Reader performance:
--------------------
constant_reader
 0:00:00.0  total
 0:00:00.0  preparing
 0:00:00.0  reading
 0:00:00.0  masking
--------------------
global_landmask
 0:00: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:02:02.6  total
 0:00:00.0  preparing
 0:02:01.5  reading
 0:00:01.6  interpolation
 0:00:00.1  interpolation_time
 0:00:01.0  rotating vectors
 0:00:00.0  masking
--------------------
Performance:
 2:22.0 total time
    4.3 configuration
    0.0 preparing main loop
      0.0 moving elements to ocean
 2:17.7 main loop
      5.9 updating elements
        0.0 oil weathering
          0.0 updating viscosities
          0.0 updating densities
        5.4 vertical mixing
    0.0 cleaning up
--------------------
===========================
Model:  OpenOil     (OpenDrift version 1.14.9)
        2998 active Oil particles  (2 deactivated, 0 scheduled)
-------------------
Environment variables:
  -----
  horizontal_diffusivity
     1) constant_reader
  -----
  land_binary_mask
     1) global_landmask
  -----
  ocean_vertical_diffusivity
  sea_floor_depth_below_sea_level
  sea_ice_area_fraction
  sea_ice_x_velocity
  sea_ice_y_velocity
  sea_surface_height
  sea_water_salinity
  sea_water_temperature
  x_sea_water_velocity
  x_wind
  y_sea_water_velocity
  y_wind
     1) https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
  -----
Readers not added for the following variables:
  ocean_mixed_layer_thickness
  sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment
  sea_surface_wave_period_at_variance_spectral_density_maximum
  sea_surface_wave_significant_height
  sea_surface_wave_stokes_drift_x_velocity
  sea_surface_wave_stokes_drift_y_velocity

Discarded readers:

Time:
        Start: 2026-03-18 08:24:09.382179 UTC
        Present: 2026-03-20 08:24:09.382179 UTC
        Calculation steps: 96 * 0:30:00 - total time: 2 days, 0:00:00
        Output steps: 49 * 1:00:00
===========================

08:26:42 INFO    opendrift:4697: Saving animation to /root/project/docs/source/gallery/animations/example_oil_ice_0.gif...
08:27:24 INFO    opendrift:3120: Time to make animation: 0:00:53.326056

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-03-18 08:24 to 2026-03-20 08:24 UTC (49 steps)'}>, <Figure size 616.885x1100 with 2 Axes>)