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:17:26 INFO    opendrift:576: OpenDriftSimulation initialised (version 1.14.8 / v1.14.8-4-g48afd36)

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:17:26 INFO    opendrift.models.openoil.openoil:1708: Oil type not specified, using default: GENERIC BUNKER C
13:17:26 INFO    opendrift.models.openoil.adios.dirjs:86: Querying ADIOS database for oil: GENERIC BUNKER C
13:17:26 INFO    opendrift.models.openoil.openoil:1717: Using density 988.1 and viscosity 0.02169233387797564 of oiltype GENERIC BUNKER C
13:17:26 INFO    opendrift.models.basemodel.environment:203: Adding a global landmask from GSHHG
13:17:30 INFO    opendrift.models.basemodel.environment:227: Fallback values will be used for the following variables which have no readers:
13:17:30 INFO    opendrift.models.basemodel.environment:230:    sea_surface_height: 0.000000
13:17:30 INFO    opendrift.models.basemodel.environment:230:    upward_sea_water_velocity: 0.000000
13:17:30 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_significant_height: 0.000000
13:17:30 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
13:17:30 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
13:17:30 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
13:17:30 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
13:17:30 INFO    opendrift.models.basemodel.environment:230:    sea_ice_area_fraction: 0.000000
13:17:30 INFO    opendrift.models.basemodel.environment:230:    sea_ice_x_velocity: 0.000000
13:17:30 INFO    opendrift.models.basemodel.environment:230:    sea_ice_y_velocity: 0.000000
13:17:30 INFO    opendrift.models.basemodel.environment:230:    sea_water_temperature: 10.000000
13:17:30 INFO    opendrift.models.basemodel.environment:230:    sea_water_salinity: 34.000000
13:17:30 INFO    opendrift.models.basemodel.environment:230:    sea_floor_depth_below_sea_level: 10000.000000
13:17:30 INFO    opendrift.models.basemodel.environment:230:    ocean_vertical_diffusivity: 0.020000
13:17:30 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:17:30 INFO    opendrift:1810: Skipping environment variable upward_sea_water_velocity because of condition ['drift:vertical_advection', 'is', False]
13:17:30 INFO    opendrift:1821: Storing previous values of element property lon because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
13:17:30 INFO    opendrift:1821: Storing previous values of element property lat because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
13:17:30 INFO    opendrift:952: Using existing reader for land_binary_mask to move elements to ocean
13:17:30 INFO    opendrift:982: All points are in ocean
13:17:30 INFO    opendrift.models.openoil.openoil:692: Oil-water surface tension is 0.035935 Nm
13:17:30 INFO    opendrift.models.openoil.openoil:705: Max water fraction not available for GENERIC BUNKER C, using default
13:17:30 INFO    opendrift:2118: 2026-01-21 13:17:26.340031 - step 1 of 96 - 3000 active elements (0 deactivated)
13:17:30 INFO    opendrift.readers:63: Opening file with xr.open_dataset
13:17:32 INFO    opendrift.readers.reader_netCDF_CF_generic:332: Detected dimensions: {'x': 'X', 'y': 'Y', 'z': 'depth', 'time': 'time'}
13:17:32 INFO    opendrift.readers.basereader:176: Variable x_sea_water_velocity will be rotated from eastward_sea_water_velocity
13:17:32 INFO    opendrift.readers.basereader:176: Variable y_sea_water_velocity will be rotated from northward_sea_water_velocity
13:17:34 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.002850502 to -0.0004373464) found for ocean_vertical_diffusivity, replacing with NaN
13:17:34 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:17:36 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.004022823 to -0.00078050414) found for ocean_vertical_diffusivity, replacing with NaN
13:17:36 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:17:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:36 INFO    opendrift:2118: 2026-01-21 13:47:26.340031 - step 2 of 96 - 3000 active elements (0 deactivated)
13:17:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:36 INFO    opendrift:2118: 2026-01-21 14:17:26.340031 - step 3 of 96 - 3000 active elements (0 deactivated)
13:17:38 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.005719042 to -0.0010995766) found for ocean_vertical_diffusivity, replacing with NaN
13:17:38 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:17:38 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:38 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:38 INFO    opendrift:2118: 2026-01-21 14:47:26.340031 - step 4 of 96 - 3000 active elements (0 deactivated)
13:17:38 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:38 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:38 INFO    opendrift:2118: 2026-01-21 15:17:26.340031 - step 5 of 96 - 3000 active elements (0 deactivated)
13:17:40 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0028588956 to -0.00050402863) found for ocean_vertical_diffusivity, replacing with NaN
13:17:40 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:17:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:41 INFO    opendrift:2118: 2026-01-21 15:47:26.340031 - step 6 of 96 - 3000 active elements (0 deactivated)
13:17:41 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:41 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:41 INFO    opendrift:2118: 2026-01-21 16:17:26.340031 - step 7 of 96 - 3000 active elements (0 deactivated)
13:17:43 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.00044311793 to -1.094205e-06) found for ocean_vertical_diffusivity, replacing with NaN
13:17:43 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:17:43 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:43 INFO    opendrift:2118: 2026-01-21 16:47:26.340031 - step 8 of 96 - 3000 active elements (0 deactivated)
13:17:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:46 INFO    opendrift:2118: 2026-01-21 17:17:26.340031 - step 9 of 96 - 2996 active elements (4 deactivated)
13:17:48 INFO    opendrift:2118: 2026-01-21 17:47:26.340031 - step 10 of 96 - 2947 active elements (53 deactivated)
13:17:49 INFO    opendrift:2118: 2026-01-21 18:17:26.340031 - step 11 of 96 - 2723 active elements (277 deactivated)
13:17:51 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0029093688 to -4.0583836e-05) found for ocean_vertical_diffusivity, replacing with NaN
13:17:51 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:17:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:51 INFO    opendrift:2118: 2026-01-21 18:47:26.340031 - step 12 of 96 - 2355 active elements (645 deactivated)
13:17:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:51 INFO    opendrift:2118: 2026-01-21 19:17:26.340031 - step 13 of 96 - 2018 active elements (982 deactivated)
13:17:53 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.007594841 to -0.0004912115) found for ocean_vertical_diffusivity, replacing with NaN
13:17:53 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:17:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:53 INFO    opendrift:2118: 2026-01-21 19:47:26.340031 - step 14 of 96 - 1859 active elements (1141 deactivated)
13:17:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:54 INFO    opendrift:2118: 2026-01-21 20:17:26.340031 - step 15 of 96 - 1796 active elements (1204 deactivated)
13:17:56 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.010429089 to -0.0015031734) found for ocean_vertical_diffusivity, replacing with NaN
13:17:56 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:17:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:56 INFO    opendrift:2118: 2026-01-21 20:47:26.340031 - step 16 of 96 - 1753 active elements (1247 deactivated)
13:17:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:56 INFO    opendrift:2118: 2026-01-21 21:17:26.340031 - step 17 of 96 - 1712 active elements (1288 deactivated)
13:17:58 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.009840391 to -0.0015006486) found for ocean_vertical_diffusivity, replacing with NaN
13:17:58 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:17:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:58 INFO    opendrift:2118: 2026-01-21 21:47:26.340031 - step 18 of 96 - 1696 active elements (1304 deactivated)
13:17:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:17:58 INFO    opendrift:2118: 2026-01-21 22:17:26.340031 - step 19 of 96 - 1691 active elements (1309 deactivated)
13:18:00 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0072454163 to -0.0012831235) found for ocean_vertical_diffusivity, replacing with NaN
13:18:00 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:01 INFO    opendrift:2118: 2026-01-21 22:47:26.340031 - step 20 of 96 - 1691 active elements (1309 deactivated)
13:18:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:01 INFO    opendrift:2118: 2026-01-21 23:17:26.340031 - step 21 of 96 - 1691 active elements (1309 deactivated)
13:18:03 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.006804044 to -0.0010382265) found for ocean_vertical_diffusivity, replacing with NaN
13:18:03 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:03 INFO    opendrift:2118: 2026-01-21 23:47:26.340031 - step 22 of 96 - 1691 active elements (1309 deactivated)
13:18:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:03 INFO    opendrift:2118: 2026-01-22 00:17:26.340031 - step 23 of 96 - 1691 active elements (1309 deactivated)
13:18:05 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0064401897 to -0.0019210458) found for ocean_vertical_diffusivity, replacing with NaN
13:18:05 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:05 INFO    opendrift:2118: 2026-01-22 00:47:26.340031 - step 24 of 96 - 1691 active elements (1309 deactivated)
13:18:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:05 INFO    opendrift:2118: 2026-01-22 01:17:26.340031 - step 25 of 96 - 1691 active elements (1309 deactivated)
13:18:08 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0073186895 to -0.0016505852) found for ocean_vertical_diffusivity, replacing with NaN
13:18:08 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:08 INFO    opendrift:2118: 2026-01-22 01:47:26.340031 - step 26 of 96 - 1691 active elements (1309 deactivated)
13:18:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:08 INFO    opendrift:2118: 2026-01-22 02:17:26.340031 - step 27 of 96 - 1691 active elements (1309 deactivated)
13:18:10 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0069580483 to -0.0013494838) found for ocean_vertical_diffusivity, replacing with NaN
13:18:10 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:10 INFO    opendrift:2118: 2026-01-22 02:47:26.340031 - step 28 of 96 - 1691 active elements (1309 deactivated)
13:18:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:10 INFO    opendrift:2118: 2026-01-22 03:17:26.340031 - step 29 of 96 - 1691 active elements (1309 deactivated)
13:18:12 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0049330816 to -0.0010556884) found for ocean_vertical_diffusivity, replacing with NaN
13:18:12 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:12 INFO    opendrift:2118: 2026-01-22 03:47:26.340031 - step 30 of 96 - 1691 active elements (1309 deactivated)
13:18:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:12 INFO    opendrift:2118: 2026-01-22 04:17:26.340031 - step 31 of 96 - 1691 active elements (1309 deactivated)
13:18:15 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0018669709 to -0.00038584156) found for ocean_vertical_diffusivity, replacing with NaN
13:18:15 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:15 INFO    opendrift:2118: 2026-01-22 04:47:26.340031 - step 32 of 96 - 1691 active elements (1309 deactivated)
13:18:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:15 INFO    opendrift:2118: 2026-01-22 05:17:26.340031 - step 33 of 96 - 1691 active elements (1309 deactivated)
13:18:17 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.00095122436 to -1.3733764e-05) found for ocean_vertical_diffusivity, replacing with NaN
13:18:17 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:17 INFO    opendrift:2118: 2026-01-22 05:47:26.340031 - step 34 of 96 - 1691 active elements (1309 deactivated)
13:18:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:17 INFO    opendrift:2118: 2026-01-22 06:17:26.340031 - step 35 of 96 - 1691 active elements (1309 deactivated)
13:18:19 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0027822782 to -0.0005362409) found for ocean_vertical_diffusivity, replacing with NaN
13:18:19 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:19 INFO    opendrift:2118: 2026-01-22 06:47:26.340031 - step 36 of 96 - 1691 active elements (1309 deactivated)
13:18:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:19 INFO    opendrift:2118: 2026-01-22 07:17:26.340031 - step 37 of 96 - 1691 active elements (1309 deactivated)
13:18:21 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.006004547 to -0.0010760586) found for ocean_vertical_diffusivity, replacing with NaN
13:18:21 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:21 INFO    opendrift:2118: 2026-01-22 07:47:26.340031 - step 38 of 96 - 1691 active elements (1309 deactivated)
13:18:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:21 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:22 INFO    opendrift:2118: 2026-01-22 08:17:26.340031 - step 39 of 96 - 1691 active elements (1309 deactivated)
13:18:24 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0070300503 to -0.0013441441) found for ocean_vertical_diffusivity, replacing with NaN
13:18:24 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:24 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:24 INFO    opendrift:2118: 2026-01-22 08:47:26.340031 - step 40 of 96 - 1691 active elements (1309 deactivated)
13:18:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:24 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:24 INFO    opendrift:2118: 2026-01-22 09:17:26.340031 - step 41 of 96 - 1691 active elements (1309 deactivated)
13:18:26 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.007115611 to -0.0015678328) found for ocean_vertical_diffusivity, replacing with NaN
13:18:26 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:26 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:26 INFO    opendrift:2118: 2026-01-22 09:47:26.340031 - step 42 of 96 - 1691 active elements (1309 deactivated)
13:18:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:26 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:26 INFO    opendrift:2118: 2026-01-22 10:17:26.340031 - step 43 of 96 - 1691 active elements (1309 deactivated)
13:18:28 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.00843652 to -0.0017494025) found for ocean_vertical_diffusivity, replacing with NaN
13:18:28 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:28 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:28 INFO    opendrift:2118: 2026-01-22 10:47:26.340031 - step 44 of 96 - 1691 active elements (1309 deactivated)
13:18:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:28 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:28 INFO    opendrift:2118: 2026-01-22 11:17:26.340031 - step 45 of 96 - 1691 active elements (1309 deactivated)
13:18:30 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0060656657 to -0.0011639112) found for ocean_vertical_diffusivity, replacing with NaN
13:18:30 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:31 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:31 INFO    opendrift:2118: 2026-01-22 11:47:26.340031 - step 46 of 96 - 1691 active elements (1309 deactivated)
13:18:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:31 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:31 INFO    opendrift:2118: 2026-01-22 12:17:26.340031 - step 47 of 96 - 1691 active elements (1309 deactivated)
13:18:33 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0038390951 to -0.0006732209) found for ocean_vertical_diffusivity, replacing with NaN
13:18:33 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:33 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:33 INFO    opendrift:2118: 2026-01-22 12:47:26.340031 - step 48 of 96 - 1691 active elements (1309 deactivated)
13:18:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:33 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:33 INFO    opendrift:2118: 2026-01-22 13:17:26.340031 - step 49 of 96 - 1691 active elements (1309 deactivated)
13:18:35 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.004030512 to -0.00039515068) found for ocean_vertical_diffusivity, replacing with NaN
13:18:35 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:35 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:35 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:35 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:35 INFO    opendrift:2118: 2026-01-22 13:47:26.340031 - step 50 of 96 - 1691 active elements (1309 deactivated)
13:18:35 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:35 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:35 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:35 INFO    opendrift:2118: 2026-01-22 14:17:26.340031 - step 51 of 96 - 1691 active elements (1309 deactivated)
13:18:37 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.004781054 to -0.00055300415) found for ocean_vertical_diffusivity, replacing with NaN
13:18:37 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:38 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:38 INFO    opendrift:2118: 2026-01-22 14:47:26.340031 - step 52 of 96 - 1691 active elements (1309 deactivated)
13:18:38 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:38 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:38 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:38 INFO    opendrift:2118: 2026-01-22 15:17:26.340031 - step 53 of 96 - 1691 active elements (1309 deactivated)
13:18:40 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0039531523 to -0.00054441835) found for ocean_vertical_diffusivity, replacing with NaN
13:18:40 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:40 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:40 INFO    opendrift:2118: 2026-01-22 15:47:26.340031 - step 54 of 96 - 1691 active elements (1309 deactivated)
13:18:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:40 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:40 INFO    opendrift:2118: 2026-01-22 16:17:26.340031 - step 55 of 96 - 1691 active elements (1309 deactivated)
13:18:42 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0013548302 to -3.1524987e-06) found for ocean_vertical_diffusivity, replacing with NaN
13:18:42 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:42 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:42 INFO    opendrift:2118: 2026-01-22 16:47:26.340031 - step 56 of 96 - 1691 active elements (1309 deactivated)
13:18:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:42 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:42 INFO    opendrift:2118: 2026-01-22 17:17:26.340031 - step 57 of 96 - 1691 active elements (1309 deactivated)
13:18:44 WARNING opendrift.readers.basereader.variables:656: Invalid values (-8.0494436e-05 to -8.0494436e-05) found for ocean_vertical_diffusivity, replacing with NaN
13:18:44 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:44 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:44 INFO    opendrift:2118: 2026-01-22 17:47:26.340031 - step 58 of 96 - 1691 active elements (1309 deactivated)
13:18:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:44 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:44 INFO    opendrift:2118: 2026-01-22 18:17:26.340031 - step 59 of 96 - 1691 active elements (1309 deactivated)
13:18:47 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0018696673 to -2.5892605e-06) found for ocean_vertical_diffusivity, replacing with NaN
13:18:47 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:47 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:47 INFO    opendrift:2118: 2026-01-22 18:47:26.340031 - step 60 of 96 - 1690 active elements (1310 deactivated)
13:18:47 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:47 INFO    opendrift:2118: 2026-01-22 19:17:26.340031 - step 61 of 96 - 1690 active elements (1310 deactivated)
13:18:49 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.005528629 to -0.0007267007) found for ocean_vertical_diffusivity, replacing with NaN
13:18:49 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:49 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:49 INFO    opendrift:2118: 2026-01-22 19:47:26.340031 - step 62 of 96 - 1690 active elements (1310 deactivated)
13:18:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:49 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:49 INFO    opendrift:2118: 2026-01-22 20:17:26.340031 - step 63 of 96 - 1690 active elements (1310 deactivated)
13:18:51 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.010215219 to -0.0013435807) found for ocean_vertical_diffusivity, replacing with NaN
13:18:51 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:51 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:51 INFO    opendrift:2118: 2026-01-22 20:47:26.340031 - step 64 of 96 - 1690 active elements (1310 deactivated)
13:18:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:52 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:52 INFO    opendrift:2118: 2026-01-22 21:17:26.340031 - step 65 of 96 - 1690 active elements (1310 deactivated)
13:18:54 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0134143345 to -0.0015085074) found for ocean_vertical_diffusivity, replacing with NaN
13:18:54 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:54 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:54 INFO    opendrift:2118: 2026-01-22 21:47:26.340031 - step 66 of 96 - 1690 active elements (1310 deactivated)
13:18:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:54 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:54 INFO    opendrift:2118: 2026-01-22 22:17:26.340031 - step 67 of 96 - 1690 active elements (1310 deactivated)
13:18:56 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.012812783 to -0.0011693143) found for ocean_vertical_diffusivity, replacing with NaN
13:18:56 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:56 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:56 INFO    opendrift:2118: 2026-01-22 22:47:26.340031 - step 68 of 96 - 1690 active elements (1310 deactivated)
13:18:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:56 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:56 INFO    opendrift:2118: 2026-01-22 23:17:26.340031 - step 69 of 96 - 1690 active elements (1310 deactivated)
13:18:58 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0103679 to -0.001030744) found for ocean_vertical_diffusivity, replacing with NaN
13:18:58 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:18:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:58 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:58 INFO    opendrift:2118: 2026-01-22 23:47:26.340031 - step 70 of 96 - 1690 active elements (1310 deactivated)
13:18:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:58 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:18:58 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:18:58 INFO    opendrift:2118: 2026-01-23 00:17:26.340031 - step 71 of 96 - 1690 active elements (1310 deactivated)
13:19:00 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0075799725 to -0.00093490136) found for ocean_vertical_diffusivity, replacing with NaN
13:19:00 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:19:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:01 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:01 INFO    opendrift:2118: 2026-01-23 00:47:26.340031 - step 72 of 96 - 1690 active elements (1310 deactivated)
13:19:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:01 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:01 INFO    opendrift:2118: 2026-01-23 01:17:26.340031 - step 73 of 96 - 1690 active elements (1310 deactivated)
13:19:03 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.004658034 to -0.0005046062) found for ocean_vertical_diffusivity, replacing with NaN
13:19:03 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:19:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:03 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:03 INFO    opendrift:2118: 2026-01-23 01:47:26.340031 - step 74 of 96 - 1690 active elements (1310 deactivated)
13:19:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:03 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:03 INFO    opendrift:2118: 2026-01-23 02:17:26.340031 - step 75 of 96 - 1690 active elements (1310 deactivated)
13:19:05 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.004039051 to -0.00071729056) found for ocean_vertical_diffusivity, replacing with NaN
13:19:05 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:19:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:05 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:05 INFO    opendrift:2118: 2026-01-23 02:47:26.340031 - step 76 of 96 - 1690 active elements (1310 deactivated)
13:19:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:05 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:05 INFO    opendrift:2118: 2026-01-23 03:17:26.340031 - step 77 of 96 - 1688 active elements (1312 deactivated)
13:19:07 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0056577623 to -0.0005264875) found for ocean_vertical_diffusivity, replacing with NaN
13:19:07 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:19:07 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:07 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:07 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:07 INFO    opendrift:2118: 2026-01-23 03:47:26.340031 - step 78 of 96 - 1688 active elements (1312 deactivated)
13:19:07 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:07 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:08 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:08 INFO    opendrift:2118: 2026-01-23 04:17:26.340031 - step 79 of 96 - 1688 active elements (1312 deactivated)
13:19:10 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0043276018 to -0.0003437034) found for ocean_vertical_diffusivity, replacing with NaN
13:19:10 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:19:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:10 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:10 INFO    opendrift:2118: 2026-01-23 04:47:26.340031 - step 80 of 96 - 1687 active elements (1313 deactivated)
13:19:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:10 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:10 INFO    opendrift:2118: 2026-01-23 05:17:26.340031 - step 81 of 96 - 1686 active elements (1314 deactivated)
13:19:12 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0022363525 to -0.00018481672) found for ocean_vertical_diffusivity, replacing with NaN
13:19:12 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:19:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:12 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:12 INFO    opendrift:2118: 2026-01-23 05:47:26.340031 - step 82 of 96 - 1686 active elements (1314 deactivated)
13:19:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:12 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:12 INFO    opendrift:2118: 2026-01-23 06:17:26.340031 - step 83 of 96 - 1686 active elements (1314 deactivated)
13:19:14 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0062386985 to -0.00048563376) found for ocean_vertical_diffusivity, replacing with NaN
13:19:14 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:19:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:15 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:15 INFO    opendrift:2118: 2026-01-23 06:47:26.340031 - step 84 of 96 - 1686 active elements (1314 deactivated)
13:19:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:15 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:15 INFO    opendrift:2118: 2026-01-23 07:17:26.340031 - step 85 of 96 - 1686 active elements (1314 deactivated)
13:19:17 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.012500524 to -0.0018049232) found for ocean_vertical_diffusivity, replacing with NaN
13:19:17 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:19:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:17 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:17 INFO    opendrift:2118: 2026-01-23 07:47:26.340031 - step 86 of 96 - 1686 active elements (1314 deactivated)
13:19:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:17 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:17 INFO    opendrift:2118: 2026-01-23 08:17:26.340031 - step 87 of 96 - 1686 active elements (1314 deactivated)
13:19:19 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.020938652 to -0.0015545163) found for ocean_vertical_diffusivity, replacing with NaN
13:19:19 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:19:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:19 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:19 INFO    opendrift:2118: 2026-01-23 08:47:26.340031 - step 88 of 96 - 1686 active elements (1314 deactivated)
13:19:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:19 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:19 INFO    opendrift:2118: 2026-01-23 09:17:26.340031 - step 89 of 96 - 1685 active elements (1315 deactivated)
13:19:21 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.032347005 to -0.001827848) found for ocean_vertical_diffusivity, replacing with NaN
13:19:21 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:19:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:21 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:21 INFO    opendrift:2118: 2026-01-23 09:47:26.340031 - step 90 of 96 - 1685 active elements (1315 deactivated)
13:19:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:21 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:21 INFO    opendrift:2118: 2026-01-23 10:17:26.340031 - step 91 of 96 - 1685 active elements (1315 deactivated)
13:19:23 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.027494216 to -0.001056209) found for ocean_vertical_diffusivity, replacing with NaN
13:19:23 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:19:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:24 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:24 INFO    opendrift:2118: 2026-01-23 10:47:26.340031 - step 92 of 96 - 1682 active elements (1318 deactivated)
13:19:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:24 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:24 INFO    opendrift:2118: 2026-01-23 11:17:26.340031 - step 93 of 96 - 1680 active elements (1320 deactivated)
13:19:26 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.027207397 to -0.0019502835) found for ocean_vertical_diffusivity, replacing with NaN
13:19:26 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:19:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:26 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:26 INFO    opendrift:2118: 2026-01-23 11:47:26.340031 - step 94 of 96 - 1672 active elements (1328 deactivated)
13:19:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:26 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:26 INFO    opendrift:2118: 2026-01-23 12:17:26.340031 - step 95 of 96 - 1659 active elements (1341 deactivated)
13:19:28 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.024707904 to -0.0017319625) found for ocean_vertical_diffusivity, replacing with NaN
13:19:28 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
13:19:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:28 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
13:19:28 INFO    opendrift:2118: 2026-01-23 12:47:26.340031 - step 96 of 96 - 1641 active elements (1359 deactivated)
13:19:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
13:19:28 INFO    opendrift.models.openoil.openoil:1189: Ice concentration above 30%, using Nordam scheme for advection in ice
<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:                                            ...


