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)

17:14:52 INFO    opendrift:513: OpenDriftSimulation initialised (version 1.14.2 / v1.14.2-88-gf09b7bd)

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))

17:14:52 INFO    opendrift.models.openoil.openoil:1708: Oil type not specified, using default: GENERIC BUNKER C
17:14:52 INFO    opendrift.models.openoil.adios.dirjs:86: Querying ADIOS database for oil: GENERIC BUNKER C
17:14:52 INFO    opendrift.models.openoil.openoil:1717: Using density 988.1 and viscosity 0.02169233387797564 of oiltype GENERIC BUNKER C
17:14:52 INFO    opendrift.models.basemodel.environment:206: Adding a global landmask from GSHHG
17:14:57 INFO    opendrift.models.basemodel.environment:229: Fallback values will be used for the following variables which have no readers:
17:14:57 INFO    opendrift.models.basemodel.environment:232:    sea_surface_height: 0.000000
17:14:57 INFO    opendrift.models.basemodel.environment:232:    upward_sea_water_velocity: 0.000000
17:14:57 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_significant_height: 0.000000
17:14:57 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
17:14:57 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
17:14:57 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
17:14:57 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
17:14:57 INFO    opendrift.models.basemodel.environment:232:    sea_ice_area_fraction: 0.000000
17:14:57 INFO    opendrift.models.basemodel.environment:232:    sea_ice_x_velocity: 0.000000
17:14:57 INFO    opendrift.models.basemodel.environment:232:    sea_ice_y_velocity: 0.000000
17:14:57 INFO    opendrift.models.basemodel.environment:232:    sea_water_temperature: 10.000000
17:14:57 INFO    opendrift.models.basemodel.environment:232:    sea_water_salinity: 34.000000
17:14:57 INFO    opendrift.models.basemodel.environment:232:    sea_floor_depth_below_sea_level: 10000.000000
17:14:57 INFO    opendrift.models.basemodel.environment:232:    ocean_vertical_diffusivity: 0.020000
17:14:57 INFO    opendrift.models.basemodel.environment:232:    ocean_mixed_layer_thickness: 50.000000

Running model

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

