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)

13:50:43 INFO    opendrift:576: OpenDriftSimulation initialised (version 1.14.7 / v1.14.7-15-g2988cd4)

Using live data from Barents 2.5 km ocean model

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

Adjusting some configuration

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

Imaginary oil spill in Hinlopen strait

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

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

13:50:48 INFO    opendrift:1803: Skipping environment variable upward_sea_water_velocity because of condition ['drift:vertical_advection', 'is', False]
13:50:48 INFO    opendrift:1814: Storing previous values of element property lon because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
13:50:48 INFO    opendrift:1814: Storing previous values of element property lat because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
13:50:48 INFO    opendrift:952: Using existing reader for land_binary_mask to move elements to ocean
13:50:48 INFO    opendrift:982: All points are in ocean
13:50:48 INFO    opendrift.models.openoil.openoil:692: Oil-water surface tension is 0.035935 Nm
13:50:48 INFO    opendrift.models.openoil.openoil:705: Max water fraction not available for GENERIC BUNKER C, using default
13:50:48 INFO    opendrift:2111: 2025-12-09 13:50:43.939787 - step 1 of 96 - 3000 active elements (0 deactivated)
13:50:48 INFO    opendrift.readers:63: Opening file with xr.open_dataset
13:50:51 INFO    opendrift.readers.reader_netCDF_CF_generic:332: Detected dimensions: {'x': 'X', 'y': 'Y', 'z': 'depth', 'time': 'time'}
13:50:51 INFO    opendrift.readers.basereader:176: Variable x_sea_water_velocity will be rotated from eastward_sea_water_velocity
13:50:51 INFO    opendrift.readers.basereader:176: Variable y_sea_water_velocity will be rotated from northward_sea_water_velocity
13:50:55 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.010476603 to -0.0023701536) found for ocean_vertical_diffusivity, replacing with NaN
13:50:55 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:50:57 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.009612147 to -0.002447895) found for ocean_vertical_diffusivity, replacing with NaN
13:50:57 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:50:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:50:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:50:57 INFO    opendrift:2111: 2025-12-09 14:20:43.939787 - step 2 of 96 - 3000 active elements (0 deactivated)
13:51:00 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.008075659 to -0.0018657603) found for ocean_vertical_diffusivity, replacing with NaN
13:51:00 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:00 INFO    opendrift:2111: 2025-12-09 14:50:43.939787 - step 3 of 96 - 3000 active elements (0 deactivated)
13:51:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:00 INFO    opendrift:2111: 2025-12-09 15:20:43.939787 - step 4 of 96 - 3000 active elements (0 deactivated)
13:51:02 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.005026889 to -0.0012187747) found for ocean_vertical_diffusivity, replacing with NaN
13:51:02 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:02 INFO    opendrift:2111: 2025-12-09 15:50:43.939787 - step 5 of 96 - 3000 active elements (0 deactivated)
13:51:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:03 INFO    opendrift:2111: 2025-12-09 16:20:43.939787 - step 6 of 96 - 3000 active elements (0 deactivated)
13:51:07 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.002544912 to -0.00047396286) found for ocean_vertical_diffusivity, replacing with NaN
13:51:07 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:07 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:07 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:09 INFO    opendrift:2111: 2025-12-09 16:50:43.939787 - step 7 of 96 - 2999 active elements (1 deactivated)
13:51:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:09 INFO    opendrift:2111: 2025-12-09 17:20:43.939787 - step 8 of 96 - 2967 active elements (33 deactivated)
13:51:11 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0013180607 to -0.00021271949) found for ocean_vertical_diffusivity, replacing with NaN
13:51:11 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:11 INFO    opendrift:2111: 2025-12-09 17:50:43.939787 - step 9 of 96 - 2696 active elements (304 deactivated)
13:51:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:12 INFO    opendrift:2111: 2025-12-09 18:20:43.939787 - step 10 of 96 - 2205 active elements (795 deactivated)
13:51:14 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0019980574 to -0.0004811797) found for ocean_vertical_diffusivity, replacing with NaN
13:51:14 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:14 INFO    opendrift:2111: 2025-12-09 18:50:43.939787 - step 11 of 96 - 1817 active elements (1183 deactivated)
13:51:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:15 INFO    opendrift:2111: 2025-12-09 19:20:43.939787 - step 12 of 96 - 1643 active elements (1357 deactivated)
13:51:17 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0042834417 to -0.0011818648) found for ocean_vertical_diffusivity, replacing with NaN
13:51:17 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:17 INFO    opendrift:2111: 2025-12-09 19:50:43.939787 - step 13 of 96 - 1589 active elements (1411 deactivated)
13:51:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:17 INFO    opendrift:2111: 2025-12-09 20:20:43.939787 - step 14 of 96 - 1565 active elements (1435 deactivated)
13:51:21 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0070529 to -0.0012001814) found for ocean_vertical_diffusivity, replacing with NaN
13:51:21 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:21 INFO    opendrift:2111: 2025-12-09 20:50:43.939787 - step 15 of 96 - 1561 active elements (1439 deactivated)
13:51:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:21 INFO    opendrift:2111: 2025-12-09 21:20:43.939787 - step 16 of 96 - 1559 active elements (1441 deactivated)
13:51:23 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.009614366 to -0.0014920095) found for ocean_vertical_diffusivity, replacing with NaN
13:51:23 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:23 INFO    opendrift:2111: 2025-12-09 21:50:43.939787 - step 17 of 96 - 1558 active elements (1442 deactivated)
13:51:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:23 INFO    opendrift:2111: 2025-12-09 22:20:43.939787 - step 18 of 96 - 1558 active elements (1442 deactivated)
13:51:26 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.010822952 to -0.0019242986) found for ocean_vertical_diffusivity, replacing with NaN
13:51:26 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:26 INFO    opendrift:2111: 2025-12-09 22:50:43.939787 - step 19 of 96 - 1558 active elements (1442 deactivated)
13:51:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:26 INFO    opendrift:2111: 2025-12-09 23:20:43.939787 - step 20 of 96 - 1558 active elements (1442 deactivated)
13:51:29 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.011456842 to -0.0025908295) found for ocean_vertical_diffusivity, replacing with NaN
13:51:29 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:29 INFO    opendrift:2111: 2025-12-09 23:50:43.939787 - step 21 of 96 - 1557 active elements (1443 deactivated)
13:51:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:29 INFO    opendrift:2111: 2025-12-10 00:20:43.939787 - step 22 of 96 - 1557 active elements (1443 deactivated)
13:51:31 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0139277885 to -0.0025465752) found for ocean_vertical_diffusivity, replacing with NaN
13:51:31 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:31 INFO    opendrift:2111: 2025-12-10 00:50:43.939787 - step 23 of 96 - 1557 active elements (1443 deactivated)
13:51:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:32 INFO    opendrift:2111: 2025-12-10 01:20:43.939787 - step 24 of 96 - 1557 active elements (1443 deactivated)
13:51:34 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.01478635 to -0.0022433798) found for ocean_vertical_diffusivity, replacing with NaN
13:51:34 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:34 INFO    opendrift:2111: 2025-12-10 01:50:43.939787 - step 25 of 96 - 1557 active elements (1443 deactivated)
13:51:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:34 INFO    opendrift:2111: 2025-12-10 02:20:43.939787 - step 26 of 96 - 1557 active elements (1443 deactivated)
13:51:36 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.013612992 to -0.0021732396) found for ocean_vertical_diffusivity, replacing with NaN
13:51:36 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:36 INFO    opendrift:2111: 2025-12-10 02:50:43.939787 - step 27 of 96 - 1557 active elements (1443 deactivated)
13:51:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:37 INFO    opendrift:2111: 2025-12-10 03:20:43.939787 - step 28 of 96 - 1557 active elements (1443 deactivated)
13:51:39 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0102647245 to -0.0018974994) found for ocean_vertical_diffusivity, replacing with NaN
13:51:39 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:39 INFO    opendrift:2111: 2025-12-10 03:50:43.939787 - step 29 of 96 - 1557 active elements (1443 deactivated)
13:51:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:39 INFO    opendrift:2111: 2025-12-10 04:20:43.939787 - step 30 of 96 - 1557 active elements (1443 deactivated)
13:51:42 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.007059658 to -0.00095628493) found for ocean_vertical_diffusivity, replacing with NaN
13:51:42 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:42 INFO    opendrift:2111: 2025-12-10 04:50:43.939787 - step 31 of 96 - 1557 active elements (1443 deactivated)
13:51:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:42 INFO    opendrift:2111: 2025-12-10 05:20:43.939787 - step 32 of 96 - 1557 active elements (1443 deactivated)
13:51:44 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0046234033 to -0.00060337887) found for ocean_vertical_diffusivity, replacing with NaN
13:51:44 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:44 INFO    opendrift:2111: 2025-12-10 05:50:43.939787 - step 33 of 96 - 1557 active elements (1443 deactivated)
13:51:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:45 INFO    opendrift:2111: 2025-12-10 06:20:43.939787 - step 34 of 96 - 1557 active elements (1443 deactivated)
13:51:47 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.005421038 to -0.00094032835) found for ocean_vertical_diffusivity, replacing with NaN
13:51:47 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:47 INFO    opendrift:2111: 2025-12-10 06:50:43.939787 - step 35 of 96 - 1557 active elements (1443 deactivated)
13:51:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:47 INFO    opendrift:2111: 2025-12-10 07:20:43.939787 - step 36 of 96 - 1557 active elements (1443 deactivated)
13:51:50 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.00680867 to -0.0008816579) found for ocean_vertical_diffusivity, replacing with NaN
13:51:50 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:50 INFO    opendrift:2111: 2025-12-10 07:50:43.939787 - step 37 of 96 - 1557 active elements (1443 deactivated)
13:51:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:50 INFO    opendrift:2111: 2025-12-10 08:20:43.939787 - step 38 of 96 - 1557 active elements (1443 deactivated)
13:51:53 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0092988545 to -0.0018833883) found for ocean_vertical_diffusivity, replacing with NaN
13:51:53 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:53 INFO    opendrift:2111: 2025-12-10 08:50:43.939787 - step 39 of 96 - 1557 active elements (1443 deactivated)
13:51:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:53 INFO    opendrift:2111: 2025-12-10 09:20:43.939787 - step 40 of 96 - 1557 active elements (1443 deactivated)
13:51:55 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.010726811 to -0.0012735456) found for ocean_vertical_diffusivity, replacing with NaN
13:51:55 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:55 INFO    opendrift:2111: 2025-12-10 09:50:43.939787 - step 41 of 96 - 1557 active elements (1443 deactivated)
13:51:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:56 INFO    opendrift:2111: 2025-12-10 10:20:43.939787 - step 42 of 96 - 1557 active elements (1443 deactivated)
13:51:58 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0120592965 to -0.00090785435) found for ocean_vertical_diffusivity, replacing with NaN
13:51:58 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:51:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:58 INFO    opendrift:2111: 2025-12-10 10:50:43.939787 - step 43 of 96 - 1557 active elements (1443 deactivated)
13:51:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:51:58 INFO    opendrift:2111: 2025-12-10 11:20:43.939787 - step 44 of 96 - 1557 active elements (1443 deactivated)
13:52:01 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.012173303 to -0.0015821628) found for ocean_vertical_diffusivity, replacing with NaN
13:52:01 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:01 INFO    opendrift:2111: 2025-12-10 11:50:43.939787 - step 45 of 96 - 1557 active elements (1443 deactivated)
13:52:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:01 INFO    opendrift:2111: 2025-12-10 12:20:43.939787 - step 46 of 96 - 1557 active elements (1443 deactivated)
13:52:03 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.012493687 to -0.0015446678) found for ocean_vertical_diffusivity, replacing with NaN
13:52:03 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:03 INFO    opendrift:2111: 2025-12-10 12:50:43.939787 - step 47 of 96 - 1557 active elements (1443 deactivated)
13:52:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:04 INFO    opendrift:2111: 2025-12-10 13:20:43.939787 - step 48 of 96 - 1557 active elements (1443 deactivated)
13:52:06 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.011995294 to -0.0015613668) found for ocean_vertical_diffusivity, replacing with NaN
13:52:06 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:06 INFO    opendrift:2111: 2025-12-10 13:50:43.939787 - step 49 of 96 - 1557 active elements (1443 deactivated)
13:52:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:06 INFO    opendrift:2111: 2025-12-10 14:20:43.939787 - step 50 of 96 - 1557 active elements (1443 deactivated)
13:52:08 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.010438122 to -0.0025474504) found for ocean_vertical_diffusivity, replacing with NaN
13:52:08 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:09 INFO    opendrift:2111: 2025-12-10 14:50:43.939787 - step 51 of 96 - 1557 active elements (1443 deactivated)
13:52:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:09 INFO    opendrift:2111: 2025-12-10 15:20:43.939787 - step 52 of 96 - 1557 active elements (1443 deactivated)
13:52:11 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.007718426 to -0.0013320418) found for ocean_vertical_diffusivity, replacing with NaN
13:52:11 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:11 INFO    opendrift:2111: 2025-12-10 15:50:43.939787 - step 53 of 96 - 1557 active elements (1443 deactivated)
13:52:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:11 INFO    opendrift:2111: 2025-12-10 16:20:43.939787 - step 54 of 96 - 1557 active elements (1443 deactivated)
13:52:14 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.004921337 to -0.0011187592) found for ocean_vertical_diffusivity, replacing with NaN
13:52:14 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:14 INFO    opendrift:2111: 2025-12-10 16:50:43.939787 - step 55 of 96 - 1557 active elements (1443 deactivated)
13:52:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:14 INFO    opendrift:2111: 2025-12-10 17:20:43.939787 - step 56 of 96 - 1557 active elements (1443 deactivated)
13:52:17 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.004156104 to -0.0006267975) found for ocean_vertical_diffusivity, replacing with NaN
13:52:17 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:17 INFO    opendrift:2111: 2025-12-10 17:50:43.939787 - step 57 of 96 - 1557 active elements (1443 deactivated)
13:52:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:17 INFO    opendrift:2111: 2025-12-10 18:20:43.939787 - step 58 of 96 - 1557 active elements (1443 deactivated)
13:52:19 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0043098032 to -0.0008154284) found for ocean_vertical_diffusivity, replacing with NaN
13:52:19 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:19 INFO    opendrift:2111: 2025-12-10 18:50:43.939787 - step 59 of 96 - 1557 active elements (1443 deactivated)
13:52:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:19 INFO    opendrift:2111: 2025-12-10 19:20:43.939787 - step 60 of 96 - 1557 active elements (1443 deactivated)
13:52:22 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0067687193 to -0.0013994476) found for ocean_vertical_diffusivity, replacing with NaN
13:52:22 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:22 INFO    opendrift:2111: 2025-12-10 19:50:43.939787 - step 61 of 96 - 1557 active elements (1443 deactivated)
13:52:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:22 INFO    opendrift:2111: 2025-12-10 20:20:43.939787 - step 62 of 96 - 1557 active elements (1443 deactivated)
13:52:24 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.009558479 to -0.0020647815) found for ocean_vertical_diffusivity, replacing with NaN
13:52:24 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:24 INFO    opendrift:2111: 2025-12-10 20:50:43.939787 - step 63 of 96 - 1557 active elements (1443 deactivated)
13:52:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:25 INFO    opendrift:2111: 2025-12-10 21:20:43.939787 - step 64 of 96 - 1556 active elements (1444 deactivated)
13:52:27 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.011648568 to -0.0024313496) found for ocean_vertical_diffusivity, replacing with NaN
13:52:27 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:27 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:27 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:27 INFO    opendrift:2111: 2025-12-10 21:50:43.939787 - step 65 of 96 - 1556 active elements (1444 deactivated)
13:52:27 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:27 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:27 INFO    opendrift:2111: 2025-12-10 22:20:43.939787 - step 66 of 96 - 1556 active elements (1444 deactivated)
13:52:30 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.012797238 to -0.0026342964) found for ocean_vertical_diffusivity, replacing with NaN
13:52:30 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:30 INFO    opendrift:2111: 2025-12-10 22:50:43.939787 - step 67 of 96 - 1554 active elements (1446 deactivated)
13:52:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:30 INFO    opendrift:2111: 2025-12-10 23:20:43.939787 - step 68 of 96 - 1554 active elements (1446 deactivated)
13:52:33 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.017103847 to -0.0027950045) found for ocean_vertical_diffusivity, replacing with NaN
13:52:33 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:33 INFO    opendrift:2111: 2025-12-10 23:50:43.939787 - step 69 of 96 - 1554 active elements (1446 deactivated)
13:52:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:33 INFO    opendrift:2111: 2025-12-11 00:20:43.939787 - step 70 of 96 - 1553 active elements (1447 deactivated)
13:52:35 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.016479056 to -0.0038252953) found for ocean_vertical_diffusivity, replacing with NaN
13:52:35 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:35 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:35 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:35 INFO    opendrift:2111: 2025-12-11 00:50:43.939787 - step 71 of 96 - 1549 active elements (1451 deactivated)
13:52:35 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:35 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:35 INFO    opendrift:2111: 2025-12-11 01:20:43.939787 - step 72 of 96 - 1546 active elements (1454 deactivated)
13:52:37 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.018067304 to -0.00323275) found for ocean_vertical_diffusivity, replacing with NaN
13:52:37 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:38 INFO    opendrift:2111: 2025-12-11 01:50:43.939787 - step 73 of 96 - 1539 active elements (1461 deactivated)
13:52:38 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:38 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:38 INFO    opendrift:2111: 2025-12-11 02:20:43.939787 - step 74 of 96 - 1533 active elements (1467 deactivated)
13:52:40 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.015525835 to -0.003287039) found for ocean_vertical_diffusivity, replacing with NaN
13:52:40 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:40 INFO    opendrift:2111: 2025-12-11 02:50:43.939787 - step 75 of 96 - 1531 active elements (1469 deactivated)
13:52:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:40 INFO    opendrift:2111: 2025-12-11 03:20:43.939787 - step 76 of 96 - 1530 active elements (1470 deactivated)
13:52:42 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.012161441 to -0.0028824056) found for ocean_vertical_diffusivity, replacing with NaN
13:52:42 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:42 INFO    opendrift:2111: 2025-12-11 03:50:43.939787 - step 77 of 96 - 1528 active elements (1472 deactivated)
13:52:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:42 INFO    opendrift:2111: 2025-12-11 04:20:43.939787 - step 78 of 96 - 1528 active elements (1472 deactivated)
13:52:45 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.009276141 to -0.0016857365) found for ocean_vertical_diffusivity, replacing with NaN
13:52:45 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:45 INFO    opendrift:2111: 2025-12-11 04:50:43.939787 - step 79 of 96 - 1528 active elements (1472 deactivated)
13:52:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:45 INFO    opendrift:2111: 2025-12-11 05:20:43.939787 - step 80 of 96 - 1528 active elements (1472 deactivated)
13:52:47 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0074663507 to -0.0009660067) found for ocean_vertical_diffusivity, replacing with NaN
13:52:47 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:48 INFO    opendrift:2111: 2025-12-11 05:50:43.939787 - step 81 of 96 - 1528 active elements (1472 deactivated)
13:52:48 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:48 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:48 INFO    opendrift:2111: 2025-12-11 06:20:43.939787 - step 82 of 96 - 1528 active elements (1472 deactivated)
13:52:50 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0072873905 to -0.0010959395) found for ocean_vertical_diffusivity, replacing with NaN
13:52:50 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:50 INFO    opendrift:2111: 2025-12-11 06:50:43.939787 - step 83 of 96 - 1528 active elements (1472 deactivated)
13:52:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:50 INFO    opendrift:2111: 2025-12-11 07:20:43.939787 - step 84 of 96 - 1527 active elements (1473 deactivated)
13:52:53 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.007519722 to -0.0007659003) found for ocean_vertical_diffusivity, replacing with NaN
13:52:53 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:53 INFO    opendrift:2111: 2025-12-11 07:50:43.939787 - step 85 of 96 - 1526 active elements (1474 deactivated)
13:52:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:53 INFO    opendrift:2111: 2025-12-11 08:20:43.939787 - step 86 of 96 - 1525 active elements (1475 deactivated)
13:52:55 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.008229474 to -0.0016921724) found for ocean_vertical_diffusivity, replacing with NaN
13:52:55 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:56 INFO    opendrift:2111: 2025-12-11 08:50:43.939787 - step 87 of 96 - 1524 active elements (1476 deactivated)
13:52:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:56 INFO    opendrift:2111: 2025-12-11 09:20:43.939787 - step 88 of 96 - 1523 active elements (1477 deactivated)
13:52:58 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.011215512 to -0.0017815083) found for ocean_vertical_diffusivity, replacing with NaN
13:52:58 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:52:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:58 INFO    opendrift:2111: 2025-12-11 09:50:43.939787 - step 89 of 96 - 1517 active elements (1483 deactivated)
13:52:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:52:58 INFO    opendrift:2111: 2025-12-11 10:20:43.939787 - step 90 of 96 - 1512 active elements (1488 deactivated)
13:53:01 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0126090385 to -0.0020051124) found for ocean_vertical_diffusivity, replacing with NaN
13:53:01 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:53:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:53:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:53:01 INFO    opendrift:2111: 2025-12-11 10:50:43.939787 - step 91 of 96 - 1503 active elements (1497 deactivated)
13:53:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:53:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:53:01 INFO    opendrift:2111: 2025-12-11 11:20:43.939787 - step 92 of 96 - 1501 active elements (1499 deactivated)
13:53:03 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.014098457 to -0.0016772362) found for ocean_vertical_diffusivity, replacing with NaN
13:53:03 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:53:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:53:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:53:04 INFO    opendrift:2111: 2025-12-11 11:50:43.939787 - step 93 of 96 - 1500 active elements (1500 deactivated)
13:53:04 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:53:04 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:53:04 INFO    opendrift:2111: 2025-12-11 12:20:43.939787 - step 94 of 96 - 1499 active elements (1501 deactivated)
13:53:06 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.014878233 to -0.00094262854) found for ocean_vertical_diffusivity, replacing with NaN
13:53:06 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:53:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:53:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:53:06 INFO    opendrift:2111: 2025-12-11 12:50:43.939787 - step 95 of 96 - 1499 active elements (1501 deactivated)
13:53:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:53:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:53:06 INFO    opendrift:2111: 2025-12-11 13:20:43.939787 - step 96 of 96 - 1499 active elements (1501 deactivated)
13:53:10 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.015567615 to -0.0010407055) found for ocean_vertical_diffusivity, replacing with NaN
13:53:10 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:53:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:53:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning

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