Print and plot results

print(o)
o.animation(background='sea_ice_area_fraction', cmap=cmocean.cm.ice,
            vmin=0, vmax=1, bgalpha=1, fast=False)
===========================
--------------------
Reader performance:
--------------------
global_landmask
 0:00:01.5  total
 0:00:00.0  preparing
 0:00:01.5  reading
 0:00:00.0  masking
--------------------
https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
 0:01:44.7  total
 0:00:00.0  preparing
 0:01:43.9  reading
 0:00:01.1  interpolation
 0:00:00.1  interpolation_time
 0:00:00.7  rotating vectors
 0:00:00.0  masking
--------------------
Performance:
 2:02.5 total time
    4.3 configuration
    0.0 preparing main loop
      0.0 moving elements to ocean
 1:58.2 main loop
      4.9 updating elements
        0.0 oil weathering
          0.0 updating viscosities
          0.0 updating densities
        4.5 vertical mixing
    0.0 cleaning up
--------------------
===========================
Model:  OpenOil     (OpenDrift version 1.14.8)
        1571 active Oil particles  (1429 deactivated, 0 scheduled)
-------------------
Environment variables:
  -----
  land_binary_mask
     1) global_landmask
  -----
  ocean_vertical_diffusivity
  sea_floor_depth_below_sea_level
  sea_ice_area_fraction
  sea_ice_x_velocity
  sea_ice_y_velocity
  sea_surface_height
  sea_water_salinity
  sea_water_temperature
  x_sea_water_velocity
  x_wind
  y_sea_water_velocity
  y_wind
     1) https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
  -----