17:14:57 INFO    opendrift:1732: Skipping environment variable upward_sea_water_velocity because of condition ['drift:vertical_advection', 'is', False]
17:14:57 INFO    opendrift:1743: Storing previous values of element property lon because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
17:14:57 INFO    opendrift:1743: Storing previous values of element property lat because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
17:14:57 INFO    opendrift:899: Using existing reader for land_binary_mask
17:14:57 INFO    opendrift:928: All points are in ocean
17:14:57 INFO    opendrift.models.openoil.openoil:692: Oil-water surface tension is 0.035935 Nm
17:14:57 INFO    opendrift.models.openoil.openoil:705: Max water fraction not available for GENERIC BUNKER C, using default
17:14:57 INFO    opendrift:2035: 2025-07-04 17:14:52.964775 - step 1 of 96 - 3000 active elements (0 deactivated)
17:14:57 INFO    opendrift.readers:61: Opening file with xr.open_dataset
17:14:58 INFO    opendrift.readers.reader_netCDF_CF_generic:332: Detected dimensions: {'x': 'X', 'y': 'Y', 'z': 'depth', 'time': 'time'}
17:14:58 INFO    opendrift.readers.basereader:176: Variable x_sea_water_velocity will be rotated from eastward_sea_water_velocity
17:14:58 INFO    opendrift.readers.basereader:176: Variable y_sea_water_velocity will be rotated from northward_sea_water_velocity
17:15:01 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.004161266 to -1.327519e-07) found for ocean_vertical_diffusivity, replacing with NaN
17:15:01 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:04 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0037316813 to -1.7585774e-07) found for ocean_vertical_diffusivity, replacing with NaN
17:15:04 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:04 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:04 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:04 INFO    opendrift:2035: 2025-07-04 17:44:52.964775 - step 2 of 96 - 3000 active elements (0 deactivated)
17:15:04 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:04 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:04 INFO    opendrift:2035: 2025-07-04 18:14:52.964775 - step 3 of 96 - 3000 active elements (0 deactivated)
17:15:07 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0030406534 to -7.926701e-08) found for ocean_vertical_diffusivity, replacing with NaN
17:15:07 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:07 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:07 INFO    opendrift:2035: 2025-07-04 18:44:52.964775 - step 4 of 96 - 3000 active elements (0 deactivated)
17:15:07 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:07 INFO    opendrift:2035: 2025-07-04 19:14:52.964775 - step 5 of 96 - 3000 active elements (0 deactivated)
17:15:10 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0019254447 to -1.09461835e-05) found for ocean_vertical_diffusivity, replacing with NaN
17:15:10 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:10 INFO    opendrift:2035: 2025-07-04 19:44:52.964775 - step 6 of 96 - 3000 active elements (0 deactivated)
17:15:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:10 INFO    opendrift:2035: 2025-07-04 20:14:52.964775 - step 7 of 96 - 3000 active elements (0 deactivated)
17:15:12 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0018324374 to -2.9753824e-05) found for ocean_vertical_diffusivity, replacing with NaN
17:15:12 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:12 INFO    opendrift:2035: 2025-07-04 20:44:52.964775 - step 8 of 96 - 3000 active elements (0 deactivated)
17:15:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:13 INFO    opendrift:2035: 2025-07-04 21:14:52.964775 - step 9 of 96 - 3000 active elements (0 deactivated)
17:15:15 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0022738765 to -6.0633494e-05) found for ocean_vertical_diffusivity, replacing with NaN
17:15:15 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:15 INFO    opendrift:2035: 2025-07-04 21:44:52.964775 - step 10 of 96 - 3000 active elements (0 deactivated)
17:15:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:15 INFO    opendrift:2035: 2025-07-04 22:14:52.964775 - step 11 of 96 - 3000 active elements (0 deactivated)
17:15:18 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0024067229 to -0.00017851584) found for ocean_vertical_diffusivity, replacing with NaN
17:15:18 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:18 INFO    opendrift:2035: 2025-07-04 22:44:52.964775 - step 12 of 96 - 3000 active elements (0 deactivated)
17:15:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:18 INFO    opendrift:2035: 2025-07-04 23:14:52.964775 - step 13 of 96 - 3000 active elements (0 deactivated)
17:15:22 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0033950808 to -0.00018556094) found for ocean_vertical_diffusivity, replacing with NaN
17:15:22 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:22 INFO    opendrift:2035: 2025-07-04 23:44:52.964775 - step 14 of 96 - 3000 active elements (0 deactivated)
17:15:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:22 INFO    opendrift:2035: 2025-07-05 00:14:52.964775 - step 15 of 96 - 3000 active elements (0 deactivated)
17:15:26 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.00418552 to -0.00012294737) found for ocean_vertical_diffusivity, replacing with NaN
17:15:26 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:28 INFO    opendrift:2035: 2025-07-05 00:44:52.964775 - step 16 of 96 - 2999 active elements (1 deactivated)
17:15:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:28 INFO    opendrift:2035: 2025-07-05 01:14:52.964775 - step 17 of 96 - 2999 active elements (1 deactivated)
17:15:30 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.004999755 to -0.00020801119) found for ocean_vertical_diffusivity, replacing with NaN
17:15:30 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:31 INFO    opendrift:2035: 2025-07-05 01:44:52.964775 - step 18 of 96 - 2999 active elements (1 deactivated)
17:15:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:31 INFO    opendrift:2035: 2025-07-05 02:14:52.964775 - step 19 of 96 - 2998 active elements (2 deactivated)
17:15:33 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.006500769 to -0.00017876513) found for ocean_vertical_diffusivity, replacing with NaN
17:15:33 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:33 INFO    opendrift:2035: 2025-07-05 02:44:52.964775 - step 20 of 96 - 2998 active elements (2 deactivated)
17:15:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:34 INFO    opendrift:2035: 2025-07-05 03:14:52.964775 - step 21 of 96 - 2998 active elements (2 deactivated)
17:15:36 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.00871235 to -3.478716e-05) found for ocean_vertical_diffusivity, replacing with NaN
17:15:36 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:36 INFO    opendrift:2035: 2025-07-05 03:44:52.964775 - step 22 of 96 - 2997 active elements (3 deactivated)
17:15:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:36 INFO    opendrift:2035: 2025-07-05 04:14:52.964775 - step 23 of 96 - 2997 active elements (3 deactivated)
17:15:39 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.009381823 to -8.409126e-05) found for ocean_vertical_diffusivity, replacing with NaN
17:15:39 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:39 INFO    opendrift:2035: 2025-07-05 04:44:52.964775 - step 24 of 96 - 2997 active elements (3 deactivated)
17:15:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:39 INFO    opendrift:2035: 2025-07-05 05:14:52.964775 - step 25 of 96 - 2997 active elements (3 deactivated)
17:15:42 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.007937315 to -0.00019404551) found for ocean_vertical_diffusivity, replacing with NaN
17:15:42 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:42 INFO    opendrift:2035: 2025-07-05 05:44:52.964775 - step 26 of 96 - 2997 active elements (3 deactivated)
17:15:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:42 INFO    opendrift:2035: 2025-07-05 06:14:52.964775 - step 27 of 96 - 2997 active elements (3 deactivated)
17:15:44 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0061478377 to -0.00014046332) found for ocean_vertical_diffusivity, replacing with NaN
17:15:44 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:45 INFO    opendrift:2035: 2025-07-05 06:44:52.964775 - step 28 of 96 - 2997 active elements (3 deactivated)
17:15:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:45 INFO    opendrift:2035: 2025-07-05 07:14:52.964775 - step 29 of 96 - 2997 active elements (3 deactivated)
17:15:47 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.004893403 to -0.00019527377) found for ocean_vertical_diffusivity, replacing with NaN
17:15:47 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:47 INFO    opendrift:2035: 2025-07-05 07:44:52.964775 - step 30 of 96 - 2997 active elements (3 deactivated)
17:15:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:48 INFO    opendrift:2035: 2025-07-05 08:14:52.964775 - step 31 of 96 - 2997 active elements (3 deactivated)
17:15:51 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.005155026 to -0.00020800867) found for ocean_vertical_diffusivity, replacing with NaN
17:15:51 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:51 INFO    opendrift:2035: 2025-07-05 08:44:52.964775 - step 32 of 96 - 2997 active elements (3 deactivated)
17:15:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:51 INFO    opendrift:2035: 2025-07-05 09:14:52.964775 - step 33 of 96 - 2997 active elements (3 deactivated)
17:15:53 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0053363782 to -0.0001888528) found for ocean_vertical_diffusivity, replacing with NaN
17:15:53 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:53 INFO    opendrift:2035: 2025-07-05 09:44:52.964775 - step 34 of 96 - 2997 active elements (3 deactivated)
17:15:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:54 INFO    opendrift:2035: 2025-07-05 10:14:52.964775 - step 35 of 96 - 2997 active elements (3 deactivated)
17:15:56 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0070892447 to -0.00012578964) found for ocean_vertical_diffusivity, replacing with NaN
17:15:56 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:56 INFO    opendrift:2035: 2025-07-05 10:44:52.964775 - step 36 of 96 - 2997 active elements (3 deactivated)
17:15:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:57 INFO    opendrift:2035: 2025-07-05 11:14:52.964775 - step 37 of 96 - 2997 active elements (3 deactivated)
17:15:59 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.009526558 to -0.000103016515) found for ocean_vertical_diffusivity, replacing with NaN
17:15:59 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:15:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:59 INFO    opendrift:2035: 2025-07-05 11:44:52.964775 - step 38 of 96 - 2997 active elements (3 deactivated)
17:15:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:15:59 INFO    opendrift:2035: 2025-07-05 12:14:52.964775 - step 39 of 96 - 2997 active elements (3 deactivated)
17:16:02 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.009563377 to -0.00040416577) found for ocean_vertical_diffusivity, replacing with NaN
17:16:02 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:02 INFO    opendrift:2035: 2025-07-05 12:44:52.964775 - step 40 of 96 - 2997 active elements (3 deactivated)
17:16:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:02 INFO    opendrift:2035: 2025-07-05 13:14:52.964775 - step 41 of 96 - 2997 active elements (3 deactivated)
17:16:05 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.010229149 to -0.00035121094) found for ocean_vertical_diffusivity, replacing with NaN
17:16:05 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:05 INFO    opendrift:2035: 2025-07-05 13:44:52.964775 - step 42 of 96 - 2997 active elements (3 deactivated)
17:16:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:05 INFO    opendrift:2035: 2025-07-05 14:14:52.964775 - step 43 of 96 - 2997 active elements (3 deactivated)
17:16:08 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0117224455 to -0.0004119614) found for ocean_vertical_diffusivity, replacing with NaN
17:16:08 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:08 INFO    opendrift:2035: 2025-07-05 14:44:52.964775 - step 44 of 96 - 2997 active elements (3 deactivated)
17:16:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:08 INFO    opendrift:2035: 2025-07-05 15:14:52.964775 - step 45 of 96 - 2997 active elements (3 deactivated)
17:16:10 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.013962276 to -0.00043455375) found for ocean_vertical_diffusivity, replacing with NaN
17:16:10 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:11 INFO    opendrift:2035: 2025-07-05 15:44:52.964775 - step 46 of 96 - 2997 active elements (3 deactivated)
17:16:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:11 INFO    opendrift:2035: 2025-07-05 16:14:52.964775 - step 47 of 96 - 2996 active elements (4 deactivated)
17:16:13 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.016501287 to -0.0005353235) found for ocean_vertical_diffusivity, replacing with NaN
17:16:13 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:14 INFO    opendrift:2035: 2025-07-05 16:44:52.964775 - step 48 of 96 - 2994 active elements (6 deactivated)
17:16:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:14 INFO    opendrift:2035: 2025-07-05 17:14:52.964775 - step 49 of 96 - 2994 active elements (6 deactivated)
17:16:16 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.01634547 to -0.0014164564) found for ocean_vertical_diffusivity, replacing with NaN
17:16:16 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:16 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:16 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:17 INFO    opendrift:2035: 2025-07-05 17:44:52.964775 - step 50 of 96 - 2989 active elements (11 deactivated)
17:16:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:17 INFO    opendrift:2035: 2025-07-05 18:14:52.964775 - step 51 of 96 - 2982 active elements (18 deactivated)
17:16:19 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.013366154 to -0.0011812702) found for ocean_vertical_diffusivity, replacing with NaN
17:16:19 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:19 INFO    opendrift:2035: 2025-07-05 18:44:52.964775 - step 52 of 96 - 2979 active elements (21 deactivated)
17:16:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:20 INFO    opendrift:2035: 2025-07-05 19:14:52.964775 - step 53 of 96 - 2974 active elements (26 deactivated)
17:16:22 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.009992571 to -0.0004895313) found for ocean_vertical_diffusivity, replacing with NaN
17:16:22 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:22 INFO    opendrift:2035: 2025-07-05 19:44:52.964775 - step 54 of 96 - 2968 active elements (32 deactivated)
17:16:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:22 INFO    opendrift:2035: 2025-07-05 20:14:52.964775 - step 55 of 96 - 2961 active elements (39 deactivated)
17:16:25 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.00808953 to -0.00033632896) found for ocean_vertical_diffusivity, replacing with NaN
17:16:25 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:25 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:25 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:25 INFO    opendrift:2035: 2025-07-05 20:44:52.964775 - step 56 of 96 - 2958 active elements (42 deactivated)
17:16:25 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:25 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:25 INFO    opendrift:2035: 2025-07-05 21:14:52.964775 - step 57 of 96 - 2954 active elements (46 deactivated)
17:16:28 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0072371485 to -0.0005341876) found for ocean_vertical_diffusivity, replacing with NaN
17:16:28 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:28 INFO    opendrift:2035: 2025-07-05 21:44:52.964775 - step 58 of 96 - 2949 active elements (51 deactivated)
17:16:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:28 INFO    opendrift:2035: 2025-07-05 22:14:52.964775 - step 59 of 96 - 2945 active elements (55 deactivated)
17:16:30 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0070027933 to -0.00045273127) found for ocean_vertical_diffusivity, replacing with NaN
17:16:30 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:30 INFO    opendrift:2035: 2025-07-05 22:44:52.964775 - step 60 of 96 - 2940 active elements (60 deactivated)
17:16:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:31 INFO    opendrift:2035: 2025-07-05 23:14:52.964775 - step 61 of 96 - 2933 active elements (67 deactivated)
17:16:34 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.00649512 to -0.00059249287) found for ocean_vertical_diffusivity, replacing with NaN
17:16:34 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:34 INFO    opendrift:2035: 2025-07-05 23:44:52.964775 - step 62 of 96 - 2921 active elements (79 deactivated)
17:16:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:34 INFO    opendrift:2035: 2025-07-06 00:14:52.964775 - step 63 of 96 - 2919 active elements (81 deactivated)
17:16:36 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0052668117 to -0.0006070284) found for ocean_vertical_diffusivity, replacing with NaN
17:16:36 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:37 INFO    opendrift:2035: 2025-07-06 00:44:52.964775 - step 64 of 96 - 2911 active elements (89 deactivated)
17:16:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:37 INFO    opendrift:2035: 2025-07-06 01:14:52.964775 - step 65 of 96 - 2900 active elements (100 deactivated)
17:16:39 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0065345853 to -0.00070319173) found for ocean_vertical_diffusivity, replacing with NaN
17:16:39 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:40 INFO    opendrift:2035: 2025-07-06 01:44:52.964775 - step 66 of 96 - 2881 active elements (119 deactivated)
17:16:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:40 INFO    opendrift:2035: 2025-07-06 02:14:52.964775 - step 67 of 96 - 2873 active elements (127 deactivated)
17:16:42 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0075991447 to -0.0007240094) found for ocean_vertical_diffusivity, replacing with NaN
17:16:42 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:43 INFO    opendrift:2035: 2025-07-06 02:44:52.964775 - step 68 of 96 - 2857 active elements (143 deactivated)
17:16:43 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:43 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:43 INFO    opendrift:2035: 2025-07-06 03:14:52.964775 - step 69 of 96 - 2839 active elements (161 deactivated)
17:16:45 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.00937609 to -0.0008662455) found for ocean_vertical_diffusivity, replacing with NaN
17:16:45 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:45 INFO    opendrift:2035: 2025-07-06 03:44:52.964775 - step 70 of 96 - 2814 active elements (186 deactivated)
17:16:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:45 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:46 INFO    opendrift:2035: 2025-07-06 04:14:52.964775 - step 71 of 96 - 2793 active elements (207 deactivated)
17:16:48 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.010860234 to -0.0006185102) found for ocean_vertical_diffusivity, replacing with NaN
17:16:48 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:48 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:48 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:48 INFO    opendrift:2035: 2025-07-06 04:44:52.964775 - step 72 of 96 - 2759 active elements (241 deactivated)
17:16:48 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:48 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:48 INFO    opendrift:2035: 2025-07-06 05:14:52.964775 - step 73 of 96 - 2737 active elements (263 deactivated)
17:16:51 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.010284586 to -0.00049652415) found for ocean_vertical_diffusivity, replacing with NaN
17:16:51 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:51 INFO    opendrift:2035: 2025-07-06 05:44:52.964775 - step 74 of 96 - 2713 active elements (287 deactivated)
17:16:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:51 INFO    opendrift:2035: 2025-07-06 06:14:52.964775 - step 75 of 96 - 2686 active elements (314 deactivated)
17:16:54 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.008839674 to -0.0008419308) found for ocean_vertical_diffusivity, replacing with NaN
17:16:54 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:54 INFO    opendrift:2035: 2025-07-06 06:44:52.964775 - step 76 of 96 - 2657 active elements (343 deactivated)
17:16:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:54 INFO    opendrift:2035: 2025-07-06 07:14:52.964775 - step 77 of 96 - 2624 active elements (376 deactivated)
17:16:57 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.00810782 to -0.000751686) found for ocean_vertical_diffusivity, replacing with NaN
17:16:57 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:16:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:57 INFO    opendrift:2035: 2025-07-06 07:44:52.964775 - step 78 of 96 - 2594 active elements (406 deactivated)
17:16:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:16:57 INFO    opendrift:2035: 2025-07-06 08:14:52.964775 - step 79 of 96 - 2566 active elements (434 deactivated)
17:17:00 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0076518347 to -0.0005042935) found for ocean_vertical_diffusivity, replacing with NaN
17:17:00 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:17:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:00 INFO    opendrift:2035: 2025-07-06 08:44:52.964775 - step 80 of 96 - 2533 active elements (467 deactivated)
17:17:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:00 INFO    opendrift:2035: 2025-07-06 09:14:52.964775 - step 81 of 96 - 2485 active elements (515 deactivated)
17:17:03 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.007651817 to -0.0003211952) found for ocean_vertical_diffusivity, replacing with NaN
17:17:03 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:17:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:03 INFO    opendrift:2035: 2025-07-06 09:44:52.964775 - step 82 of 96 - 2440 active elements (560 deactivated)
17:17:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:03 INFO    opendrift:2035: 2025-07-06 10:14:52.964775 - step 83 of 96 - 2399 active elements (601 deactivated)
17:17:06 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.007349758 to -0.0003017104) found for ocean_vertical_diffusivity, replacing with NaN
17:17:06 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:17:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:06 INFO    opendrift:2035: 2025-07-06 10:44:52.964775 - step 84 of 96 - 2337 active elements (663 deactivated)
17:17:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:06 INFO    opendrift:2035: 2025-07-06 11:14:52.964775 - step 85 of 96 - 2267 active elements (733 deactivated)
17:17:09 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0061893784 to -0.00019114988) found for ocean_vertical_diffusivity, replacing with NaN
17:17:09 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:17:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:09 INFO    opendrift:2035: 2025-07-06 11:44:52.964775 - step 86 of 96 - 2194 active elements (806 deactivated)
17:17:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:09 INFO    opendrift:2035: 2025-07-06 12:14:52.964775 - step 87 of 96 - 2093 active elements (907 deactivated)
17:17:11 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.006584637 to -0.0007401208) found for ocean_vertical_diffusivity, replacing with NaN
17:17:11 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:17:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:12 INFO    opendrift:2035: 2025-07-06 12:44:52.964775 - step 88 of 96 - 1984 active elements (1016 deactivated)
17:17:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:12 INFO    opendrift:2035: 2025-07-06 13:14:52.964775 - step 89 of 96 - 1878 active elements (1122 deactivated)
17:17:14 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.0075032213 to -0.001109632) found for ocean_vertical_diffusivity, replacing with NaN
17:17:14 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:17:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:15 INFO    opendrift:2035: 2025-07-06 13:44:52.964775 - step 90 of 96 - 1730 active elements (1270 deactivated)
17:17:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:15 INFO    opendrift:2035: 2025-07-06 14:14:52.964775 - step 91 of 96 - 1575 active elements (1425 deactivated)
17:17:17 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.008944918 to -0.0013201785) found for ocean_vertical_diffusivity, replacing with NaN
17:17:17 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:17:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:18 INFO    opendrift:2035: 2025-07-06 14:44:52.964775 - step 92 of 96 - 1430 active elements (1570 deactivated)
17:17:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:18 INFO    opendrift:2035: 2025-07-06 15:14:52.964775 - step 93 of 96 - 1280 active elements (1720 deactivated)
17:17:20 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.010414142 to -0.0018792652) found for ocean_vertical_diffusivity, replacing with NaN
17:17:20 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:17:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:20 INFO    opendrift:2035: 2025-07-06 15:44:52.964775 - step 94 of 96 - 1112 active elements (1888 deactivated)
17:17:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:21 INFO    opendrift:2035: 2025-07-06 16:14:52.964775 - step 95 of 96 - 997 active elements (2003 deactivated)
17:17:23 WARNING opendrift.readers.basereader.variables:656: Invalid values (-0.012742383 to -0.0012981863) found for ocean_vertical_diffusivity, replacing with NaN
17:17:23 WARNING opendrift.readers.basereader.variables:659: (allowed range: [0, 1])
17:17:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:23 INFO    opendrift:2035: 2025-07-06 16:44:52.964775 - step 96 of 96 - 872 active elements (2128 deactivated)
17:17:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
17:17:23 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-07-04T17:14:52.964775 ... 2...

    standard_name :

    time

    long_name :

    time

    dtype :

    <class 'numpy.float64'>

    array(['2025-07-04T17:14:52.964775000', '2025-07-04T18:14:52.964775000',
           '2025-07-04T19:14:52.964775000', '2025-07-04T20:14:52.964775000',
           '2025-07-04T21:14:52.964775000', '2025-07-04T22:14:52.964775000',
           '2025-07-04T23:14:52.964775000', '2025-07-05T00:14:52.964775000',
           '2025-07-05T01:14:52.964775000', '2025-07-05T02:14:52.964775000',
           '2025-07-05T03:14:52.964775000', '2025-07-05T04:14:52.964775000',
           '2025-07-05T05:14:52.964775000', '2025-07-05T06:14:52.964775000',
           '2025-07-05T07:14:52.964775000', '2025-07-05T08:14:52.964775000',
           '2025-07-05T09:14:52.964775000', '2025-07-05T10:14:52.964775000',
           '2025-07-05T11:14:52.964775000', '2025-07-05T12:14:52.964775000',
           '2025-07-05T13:14:52.964775000', '2025-07-05T14:14:52.964775000',
           '2025-07-05T15:14:52.964775000', '2025-07-05T16:14:52.964775000',
           '2025-07-05T17:14:52.964775000', '2025-07-05T18:14:52.964775000',
           '2025-07-05T19:14:52.964775000', '2025-07-05T20:14:52.964775000',
           '2025-07-05T21:14:52.964775000', '2025-07-05T22:14:52.964775000',
           '2025-07-05T23:14:52.964775000', '2025-07-06T00:14:52.964775000',
           '2025-07-06T01:14:52.964775000', '2025-07-06T02:14:52.964775000',
           '2025-07-06T03:14:52.964775000', '2025-07-06T04:14:52.964775000',
           '2025-07-06T05:14:52.964775000', '2025-07-06T06:14:52.964775000',
           '2025-07-06T07:14:52.964775000', '2025-07-06T08:14:52.964775000',
           '2025-07-06T09:14:52.964775000', '2025-07-06T10:14:52.964775000',
           '2025-07-06T11:14:52.964775000', '2025-07-06T12:14:52.964775000',
           '2025-07-06T13:14:52.964775000', '2025-07-06T14:14:52.964775000',
           '2025-07-06T15:14:52.964775000', '2025-07-06T16:14:52.964775000',
           '2025-07-06T17:14:52.964775000'], 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., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  1.],
           [ 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., ..., nan, nan, nan],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           ...,
           [ 1.,  1.,  1., ...,  1.,  1.,  0.],
           [ 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., ...,     nan,     nan,     nan],
           [     0.,   3600.,   7200., ..., 165600., 169200., 172800.],
           ...,
           [     0.,   3600.,   7200., ..., 165600., 169200., 171000.],
           [     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., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan]],
          shape=(3000, 49), dtype=float32)

  • lon
    (trajectory, time)
    float32
    19.2 19.15 19.11 ... nan nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    axis :

    X

    minval :

    18.822432

    maxval :

    19.728199

    array([[19.195274, 19.148632, 19.111015, ...,       nan,       nan,
                  nan],
           [19.191893, 19.142187, 19.1418  , ...,       nan,       nan,
                  nan],
           [19.193327, 19.166204, 19.15524 , ..., 19.29427 , 19.290596,
            19.25856 ],
           ...,
           [19.190422, 19.168274, 19.160885, ..., 19.24517 , 19.234701,
            19.23248 ],
           [19.187893, 19.123146, 19.109854, ...,       nan,       nan,
                  nan],
           [19.1907  , 19.161379, 19.140356, ...,       nan,       nan,
                  nan]], shape=(3000, 49), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    axis :

    Y

    minval :

    79.572624

    maxval :

    79.76473

    array([[79.59931 , 79.60405 , 79.60682 , ...,       nan,       nan,
                  nan],
           [79.598854, 79.60169 , 79.60677 , ...,       nan,       nan,
                  nan],
           [79.59855 , 79.60331 , 79.60967 , ..., 79.73504 , 79.74171 ,
            79.74606 ],
           ...,
           [79.59836 , 79.59996 , 79.60392 , ..., 79.74531 , 79.75156 ,
            79.75524 ],
           [79.59878 , 79.60379 , 79.607605, ...,       nan,       nan,
                  nan],
           [79.59873 , 79.60167 , 79.60656 , ...,       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.45525372

    maxval :

    0.0

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

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

    dtype :

    <class 'numpy.float32'>

    units :

    m2/s

    description :

    Kinematic viscosity of oil (emulsion)

    minval :

    0.021692334

    maxval :

    0.11361812

    array([[0.02169233, 0.11047644, 0.10998931, ...,        nan,        nan,
                   nan],
           [0.02169233, 0.11007517, 0.11037365, ...,        nan,        nan,
                   nan],
           [0.02169233, 0.11022793, 0.11042295, ..., 0.09061252, 0.09086485,
            0.09082169],
           ...,
           [0.02169233, 0.11031399, 0.11032799, ..., 0.09142069, 0.0913781 ,
            0.0911952 ],
           [0.02169233, 0.1100581 , 0.10993555, ...,        nan,        nan,
                   nan],
           [0.02169233, 0.11042711, 0.11040375, ...,        nan,        nan,
                   nan]], shape=(3000, 49), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    kg/m^3

    minval :

    988.1

    maxval :

    995.5477

    array([[988.1    , 995.42645, 995.40735, ...,       nan,       nan,
                  nan],
           [988.1    , 995.4107 , 995.4224 , ...,       nan,       nan,
                  nan],
           [988.1    , 995.41675, 995.4244 , ..., 994.56476, 994.5769 ,
            994.5748 ],
           ...,
           [988.1    , 995.4201 , 995.42065, ..., 994.6035 , 994.6015 ,
            994.5927 ],
           [988.1    , 995.41003, 995.4052 , ...,       nan,       nan,
                  nan],
           [988.1    , 995.42456, 995.42365, ...,       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., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 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., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 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., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 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., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 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., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 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., ..., nan, nan, nan],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           ...,
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 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., ..., nan, nan, nan],
           [ 3.,  3.,  3., ...,  3.,  3.,  3.],
           ...,
           [ 3.,  3.,  3., ...,  3.,  3.,  3.],
           [ 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., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 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., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 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, ...,   nan,   nan,   nan],
           [0.001, 0.001, 0.001, ..., 0.001, 0.001, 0.001],
           ...,
           [0.001, 0.001, 0.001, ..., 0.001, 0.001, 0.001],
           [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.0021178133

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

    minval :

    -0.31912243

    maxval :

    0.42718613

    array([[-0.05641371, -0.14377157,  0.08761997, ...,         nan,
                    nan,         nan],
           [-0.2064487 ,  0.00562726,  0.0361237 , ...,         nan,
                    nan,         nan],
           [-0.15561533, -0.06263353,  0.08226167, ...,  0.04571818,
            -0.07934075, -0.07471409],
           ...,
           [-0.20295441, -0.0335131 ,  0.07111118, ...,  0.01753047,
            -0.13615668, -0.02868745],
           [-0.2152395 ,  0.01352025,  0.12318778, ...,         nan,
                    nan,         nan],
           [-0.12681845, -0.03649944, -0.02313464, ...,         nan,
                    nan,         nan]], shape=(3000, 49), dtype=float32)

  • y_sea_water_velocity
    (trajectory, time)
    float32
    0.2543 0.04278 -0.05633 ... nan nan

    minval :

    -0.2850785

    maxval :

    0.44436228

    array([[ 0.2543064 ,  0.04277782, -0.05633118, ...,         nan,
                    nan,         nan],
           [ 0.20843747,  0.14740714, -0.05484227, ...,         nan,
                    nan,         nan],
           [ 0.1819615 ,  0.08718619, -0.11116382, ..., -0.00584967,
             0.02482928,  0.04756401],
           ...,
           [ 0.08803403,  0.11434535,  0.01940844, ...,  0.10956893,
             0.03122267,  0.03919407],
           [ 0.18729106,  0.14621145, -0.01700595, ...,         nan,
                    nan,         nan],
           [ 0.1983785 ,  0.12352369, -0.04598727, ...,         nan,
                    nan,         nan]], shape=(3000, 49), dtype=float32)

  • x_wind
    (trajectory, time)
    float32
    -1.84 -2.031 -1.989 ... nan nan nan

    minval :

    -4.64724

    maxval :

    3.700867

    array([[-1.8401947 , -2.0307891 , -1.9892282 , ...,         nan,
                    nan,         nan],
           [-2.0833268 , -2.0281014 , -2.5351765 , ...,         nan,
                    nan,         nan],
           [-2.8643618 , -3.216069  , -3.1997237 , ..., -0.65315324,
            -1.2861427 , -0.66422355],
           ...,
           [-1.8948915 , -1.7894654 , -1.9905291 , ..., -0.45738816,
            -1.7658888 , -0.06998288],
           [-2.587947  , -2.3767173 , -2.55124   , ...,         nan,
                    nan,         nan],
           [-2.389606  , -3.4056776 , -2.0837898 , ...,         nan,
                    nan,         nan]], shape=(3000, 49), dtype=float32)

  • y_wind
    (trajectory, time)
    float32
    0.6632 -1.039 0.811 ... nan nan nan

    minval :

    -2.1853464

    maxval :

    6.943157

    array([[ 0.6632208 , -1.0394273 ,  0.8109766 , ...,         nan,
                    nan,         nan],
           [-0.26719078,  0.5355398 , -0.33488005, ...,         nan,
                    nan,         nan],
           [-0.7520711 , -0.04717249, -0.01847894, ...,  2.719156  ,
             4.2494454 ,  3.1935933 ],
           ...,
           [-0.55541074, -0.37615138,  1.0003211 , ...,  3.9340985 ,
             2.7919319 ,  3.1850066 ],
           [-0.6320049 , -0.6114157 , -0.24604328, ...,         nan,
                    nan,         nan],
           [-0.5298153 ,  0.59236324, -0.4340918 , ...,         nan,
                    nan,         nan]], shape=(3000, 49), dtype=float32)

  • sea_surface_height
    (trajectory, time)
    float32
    -0.8377 -0.7075 -0.5147 ... nan nan

    minval :

    -0.83941036

    maxval :

    -0.31428507

    array([[-0.83766615, -0.70748925, -0.5146781 , ...,         nan,
                    nan,         nan],
           [-0.83781374, -0.7083188 , -0.51605904, ...,         nan,
                    nan,         nan],
           [-0.8380074 , -0.70861864, -0.51627207, ..., -0.61341417,
            -0.7184755 , -0.7853731 ],
           ...,
           [-0.83803046, -0.7101185 , -0.5175551 , ..., -0.61375034,
            -0.7187615 , -0.78538036],
           [-0.83774817, -0.70662206, -0.5144497 , ...,         nan,
                    nan,         nan],
           [-0.8378454 , -0.7090695 , -0.516019  , ...,         nan,
                    nan,         nan]], shape=(3000, 49), dtype=float32)

  • sea_surface_wave_significant_height
    (trajectory, time)
    float32
    0.09412 0.128 0.1135 ... nan nan

    minval :

    3.4868182e-05

    maxval :

    1.1999053

    array([[0.09412398, 0.12803105, 0.11352192, ...,        nan,        nan,
                   nan],
           [0.10852639, 0.10823994, 0.16086589, ...,        nan,        nan,
                   nan],
           [0.21574637, 0.25449497, 0.2518689 , ..., 0.19238229, 0.48491392,
            0.26174968],
           ...,
           [0.09591773, 0.08225444, 0.12208606, ..., 0.3858838 , 0.26846588,
            0.24966948],
           [0.17458373, 0.14815632, 0.1616063 , ...,        nan,        nan,
                   nan],
           [0.14737666, 0.29395851, 0.11145315, ...,        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., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 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., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 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., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 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.671 1.949 1.836 ... nan nan nan

    minval :

    0.032169495

    maxval :

    5.9676456

    array([[1.6713973, 1.9493384, 1.8355637, ...,       nan,       nan,
                  nan],
           [1.7947223, 1.7923522, 2.185051 , ...,       nan,       nan,
                  nan],
           [2.5304708, 2.7483318, 2.7341154, ..., 2.3895278, 3.7936907,
            2.7872293],
           ...,
           [1.6872483, 1.5624614, 1.9035425, ..., 3.3842144, 2.8227613,
            2.7221515],
           [2.2763104, 2.0969577, 2.1900735, ...,       nan,       nan,
                  nan],
           [2.0914328, 2.9537425, 1.8187615, ...,       nan,       nan,
                  nan]], shape=(3000, 49), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.00031643617

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

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

    minval :

    -0.0084534995

    maxval :

    0.0

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

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

    minval :

    0.0

    maxval :

    0.007534966

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

  • sea_water_temperature
    (trajectory, time)
    float32
    0.8571 0.8745 0.8869 ... nan nan

    minval :

    0.68453664

    maxval :

    275.76245

    array([[  0.8571021 ,   0.8745475 ,   0.88691986, ...,          nan,
                     nan,          nan],
           [  0.8595025 ,   0.8885725 ,   0.8659982 , ...,          nan,
                     nan,          nan],
           [  0.8613846 ,   0.86736465,   0.8503495 , ...,   2.1312017 ,
              2.135333  , 275.28778   ],
           ...,
           [  0.86227113,   0.87798876,   0.8645746 , ...,   2.0887866 ,
              2.1109083 ,   2.1235008 ],
           [  0.8596623 ,   0.892982  ,   0.8843943 , ...,          nan,
                     nan,          nan],
           [  0.860124  ,   0.87566704,   0.8673824 , ...,          nan,
                     nan,          nan]], shape=(3000, 49), dtype=float32)

  • sea_water_salinity
    (trajectory, time)
    float32
    33.78 33.78 33.78 ... nan nan nan

    minval :

    33.44262

    maxval :

    33.82684

    array([[33.77563 , 33.776   , 33.776436, ...,       nan,       nan,
                  nan],
           [33.775543, 33.77684 , 33.776283, ...,       nan,       nan,
                  nan],
           [33.77572 , 33.77651 , 33.77558 , ..., 33.694443, 33.68427 ,
            33.666325],
           ...,
           [33.775593, 33.77746 , 33.77712 , ..., 33.627735, 33.672783,
            33.687046],
           [33.775322, 33.7761  , 33.77619 , ...,       nan,       nan,
                  nan],
           [33.77551 , 33.776848, 33.776333, ...,       nan,       nan,
                  nan]], shape=(3000, 49), dtype=float32)

  • sea_floor_depth_below_sea_level
    (trajectory, time)
    float32
    170.4 186.0 196.0 ... nan nan nan

    minval :

    10.0

    maxval :

    247.74042

    array([[170.4251  , 185.95155 , 195.95549 , ...,        nan,        nan,
                   nan],
           [170.93388 , 183.45432 , 190.0671  , ...,        nan,        nan,
                   nan],
           [170.35571 , 180.93307 , 187.62726 , ...,  22.598293,  18.373354,
             17.463322],
           ...,
           [170.9374  , 176.95679 , 182.81964 , ...,  16.654173,  14.02806 ,
             12.693246],
           [171.87785 , 189.58444 , 197.26588 , ...,        nan,        nan,
                   nan],
           [171.14177 , 180.53703 , 190.38316 , ...,        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, ...,  nan,  nan,  nan],
           [0.02, 0.02, 0.02, ..., 0.02, 0.02, 0.02],
           ...,
           [0.02, 0.02, 0.02, ..., 0.02, 0.02, 0.02],
           [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., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 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., ..., nan, nan, nan],
           [50., 50., 50., ..., 50., 50., 50.],
           ...,
           [50., 50., 50., ..., 50., 50., 50.],
           [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-07-04 17:14:52.964775', '2025-07-04 18:14:52.964775',
                   '2025-07-04 19:14:52.964775', '2025-07-04 20:14:52.964775',
                   '2025-07-04 21:14:52.964775', '2025-07-04 22:14:52.964775',
                   '2025-07-04 23:14:52.964775', '2025-07-05 00:14:52.964775',
                   '2025-07-05 01:14:52.964775', '2025-07-05 02:14:52.964775',
                   '2025-07-05 03:14:52.964775', '2025-07-05 04:14:52.964775',
                   '2025-07-05 05:14:52.964775', '2025-07-05 06:14:52.964775',
                   '2025-07-05 07:14:52.964775', '2025-07-05 08:14:52.964775',
                   '2025-07-05 09:14:52.964775', '2025-07-05 10:14:52.964775',
                   '2025-07-05 11:14:52.964775', '2025-07-05 12:14:52.964775',
                   '2025-07-05 13:14:52.964775', '2025-07-05 14:14:52.964775',
                   '2025-07-05 15:14:52.964775', '2025-07-05 16:14:52.964775',
                   '2025-07-05 17:14:52.964775', '2025-07-05 18:14:52.964775',
                   '2025-07-05 19:14:52.964775', '2025-07-05 20:14:52.964775',
                   '2025-07-05 21:14:52.964775', '2025-07-05 22:14:52.964775',
                   '2025-07-05 23:14:52.964775', '2025-07-06 00:14:52.964775',
                   '2025-07-06 01:14:52.964775', '2025-07-06 02:14:52.964775',
                   '2025-07-06 03:14:52.964775', '2025-07-06 04:14:52.964775',
                   '2025-07-06 05:14:52.964775', '2025-07-06 06:14:52.964775',
                   '2025-07-06 07:14:52.964775', '2025-07-06 08:14:52.964775',
                   '2025-07-06 09:14:52.964775', '2025-07-06 10:14:52.964775',
                   '2025-07-06 11:14:52.964775', '2025-07-06 12:14:52.964775',
                   '2025-07-06 13:14:52.964775', '2025-07-06 14:14:52.964775',
                   '2025-07-06 15:14:52.964775', '2025-07-06 16:14:52.964775',
                   '2025-07-06 17:14:52.964775'],
                  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-07-11 17:14:57.201009

  date_created :

  2025-07-11T17:14:57.201026

  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-07-04 17:14:52.964775

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

  seed_lon :

  19.1909

  seed_lat :

  79.5986

  seed_radius :

  50.0

  seed_time :

  2025-07-04 17:14:52.964775

  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-07-11 17:14:57.336788

  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-07-06 17:14:52.964775

  time_coverage_duration :

  P2DT0H0M0S

  time_coverage_resolution :

  P0DT0H30M0S

  performance :

  -------------------- Reader performance: -------------------- global_landmask 0:00:01.5 total 0:00:00.0 preparing 0:00:01.5 reading 0:00:00.0 masking -------------------- https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be 0:02:10.5 total 0:00:00.0 preparing 0:02:09.4 reading 0:00:01.5 interpolation 0:00:00.1 interpolation_time 0:00:01.0 rotating vectors 0:00:00.0 masking -------------------- Performance: 2:31.1 total time 4.2 configuration 0.1 preparing main loop 0.0 moving elements to ocean 2:26.7 main loop 6.4 updating elements 0.0 oil weathering 0.0 updating viscosities 0.0 updating densities 5.8 vertical mixing 0.0 cleaning up --------------------

  geospatial_bounds_crs :

  EPSG:4326

  geospatial_bounds_vertical_crs :

  EPSG:5831

  geospatial_lat_min :

  79.572624

  geospatial_lat_max :

  79.76473

  geospatial_lat_units :

  degrees_north

  geospatial_lat_resolution :

  point

  geospatial_lon_min :

  18.822432

  geospatial_lon_max :

  19.728199

  geospatial_lon_units :

  degrees_east

  geospatial_lon_resolution :

  point

  runtime :

  0:02:31.192135

  geospatial_vertical_min :

  -0.45525372

  geospatial_vertical_max :

  0.0

  geospatial_vertical_positive :

  up

Print and plot results

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

===========================
--------------------
Reader performance:
--------------------
global_landmask
 0:00:01.5  total
 0:00:00.0  preparing
 0:00:01.5  reading
 0:00:00.0  masking
--------------------
https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
 0:02:10.5  total
 0:00:00.0  preparing
 0:02:09.4  reading
 0:00:01.5  interpolation
 0:00:00.1  interpolation_time
 0:00:01.0  rotating vectors
 0:00:00.0  masking
--------------------
Performance:
 2:31.1 total time
    4.2 configuration
    0.1 preparing main loop
      0.0 moving elements to ocean
 2:26.7 main loop
      6.4 updating elements
        0.0 oil weathering
          0.0 updating viscosities
          0.0 updating densities
        5.8 vertical mixing
    0.0 cleaning up
--------------------
===========================
Model:  OpenOil     (OpenDrift version 1.14.2)
        678 active Oil particles  (2322 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-07-04 17:14:52.964775 UTC
        Present: 2025-07-06 17:14:52.964775 UTC
        Calculation steps: 96 * 0:30:00 - total time: 2 days, 0:00:00
        Output steps: 49 * 1:00:00
===========================

/opt/conda/envs/opendrift/lib/python3.13/site-packages/cartopy/mpl/geoaxes.py:1692: UserWarning: No data for colormapping provided via 'c'. Parameters 'vmin', 'vmax' will be ignored
  result = super().scatter(*args, **kwargs)
17:17:35 INFO    opendrift:4553: Saving animation to /root/project/docs/source/gallery/animations/example_oil_ice_0.gif...
17:18:15 INFO    opendrift:2994: Time to make animation: 0:00:51.360257

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-07-04 17:14 to 2025-07-06 17:14 UTC (49 steps)'}>, <Figure size 568.32x1100 with 2 Axes>)