xarray.Dataset

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

    cf_role :

    trajectory_id

    dtype :

    <class 'numpy.int32'>

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

  • time
    (time)
    datetime64[ns]
    2025-12-09T13:50:43.939787 ... 2...

    standard_name :

    time

    long_name :

    time

    dtype :

    <class 'numpy.float64'>

    array(['2025-12-09T13:50:43.939787000', '2025-12-09T14:50:43.939787000',
           '2025-12-09T15:50:43.939787000', '2025-12-09T16:50:43.939787000',
           '2025-12-09T17:50:43.939787000', '2025-12-09T18:50:43.939787000',
           '2025-12-09T19:50:43.939787000', '2025-12-09T20:50:43.939787000',
           '2025-12-09T21:50:43.939787000', '2025-12-09T22:50:43.939787000',
           '2025-12-09T23:50:43.939787000', '2025-12-10T00:50:43.939787000',
           '2025-12-10T01:50:43.939787000', '2025-12-10T02:50:43.939787000',
           '2025-12-10T03:50:43.939787000', '2025-12-10T04:50:43.939787000',
           '2025-12-10T05:50:43.939787000', '2025-12-10T06:50:43.939787000',
           '2025-12-10T07:50:43.939787000', '2025-12-10T08:50:43.939787000',
           '2025-12-10T09:50:43.939787000', '2025-12-10T10:50:43.939787000',
           '2025-12-10T11:50:43.939787000', '2025-12-10T12:50:43.939787000',
           '2025-12-10T13:50:43.939787000', '2025-12-10T14:50:43.939787000',
           '2025-12-10T15:50:43.939787000', '2025-12-10T16:50:43.939787000',
           '2025-12-10T17:50:43.939787000', '2025-12-10T18:50:43.939787000',
           '2025-12-10T19:50:43.939787000', '2025-12-10T20:50:43.939787000',
           '2025-12-10T21:50:43.939787000', '2025-12-10T22:50:43.939787000',
           '2025-12-10T23:50:43.939787000', '2025-12-11T00:50:43.939787000',
           '2025-12-11T01:50:43.939787000', '2025-12-11T02:50:43.939787000',
           '2025-12-11T03:50:43.939787000', '2025-12-11T04:50:43.939787000',
           '2025-12-11T05:50:43.939787000', '2025-12-11T06:50:43.939787000',
           '2025-12-11T07:50:43.939787000', '2025-12-11T08:50:43.939787000',
           '2025-12-11T09:50:43.939787000', '2025-12-11T10:50:43.939787000',
           '2025-12-11T11:50:43.939787000', '2025-12-11T12:50:43.939787000',
           '2025-12-11T13:50:43.939787000'], dtype='datetime64[ns]')

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

    dtype :

    <class 'numpy.int32'>

    valid_range :

    [0 1]

    flag_values :

    [0 1]

    flag_meanings :

    active stranded

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

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

    dtype :

    <class 'numpy.int32'>

    minval :

    0

    maxval :

    1

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

  • age_seconds
    (trajectory, time)
    float32
    0.0 3.6e+03 7.2e+03 ... nan nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    s

    minval :

    0.0

    maxval :

    172800.0

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

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

    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., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan]],
          shape=(3000, 49), dtype=float32)

  • lon
    (trajectory, time)
    float32
    19.2 19.17 19.14 ... nan nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    axis :

    X

    minval :

    18.18124

    maxval :

    19.3017

    array([[19.195274, 19.168226, 19.140812, ...,       nan,       nan,
                  nan],
           [19.191893, 19.161783, 19.171453, ..., 18.378195, 18.292667,
            18.228619],
           [19.193327, 19.185986, 19.184488, ...,       nan,       nan,
                  nan],
           ...,
           [19.190422, 19.187866, 19.191723, ...,       nan,       nan,
                  nan],
           [19.187893, 19.14277 , 19.140345, ...,       nan,       nan,
                  nan],
           [19.1907  , 19.18075 , 19.169445, ...,       nan,       nan,
                  nan]], shape=(3000, 49), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    axis :

    Y

    minval :

    79.51721

    maxval :

    79.8118

    array([[79.59931 , 79.600464, 79.59671 , ...,       nan,       nan,
                  nan],
           [79.598854, 79.59687 , 79.595695, ..., 79.764725, 79.774956,
            79.778595],
           [79.59855 , 79.59911 , 79.6003  , ...,       nan,       nan,
                  nan],
           ...,
           [79.59836 , 79.59576 , 79.59479 , ...,       nan,       nan,
                  nan],
           [79.59878 , 79.59881 , 79.595024, ...,       nan,       nan,
                  nan],
           [79.59873 , 79.59765 , 79.5969  , ...,       nan,       nan,
                  nan]], shape=(3000, 49), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    m

    standard_name :

    z

    long_name :

    vertical position

    axis :

    Z

    positive :

    up

    minval :

    0.0

    maxval :

    0.0

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

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

    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, ...,  nan,  nan,  nan],
           [0.03, 0.03, 0.03, ..., 0.03, 0.03, 0.03],
           [0.03, 0.03, 0.03, ...,  nan,  nan,  nan],
           ...,
           [0.03, 0.03, 0.03, ...,  nan,  nan,  nan],
           [0.03, 0.03, 0.03, ...,  nan,  nan,  nan],
           [0.03, 0.03, 0.03, ...,  nan,  nan,  nan]],
          shape=(3000, 49), dtype=float32)

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

    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., ..., nan, nan, nan],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ..., nan, nan, nan],
           ...,
           [ 1.,  1.,  1., ..., nan, nan, nan],
           [ 1.,  1.,  1., ..., nan, nan, nan],
           [ 1.,  1.,  1., ..., nan, nan, nan]],
          shape=(3000, 49), dtype=float32)

  • terminal_velocity
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... nan nan nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    m/s

    description :

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

    minval :

    0.0

    maxval :

    0.048275106

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

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

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.32936665

    maxval :

    0.32936665

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

  • viscosity
    (trajectory, time)
    float32
    0.02169 0.1366 0.1366 ... nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    m2/s

    description :

    Kinematic viscosity of oil (emulsion)

    minval :

    0.021692334

    maxval :

    0.14919476

    array([[0.02169233, 0.13659516, 0.136553  , ...,        nan,        nan,
                   nan],
           [0.02169233, 0.13661721, 0.13665031, ..., 0.14647959, 0.1460538 ,
            0.14576118],
           [0.02169233, 0.13670938, 0.1366594 , ...,        nan,        nan,
                   nan],
           ...,
           [0.02169233, 0.13667758, 0.13672887, ...,        nan,        nan,
                   nan],
           [0.02169233, 0.13662696, 0.13661203, ...,        nan,        nan,
                   nan],
           [0.02169233, 0.1366289 , 0.13648881, ...,        nan,        nan,
                   nan]], shape=(3000, 49), dtype=float32)

  • density
    (trajectory, time)
    float32
    988.1 996.3 996.3 ... nan nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    kg/m^3

    minval :

    988.1

    maxval :

    996.7174

    array([[988.1    , 996.3402 , 996.33887, ...,       nan,       nan,
                  nan],
           [988.1    , 996.3409 , 996.3419 , ..., 996.639  , 996.6266 ,
            996.618  ],
           [988.1    , 996.34375, 996.3422 , ...,       nan,       nan,
                  nan],
           ...,
           [988.1    , 996.3428 , 996.34436, ...,       nan,       nan,
                  nan],
           [988.1    , 996.3412 , 996.3407 , ...,       nan,       nan,
                  nan],
           [988.1    , 996.34125, 996.33685, ...,       nan,       nan,
                  nan]], shape=(3000, 49), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    s

    minval :

    0.0

    maxval :

    0.0

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

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

    dtype :

    <class 'numpy.float32'>

    units :

    m2

    minval :

    0.0

    maxval :

    0.0

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

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

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.0

    maxval :

    0.0

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

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

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.0

    maxval :

    0.0

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

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

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.0

    maxval :

    0.0

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

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

    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., ..., nan, nan, nan],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ..., nan, nan, nan],
           ...,
           [ 1.,  1.,  1., ..., nan, nan, nan],
           [ 1.,  1.,  1., ..., nan, nan, nan],
           [ 1.,  1.,  1., ..., nan, nan, nan]],
          shape=(3000, 49), dtype=float32)

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

    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., ..., nan, nan, nan],
           [ 3.,  3.,  3., ...,  3.,  3.,  3.],
           [ 3.,  3.,  3., ..., nan, nan, nan],
           ...,
           [ 3.,  3.,  3., ..., nan, nan, nan],
           [ 3.,  3.,  3., ..., nan, nan, nan],
           [ 3.,  3.,  3., ..., nan, nan, nan]],
          shape=(3000, 49), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    %

    minval :

    0.0

    maxval :

    0.0

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

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

    dtype :

    <class 'numpy.float32'>

    units :

    %

    minval :

    0.0

    maxval :

    0.0

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

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

    dtype :

    <class 'numpy.float32'>

    units :

    m

    minval :

    0.001

    maxval :

    0.001

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

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

    dtype :

    <class 'numpy.float32'>

    units :

    m

    minval :

    0.0

    maxval :

    0.0029642938

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

  • x_sea_water_velocity
    (trajectory, time)
    float32
    0.04882 -0.1049 0.07916 ... nan nan

    minval :

    -0.48157084

    maxval :

    0.30787897

    array([[ 0.04881758, -0.10493048,  0.07916058, ...,         nan,
                    nan,         nan],
           [-0.10048569,  0.05191835,  0.02496086, ..., -0.1174022 ,
            -0.11067627, -0.13091078],
           [-0.04858498, -0.02157702,  0.06072359, ...,         nan,
                    nan,         nan],
           ...,
           [-0.0958524 ,  0.01991254,  0.07036792, ...,         nan,
                    nan,         nan],
           [-0.10970891,  0.05634799,  0.11450619, ...,         nan,
                    nan,         nan],
           [-0.02070761,  0.00901726, -0.03566274, ...,         nan,
                    nan,         nan]], shape=(3000, 49), dtype=float32)

  • y_sea_water_velocity
    (trajectory, time)
    float32
    0.219 -0.0723 -0.2494 ... nan nan

    minval :

    -0.5088361

    maxval :

    0.56098723

    array([[ 0.2190409 , -0.0723011 , -0.24943009, ...,         nan,
                    nan,         nan],
           [ 0.16612974,  0.0147258 , -0.20532297, ...,  0.03142954,
             0.00863571,  0.07884721],
           [ 0.14190815,  0.00128305, -0.2331111 , ...,         nan,
                    nan,         nan],
           ...,
           [ 0.04224634,  0.0293895 , -0.09982891, ...,         nan,
                    nan,         nan],
           [ 0.13732214, -0.01681053, -0.21742097, ...,         nan,
                    nan,         nan],
           [ 0.153649  ,  0.02692819, -0.19539833, ...,         nan,
                    nan,         nan]], shape=(3000, 49), dtype=float32)

  • x_wind
    (trajectory, time)
    float32
    -1.327 -0.9793 -0.9838 ... nan nan

    minval :

    -9.529897

    maxval :

    4.424456

    array([[-1.3270767 , -0.9792882 , -0.983787  , ...,         nan,
                    nan,         nan],
           [-1.5777341 , -1.0700068 , -1.6331271 , ..., -8.570448  ,
            -7.9184985 , -7.6711626 ],
           [-2.3780484 , -2.2400208 , -2.213462  , ...,         nan,
                    nan,         nan],
           ...,
           [-1.405965  , -0.95417625, -1.2871306 , ...,         nan,
                    nan,         nan],
           [-2.0705621 , -1.2889957 , -1.5004244 , ...,         nan,
                    nan,         nan],
           [-1.8849431 , -2.4844568 , -1.1905462 , ...,         nan,
                    nan,         nan]], shape=(3000, 49), dtype=float32)

  • y_wind
    (trajectory, time)
    float32
    -1.822 -3.21 -1.714 ... nan nan nan

    minval :

    -4.690834

    maxval :

    7.4608884

    array([[-1.8215622, -3.209966 , -1.7141391, ...,        nan,        nan,
                   nan],
           [-2.7617633, -1.6889145, -2.8539894, ...,  0.893284 ,  2.1654527,
             2.1839144],
           [-3.2614346, -2.2228413, -2.4578145, ...,        nan,        nan,
                   nan],
           ...,
           [-3.0655236, -2.6403558, -1.5813134, ...,        nan,        nan,
                   nan],
           [-3.1205268, -2.8124583, -2.7711535, ...,        nan,        nan,
                   nan],
           [-3.0263567, -1.6200793, -2.9443917, ...,        nan,        nan,
                   nan]], shape=(3000, 49), dtype=float32)

  • sea_surface_height
    (trajectory, time)
    float32
    -0.8181 -0.6761 -0.4571 ... nan nan

    minval :

    -1.020889

    maxval :

    -0.12343603

    array([[-0.81807995, -0.6760642 , -0.4571412 , ...,         nan,
                    nan,         nan],
           [-0.8183964 , -0.6770083 , -0.45995104, ..., -0.6996801 ,
            -0.7632353 , -0.7825239 ],
           [-0.8184969 , -0.67727786, -0.46001154, ...,         nan,
                    nan,         nan],
           ...,
           [-0.8186732 , -0.67859036, -0.46185327, ...,         nan,
                    nan,         nan],
           [-0.8185507 , -0.67531216, -0.45738584, ...,         nan,
                    nan,         nan],
           [-0.81848913, -0.6776465 , -0.45956215, ...,         nan,
                    nan,         nan]], shape=(3000, 49), dtype=float32)

  • sea_surface_wave_significant_height
    (trajectory, time)
    float32
    0.1249 0.2771 0.09609 ... nan nan

    minval :

    2.0024085e-05

    maxval :

    2.2547238

    array([[0.12494883, 0.27706704, 0.0960903 , ...,        nan,        nan,
                   nan],
           [0.24886788, 0.09833473, 0.26598403, ..., 1.8265629 , 1.657838  ,
            1.5649588 ],
           [0.4007849 , 0.24498442, 0.26913056, ...,        nan,        nan,
                   nan],
           ...,
           [0.27980462, 0.1938955 , 0.10226851, ...,        nan,        nan,
                   nan],
           [0.3450129 , 0.23545723, 0.24429189, ...,        nan,        nan,
                   nan],
           [0.3127114 , 0.21641067, 0.2481363 , ...,        nan,        nan,
                   nan]], shape=(3000, 49), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

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

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

    minval :

    0.0

    maxval :

    0.0

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

    minval :

    0.0

    maxval :

    0.0

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

  • sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment
    (trajectory, time)
    float32
    1.926 2.868 1.689 ... nan nan nan

    minval :

    0.024378441

    maxval :

    8.180432

    array([[1.9257312, 2.8676226, 1.6887655, ...,       nan,       nan,
                  nan],
           [2.7177782, 1.7083743, 2.809683 , ..., 7.3628664, 7.0145636,
            6.8152385],
           [3.448937 , 2.69649  , 2.8262534, ...,       nan,       nan,
                  nan],
           ...,
           [2.881755 , 2.3989072, 1.7422101, ...,       nan,       nan,
                  nan],
           [3.19998  , 2.6435382, 2.692676 , ...,       nan,       nan,
                  nan],
           [3.0465016, 2.5343635, 2.7137806, ...,       nan,       nan,
                  nan]], shape=(3000, 49), dtype=float32)

  • sea_ice_area_fraction
    (trajectory, time)
    float32
    4.797e-06 3.542e-06 ... nan nan

    minval :

    0.0

    maxval :

    0.0724162

    array([[4.7973358e-06, 3.5416717e-06, 3.4201535e-06, ...,           nan,
                      nan,           nan],
           [2.8943782e-06, 1.4170680e-06, 0.0000000e+00, ..., 2.3907527e-02,
            3.9272744e-02, 4.4297107e-02],
           [2.7323397e-06, 1.0899439e-06, 3.1416157e-06, ...,           nan,
                      nan,           nan],
           ...,
           [1.4935536e-06, 0.0000000e+00, 0.0000000e+00, ...,           nan,
                      nan,           nan],
           [1.5145941e-06, 4.9005826e-06, 2.3068442e-06, ...,           nan,
                      nan,           nan],
           [2.2959443e-06, 0.0000000e+00, 4.7271254e-07, ...,           nan,
                      nan,           nan]], shape=(3000, 49), dtype=float32)

  • sea_ice_x_velocity
    (trajectory, time)
    float32
    -0.00159 0.0 0.0 ... nan nan nan

    minval :

    -0.27889296

    maxval :

    0.09194609

    array([[-1.5897190e-03,  0.0000000e+00,  0.0000000e+00, ...,
                       nan,            nan,            nan],
           [-9.6004282e-04,  0.0000000e+00,  0.0000000e+00, ...,
            -7.9935163e-02, -1.5417255e-01, -1.7461710e-01],
           [-9.0401934e-04,  0.0000000e+00,  7.0994000e-05, ...,
                       nan,            nan,            nan],
           ...,
           [-4.9498910e-04,  0.0000000e+00,  0.0000000e+00, ...,
                       nan,            nan,            nan],
           [-5.0388253e-04,  0.0000000e+00,  0.0000000e+00, ...,
                       nan,            nan,            nan],
           [-7.6192286e-04,  0.0000000e+00,  0.0000000e+00, ...,
                       nan,            nan,            nan]],
          shape=(3000, 49), dtype=float32)

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

    minval :

    -0.12241626

    maxval :

    0.29631364

    array([[0.00395192, 0.        , 0.        , ...,        nan,        nan,
                   nan],
           [0.0023862 , 0.        , 0.        , ..., 0.02869869, 0.06252968,
            0.09323971],
           [0.00224711, 0.        , 0.00011558, ...,        nan,        nan,
                   nan],
           ...,
           [0.00123021, 0.        , 0.        , ...,        nan,        nan,
                   nan],
           [0.00125216, 0.        , 0.        , ...,        nan,        nan,
                   nan],
           [0.00189366, 0.        , 0.        , ...,        nan,        nan,
                   nan]], shape=(3000, 49), dtype=float32)

  • sea_water_temperature
    (trajectory, time)
    float32
    -0.5062 -0.5078 -0.503 ... nan nan

    minval :

    -1.0801514

    maxval :

    272.4563

    array([[ -0.50621456,  -0.50777775,  -0.50301725, ...,          nan,
                     nan,          nan],
           [ -0.507483  ,  -0.5179258 ,  -0.5145264 , ...,  -0.9531509 ,
             -0.9327376 , 272.223     ],
           [ -0.5074593 ,  -0.51556844,  -0.5062683 , ...,          nan,
                     nan,          nan],
           ...,
           [ -0.50833815,  -0.51850545,  -0.5149657 , ...,          nan,
                     nan,          nan],
           [ -0.5085174 ,  -0.506661  ,  -0.50796974, ...,          nan,
                     nan,          nan],
           [ -0.507891  ,  -0.5180763 ,  -0.5118096 , ...,          nan,
                     nan,          nan]], shape=(3000, 49), dtype=float32)

  • sea_water_salinity
    (trajectory, time)
    float32
    33.89 33.9 33.91 ... nan nan nan

    minval :

    33.80649

    maxval :

    33.932465

    array([[33.894478, 33.901424, 33.906044, ...,       nan,       nan,
                  nan],
           [33.89592 , 33.903748, 33.903168, ..., 33.820686, 33.826946,
            33.830246],
           [33.895584, 33.898693, 33.899082, ...,       nan,       nan,
                  nan],
           ...,
           [33.89671 , 33.89922 , 33.899963, ...,       nan,       nan,
                  nan],
           [33.897354, 33.905228, 33.906784, ...,       nan,       nan,
                  nan],
           [33.896404, 33.900448, 33.90309 , ...,       nan,       nan,
                  nan]], shape=(3000, 49), dtype=float32)

  • sea_floor_depth_below_sea_level
    (trajectory, time)
    float32
    170.4 177.7 176.1 ... nan nan nan

    minval :

    72.72264

    maxval :

    293.2448

    array([[170.42511, 177.73451, 176.05069, ...,       nan,       nan,
                  nan],
           [170.93388, 173.12967, 169.95255, ..., 225.60559, 254.94885,
            266.51163],
           [170.35571, 172.60345, 173.88156, ...,       nan,       nan,
                  nan],
           ...,
           [170.93738, 167.96635, 166.02603, ...,       nan,       nan,
                  nan],
           [171.87787, 179.0039 , 173.47043, ...,       nan,       nan,
                  nan],
           [171.14177, 171.68904, 172.06404, ...,       nan,       nan,
                  nan]], shape=(3000, 49), dtype=float32)

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

    minval :

    0.02

    maxval :

    0.02

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

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

    minval :

    0.0

    maxval :

    0.0

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

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

    minval :

    50.0

    maxval :

    50.0

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

• Attributes: (165)

  Conventions :

  CF-1.11, ACDD-1.3

  standard_name_vocabulary :

  CF Standard Name Table v85

  featureType :

  trajectory

  title :

  OpenDrift trajectory simulation

  summary :

  Output from simulation with OpenDrift framework

  keywords :

  trajectory, drift, lagrangian, simulation

  history :

  Created 2025-12-16 13:50:48.848105

  date_created :

  2025-12-16T13:50:48.848125

  source :

  Output from simulation with OpenDrift

  model_url :

  https://github.com/OpenDrift/opendrift

  opendrift_class :

  OpenOil

  opendrift_module :

  opendrift.models.openoil.openoil

  readers :

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

  time_coverage_start :

  2025-12-09 13:50:43.939787

  time_step_calculation :

  0:30:00

  time_step_output :

  1:00:00

  config_environment:constant:x_sea_water_velocity :

  None

  config_environment:fallback:x_sea_water_velocity :

  None

  config_environment:constant:y_sea_water_velocity :

  None

  config_environment:fallback:y_sea_water_velocity :

  None

  config_environment:constant:x_wind :

  None

  config_environment:fallback:x_wind :

  None

  config_environment:constant:y_wind :

  None

  config_environment:fallback:y_wind :

  None

  config_environment:constant:sea_surface_height :

  None

  config_environment:fallback:sea_surface_height :

  0

  config_environment:constant:upward_sea_water_velocity :

  None

  config_environment:fallback:upward_sea_water_velocity :

  0

  config_environment:constant:sea_surface_wave_significant_height :

  None

  config_environment:fallback:sea_surface_wave_significant_height :

  0

  config_environment:constant:sea_surface_wave_stokes_drift_x_velocity :

  None

  config_environment:fallback:sea_surface_wave_stokes_drift_x_velocity :

  0

  config_environment:constant:sea_surface_wave_stokes_drift_y_velocity :

  None

  config_environment:fallback:sea_surface_wave_stokes_drift_y_velocity :

  0

  config_environment:constant:sea_surface_wave_period_at_variance_spectral_density_maximum :

  None

  config_environment:fallback:sea_surface_wave_period_at_variance_spectral_density_maximum :

  0

  config_environment:constant:sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment :

  None

  config_environment:fallback:sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment :

  0

  config_environment:constant:sea_ice_area_fraction :

  None

  config_environment:fallback:sea_ice_area_fraction :

  0

  config_environment:constant:sea_ice_x_velocity :

  None

  config_environment:fallback:sea_ice_x_velocity :

  0

  config_environment:constant:sea_ice_y_velocity :

  None

  config_environment:fallback:sea_ice_y_velocity :

  0

  config_environment:constant:sea_water_temperature :

  None

  config_environment:fallback:sea_water_temperature :

  10

  config_environment:constant:sea_water_salinity :

  None

  config_environment:fallback:sea_water_salinity :

  34

  config_environment:constant:sea_floor_depth_below_sea_level :

  None

  config_environment:fallback:sea_floor_depth_below_sea_level :

  10000

  config_environment:constant:ocean_vertical_diffusivity :

  None

  config_environment:fallback:ocean_vertical_diffusivity :

  0.02

  config_environment:constant:land_binary_mask :

  None

  config_environment:fallback:land_binary_mask :

  None

  config_environment:constant:ocean_mixed_layer_thickness :

  None

  config_environment:fallback:ocean_mixed_layer_thickness :

  50

  config_general:use_auto_landmask :

  True

  config_drift:current_uncertainty :

  0.05

  config_drift:current_uncertainty_uniform :

  0

  config_drift:max_speed :

  1.3

  config_readers:max_number_of_fails :

  1

  config_general:simulation_name :

  config_general:coastline_action :

  stranding

  config_general:coastline_approximation_precision :

  0.001

  config_general:time_step_minutes :

  60

  config_general:time_step_output_minutes :

  None

  config_seed:ocean_only :

  True

  config_seed:number :

  1

  config_drift:max_age_seconds :

  None

  config_drift:advection_scheme :

  euler

  config_drift:horizontal_diffusivity :

  10.0

  config_drift:profiles_depth :

  50

  config_drift:wind_uncertainty :

  0.5

  config_drift:relative_wind :

  False

  config_drift:deactivate_north_of :

  None

  config_drift:deactivate_south_of :

  None

  config_drift:deactivate_east_of :

  None

  config_drift:deactivate_west_of :

  None

  config_seed:origin_marker :

  0

  config_seed:z :

  0

  config_seed:wind_drift_factor :

  0.03

  config_seed:current_drift_factor :

  1

  config_seed:terminal_velocity :

  0.0

  config_seed:oil_film_thickness :

  0.001

  config_drift:vertical_advection :

  False

  config_drift:vertical_advection_at_surface :

  False

  config_drift:water_column_stretching :

  False

  config_drift:vertical_advection_correction :

  False

  config_drift:vertical_mixing :

  True

  config_drift:vertical_mixing_at_surface :

  False

  config_vertical_mixing:timestep :

  60

  config_vertical_mixing:diffusivitymodel :

  environment

  config_vertical_mixing:background_diffusivity :

  1.2e-05

  config_vertical_mixing:TSprofiles :

  False

  config_drift:wind_drift_depth :

  0.1

  config_drift:stokes_drift :

  True

  config_drift:stokes_drift_profile :

  Phillips

  config_drift:use_tabularised_stokes_drift :

  False

  config_drift:tabularised_stokes_drift_fetch :

  25000

  config_general:seafloor_action :

  lift_to_seafloor

  config_drift:truncate_ocean_model_below_m :

  3.0

  config_seed:seafloor :

  False

  config_seed:m3_per_hour :

  1

  config_seed:droplet_size_distribution :

  uniform

  config_seed:droplet_diameter_mu :

  0.001

  config_seed:droplet_diameter_sigma :

  0.0005

  config_seed:droplet_diameter_min_subsea :

  0.0005

  config_seed:droplet_diameter_max_subsea :

  0.005

  config_processes:dispersion :

  False

  config_processes:evaporation :

  False

  config_processes:emulsification :

  False

  config_processes:biodegradation :

  False

  config_biodegradation:method :

  Adcroft

  config_processes:update_oilfilm_thickness :

  False

  config_wave_entrainment:droplet_size_distribution :

  Johansen et al. (2015)

  config_wave_entrainment:entrainment_rate :

  Li et al. (2017)

  config_seed:oil_type :

  GENERIC BUNKER C

  opendrift_version :

  1.14.7

  seed_lon :

  19.1909

  seed_lat :

  79.5986

  seed_radius :

  50.0

  seed_time :

  2025-12-09 13:50:43.939787

  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 :

  2025-12-16 13:50:48.861006

  reader_x_sea_water_velocity :

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

  reader_y_sea_water_velocity :

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

  reader_x_wind :

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

  reader_y_wind :

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

  reader_sea_surface_height :

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

  reader_sea_surface_wave_significant_height :

  0

  reader_sea_surface_wave_stokes_drift_x_velocity :

  0

  reader_sea_surface_wave_stokes_drift_y_velocity :

  0

  reader_sea_surface_wave_period_at_variance_spectral_density_maximum :

  0

  reader_sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment :

  0

  reader_sea_ice_area_fraction :

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

  reader_sea_ice_x_velocity :

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

  reader_sea_ice_y_velocity :

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

  reader_sea_water_temperature :

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

  reader_sea_water_salinity :

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

  reader_sea_floor_depth_below_sea_level :

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

  reader_ocean_vertical_diffusivity :

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

  reader_land_binary_mask :

  ['global_landmask']

  reader_ocean_mixed_layer_thickness :

  50

  time_coverage_end :

  2025-12-11 13:50:43.939787

  time_coverage_duration :

  P2DT0H0M0S

  time_coverage_resolution :

  P0DT0H30M0S

  performance :

  -------------------- Reader performance: -------------------- global_landmask 0:00:01.7 total 0:00:00.0 preparing 0:00:01.7 reading 0:00:00.0 masking -------------------- https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be 0:02:05.1 total 0:00:00.0 preparing 0:02:04.3 reading 0:00:01.3 interpolation 0:00:00.1 interpolation_time 0:00:00.8 rotating vectors 0:00:00.0 masking -------------------- Performance: 2:26.3 total time 4.9 configuration 0.0 preparing main loop 0.0 moving elements to ocean 2:21.3 main loop 5.8 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.51721

  geospatial_lat_max :

  79.8118

  geospatial_lat_units :

  degrees_north

  geospatial_lat_resolution :

  point

  geospatial_lon_min :

  18.18124

  geospatial_lon_max :

  19.3017

  geospatial_lon_units :

  degrees_east

  geospatial_lon_resolution :

  point

  runtime :

  0:02:26.342717

  geospatial_vertical_min :

  0.0

  geospatial_vertical_max :

  0.0

  geospatial_vertical_positive :

  up

Print and plot results

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

===========================
--------------------
Reader performance:
--------------------
global_landmask
 0:00:01.7  total
 0:00:00.0  preparing
 0:00:01.7  reading
 0:00:00.0  masking
--------------------
https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
 0:02:05.1  total
 0:00:00.0  preparing
 0:02:04.3  reading
 0:00:01.3  interpolation
 0:00:00.1  interpolation_time
 0:00:00.8  rotating vectors
 0:00:00.0  masking
--------------------
Performance:
 2:26.3 total time
    4.9 configuration
    0.0 preparing main loop
      0.0 moving elements to ocean
 2:21.3 main loop
      5.8 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.7)
        1499 active Oil particles  (1501 deactivated, 0 scheduled)
