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)

09:55:37 INFO    opendrift:509: OpenDriftSimulation initialised (version 1.14.2 / v1.14.2-59-g60ee835)

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.utcnow() - timedelta(days=7))

/root/project/examples/example_oil_ice.py:29: DeprecationWarning: datetime.datetime.utcnow() is deprecated and scheduled for removal in a future version. Use timezone-aware objects to represent datetimes in UTC: datetime.datetime.now(datetime.UTC).
  number=3000, time=datetime.utcnow() - timedelta(days=7))
09:55:37 INFO    opendrift.models.openoil.openoil:1703: Oil type not specified, using default: GENERIC BUNKER C
09:55:37 INFO    opendrift.models.openoil.adios.dirjs:86: Querying ADIOS database for oil: GENERIC BUNKER C
09:55:37 INFO    opendrift.models.openoil.openoil:1712: Using density 988.1 and viscosity 0.02169233387797564 of oiltype GENERIC BUNKER C
09:55:37 INFO    opendrift.models.basemodel.environment:206: Adding a global landmask from GSHHG
09:55:42 INFO    opendrift.models.basemodel.environment:229: Fallback values will be used for the following variables which have no readers:
09:55:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_height: 0.000000
09:55:42 INFO    opendrift.models.basemodel.environment:232:    upward_sea_water_velocity: 0.000000
09:55:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_significant_height: 0.000000
09:55:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
09:55:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
09:55:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
09:55:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
09:55:42 INFO    opendrift.models.basemodel.environment:232:    sea_ice_area_fraction: 0.000000
09:55:42 INFO    opendrift.models.basemodel.environment:232:    sea_ice_x_velocity: 0.000000
09:55:42 INFO    opendrift.models.basemodel.environment:232:    sea_ice_y_velocity: 0.000000
09:55:42 INFO    opendrift.models.basemodel.environment:232:    sea_water_temperature: 10.000000
09:55:42 INFO    opendrift.models.basemodel.environment:232:    sea_water_salinity: 34.000000
09:55:42 INFO    opendrift.models.basemodel.environment:232:    sea_floor_depth_below_sea_level: 10000.000000
09:55:42 INFO    opendrift.models.basemodel.environment:232:    ocean_vertical_diffusivity: 0.020000
09:55:42 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)