Readers not added for the following variables:
  ocean_mixed_layer_thickness
  sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment
  sea_surface_wave_period_at_variance_spectral_density_maximum
  sea_surface_wave_significant_height
  sea_surface_wave_stokes_drift_x_velocity
  sea_surface_wave_stokes_drift_y_velocity

Discarded readers:

Time:
        Start: 2026-01-21 13:17:26.340031 UTC
        Present: 2026-01-23 13:17:26.340031 UTC
        Calculation steps: 96 * 0:30:00 - total time: 2 days, 0:00:00
        Output steps: 49 * 1:00:00
===========================

13:19:40 INFO    opendrift:4675: Saving animation to /root/project/docs/source/gallery/animations/example_oil_ice_0.gif...
13:20:25 INFO    opendrift:3101: Time to make animation: 0:00:56.070799
../_images/example_oil_ice_0.gif
o.plot(background='sea_ice_area_fraction', cmap=cmocean.cm.ice,
       vmin=0, vmax=1, bgalpha=1, fast=False)
OpenDrift - OpenOil (GENERIC BUNKER C) 2026-01-21 13:17 to 2026-01-23 13:17 UTC (49 steps)
(<GeoAxes: title={'center': 'OpenDrift - OpenOil (GENERIC BUNKER C)\n2026-01-21 13:17 to 2026-01-23 13:17 UTC (49 steps)'}>, <Figure size 683.421x1100 with 2 Axes>)

Total running time of the script: (3 minutes 15.348 seconds)

Gallery generated by Sphinx-Gallery