-------------------
Environment variables:
  -----
  land_binary_mask
     1) global_landmask
  -----
  ocean_vertical_diffusivity
  sea_floor_depth_below_sea_level
  sea_ice_area_fraction
  sea_ice_x_velocity
  sea_ice_y_velocity
  sea_surface_height
  sea_water_salinity
  sea_water_temperature
  x_sea_water_velocity
  x_wind
  y_sea_water_velocity
  y_wind
     1) https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
  -----
Readers not added for the following variables:
  ocean_mixed_layer_thickness
  sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment
  sea_surface_wave_period_at_variance_spectral_density_maximum
  sea_surface_wave_significant_height
  sea_surface_wave_stokes_drift_x_velocity
  sea_surface_wave_stokes_drift_y_velocity

Discarded readers:

Time:
        Start: 2025-12-09 13:50:43.939787 UTC
        Present: 2025-12-11 13:50:43.939787 UTC
        Calculation steps: 96 * 0:30:00 - total time: 2 days, 0:00:00
        Output steps: 49 * 1:00:00
===========================

13:53:21 INFO    opendrift:4660: Saving animation to /root/project/docs/source/gallery/animations/example_oil_ice_0.gif...
13:54:05 INFO    opendrift:3094: Time to make animation: 0:00:55.223169

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)\n2025-12-09 13:50 to 2025-12-11 13:50 UTC (49 steps)'}>, <Figure size 560.947x1100 with 2 Axes>)