09:55:42 INFO    opendrift:885: Using existing reader for land_binary_mask
09:55:42 INFO    opendrift:914: All points are in ocean
09:55:42 INFO    opendrift:2035: Storing previous position of elements for coastline interaction
09:55:42 INFO    opendrift.models.openoil.openoil:687: Oil-water surface tension is 0.035935 Nm
09:55:42 INFO    opendrift.models.openoil.openoil:700: Max water fraction not available for GENERIC BUNKER C, using default
09:55:42 INFO    opendrift:2096: 2025-06-19 09:55:37.704451 - step 1 of 96 - 3000 active elements (0 deactivated)
09:55:42 INFO    opendrift.readers:61: Opening file with xr.open_dataset
09:55:44 INFO    opendrift.readers.reader_netCDF_CF_generic:332: Detected dimensions: {'x': 'X', 'y': 'Y', 'z': 'depth', 'time': 'time'}
09:55:44 INFO    opendrift.readers.basereader:176: Variable x_sea_water_velocity will be rotated from eastward_sea_water_velocity
09:55:44 INFO    opendrift.readers.basereader:176: Variable y_sea_water_velocity will be rotated from northward_sea_water_velocity
09:55:47 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0038807474 to -0.00063441036) found for ocean_vertical_diffusivity, replacing with NaN
09:55:47 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:55:49 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0049419613 to -0.0002660835) found for ocean_vertical_diffusivity, replacing with NaN
09:55:49 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:55:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:55:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:55:49 INFO    opendrift:2096: 2025-06-19 10:25:37.704451 - step 2 of 96 - 3000 active elements (0 deactivated)
09:55:52 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.006478963 to -0.00068334315) found for ocean_vertical_diffusivity, replacing with NaN
09:55:52 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:55:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:55:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:55:52 INFO    opendrift:2096: 2025-06-19 10:55:37.704451 - step 3 of 96 - 3000 active elements (0 deactivated)
09:55:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:55:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:55:52 INFO    opendrift:2096: 2025-06-19 11:25:37.704451 - step 4 of 96 - 3000 active elements (0 deactivated)
09:55:56 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.008040514 to -0.0008865317) found for ocean_vertical_diffusivity, replacing with NaN
09:55:56 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:55:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:55:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:55:56 INFO    opendrift:2096: 2025-06-19 11:55:37.704451 - step 5 of 96 - 3000 active elements (0 deactivated)
09:55:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:55:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:55:58 INFO    opendrift:2096: 2025-06-19 12:25:37.704451 - step 6 of 96 - 2996 active elements (4 deactivated)
09:56:00 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0107278945 to -0.00066464697) found for ocean_vertical_diffusivity, replacing with NaN
09:56:00 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:00 INFO    opendrift:2096: 2025-06-19 12:55:37.704451 - step 7 of 96 - 2973 active elements (27 deactivated)
09:56:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:01 INFO    opendrift:2096: 2025-06-19 13:25:37.704451 - step 8 of 96 - 2923 active elements (77 deactivated)
09:56:03 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.016850984 to -0.0017077824) found for ocean_vertical_diffusivity, replacing with NaN
09:56:03 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:03 INFO    opendrift:2096: 2025-06-19 13:55:37.704451 - step 9 of 96 - 2869 active elements (131 deactivated)
09:56:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:03 INFO    opendrift:2096: 2025-06-19 14:25:37.704451 - step 10 of 96 - 2809 active elements (191 deactivated)
09:56:06 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.02185658 to -0.002996402) found for ocean_vertical_diffusivity, replacing with NaN
09:56:06 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:06 INFO    opendrift:2096: 2025-06-19 14:55:37.704451 - step 11 of 96 - 2773 active elements (227 deactivated)
09:56:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:06 INFO    opendrift:2096: 2025-06-19 15:25:37.704451 - step 12 of 96 - 2756 active elements (244 deactivated)
09:56:09 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.020658039 to -0.0026508337) found for ocean_vertical_diffusivity, replacing with NaN
09:56:09 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:09 INFO    opendrift:2096: 2025-06-19 15:55:37.704451 - step 13 of 96 - 2752 active elements (248 deactivated)
09:56:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:09 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:09 INFO    opendrift:2096: 2025-06-19 16:25:37.704451 - step 14 of 96 - 2749 active elements (251 deactivated)
09:56:11 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.01468062 to -0.0011775766) found for ocean_vertical_diffusivity, replacing with NaN
09:56:11 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:12 INFO    opendrift:2096: 2025-06-19 16:55:37.704451 - step 15 of 96 - 2746 active elements (254 deactivated)
09:56:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:12 INFO    opendrift:2096: 2025-06-19 17:25:37.704451 - step 16 of 96 - 2745 active elements (255 deactivated)
09:56:14 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.011350506 to -0.0010385409) found for ocean_vertical_diffusivity, replacing with NaN
09:56:14 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:14 INFO    opendrift:2096: 2025-06-19 17:55:37.704451 - step 17 of 96 - 2741 active elements (259 deactivated)
09:56:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:14 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:14 INFO    opendrift:2096: 2025-06-19 18:25:37.704451 - step 18 of 96 - 2728 active elements (272 deactivated)
09:56:17 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.008115157 to -0.00028125287) found for ocean_vertical_diffusivity, replacing with NaN
09:56:17 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:17 INFO    opendrift:2096: 2025-06-19 18:55:37.704451 - step 19 of 96 - 2695 active elements (305 deactivated)
09:56:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:17 INFO    opendrift:2096: 2025-06-19 19:25:37.704451 - step 20 of 96 - 2648 active elements (352 deactivated)
09:56:20 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.006423709 to -0.00010702345) found for ocean_vertical_diffusivity, replacing with NaN
09:56:20 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:20 INFO    opendrift:2096: 2025-06-19 19:55:37.704451 - step 21 of 96 - 2568 active elements (432 deactivated)
09:56:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:20 INFO    opendrift:2096: 2025-06-19 20:25:37.704451 - step 22 of 96 - 2479 active elements (521 deactivated)
09:56:22 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0056433515 to -5.2381136e-05) found for ocean_vertical_diffusivity, replacing with NaN
09:56:22 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:23 INFO    opendrift:2096: 2025-06-19 20:55:37.704451 - step 23 of 96 - 2377 active elements (623 deactivated)
09:56:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:23 INFO    opendrift:2096: 2025-06-19 21:25:37.704451 - step 24 of 96 - 2296 active elements (704 deactivated)
09:56:25 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0051176753 to -0.00012112187) found for ocean_vertical_diffusivity, replacing with NaN
09:56:25 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:25 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:25 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:25 INFO    opendrift:2096: 2025-06-19 21:55:37.704451 - step 25 of 96 - 2203 active elements (797 deactivated)
09:56:25 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:25 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:26 INFO    opendrift:2096: 2025-06-19 22:25:37.704451 - step 26 of 96 - 2140 active elements (860 deactivated)
09:56:28 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0044715037 to -0.00021881457) found for ocean_vertical_diffusivity, replacing with NaN
09:56:28 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:28 INFO    opendrift:2096: 2025-06-19 22:55:37.704451 - step 27 of 96 - 2076 active elements (924 deactivated)
09:56:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:28 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:28 INFO    opendrift:2096: 2025-06-19 23:25:37.704451 - step 28 of 96 - 2009 active elements (991 deactivated)
09:56:31 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.004009147 to -0.00026988244) found for ocean_vertical_diffusivity, replacing with NaN
09:56:31 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:31 INFO    opendrift:2096: 2025-06-19 23:55:37.704451 - step 29 of 96 - 1968 active elements (1032 deactivated)
09:56:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:31 INFO    opendrift:2096: 2025-06-20 00:25:37.704451 - step 30 of 96 - 1935 active elements (1065 deactivated)
09:56:34 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0052477447 to -0.00079089846) found for ocean_vertical_diffusivity, replacing with NaN
09:56:34 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:34 INFO    opendrift:2096: 2025-06-20 00:55:37.704451 - step 31 of 96 - 1911 active elements (1089 deactivated)
09:56:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:34 INFO    opendrift:2096: 2025-06-20 01:25:37.704451 - step 32 of 96 - 1896 active elements (1104 deactivated)
09:56:36 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0060557453 to -0.0012976914) found for ocean_vertical_diffusivity, replacing with NaN
09:56:36 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:36 INFO    opendrift:2096: 2025-06-20 01:55:37.704451 - step 33 of 96 - 1876 active elements (1124 deactivated)
09:56:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:36 INFO    opendrift:2096: 2025-06-20 02:25:37.704451 - step 34 of 96 - 1862 active elements (1138 deactivated)
09:56:39 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.013219904 to -0.0014358633) found for ocean_vertical_diffusivity, replacing with NaN
09:56:39 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:39 INFO    opendrift:2096: 2025-06-20 02:55:37.704451 - step 35 of 96 - 1852 active elements (1148 deactivated)
09:56:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:39 INFO    opendrift:2096: 2025-06-20 03:25:37.704451 - step 36 of 96 - 1837 active elements (1163 deactivated)
09:56:41 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.015195384 to -0.0009940694) found for ocean_vertical_diffusivity, replacing with NaN
09:56:41 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:41 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:41 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:41 INFO    opendrift:2096: 2025-06-20 03:55:37.704451 - step 37 of 96 - 1808 active elements (1192 deactivated)
09:56:41 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:41 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:42 INFO    opendrift:2096: 2025-06-20 04:25:37.704451 - step 38 of 96 - 1770 active elements (1230 deactivated)
09:56:44 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.007843118 to -0.00037491985) found for ocean_vertical_diffusivity, replacing with NaN
09:56:44 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:44 INFO    opendrift:2096: 2025-06-20 04:55:37.704451 - step 39 of 96 - 1722 active elements (1278 deactivated)
09:56:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:44 INFO    opendrift:2096: 2025-06-20 05:25:37.704451 - step 40 of 96 - 1672 active elements (1328 deactivated)
09:56:47 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.005326221 to -0.0005178188) found for ocean_vertical_diffusivity, replacing with NaN
09:56:47 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:47 INFO    opendrift:2096: 2025-06-20 05:55:37.704451 - step 41 of 96 - 1637 active elements (1363 deactivated)
09:56:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:47 INFO    opendrift:2096: 2025-06-20 06:25:37.704451 - step 42 of 96 - 1608 active elements (1392 deactivated)
09:56:49 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0071207583 to -0.000338716) found for ocean_vertical_diffusivity, replacing with NaN
09:56:49 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:49 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:50 INFO    opendrift:2096: 2025-06-20 06:55:37.704451 - step 43 of 96 - 1580 active elements (1420 deactivated)
09:56:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:50 INFO    opendrift:2096: 2025-06-20 07:25:37.704451 - step 44 of 96 - 1557 active elements (1443 deactivated)
09:56:52 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0050919075 to -4.1905372e-05) found for ocean_vertical_diffusivity, replacing with NaN
09:56:52 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:52 INFO    opendrift:2096: 2025-06-20 07:55:37.704451 - step 45 of 96 - 1533 active elements (1467 deactivated)
09:56:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:52 INFO    opendrift:2096: 2025-06-20 08:25:37.704451 - step 46 of 96 - 1503 active elements (1497 deactivated)
09:56:55 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0031560864 to -6.384925e-05) found for ocean_vertical_diffusivity, replacing with NaN
09:56:55 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:55 INFO    opendrift:2096: 2025-06-20 08:55:37.704451 - step 47 of 96 - 1477 active elements (1523 deactivated)
09:56:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:55 INFO    opendrift:2096: 2025-06-20 09:25:37.704451 - step 48 of 96 - 1461 active elements (1539 deactivated)
09:56:57 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0032376281 to -1.6135262e-05) found for ocean_vertical_diffusivity, replacing with NaN
09:56:57 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:56:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:57 INFO    opendrift:2096: 2025-06-20 09:55:37.704451 - step 49 of 96 - 1440 active elements (1560 deactivated)
09:56:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:56:58 INFO    opendrift:2096: 2025-06-20 10:25:37.704451 - step 50 of 96 - 1429 active elements (1571 deactivated)
09:57:00 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.004359997 to -8.4807005e-05) found for ocean_vertical_diffusivity, replacing with NaN
09:57:00 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:00 INFO    opendrift:2096: 2025-06-20 10:55:37.704451 - step 51 of 96 - 1417 active elements (1583 deactivated)
09:57:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:00 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:00 INFO    opendrift:2096: 2025-06-20 11:25:37.704451 - step 52 of 96 - 1401 active elements (1599 deactivated)
09:57:03 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0038768533 to -0.00012425931) found for ocean_vertical_diffusivity, replacing with NaN
09:57:03 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:03 INFO    opendrift:2096: 2025-06-20 11:55:37.704451 - step 53 of 96 - 1398 active elements (1602 deactivated)
09:57:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:03 INFO    opendrift:2096: 2025-06-20 12:25:37.704451 - step 54 of 96 - 1392 active elements (1608 deactivated)
09:57:05 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0024419473 to -8.36557e-05) found for ocean_vertical_diffusivity, replacing with NaN
09:57:05 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:05 INFO    opendrift:2096: 2025-06-20 12:55:37.704451 - step 55 of 96 - 1386 active elements (1614 deactivated)
09:57:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:05 INFO    opendrift:2096: 2025-06-20 13:25:37.704451 - step 56 of 96 - 1385 active elements (1615 deactivated)
09:57:08 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0066331998 to -0.00016302927) found for ocean_vertical_diffusivity, replacing with NaN
09:57:08 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:08 INFO    opendrift:2096: 2025-06-20 13:55:37.704451 - step 57 of 96 - 1382 active elements (1618 deactivated)
09:57:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:08 INFO    opendrift:2096: 2025-06-20 14:25:37.704451 - step 58 of 96 - 1381 active elements (1619 deactivated)
09:57:10 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.024608716 to -0.0009993697) found for ocean_vertical_diffusivity, replacing with NaN
09:57:10 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:11 INFO    opendrift:2096: 2025-06-20 14:55:37.704451 - step 59 of 96 - 1381 active elements (1619 deactivated)
09:57:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:11 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:11 INFO    opendrift:2096: 2025-06-20 15:25:37.704451 - step 60 of 96 - 1381 active elements (1619 deactivated)
09:57:13 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.041307863 to -0.0012453041) found for ocean_vertical_diffusivity, replacing with NaN
09:57:13 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:13 INFO    opendrift:2096: 2025-06-20 15:55:37.704451 - step 61 of 96 - 1376 active elements (1624 deactivated)
09:57:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:13 INFO    opendrift:2096: 2025-06-20 16:25:37.704451 - step 62 of 96 - 1375 active elements (1625 deactivated)
09:57:16 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.025373979 to -0.00069156836) found for ocean_vertical_diffusivity, replacing with NaN
09:57:16 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:16 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:16 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:16 INFO    opendrift:2096: 2025-06-20 16:55:37.704451 - step 63 of 96 - 1369 active elements (1631 deactivated)
09:57:16 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:16 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:16 INFO    opendrift:2096: 2025-06-20 17:25:37.704451 - step 64 of 96 - 1362 active elements (1638 deactivated)
09:57:18 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.006251681 to -8.535373e-05) found for ocean_vertical_diffusivity, replacing with NaN
09:57:18 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:18 INFO    opendrift:2096: 2025-06-20 17:55:37.704451 - step 65 of 96 - 1356 active elements (1644 deactivated)
09:57:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:18 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:18 INFO    opendrift:2096: 2025-06-20 18:25:37.704451 - step 66 of 96 - 1346 active elements (1654 deactivated)
09:57:21 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.013048751 to -4.221712e-05) found for ocean_vertical_diffusivity, replacing with NaN
09:57:21 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:21 INFO    opendrift:2096: 2025-06-20 18:55:37.704451 - step 67 of 96 - 1342 active elements (1658 deactivated)
09:57:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:21 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:21 INFO    opendrift:2096: 2025-06-20 19:25:37.704451 - step 68 of 96 - 1339 active elements (1661 deactivated)
09:57:24 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.010597968 to -1.050836e-05) found for ocean_vertical_diffusivity, replacing with NaN
09:57:24 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:24 INFO    opendrift:2096: 2025-06-20 19:55:37.704451 - step 69 of 96 - 1338 active elements (1662 deactivated)
09:57:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:24 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:24 INFO    opendrift:2096: 2025-06-20 20:25:37.704451 - step 70 of 96 - 1337 active elements (1663 deactivated)
09:57:26 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.008593804 to -7.0348565e-06) found for ocean_vertical_diffusivity, replacing with NaN
09:57:26 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:26 INFO    opendrift:2096: 2025-06-20 20:55:37.704451 - step 71 of 96 - 1337 active elements (1663 deactivated)
09:57:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:26 INFO    opendrift:2096: 2025-06-20 21:25:37.704451 - step 72 of 96 - 1337 active elements (1663 deactivated)
09:57:29 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0092997225 to -0.00010527169) found for ocean_vertical_diffusivity, replacing with NaN
09:57:29 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:29 INFO    opendrift:2096: 2025-06-20 21:55:37.704451 - step 73 of 96 - 1335 active elements (1665 deactivated)
09:57:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:29 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:29 INFO    opendrift:2096: 2025-06-20 22:25:37.704451 - step 74 of 96 - 1332 active elements (1668 deactivated)
09:57:31 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.012323378 to -3.76049e-05) found for ocean_vertical_diffusivity, replacing with NaN
09:57:31 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:31 INFO    opendrift:2096: 2025-06-20 22:55:37.704451 - step 75 of 96 - 1332 active elements (1668 deactivated)
09:57:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:31 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:32 INFO    opendrift:2096: 2025-06-20 23:25:37.704451 - step 76 of 96 - 1332 active elements (1668 deactivated)
09:57:34 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.014206214 to -1.326749e-05) found for ocean_vertical_diffusivity, replacing with NaN
09:57:34 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:34 INFO    opendrift:2096: 2025-06-20 23:55:37.704451 - step 77 of 96 - 1332 active elements (1668 deactivated)
09:57:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:34 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:34 INFO    opendrift:2096: 2025-06-21 00:25:37.704451 - step 78 of 96 - 1331 active elements (1669 deactivated)
09:57:36 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.008003669 to -7.797356e-05) found for ocean_vertical_diffusivity, replacing with NaN
09:57:36 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:36 INFO    opendrift:2096: 2025-06-21 00:55:37.704451 - step 79 of 96 - 1331 active elements (1669 deactivated)
09:57:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:36 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:37 INFO    opendrift:2096: 2025-06-21 01:25:37.704451 - step 80 of 96 - 1331 active elements (1669 deactivated)
09:57:39 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0024179 to -0.00023085844) found for ocean_vertical_diffusivity, replacing with NaN
09:57:39 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:39 INFO    opendrift:2096: 2025-06-21 01:55:37.704451 - step 81 of 96 - 1330 active elements (1670 deactivated)
09:57:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:39 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:39 INFO    opendrift:2096: 2025-06-21 02:25:37.704451 - step 82 of 96 - 1329 active elements (1671 deactivated)
09:57:42 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.008795474 to -0.00028054838) found for ocean_vertical_diffusivity, replacing with NaN
09:57:42 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:42 INFO    opendrift:2096: 2025-06-21 02:55:37.704451 - step 83 of 96 - 1328 active elements (1672 deactivated)
09:57:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:42 INFO    opendrift:2096: 2025-06-21 03:25:37.704451 - step 84 of 96 - 1328 active elements (1672 deactivated)
09:57:44 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.026724093 to -0.00069174176) found for ocean_vertical_diffusivity, replacing with NaN
09:57:44 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:44 INFO    opendrift:2096: 2025-06-21 03:55:37.704451 - step 85 of 96 - 1328 active elements (1672 deactivated)
09:57:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:44 INFO    opendrift:2096: 2025-06-21 04:25:37.704451 - step 86 of 96 - 1328 active elements (1672 deactivated)
09:57:47 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.028631026 to -0.0012786736) found for ocean_vertical_diffusivity, replacing with NaN
09:57:47 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:47 INFO    opendrift:2096: 2025-06-21 04:55:37.704451 - step 87 of 96 - 1327 active elements (1673 deactivated)
09:57:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:47 INFO    opendrift:2096: 2025-06-21 05:25:37.704451 - step 88 of 96 - 1325 active elements (1675 deactivated)
09:57:50 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.008696288 to -0.0014969428) found for ocean_vertical_diffusivity, replacing with NaN
09:57:50 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:50 INFO    opendrift:2096: 2025-06-21 05:55:37.704451 - step 89 of 96 - 1321 active elements (1679 deactivated)
09:57:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:50 INFO    opendrift:2096: 2025-06-21 06:25:37.704451 - step 90 of 96 - 1321 active elements (1679 deactivated)
09:57:53 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.01108552 to -0.0002603525) found for ocean_vertical_diffusivity, replacing with NaN
09:57:53 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:53 INFO    opendrift:2096: 2025-06-21 06:55:37.704451 - step 91 of 96 - 1321 active elements (1679 deactivated)
09:57:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:53 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:53 INFO    opendrift:2096: 2025-06-21 07:25:37.704451 - step 92 of 96 - 1321 active elements (1679 deactivated)
09:57:56 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.008484513 to -0.00030180786) found for ocean_vertical_diffusivity, replacing with NaN
09:57:56 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:56 INFO    opendrift:2096: 2025-06-21 07:55:37.704451 - step 93 of 96 - 1321 active elements (1679 deactivated)
09:57:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:56 INFO    opendrift:2096: 2025-06-21 08:25:37.704451 - step 94 of 96 - 1321 active elements (1679 deactivated)
09:57:59 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0043270783 to -0.00032995394) found for ocean_vertical_diffusivity, replacing with NaN
09:57:59 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:57:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:59 INFO    opendrift:2096: 2025-06-21 08:55:37.704451 - step 95 of 96 - 1320 active elements (1680 deactivated)
09:57:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:57:59 INFO    opendrift:2096: 2025-06-21 09:25:37.704451 - step 96 of 96 - 1320 active elements (1680 deactivated)
09:58:02 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0036665183 to -0.00047354965) found for ocean_vertical_diffusivity, replacing with NaN
09:58:02 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
09:58:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
09:58:02 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning

<xarray.Dataset> Size: 26MB
Dimensions:                                                                              (
                                                                                          trajectory: 3000,
                                                                                          time: 49)
Coordinates:
  * trajectory                                                                           (trajectory) int64 24kB ...
  * time                                                                                 (time) datetime64[ns] 392B ...
Data variables: (12/44)
    status                                                                               (trajectory, time) float32 588kB ...
    moving                                                                               (trajectory, time) float32 588kB ...
    age_seconds                                                                          (trajectory, time) float32 588kB ...
    origin_marker                                                                        (trajectory, time) float32 588kB ...
    lon                                                                                  (trajectory, time) float32 588kB ...
    lat                                                                                  (trajectory, time) float32 588kB ...
    ...                                                                                   ...
    sea_water_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/164)
    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-06-19T09:55:37.704451 ... 2...

    standard_name :

    time

    long_name :

    time

    dtype :

    <class 'numpy.float64'>

    array(['2025-06-19T09:55:37.704451000', '2025-06-19T10:55:37.704451000',
           '2025-06-19T11:55:37.704451000', '2025-06-19T12:55:37.704451000',
           '2025-06-19T13:55:37.704451000', '2025-06-19T14:55:37.704451000',
           '2025-06-19T15:55:37.704451000', '2025-06-19T16:55:37.704451000',
           '2025-06-19T17:55:37.704451000', '2025-06-19T18:55:37.704451000',
           '2025-06-19T19:55:37.704451000', '2025-06-19T20:55:37.704451000',
           '2025-06-19T21:55:37.704451000', '2025-06-19T22:55:37.704451000',
           '2025-06-19T23:55:37.704451000', '2025-06-20T00:55:37.704451000',
           '2025-06-20T01:55:37.704451000', '2025-06-20T02:55:37.704451000',
           '2025-06-20T03:55:37.704451000', '2025-06-20T04:55:37.704451000',
           '2025-06-20T05:55:37.704451000', '2025-06-20T06:55:37.704451000',
           '2025-06-20T07:55:37.704451000', '2025-06-20T08:55:37.704451000',
           '2025-06-20T09:55:37.704451000', '2025-06-20T10:55:37.704451000',
           '2025-06-20T11:55:37.704451000', '2025-06-20T12:55:37.704451000',
           '2025-06-20T13:55:37.704451000', '2025-06-20T14:55:37.704451000',
           '2025-06-20T15:55:37.704451000', '2025-06-20T16:55:37.704451000',
           '2025-06-20T17:55:37.704451000', '2025-06-20T18:55:37.704451000',
           '2025-06-20T19:55:37.704451000', '2025-06-20T20:55:37.704451000',
           '2025-06-20T21:55:37.704451000', '2025-06-20T22:55:37.704451000',
           '2025-06-20T23:55:37.704451000', '2025-06-21T00:55:37.704451000',
           '2025-06-21T01:55:37.704451000', '2025-06-21T02:55:37.704451000',
           '2025-06-21T03:55:37.704451000', '2025-06-21T04:55:37.704451000',
           '2025-06-21T05:55:37.704451000', '2025-06-21T06:55:37.704451000',
           '2025-06-21T07:55:37.704451000', '2025-06-21T08:55:37.704451000',
           '2025-06-21T09:55:37.704451000'], dtype='datetime64[ns]')

• Data variables: (44)
  • 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., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan]],
          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., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ..., nan, nan, nan],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           ...,
           [ 1.,  1.,  1., ..., nan, nan, nan],
           [ 1.,  1.,  1., ..., nan, nan, nan],
           [ 1.,  1.,  1., ..., nan, nan, nan]],
          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., ..., 165600., 169200., 172800.],
           [     0.,   3600.,   7200., ...,     nan,     nan,     nan],
           [     0.,   3600.,   7200., ..., 165600., 169200., 172800.],
           ...,
           [     0.,   3600.,   7200., ...,     nan,     nan,     nan],
           [     0.,   3600.,   7200., ...,     nan,     nan,     nan],
           [     0.,   3600.,   7200., ...,     nan,     nan,     nan]],
          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., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan]],
          shape=(3000, 49), dtype=float32)

  • lon
    (trajectory, time)
    float32
    19.2 19.21 19.25 ... nan nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    axis :

    X

    minval :

    19.011507

    maxval :

    20.040127

    array([[19.195274, 19.214897, 19.247133, ..., 19.567202, 19.565477,
            19.564339],
           [19.191893, 19.2157  , 19.232859, ...,       nan,       nan,
                  nan],
           [19.193327, 19.230707, 19.276842, ..., 19.55558 , 19.558638,
            19.560644],
           ...,
           [19.190422, 19.19865 , 19.21901 , ...,       nan,       nan,
                  nan],
           [19.187893, 19.22117 , 19.25207 , ...,       nan,       nan,
                  nan],
           [19.1907  , 19.229303, 19.229008, ...,       nan,       nan,
                  nan]], shape=(3000, 49), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    axis :

    Y

    minval :

    79.457535

    maxval :

    79.654564

    array([[79.59931 , 79.597626, 79.59428 , ..., 79.50743 , 79.50375 ,
            79.49495 ],
           [79.598854, 79.59528 , 79.58706 , ...,       nan,       nan,
                  nan],
           [79.59855 , 79.58983 , 79.586365, ..., 79.49775 , 79.49136 ,
            79.48704 ],
           ...,
           [79.59836 , 79.58963 , 79.58465 , ...,       nan,       nan,
                  nan],
           [79.59878 , 79.59089 , 79.582504, ...,       nan,       nan,
                  nan],
           [79.59873 , 79.595764, 79.589005, ...,       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 :

    -6.623588

    maxval :

    0.0

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

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

    dtype :

    <class 'numpy.float32'>

    units :

    1

    description :

    Elements are moved with this fraction of the current vector, in addition to currents and Stokes drift

    minval :

    1.0

    maxval :

    1.0

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

    array([[0.        , 0.        , 0.        , ..., 0.        , 0.        ,
            0.        ],
           [0.        , 0.        , 0.        , ...,        nan,        nan,
                   nan],
           [0.        , 0.        , 0.00571372, ..., 0.00572256, 0.00571701,
            0.0057243 ],
           ...,
           [0.        , 0.        , 0.        , ...,        nan,        nan,
                   nan],
           [0.        , 0.        , 0.        , ...,        nan,        nan,
                   nan],
           [0.        , 0.        , 0.        , ...,        nan,        nan,
                   nan]], shape=(3000, 49), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.32936665

    maxval :

    0.32936665

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

  • viscosity
    (trajectory, time)
    float32
    0.02169 0.1305 0.1308 ... nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    m2/s

    description :

    Kinematic viscosity of oil (emulsion)

    minval :

    0.021692334

    maxval :

    0.1376163

    array([[0.02169233, 0.13052559, 0.13082346, ..., 0.1276877 , 0.12641472,
            0.12663804],
           [0.02169233, 0.13028891, 0.13063821, ...,        nan,        nan,
                   nan],
           [0.02169233, 0.13073793, 0.13162836, ..., 0.12898938, 0.12938222,
            0.1288819 ],
           ...,
           [0.02169233, 0.13015871, 0.1310593 , ...,        nan,        nan,
                   nan],
           [0.02169233, 0.13021979, 0.1314924 , ...,        nan,        nan,
                   nan],
           [0.02169233, 0.13015749, 0.13090535, ...,        nan,        nan,
                   nan]], shape=(3000, 49), dtype=float32)

  • density
    (trajectory, time)
    float32
    988.1 996.1 996.2 ... nan nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    kg/m^3

    minval :

    988.1

    maxval :

    996.3721

    array([[988.1    , 996.14526, 996.155  , ..., 996.05084, 996.00775,
            996.0154 ],
           [988.1    , 996.13745, 996.149  , ...,       nan,       nan,
                  nan],
           [988.1    , 996.1522 , 996.1814 , ..., 996.0944 , 996.1075 ,
            996.0908 ],
           ...,
           [988.1    , 996.1332 , 996.1628 , ...,       nan,       nan,
                  nan],
           [988.1    , 996.1352 , 996.17694, ...,       nan,       nan,
                  nan],
           [988.1    , 996.1331 , 996.1577 , ...,       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., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan]],
          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., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan]],
          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., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan]],
          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., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan]],
          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., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan]],
          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., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ..., nan, nan, nan],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           ...,
           [ 1.,  1.,  1., ..., nan, nan, nan],
           [ 1.,  1.,  1., ..., nan, nan, nan],
           [ 1.,  1.,  1., ..., nan, nan, nan]],
          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., ...,  3.,  3.,  3.],
           [ 3.,  3.,  3., ..., nan, nan, nan],
           [ 3.,  3.,  3., ...,  3.,  3.,  3.],
           ...,
           [ 3.,  3.,  3., ..., nan, nan, nan],
           [ 3.,  3.,  3., ..., nan, nan, nan],
           [ 3.,  3.,  3., ..., nan, nan, nan]],
          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., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan]],
          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., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan]],
          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, ..., 0.001, 0.001, 0.001],
           [0.001, 0.001, 0.001, ...,   nan,   nan,   nan],
           [0.001, 0.001, 0.001, ..., 0.001, 0.001, 0.001],
           ...,
           [0.001, 0.001, 0.001, ...,   nan,   nan,   nan],
           [0.001, 0.001, 0.001, ...,   nan,   nan,   nan],
           [0.001, 0.001, 0.001, ...,   nan,   nan,   nan]],
          shape=(3000, 49), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    m

    minval :

    0.0

    maxval :

    0.0029854549

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

  • x_sea_water_velocity
    (trajectory, time)
    float32
    0.1118 0.03219 -0.04996 ... nan nan

    minval :

    -0.45580676

    maxval :

    0.39298648

    array([[ 0.11175027,  0.03218732, -0.04996448, ..., -0.01244314,
             0.01776481,  0.00437073],
           [-0.03608771, -0.04968661, -0.01910808, ...,         nan,
                    nan,         nan],
           [ 0.01541927, -0.02013632,  0.03673995, ..., -0.11148575,
            -0.02108256, -0.04441232],
           ...,
           [-0.03068552, -0.06208014,  0.01296356, ...,         nan,
                    nan,         nan],
           [-0.04378641,  0.01207774, -0.04730261, ...,         nan,
                    nan,         nan],
           [ 0.04418846,  0.00046961, -0.01643308, ...,         nan,
                    nan,         nan]], shape=(3000, 49), dtype=float32)

  • y_sea_water_velocity
    (trajectory, time)
    float32
    0.002295 0.04216 ... nan nan

    minval :

    -0.4602833

    maxval :

    0.5444253

    array([[ 0.00229533,  0.04216388, -0.00903554, ..., -0.11449368,
            -0.1842655 , -0.15512033],
           [-0.0466179 , -0.03893718,  0.07403764, ...,         nan,
                    nan,         nan],
           [-0.07291572,  0.03731424,  0.05009178, ..., -0.21274993,
            -0.26648   , -0.13118546],
           ...,
           [-0.16900107,  0.05004745,  0.00547197, ...,         nan,
                    nan,         nan],
           [-0.07035927,  0.01788829,  0.0106231 , ...,         nan,
                    nan,         nan],
           [-0.05766536, -0.03645582, -0.05158135, ...,         nan,
                    nan,         nan]], shape=(3000, 49), dtype=float32)

  • x_wind
    (trajectory, time)
    float32
    6.404 5.771 5.061 ... nan nan nan

    minval :

    -3.3859003

    maxval :

    7.7600255

    array([[6.404282  , 5.7709804 , 5.060549  , ..., 0.47823152, 2.2747216 ,
            2.2215564 ],
           [6.1761866 , 6.5879245 , 5.227361  , ...,        nan,        nan,
                   nan],
           [5.3795953 , 6.454175  , 5.454363  , ..., 2.0257015 , 1.8175057 ,
            1.923351  ],
           ...,
           [6.365642  , 5.8570457 , 4.4433804 , ...,        nan,        nan,
                   nan],
           [5.6977944 , 5.4843783 , 5.8889003 , ...,        nan,        nan,
                   nan],
           [5.8758965 , 5.106558  , 5.2857914 , ...,        nan,        nan,
                   nan]], shape=(3000, 49), dtype=float32)

  • y_wind
    (trajectory, time)
    float32
    -4.008 -5.707 -5.932 ... nan nan

    minval :

    -6.948589

    maxval :

    2.0303602

    array([[-4.007551  , -5.7074313 , -5.9323435 , ..., -0.6758151 ,
            -0.42365012,  0.28938526],
           [-4.988847  , -5.077519  , -5.1502013 , ...,         nan,
                    nan,         nan],
           [-5.485705  , -4.680029  , -4.1079016 , ..., -0.21449064,
            -0.42132393,  0.41449782],
           ...,
           [-5.319088  , -5.4226656 , -5.0266924 , ...,         nan,
                    nan,         nan],
           [-5.382487  , -5.744465  , -4.6711826 , ...,         nan,
                    nan,         nan],
           [-5.266737  , -4.938289  , -4.4471297 , ...,         nan,
                    nan,         nan]], shape=(3000, 49), dtype=float32)

  • sea_surface_height
    (trajectory, time)
    float32
    -0.2822 -0.3697 -0.4651 ... nan nan

    minval :

    -1.0802746

    maxval :

    -0.26634133

    array([[-0.28224668, -0.3697369 , -0.4651174 , ..., -0.39807087,
            -0.31386065, -0.34074962],
           [-0.28222752, -0.36936644, -0.46422958, ...,         nan,
                    nan,         nan],
           [-0.28217003, -0.36820462, -0.46197373, ..., -0.397043  ,
            -0.31042248, -0.33842304],
           ...,
           [-0.28218007, -0.36900187, -0.46446502, ...,         nan,
                    nan,         nan],
           [-0.2822666 , -0.36859223, -0.4622909 , ...,         nan,
                    nan,         nan],
           [-0.282226  , -0.36915126, -0.46482435, ...,         nan,
                    nan,         nan]], shape=(3000, 49), dtype=float32)

  • upward_sea_water_velocity
    (trajectory, time)
    float32
    -2.447e-05 -2.605e-05 ... nan nan

    minval :

    -8.7226355e-05

    maxval :

    9.560622e-05

    array([[-2.4473398e-05, -2.6049287e-05, -2.7191096e-05, ...,
             4.0343704e-05, -4.8339402e-06, -1.4261548e-05],
           [-2.4520787e-05, -2.5572150e-05, -2.5880741e-05, ...,
                       nan,            nan,            nan],
           [-2.4440142e-05, -2.3937479e-05, -2.4612675e-05, ...,
             4.4776549e-05, -2.0768236e-06, -1.2405471e-05],
           ...,
           [-2.4505112e-05, -2.5170983e-05, -2.6077720e-05, ...,
                       nan,            nan,            nan],
           [-2.4632951e-05, -2.4521760e-05, -2.4558069e-05, ...,
                       nan,            nan,            nan],
           [-2.4542131e-05, -2.5122821e-05, -2.6300988e-05, ...,
                       nan,            nan,            nan]],
          shape=(3000, 49), dtype=float32)

  • sea_surface_wave_significant_height
    (trajectory, time)
    float32
    1.404 1.621 1.496 ... nan nan nan

    minval :

    5.2378186e-06

    maxval :

    2.3443701

    array([[1.4040521 , 1.6206229 , 1.4957256 , ..., 0.01686161, 0.1317044 ,
            0.1234688 ],
           [1.5506332 , 1.701876  , 1.324707  , ...,        nan,        nan,
                   nan],
           [1.4522119 , 1.5635524 , 1.1469733 , ..., 0.10207702, 0.08562868,
            0.09522875],
           ...,
           [1.6928269 , 1.567273  , 1.1072773 , ...,        nan,        nan,
                   nan],
           [1.5113262 , 1.5517012 , 1.3898777 , ...,        nan,        nan,
                   nan],
           [1.5317109 , 1.2414054 , 1.173827  , ...,        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., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan]],
          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., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan]],
          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., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           ...,
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           [ 0.,  0.,  0., ..., nan, nan, nan]],
          shape=(3000, 49), dtype=float32)

  • sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment
    (trajectory, time)
    float32
    6.455 6.935 6.663 ... nan nan nan

    minval :

    0.012468231

    maxval :

    8.341471

    array([[6.455371 , 6.9353848, 6.6627812, ..., 0.7074229, 1.9771049,
            1.914292 ],
           [6.7839727, 7.1071177, 6.2703166, ...,       nan,       nan,
                  nan],
           [6.5651484, 6.8121753, 5.8345346, ..., 1.7405783, 1.5941869,
            1.6811777],
           ...,
           [7.088198 , 6.820276 , 5.732681 , ...,       nan,       nan,
                  nan],
           [6.6974382, 6.786309 , 6.422703 , ...,       nan,       nan,
                  nan],
           [6.742454 , 6.0699673, 5.9024405, ...,       nan,       nan,
                  nan]], shape=(3000, 49), dtype=float32)

  • sea_ice_area_fraction
    (trajectory, time)
    float32
    0.05874 0.06037 0.0648 ... nan nan

    minval :

    0.011770496

    maxval :

    0.2377118

    array([[0.05873717, 0.06037185, 0.06480418, ..., 0.09214708, 0.09782857,
            0.10612968],
           [0.0582911 , 0.06220036, 0.06576943, ...,        nan,        nan,
                   nan],
           [0.05882456, 0.06906744, 0.07643416, ..., 0.11038814, 0.12183934,
            0.12654322],
           ...,
           [0.05830391, 0.05915367, 0.06214987, ...,        nan,        nan,
                   nan],
           [0.05743437, 0.06591115, 0.07261294, ...,        nan,        nan,
                   nan],
           [0.05810613, 0.06455012, 0.06418598, ...,        nan,        nan,
                   nan]], shape=(3000, 49), dtype=float32)

  • sea_ice_x_velocity
    (trajectory, time)
    float32
    0.03905 0.04502 0.045 ... nan nan

    minval :

    -0.15717354

    maxval :

    0.10390354

    array([[ 0.0390547 ,  0.04502173,  0.04500297, ...,  0.00650267,
             0.01160827,  0.00844864],
           [ 0.03878018,  0.0438885 ,  0.04360198, ...,         nan,
                    nan,         nan],
           [ 0.03856081,  0.04060725,  0.03981784, ..., -0.00810718,
            -0.00552744, -0.00640063],
           ...,
           [ 0.03846107,  0.04237122,  0.04309984, ...,         nan,
                    nan,         nan],
           [ 0.03876386,  0.04164328,  0.04048008, ...,         nan,
                    nan,         nan],
           [ 0.03870915,  0.04314725,  0.0447427 , ...,         nan,
                    nan,         nan]], shape=(3000, 49), dtype=float32)

  • sea_ice_y_velocity
    (trajectory, time)
    float32
    -0.02023 -0.03676 ... nan nan

    minval :

    -0.124134034

    maxval :

    0.17111596

    array([[-0.02023249, -0.03675835, -0.05269475, ...,  0.00481401,
            -0.01084783, -0.00827537],
           [-0.01979725, -0.03880931, -0.05526684, ...,         nan,
                    nan,         nan],
           [-0.02036427, -0.04656222, -0.0657973 , ...,  0.01369401,
            -0.00234134, -0.00455575],
           ...,
           [-0.01983754, -0.04328741, -0.05584898, ...,         nan,
                    nan,         nan],
           [-0.0189441 , -0.0440055 , -0.06232903, ...,         nan,
                    nan,         nan],
           [-0.01961701, -0.04196801, -0.05265203, ...,         nan,
                    nan,         nan]], shape=(3000, 49), dtype=float32)

  • sea_water_temperature
    (trajectory, time)
    float32
    -0.1861 -0.2096 -0.2504 ... nan nan

    minval :

    -0.5648154

    maxval :

    273.8252

    array([[-1.8605775e-01, -2.0964758e-01, -2.5039744e-01, ...,
            -3.6159746e-02, -6.5494897e-03,  2.7312695e+02],
           [-1.8328544e-01, -2.1933509e-01, -2.6008898e-01, ...,
                       nan,            nan,            nan],
           [-1.8600745e-01, -2.5387120e-01, -3.0623031e-01, ...,
            -1.4342694e-01, -1.3954557e-01,  2.7301373e+02],
           ...,
           [-1.8301073e-01, -1.9866085e-01, -2.4111220e-01, ...,
                       nan,            nan,            nan],
           [-1.7852506e-01, -2.3739566e-01, -2.9226416e-01, ...,
                       nan,            nan,            nan],
           [-1.8218377e-01, -2.3035160e-01, -2.5311884e-01, ...,
                       nan,            nan,            nan]],
          shape=(3000, 49), dtype=float32)

  • sea_water_salinity
    (trajectory, time)
    float32
    33.61 33.6 33.58 ... nan nan nan

    minval :

    33.004204

    maxval :

    33.64183

    array([[33.61051 , 33.603012, 33.583435, ..., 33.393116, 33.39457 ,
            33.391407],
           [33.612663, 33.602425, 33.589764, ...,       nan,       nan,
                  nan],
           [33.61175 , 33.592052, 33.562107, ..., 33.388653, 33.388905,
            33.388405],
           ...,
           [33.61359 , 33.612495, 33.59734 , ...,       nan,       nan,
                  nan],
           [33.615215, 33.5983  , 33.574894, ...,       nan,       nan,
                  nan],
           [33.61342 , 33.59389 , 33.593224, ...,       nan,       nan,
                  nan]], shape=(3000, 49), dtype=float32)

  • sea_floor_depth_below_sea_level
    (trajectory, time)
    float32
    170.4 164.5 154.6 ... nan nan nan

    minval :

    31.612967

    maxval :

    196.33914

    array([[170.42508 , 164.4656  , 154.57133 , ...,  93.585396,  91.83893 ,
             88.9146  ],
           [170.93388 , 162.4598  , 149.1992  , ...,        nan,        nan,
                   nan],
           [170.35571 , 153.75998 , 141.26956 , ...,  86.73605 ,  81.57847 ,
             79.50233 ],
           ...,
           [170.93738 , 157.27959 , 147.03888 , ...,        nan,        nan,
                   nan],
           [171.87785 , 156.52078 , 139.81734 , ...,        nan,        nan,
                   nan],
           [171.14175 , 159.7099  , 152.68683 , ...,        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.024936968

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

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

    minval :

    0.0

    maxval :

    0.0

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

• Attributes: (164)

  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-06-26 09:55:42.110216

  date_created :

  2025-06-26T09:55:42.110233

  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-06-19 09:55:37.704451

  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: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-06-19 09:55:37.704451

  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-06-26 09:55:42.473813

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

  ['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-06-21 09:55:37.704451

  time_coverage_duration :

  P2DT0H0M0S

  time_coverage_resolution :

  P0DT0H30M0S

  performance :

  -------------------- Reader performance: -------------------- global_landmask 0:00:01.6 total 0:00:00.0 preparing 0:00:01.6 reading 0:00:00.0 masking -------------------- https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be 0:02:03.6 total 0:00:00.0 preparing 0:02:02.8 reading 0:00:01.3 interpolation 0:00:00.1 interpolation_time 0:00:00.7 rotating vectors 0:00:00.0 masking -------------------- Performance: 2:25.1 total time 4.6 configuration 0.3 preparing main loop 0.0 moving elements to ocean 2:20.0 main loop 5.7 updating elements 0.0 oil weathering 0.0 updating viscosities 0.0 updating densities 5.3 vertical mixing 0.0 cleaning up --------------------

  geospatial_bounds_crs :

  EPSG:4326

  geospatial_bounds_vertical_crs :

  EPSG:5831

  geospatial_lat_min :

  79.457535

  geospatial_lat_max :

  79.654564

  geospatial_lat_units :

  degrees_north

  geospatial_lat_resolution :

  point

  geospatial_lon_min :

  19.011507

  geospatial_lon_max :

  20.040127

  geospatial_lon_units :

  degrees_east

  geospatial_lon_resolution :

  point

  runtime :

  0:02:25.125074

  geospatial_vertical_min :

  -6.623588

  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.6  total
 0:00:00.0  preparing
 0:00:01.6  reading
 0:00:00.0  masking
--------------------
https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
 0:02:03.6  total
 0:00:00.0  preparing
 0:02:02.8  reading
 0:00:01.3  interpolation
 0:00:00.1  interpolation_time
 0:00:00.7  rotating vectors
 0:00:00.0  masking
--------------------
Performance:
 2:25.1 total time
    4.6 configuration
    0.3 preparing main loop
      0.0 moving elements to ocean
 2:20.0 main loop
      5.7 updating elements
        0.0 oil weathering
          0.0 updating viscosities
          0.0 updating densities
        5.3 vertical mixing
    0.0 cleaning up
--------------------
===========================
Model:  OpenOil     (OpenDrift version 1.14.2)
        1318 active Oil particles  (1682 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
  upward_sea_water_velocity
  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-06-19 09:55:37.704451 UTC
        Present: 2025-06-21 09:55:37.704451 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)
09:58:13 INFO    opendrift:4607: Saving animation to /root/project/docs/source/gallery/animations/example_oil_ice_0.gif...
09:58:53 INFO    opendrift:3055: Time to make animation: 0:00:50.288175

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-06-19 09:55 to 2025-06-21 09:55 UTC (49 steps)'}>, <Figure size 610.745x1100 with 2 Axes>)