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)

15:18:17 INFO    opendrift:509: OpenDriftSimulation initialised (version 1.13.1 / v1.13.1-12-gca7f755)

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

15:18:17 INFO    opendrift.models.openoil.openoil:1701: Oil type not specified, using default: GENERIC BUNKER C
15:18:17 INFO    opendrift.models.openoil.adios.dirjs:86: Querying ADIOS database for oil: GENERIC BUNKER C
15:18:17 INFO    opendrift.models.openoil.openoil:1710: Using density 988.1 and viscosity 0.021692333877975645 of oiltype GENERIC BUNKER C
15:18:17 INFO    opendrift.models.basemodel.environment:206: Adding a global landmask from GSHHG
15:18:21 INFO    opendrift.models.basemodel.environment:229: Fallback values will be used for the following variables which have no readers:
15:18:21 INFO    opendrift.models.basemodel.environment:232:    sea_surface_height: 0.000000
15:18:21 INFO    opendrift.models.basemodel.environment:232:    upward_sea_water_velocity: 0.000000
15:18:21 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_significant_height: 0.000000
15:18:21 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
15:18:21 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
15:18:21 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
15:18:21 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
15:18:21 INFO    opendrift.models.basemodel.environment:232:    sea_ice_area_fraction: 0.000000
15:18:21 INFO    opendrift.models.basemodel.environment:232:    sea_ice_x_velocity: 0.000000
15:18:21 INFO    opendrift.models.basemodel.environment:232:    sea_ice_y_velocity: 0.000000
15:18:21 INFO    opendrift.models.basemodel.environment:232:    sea_water_temperature: 10.000000
15:18:21 INFO    opendrift.models.basemodel.environment:232:    sea_water_salinity: 34.000000
15:18:21 INFO    opendrift.models.basemodel.environment:232:    sea_floor_depth_below_sea_level: 10000.000000
15:18:21 INFO    opendrift.models.basemodel.environment:232:    ocean_vertical_diffusivity: 0.020000
15:18:21 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)

15:18:21 INFO    opendrift:919: Using existing reader for land_binary_mask
15:18:21 INFO    opendrift:930: All points are in ocean
15:18:21 INFO    opendrift.models.openoil.openoil:685: Oil-water surface tension is 0.035935 Nm
15:18:21 INFO    opendrift.models.openoil.openoil:698: Max water fraction not available for GENERIC BUNKER C, using default
15:18:21 INFO    opendrift:2056: 2025-03-27 15:18:17.209134 - step 1 of 96 - 3000 active elements (0 deactivated)
15:18:21 INFO    opendrift.readers:61: Opening file with xr.open_dataset
15:20:23 INFO    opendrift.readers.reader_netCDF_CF_generic:332: Detected dimensions: {'x': 'X', 'y': 'Y', 'z': 'depth', 'time': 'time'}
15:20:23 INFO    opendrift.readers.basereader:176: Variable x_sea_water_velocity will be rotated from eastward_sea_water_velocity
15:20:23 INFO    opendrift.readers.basereader:176: Variable y_sea_water_velocity will be rotated from northward_sea_water_velocity
15:20:27 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0003610238 to -3.6268227e-07) found for ocean_vertical_diffusivity, replacing with NaN
15:20:27 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:20:30 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0027372167 to -0.00033993289) found for ocean_vertical_diffusivity, replacing with NaN
15:20:30 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:20:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:30 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:30 INFO    opendrift:2056: 2025-03-27 15:48:17.209134 - step 2 of 96 - 3000 active elements (0 deactivated)
15:20:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:30 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:30 INFO    opendrift:2056: 2025-03-27 16:18:17.209134 - step 3 of 96 - 3000 active elements (0 deactivated)
15:20:33 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.02604523 to -0.0030932112) found for ocean_vertical_diffusivity, replacing with NaN
15:20:33 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:20:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:33 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:33 INFO    opendrift:2056: 2025-03-27 16:48:17.209134 - step 4 of 96 - 3000 active elements (0 deactivated)
15:20:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:33 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:33 INFO    opendrift:2056: 2025-03-27 17:18:17.209134 - step 5 of 96 - 3000 active elements (0 deactivated)
15:20:37 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.05120443 to -0.006118676) found for ocean_vertical_diffusivity, replacing with NaN
15:20:37 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:20:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:37 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:37 INFO    opendrift:2056: 2025-03-27 17:48:17.209134 - step 6 of 96 - 3000 active elements (0 deactivated)
15:20:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:37 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:37 INFO    opendrift:2056: 2025-03-27 18:18:17.209134 - step 7 of 96 - 3000 active elements (0 deactivated)
15:20:40 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.039279774 to -0.0064136153) found for ocean_vertical_diffusivity, replacing with NaN
15:20:40 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:20:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:40 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:40 INFO    opendrift:2056: 2025-03-27 18:48:17.209134 - step 8 of 96 - 3000 active elements (0 deactivated)
15:20:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:42 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:43 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:43 INFO    opendrift:2056: 2025-03-27 19:18:17.209134 - step 9 of 96 - 2999 active elements (1 deactivated)
15:20:47 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.016999327 to -0.004190925) found for ocean_vertical_diffusivity, replacing with NaN
15:20:47 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:20:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:47 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:47 INFO    opendrift:2056: 2025-03-27 19:48:17.209134 - step 10 of 96 - 2999 active elements (1 deactivated)
15:20:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:47 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:47 INFO    opendrift:2056: 2025-03-27 20:18:17.209134 - step 11 of 96 - 2997 active elements (3 deactivated)
15:20:50 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.009858792 to -0.003799685) found for ocean_vertical_diffusivity, replacing with NaN
15:20:50 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:20:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:50 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:50 INFO    opendrift:2056: 2025-03-27 20:48:17.209134 - step 12 of 96 - 2995 active elements (5 deactivated)
15:20:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:50 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:50 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:50 INFO    opendrift:2056: 2025-03-27 21:18:17.209134 - step 13 of 96 - 2991 active elements (9 deactivated)
15:20:54 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.008204452 to -0.0034688797) found for ocean_vertical_diffusivity, replacing with NaN
15:20:54 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:20:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:54 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:54 INFO    opendrift:2056: 2025-03-27 21:48:17.209134 - step 14 of 96 - 2979 active elements (21 deactivated)
15:20:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:54 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:54 INFO    opendrift:2056: 2025-03-27 22:18:17.209134 - step 15 of 96 - 2972 active elements (28 deactivated)
15:20:57 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.012582601 to -0.0042569116) found for ocean_vertical_diffusivity, replacing with NaN
15:20:57 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:20:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:57 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:57 INFO    opendrift:2056: 2025-03-27 22:48:17.209134 - step 16 of 96 - 2961 active elements (39 deactivated)
15:20:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:57 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:20:58 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:20:58 INFO    opendrift:2056: 2025-03-27 23:18:17.209134 - step 17 of 96 - 2944 active elements (56 deactivated)
15:21:01 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.01044945 to -0.0026618347) found for ocean_vertical_diffusivity, replacing with NaN
15:21:01 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:01 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:01 INFO    opendrift:2056: 2025-03-27 23:48:17.209134 - step 18 of 96 - 2938 active elements (62 deactivated)
15:21:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:01 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:01 INFO    opendrift:2056: 2025-03-28 00:18:17.209134 - step 19 of 96 - 2929 active elements (71 deactivated)
15:21:05 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.009848064 to -0.0016612484) found for ocean_vertical_diffusivity, replacing with NaN
15:21:05 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:05 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:05 INFO    opendrift:2056: 2025-03-28 00:48:17.209134 - step 20 of 96 - 2918 active elements (82 deactivated)
15:21:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:05 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:05 INFO    opendrift:2056: 2025-03-28 01:18:17.209134 - step 21 of 96 - 2915 active elements (85 deactivated)
15:21:08 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.009609729 to -0.0006919717) found for ocean_vertical_diffusivity, replacing with NaN
15:21:08 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:08 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:08 INFO    opendrift:2056: 2025-03-28 01:48:17.209134 - step 22 of 96 - 2911 active elements (89 deactivated)
15:21:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:08 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:08 INFO    opendrift:2056: 2025-03-28 02:18:17.209134 - step 23 of 96 - 2904 active elements (96 deactivated)
15:21:12 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.002400813 to -3.9718412e-05) found for ocean_vertical_diffusivity, replacing with NaN
15:21:12 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:12 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:12 INFO    opendrift:2056: 2025-03-28 02:48:17.209134 - step 24 of 96 - 2896 active elements (104 deactivated)
15:21:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:12 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:12 INFO    opendrift:2056: 2025-03-28 03:18:17.209134 - step 25 of 96 - 2892 active elements (108 deactivated)
15:21:15 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0027180186 to -3.6463443e-06) found for ocean_vertical_diffusivity, replacing with NaN
15:21:15 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:15 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:15 INFO    opendrift:2056: 2025-03-28 03:48:17.209134 - step 26 of 96 - 2879 active elements (121 deactivated)
15:21:15 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:15 INFO    opendrift:2056: 2025-03-28 04:18:17.209134 - step 27 of 96 - 2874 active elements (126 deactivated)
15:21:19 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.017467154 to -0.0029561822) found for ocean_vertical_diffusivity, replacing with NaN
15:21:19 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:19 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:19 INFO    opendrift:2056: 2025-03-28 04:48:17.209134 - step 28 of 96 - 2862 active elements (138 deactivated)
15:21:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:19 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:19 INFO    opendrift:2056: 2025-03-28 05:18:17.209134 - step 29 of 96 - 2857 active elements (143 deactivated)
15:21:22 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.047463994 to -0.0066782907) found for ocean_vertical_diffusivity, replacing with NaN
15:21:22 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:22 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:22 INFO    opendrift:2056: 2025-03-28 05:48:17.209134 - step 30 of 96 - 2850 active elements (150 deactivated)
15:21:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:23 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:23 INFO    opendrift:2056: 2025-03-28 06:18:17.209134 - step 31 of 96 - 2846 active elements (154 deactivated)
15:21:26 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.04395428 to -0.007816557) found for ocean_vertical_diffusivity, replacing with NaN
15:21:26 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:26 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:26 INFO    opendrift:2056: 2025-03-28 06:48:17.209134 - step 32 of 96 - 2841 active elements (159 deactivated)
15:21:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:26 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:26 INFO    opendrift:2056: 2025-03-28 07:18:17.209134 - step 33 of 96 - 2837 active elements (163 deactivated)
15:21:30 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.019060358 to -0.002278598) found for ocean_vertical_diffusivity, replacing with NaN
15:21:30 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:30 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:30 INFO    opendrift:2056: 2025-03-28 07:48:17.209134 - step 34 of 96 - 2835 active elements (165 deactivated)
15:21:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:30 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:30 INFO    opendrift:2056: 2025-03-28 08:18:17.209134 - step 35 of 96 - 2835 active elements (165 deactivated)
15:21:33 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.013104989 to -0.0017008776) found for ocean_vertical_diffusivity, replacing with NaN
15:21:33 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:33 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:33 INFO    opendrift:2056: 2025-03-28 08:48:17.209134 - step 36 of 96 - 2835 active elements (165 deactivated)
15:21:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:33 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:34 INFO    opendrift:2056: 2025-03-28 09:18:17.209134 - step 37 of 96 - 2834 active elements (166 deactivated)
15:21:37 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.011059331 to -0.0012137312) found for ocean_vertical_diffusivity, replacing with NaN
15:21:37 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:37 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:37 INFO    opendrift:2056: 2025-03-28 09:48:17.209134 - step 38 of 96 - 2833 active elements (167 deactivated)
15:21:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:37 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:37 INFO    opendrift:2056: 2025-03-28 10:18:17.209134 - step 39 of 96 - 2831 active elements (169 deactivated)
15:21:40 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.011342908 to -0.0014484846) found for ocean_vertical_diffusivity, replacing with NaN
15:21:40 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:40 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:41 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:41 INFO    opendrift:2056: 2025-03-28 10:48:17.209134 - step 40 of 96 - 2828 active elements (172 deactivated)
15:21:41 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:41 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:41 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:41 INFO    opendrift:2056: 2025-03-28 11:18:17.209134 - step 41 of 96 - 2827 active elements (173 deactivated)
15:21:44 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.011752667 to -0.0010711044) found for ocean_vertical_diffusivity, replacing with NaN
15:21:44 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:44 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:44 INFO    opendrift:2056: 2025-03-28 11:48:17.209134 - step 42 of 96 - 2825 active elements (175 deactivated)
15:21:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:44 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:44 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:44 INFO    opendrift:2056: 2025-03-28 12:18:17.209134 - step 43 of 96 - 2823 active elements (177 deactivated)
15:21:47 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.01379464 to -0.0011370624) found for ocean_vertical_diffusivity, replacing with NaN
15:21:47 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:47 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:47 INFO    opendrift:2056: 2025-03-28 12:48:17.209134 - step 44 of 96 - 2821 active elements (179 deactivated)
15:21:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:47 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:48 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:48 INFO    opendrift:2056: 2025-03-28 13:18:17.209134 - step 45 of 96 - 2815 active elements (185 deactivated)
15:21:51 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.012695961 to -0.00019790627) found for ocean_vertical_diffusivity, replacing with NaN
15:21:51 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:51 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:51 INFO    opendrift:2056: 2025-03-28 13:48:17.209134 - step 46 of 96 - 2813 active elements (187 deactivated)
15:21:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:51 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:51 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:51 INFO    opendrift:2056: 2025-03-28 14:18:17.209134 - step 47 of 96 - 2807 active elements (193 deactivated)
15:21:54 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0015077526 to -1.3774286e-05) found for ocean_vertical_diffusivity, replacing with NaN
15:21:54 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:54 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:54 INFO    opendrift:2056: 2025-03-28 14:48:17.209134 - step 48 of 96 - 2801 active elements (199 deactivated)
15:21:54 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:21:55 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:55 INFO    opendrift:2056: 2025-03-28 15:18:17.209134 - step 49 of 96 - 2794 active elements (206 deactivated)
15:21:58 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0022711365 to -2.1834626e-06) found for ocean_vertical_diffusivity, replacing with NaN
15:21:58 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:21:58 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:58 INFO    opendrift:2056: 2025-03-28 15:48:17.209134 - step 50 of 96 - 2781 active elements (219 deactivated)
15:21:58 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:21:58 INFO    opendrift:2056: 2025-03-28 16:18:17.209134 - step 51 of 96 - 2764 active elements (236 deactivated)
15:22:01 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.012082136 to -0.0015956464) found for ocean_vertical_diffusivity, replacing with NaN
15:22:01 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:01 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:01 INFO    opendrift:2056: 2025-03-28 16:48:17.209134 - step 52 of 96 - 2749 active elements (251 deactivated)
15:22:01 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:01 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:01 INFO    opendrift:2056: 2025-03-28 17:18:17.209134 - step 53 of 96 - 2735 active elements (265 deactivated)
15:22:04 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.05509226 to -0.006366697) found for ocean_vertical_diffusivity, replacing with NaN
15:22:04 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:04 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:04 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:05 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:05 INFO    opendrift:2056: 2025-03-28 17:48:17.209134 - step 54 of 96 - 2707 active elements (293 deactivated)
15:22:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:05 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:05 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:05 INFO    opendrift:2056: 2025-03-28 18:18:17.209134 - step 55 of 96 - 2693 active elements (307 deactivated)
15:22:08 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.07306447 to -0.0025496795) found for ocean_vertical_diffusivity, replacing with NaN
15:22:08 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:08 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:08 INFO    opendrift:2056: 2025-03-28 18:48:17.209134 - step 56 of 96 - 2672 active elements (328 deactivated)
15:22:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:08 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:08 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:08 INFO    opendrift:2056: 2025-03-28 19:18:17.209134 - step 57 of 96 - 2662 active elements (338 deactivated)
15:22:12 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.04268506 to -0.0014226226) found for ocean_vertical_diffusivity, replacing with NaN
15:22:12 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:12 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:12 INFO    opendrift:2056: 2025-03-28 19:48:17.209134 - step 58 of 96 - 2640 active elements (360 deactivated)
15:22:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:12 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:12 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:12 INFO    opendrift:2056: 2025-03-28 20:18:17.209134 - step 59 of 96 - 2614 active elements (386 deactivated)
15:22:15 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.025209092 to -0.0035906238) found for ocean_vertical_diffusivity, replacing with NaN
15:22:15 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:15 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:15 INFO    opendrift:2056: 2025-03-28 20:48:17.209134 - step 60 of 96 - 2597 active elements (403 deactivated)
15:22:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:15 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:16 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:16 INFO    opendrift:2056: 2025-03-28 21:18:17.209134 - step 61 of 96 - 2572 active elements (428 deactivated)
15:22:19 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.01700509 to -0.0012391782) found for ocean_vertical_diffusivity, replacing with NaN
15:22:19 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:19 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:19 INFO    opendrift:2056: 2025-03-28 21:48:17.209134 - step 62 of 96 - 2542 active elements (458 deactivated)
15:22:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:19 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:19 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:19 INFO    opendrift:2056: 2025-03-28 22:18:17.209134 - step 63 of 96 - 2514 active elements (486 deactivated)
15:22:22 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.014875397 to -0.00093947904) found for ocean_vertical_diffusivity, replacing with NaN
15:22:22 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:22 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:23 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:23 INFO    opendrift:2056: 2025-03-28 22:48:17.209134 - step 64 of 96 - 2491 active elements (509 deactivated)
15:22:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:23 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:23 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:23 INFO    opendrift:2056: 2025-03-28 23:18:17.209134 - step 65 of 96 - 2460 active elements (540 deactivated)
15:22:26 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.015767789 to -0.0014557505) found for ocean_vertical_diffusivity, replacing with NaN
15:22:26 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:26 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:26 INFO    opendrift:2056: 2025-03-28 23:48:17.209134 - step 66 of 96 - 2423 active elements (577 deactivated)
15:22:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:26 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:27 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:27 INFO    opendrift:2056: 2025-03-29 00:18:17.209134 - step 67 of 96 - 2397 active elements (603 deactivated)
15:22:30 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.013832853 to -0.0017330198) found for ocean_vertical_diffusivity, replacing with NaN
15:22:30 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:30 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:30 INFO    opendrift:2056: 2025-03-29 00:48:17.209134 - step 68 of 96 - 2364 active elements (636 deactivated)
15:22:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:30 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:30 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:30 INFO    opendrift:2056: 2025-03-29 01:18:17.209134 - step 69 of 96 - 2336 active elements (664 deactivated)
15:22:33 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.02184638 to -0.001083417) found for ocean_vertical_diffusivity, replacing with NaN
15:22:33 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:33 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:33 INFO    opendrift:2056: 2025-03-29 01:48:17.209134 - step 70 of 96 - 2316 active elements (684 deactivated)
15:22:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:33 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:34 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:34 INFO    opendrift:2056: 2025-03-29 02:18:17.209134 - step 71 of 96 - 2296 active elements (704 deactivated)
15:22:37 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.010133364 to -1.1579495e-05) found for ocean_vertical_diffusivity, replacing with NaN
15:22:37 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:37 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:37 INFO    opendrift:2056: 2025-03-29 02:48:17.209134 - step 72 of 96 - 2268 active elements (732 deactivated)
15:22:37 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:37 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:37 INFO    opendrift:2056: 2025-03-29 03:18:17.209134 - step 73 of 96 - 2246 active elements (754 deactivated)
15:22:41 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0010443755 to -1.7887218e-06) found for ocean_vertical_diffusivity, replacing with NaN
15:22:41 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:41 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:41 INFO    opendrift:2056: 2025-03-29 03:48:17.209134 - step 74 of 96 - 2233 active elements (767 deactivated)
15:22:41 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:41 INFO    opendrift:2056: 2025-03-29 04:18:17.209134 - step 75 of 96 - 2214 active elements (786 deactivated)
15:22:45 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.006274988 to -1.791827e-05) found for ocean_vertical_diffusivity, replacing with NaN
15:22:45 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:45 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:45 INFO    opendrift:2056: 2025-03-29 04:48:17.209134 - step 76 of 96 - 2203 active elements (797 deactivated)
15:22:45 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:45 INFO    opendrift:2056: 2025-03-29 05:18:17.209134 - step 77 of 96 - 2188 active elements (812 deactivated)
15:22:48 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.040637493 to -0.0057002828) found for ocean_vertical_diffusivity, replacing with NaN
15:22:48 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:48 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:48 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:48 INFO    opendrift:2056: 2025-03-29 05:48:17.209134 - step 78 of 96 - 2179 active elements (821 deactivated)
15:22:48 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:48 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:48 INFO    opendrift:2056: 2025-03-29 06:18:17.209134 - step 79 of 96 - 2167 active elements (833 deactivated)
15:22:52 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.084265284 to -0.002408678) found for ocean_vertical_diffusivity, replacing with NaN
15:22:52 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:52 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:52 INFO    opendrift:2056: 2025-03-29 06:48:17.209134 - step 80 of 96 - 2160 active elements (840 deactivated)
15:22:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:52 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:52 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:52 INFO    opendrift:2056: 2025-03-29 07:18:17.209134 - step 81 of 96 - 2155 active elements (845 deactivated)
15:22:55 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.0764121 to -0.0012011089) found for ocean_vertical_diffusivity, replacing with NaN
15:22:55 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:55 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:56 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:56 INFO    opendrift:2056: 2025-03-29 07:48:17.209134 - step 82 of 96 - 2154 active elements (846 deactivated)
15:22:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:56 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:56 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:56 INFO    opendrift:2056: 2025-03-29 08:18:17.209134 - step 83 of 96 - 2150 active elements (850 deactivated)
15:22:59 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.026207607 to -0.0046728775) found for ocean_vertical_diffusivity, replacing with NaN
15:22:59 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:22:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:59 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:59 INFO    opendrift:2056: 2025-03-29 08:48:17.209134 - step 84 of 96 - 2147 active elements (853 deactivated)
15:22:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:59 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:22:59 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:22:59 INFO    opendrift:2056: 2025-03-29 09:18:17.209134 - step 85 of 96 - 2145 active elements (855 deactivated)
15:23:03 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.01353099 to -0.0008871437) found for ocean_vertical_diffusivity, replacing with NaN
15:23:03 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:23:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:03 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:23:03 INFO    opendrift:2056: 2025-03-29 09:48:17.209134 - step 86 of 96 - 2144 active elements (856 deactivated)
15:23:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:03 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:03 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:23:03 INFO    opendrift:2056: 2025-03-29 10:18:17.209134 - step 87 of 96 - 2144 active elements (856 deactivated)
15:23:06 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.011299693 to -0.0003592298) found for ocean_vertical_diffusivity, replacing with NaN
15:23:06 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:23:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:06 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:23:06 INFO    opendrift:2056: 2025-03-29 10:48:17.209134 - step 88 of 96 - 2139 active elements (861 deactivated)
15:23:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:06 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:06 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:23:06 INFO    opendrift:2056: 2025-03-29 11:18:17.209134 - step 89 of 96 - 2139 active elements (861 deactivated)
15:23:10 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.010986348 to -0.00032930812) found for ocean_vertical_diffusivity, replacing with NaN
15:23:10 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:23:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:10 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:23:10 INFO    opendrift:2056: 2025-03-29 11:48:17.209134 - step 90 of 96 - 2138 active elements (862 deactivated)
15:23:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:10 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:10 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:23:10 INFO    opendrift:2056: 2025-03-29 12:18:17.209134 - step 91 of 96 - 2137 active elements (863 deactivated)
15:23:13 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.010658032 to -0.0007045887) found for ocean_vertical_diffusivity, replacing with NaN
15:23:13 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:23:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:13 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:23:13 INFO    opendrift:2056: 2025-03-29 12:48:17.209134 - step 92 of 96 - 2137 active elements (863 deactivated)
15:23:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:13 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:13 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:23:13 INFO    opendrift:2056: 2025-03-29 13:18:17.209134 - step 93 of 96 - 2137 active elements (863 deactivated)
15:23:17 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.030265233 to -0.00043276537) found for ocean_vertical_diffusivity, replacing with NaN
15:23:17 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:23:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:17 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:23:17 INFO    opendrift:2056: 2025-03-29 13:48:17.209134 - step 94 of 96 - 2137 active elements (863 deactivated)
15:23:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:17 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:17 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:23:17 INFO    opendrift:2056: 2025-03-29 14:18:17.209134 - step 95 of 96 - 2136 active elements (864 deactivated)
15:23:20 WARNING opendrift.readers.basereader.variables:654: Invalid values (-0.023258379 to -2.6015156e-05) found for ocean_vertical_diffusivity, replacing with NaN
15:23:20 WARNING opendrift.readers.basereader.variables:657: (allowed range: [0, 1])
15:23:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:20 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice
15:23:20 INFO    opendrift:2056: 2025-03-29 14:48:17.209134 - step 96 of 96 - 2135 active elements (865 deactivated)
15:23:20 WARNING opendrift.readers.interpolation.interpolators:118: Only NaNs input to linearNDFast - returning
15:23:20 INFO    opendrift.models.openoil.openoil:1182: Ice concentration above 30%, using Nordam scheme for advection in ice

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

  • time
    (time)
    datetime64[ns]
    2025-03-27T15:18:17.209134 ... 2...

    standard_name :

    time

    long_name :

    time

    dtype :

    <class 'numpy.float64'>

    array(['2025-03-27T15:18:17.209134000', '2025-03-27T16:18:17.209134000',
           '2025-03-27T17:18:17.209134000', '2025-03-27T18:18:17.209134000',
           '2025-03-27T19:18:17.209134000', '2025-03-27T20:18:17.209134000',
           '2025-03-27T21:18:17.209134000', '2025-03-27T22:18:17.209134000',
           '2025-03-27T23:18:17.209134000', '2025-03-28T00:18:17.209134000',
           '2025-03-28T01:18:17.209134000', '2025-03-28T02:18:17.209134000',
           '2025-03-28T03:18:17.209134000', '2025-03-28T04:18:17.209134000',
           '2025-03-28T05:18:17.209134000', '2025-03-28T06:18:17.209134000',
           '2025-03-28T07:18:17.209134000', '2025-03-28T08:18:17.209134000',
           '2025-03-28T09:18:17.209134000', '2025-03-28T10:18:17.209134000',
           '2025-03-28T11:18:17.209134000', '2025-03-28T12:18:17.209134000',
           '2025-03-28T13:18:17.209134000', '2025-03-28T14:18:17.209134000',
           '2025-03-28T15:18:17.209134000', '2025-03-28T16:18:17.209134000',
           '2025-03-28T17:18:17.209134000', '2025-03-28T18:18:17.209134000',
           '2025-03-28T19:18:17.209134000', '2025-03-28T20:18:17.209134000',
           '2025-03-28T21:18:17.209134000', '2025-03-28T22:18:17.209134000',
           '2025-03-28T23:18:17.209134000', '2025-03-29T00:18:17.209134000',
           '2025-03-29T01:18:17.209134000', '2025-03-29T02:18:17.209134000',
           '2025-03-29T03:18:17.209134000', '2025-03-29T04:18:17.209134000',
           '2025-03-29T05:18:17.209134000', '2025-03-29T06:18:17.209134000',
           '2025-03-29T07:18:17.209134000', '2025-03-29T08:18:17.209134000',
           '2025-03-29T09:18:17.209134000', '2025-03-29T10:18:17.209134000',
           '2025-03-29T11:18:17.209134000', '2025-03-29T12:18:17.209134000',
           '2025-03-29T13:18:17.209134000', '2025-03-29T14:18:17.209134000',
           '2025-03-29T15:18:17.209134000'], 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., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan]], 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., ..., nan, nan, nan],
           ...,
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ..., nan, nan, nan]], 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., ...,     nan,     nan,     nan],
           ...,
           [     0.,   3600.,   7200., ..., 165600., 169200., 172800.],
           [     0.,   3600.,   7200., ..., 165600., 169200., 172800.],
           [     0.,   3600.,   7200., ...,     nan,     nan,     nan]],
          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., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan]], dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    axis :

    X

    minval :

    18.717579

    maxval :

    19.329117

    array([[19.195274, 19.173471, 19.166887, ..., 19.002123, 19.001467,
            19.015726],
           [19.191893, 19.177912, 19.182077, ...,       nan,       nan,
                  nan],
           [19.193327, 19.202751, 19.198984, ...,       nan,       nan,
                  nan],
           ...,
           [19.190422, 19.193575, 19.202   , ..., 19.021955, 19.03009 ,
            19.051191],
           [19.187893, 19.165068, 19.167374, ..., 19.03456 , 19.04575 ,
            19.04256 ],
           [19.1907  , 19.192673, 19.186674, ...,       nan,       nan,
                  nan]], dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    axis :

    Y

    minval :

    79.419716

    maxval :

    79.61108

    array([[79.59931 , 79.59458 , 79.59429 , ..., 79.516525, 79.50844 ,
            79.504105],
           [79.598854, 79.59388 , 79.59525 , ...,       nan,       nan,
                  nan],
           [79.59855 , 79.59485 , 79.5905  , ...,       nan,       nan,
                  nan],
           ...,
           [79.59836 , 79.59492 , 79.59524 , ..., 79.51213 , 79.50883 ,
            79.50658 ],
           [79.59878 , 79.59678 , 79.595825, ..., 79.47051 , 79.463806,
            79.4529  ],
           [79.59873 , 79.5942  , 79.595215, ...,       nan,       nan,
                  nan]], 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 :

    -121.51813

    maxval :

    0.0

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

    array([[0.        , 0.        , 0.        , ..., 0.04654489, 0.04654351,
            0.0465438 ],
           [0.        , 0.        , 0.        , ...,        nan,        nan,
                   nan],
           [0.        , 0.        , 0.        , ...,        nan,        nan,
                   nan],
           ...,
           [0.        , 0.        , 0.        , ..., 0.01639567, 0.01640037,
            0.01639698],
           [0.        , 0.        , 0.        , ..., 0.04259444, 0.0425945 ,
            0.0426034 ],
           [0.        , 0.        , 0.        , ...,        nan,        nan,
                   nan]], 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, ...,        nan,        nan,
                   nan],
           ...,
           [0.32936665, 0.32936665, 0.32936665, ..., 0.32936665, 0.32936665,
            0.32936665],
           [0.32936665, 0.32936665, 0.32936665, ..., 0.32936665, 0.32936665,
            0.32936665],
           [0.32936665, 0.32936665, 0.32936665, ...,        nan,        nan,
                   nan]], dtype=float32)

  • viscosity
    (trajectory, time)
    float32
    0.02169 0.1651 0.1649 ... nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    m2/s

    description :

    Kinematic viscosity of oil (emulsion)

    minval :

    0.021692334

    maxval :

    0.16646668

    array([[0.02169233, 0.1650736 , 0.16488479, ..., 0.1597472 , 0.15978074,
            0.15975209],
           [0.02169233, 0.16490114, 0.16497162, ...,        nan,        nan,
                   nan],
           [0.02169233, 0.16504842, 0.16511957, ...,        nan,        nan,
                   nan],
           ...,
           [0.02169233, 0.16503204, 0.1651158 , ..., 0.15987523, 0.1597539 ,
            0.15983255],
           [0.02169233, 0.16489422, 0.16482379, ..., 0.16333601, 0.16330485,
            0.16291037],
           [0.02169233, 0.16501063, 0.16493008, ...,        nan,        nan,
                   nan]], dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    kg/m^3

    minval :

    988.1

    maxval :

    997.18353

    array([[988.1    , 997.1478 , 997.143  , ..., 997.0085 , 997.00934,
            997.0086 ],
           [988.1    , 997.14343, 997.1452 , ...,       nan,       nan,
                  nan],
           [988.1    , 997.1472 , 997.14905, ...,       nan,       nan,
                  nan],
           ...,
           [988.1    , 997.1468 , 997.1489 , ..., 997.0119 , 997.00867,
            997.01074],
           [988.1    , 997.14325, 997.1414 , ..., 997.1029 , 997.1021 ,
            997.0918 ],
           [988.1    , 997.14624, 997.14417, ...,       nan,       nan,
                  nan]], 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., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan]], 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., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan]], 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., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan]], 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., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan]], 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., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan]], 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., ..., nan, nan, nan],
           ...,
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ..., nan, nan, nan]], 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., ..., nan, nan, nan],
           ...,
           [ 3.,  3.,  3., ...,  3.,  3.,  3.],
           [ 3.,  3.,  3., ...,  3.,  3.,  3.],
           [ 3.,  3.,  3., ..., nan, nan, nan]], 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., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan]], 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., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan]], 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, ...,   nan,   nan,   nan],
           ...,
           [0.001, 0.001, 0.001, ..., 0.001, 0.001, 0.001],
           [0.001, 0.001, 0.001, ..., 0.001, 0.001, 0.001],
           [0.001, 0.001, 0.001, ...,   nan,   nan,   nan]], 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.0029986985

    array([[0.        , 0.        , 0.        , ..., 0.00279284, 0.00279284,
            0.00279284],
           [0.        , 0.        , 0.        , ...,        nan,        nan,
                   nan],
           [0.        , 0.        , 0.        , ...,        nan,        nan,
                   nan],
           ...,
           [0.        , 0.        , 0.        , ..., 0.00133724, 0.00133724,
            0.00133724],
           [0.        , 0.        , 0.        , ..., 0.00234738, 0.00234738,
            0.00234738],
           [0.        , 0.        , 0.        , ...,        nan,        nan,
                   nan]], dtype=float32)

  • x_sea_water_velocity
    (trajectory, time)
    float32
    0.2167 -0.1164 -0.0507 ... nan nan

    minval :

    -0.37956378

    maxval :

    0.35684502

    array([[ 0.21667446, -0.11635908, -0.05070413, ...,  0.03812661,
             0.21226154,  0.18440679],
           [ 0.06860436,  0.03710322, -0.01434156, ...,         nan,
                    nan,         nan],
           [ 0.12029263, -0.02979132, -0.15389013, ...,         nan,
                    nan,         nan],
           ...,
           [ 0.0739558 ,  0.00544716, -0.14188348, ...,  0.12723124,
             0.17593452,  0.20391235],
           [ 0.0605474 ,  0.04315057, -0.17664507, ...,  0.08169084,
             0.17661567,  0.05816387],
           [ 0.14879149, -0.00200194, -0.10342697, ...,         nan,
                    nan,         nan]], dtype=float32)

  • y_sea_water_velocity
    (trajectory, time)
    float32
    -0.0871 -0.09177 ... nan nan

    minval :

    -0.6841104

    maxval :

    0.55537975

    array([[-0.08709584, -0.09177262,  0.09123033, ..., -0.13535647,
            -0.41412652, -0.50803626],
           [-0.13704118,  0.0080362 ,  0.08246801, ...,         nan,
                    nan,         nan],
           [-0.16417277, -0.03865303,  0.08343169, ...,         nan,
                    nan,         nan],
           ...,
           [-0.26062942, -0.00560117,  0.01283811, ..., -0.22351435,
            -0.38485157, -0.48117328],
           [-0.16079438, -0.00672221,  0.0378859 , ..., -0.15374841,
            -0.4482141 , -0.3975203 ],
           [-0.14834888, -0.00333169,  0.03995663, ...,         nan,
                    nan,         nan]], dtype=float32)

  • x_wind
    (trajectory, time)
    float32
    2.431 2.514 2.29 ... nan nan nan

    minval :

    -4.2348137

    maxval :

    8.95596

    array([[2.4306226 , 2.5138652 , 2.2899826 , ..., 3.4242523 , 3.6334057 ,
            5.755796  ],
           [2.1922438 , 2.4164867 , 3.3013854 , ...,        nan,        nan,
                   nan],
           [1.3934088 , 1.2315664 , 2.5658085 , ...,        nan,        nan,
                   nan],
           ...,
           [2.373311  , 2.5214512 , 2.8318386 , ..., 4.561341  , 4.4757433 ,
            5.179099  ],
           [1.7078764 , 2.1497426 , 2.6585393 , ..., 5.050564  , 4.211967  ,
            3.3506088 ],
           [1.8888566 , 0.99848884, 3.4338706 , ...,        nan,        nan,
                   nan]], dtype=float32)

  • y_wind
    (trajectory, time)
    float32
    -2.788 -3.952 -3.902 ... nan nan

    minval :

    -12.43363

    maxval :

    -1.4763027

    array([[-2.7879066, -3.9516175, -3.9016013, ..., -8.306592 , -8.870035 ,
            -9.598288 ],
           [-3.7420728, -2.4255493, -3.2575674, ...,        nan,        nan,
                   nan],
           [-4.239407 , -2.927578 , -3.1629472, ...,        nan,        nan,
                   nan],
           ...,
           [-4.0540075, -3.292201 , -3.9553607, ..., -8.525614 , -8.070749 ,
            -9.427613 ],
           [-4.11392  , -3.5078433, -2.8367844, ..., -9.400041 , -8.234567 ,
            -8.534998 ],
           [-4.0112844, -2.3327153, -4.1284866, ...,        nan,        nan,
                   nan]], dtype=float32)

  • sea_surface_height
    (trajectory, time)
    float32
    -0.0528 -0.3207 -0.754 ... nan nan

    minval :

    -1.2267007

    maxval :

    0.054150004

    array([[-0.05280218, -0.32071897, -0.7539749 , ..., -0.51393884,
            -0.28479853, -0.12317982],
           [-0.05260313, -0.31993514, -0.7525424 , ...,         nan,
                    nan,         nan],
           [-0.05246046, -0.318341  , -0.7499236 , ...,         nan,
                    nan,         nan],
           ...,
           [-0.05238155, -0.31916124, -0.7504521 , ..., -0.5182542 ,
            -0.28905225, -0.12837932],
           [-0.05257996, -0.32261908, -0.754192  , ..., -0.54277194,
            -0.31813744, -0.1455281 ],
           [-0.05255093, -0.31883338, -0.752053  , ...,         nan,
                    nan,         nan]], dtype=float32)

  • upward_sea_water_velocity
    (trajectory, time)
    float32
    -3.46e-05 -0.00011 ... nan nan

    minval :

    -0.00022860074

    maxval :

    0.00021326213

    array([[-3.45965564e-05, -1.10030822e-04, -1.12913243e-04, ...,
             3.36888188e-05,  7.56681038e-05,  8.48517884e-05],
           [-3.46294801e-05, -1.09989604e-04, -1.13875343e-04, ...,
                        nan,             nan,             nan],
           [-3.45566295e-05, -1.10032161e-04, -1.13790149e-04, ...,
                        nan,             nan,             nan],
           ...,
           [-3.46085508e-05, -1.10044122e-04, -1.14875365e-04, ...,
             3.43569154e-05,  7.57786547e-05,  8.42711088e-05],
           [-3.47238383e-05, -1.10070403e-04, -1.13243557e-04, ...,
             3.57021745e-05,  8.27945260e-05,  9.65494037e-05],
           [-3.46457346e-05, -1.10002373e-04, -1.14100258e-04, ...,
                        nan,             nan,             nan]], dtype=float32)

  • sea_surface_wave_significant_height
    (trajectory, time)
    float32
    0.3365 0.5396 0.5035 ... nan nan

    minval :

    0.0747244

    maxval :

    4.598056

    array([[0.3365366 , 0.539596  , 0.5034762 , ..., 1.9858344 , 2.260227  ,
            3.081305  ],
           [0.46270242, 0.28837836, 0.52916783, ...,        nan,        nan,
                   nan],
           [0.48988828, 0.2481517 , 0.40805516, ...,        nan,        nan,
                   nan],
           ...,
           [0.5428624 , 0.42302904, 0.58213896, ..., 2.299901  , 2.095164  ,
            2.8462927 ],
           [0.48809308, 0.41638833, 0.37183332, ..., 2.801176  , 2.1044993 ,
            2.0681903 ],
           [0.48359126, 0.1583881 , 0.7093623 , ...,        nan,        nan,
                   nan]], 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., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan]], 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., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan]], 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., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan]], dtype=float32)

  • sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment
    (trajectory, time)
    float32
    3.16 4.002 3.866 ... nan nan nan

    minval :

    1.4892265

    maxval :

    11.681987

    array([[3.1604266, 4.0018787, 3.865619 , ..., 7.6771684, 8.19041  ,
            9.563061 ],
           [3.7057872, 2.9255729, 3.9630203, ...,       nan,       nan,
                  nan],
           [3.8130994, 2.7138648, 3.4800782, ...,       nan,       nan,
                  nan],
           ...,
           [4.013973 , 3.5433552, 4.156645 , ..., 8.26198  , 7.8856697,
            9.191139 ],
           [3.8061063, 3.5154335, 3.3220315, ..., 9.118004 , 7.9032183,
            7.834744 ],
           [3.7885132, 2.1681576, 4.588425 , ...,       nan,       nan,
                  nan]], dtype=float32)

  • sea_ice_area_fraction
    (trajectory, time)
    float32
    0.8619 0.7942 0.7692 ... nan nan

    minval :

    0.4407441

    maxval :

    0.925149

    array([[0.8618688 , 0.79419446, 0.76918447, ..., 0.83293056, 0.81508714,
            0.8090444 ],
           [0.8635818 , 0.7940279 , 0.78978354, ...,        nan,        nan,
                   nan],
           [0.8628792 , 0.8260913 , 0.7713555 , ...,        nan,        nan,
                   nan],
           ...,
           [0.86433625, 0.8188036 , 0.8030685 , ..., 0.814806  , 0.8130601 ,
            0.8191261 ],
           [0.86560094, 0.79473794, 0.7809452 , ..., 0.68843   , 0.70959884,
            0.70613384],
           [0.8641846 , 0.8121934 , 0.79305756, ...,        nan,        nan,
                   nan]], dtype=float32)

  • sea_ice_x_velocity
    (trajectory, time)
    float32
    -0.01443 -0.01141 ... nan nan

    minval :

    -0.10480787

    maxval :

    0.14134997

    array([[-0.01442876, -0.01141203, -0.02034282, ...,  0.03608811,
             0.05135135,  0.06236253],
           [-0.01445998, -0.01101573, -0.01927719, ...,         nan,
                    nan,         nan],
           [-0.01435197, -0.00976923, -0.01847446, ...,         nan,
                    nan,         nan],
           ...,
           [-0.01441683, -0.01025101, -0.01818361, ...,  0.03372276,
             0.05234284,  0.06745415],
           [-0.01458928, -0.01152821, -0.02018995, ...,  0.01728074,
             0.02973689,  0.03581415],
           [-0.01447895, -0.0101911 , -0.01900607, ...,         nan,
                    nan,         nan]], dtype=float32)

  • sea_ice_y_velocity
    (trajectory, time)
    float32
    -0.06275 -0.07013 ... nan nan

    minval :

    -0.30817115

    maxval :

    0.03194577

    array([[-0.06274709, -0.07013223, -0.06176553, ..., -0.09728571,
            -0.13141827, -0.15652607],
           [-0.06373025, -0.07038913, -0.05893559, ...,         nan,
                    nan,         nan],
           [-0.06383684, -0.06610526, -0.0610414 , ...,         nan,
                    nan,         nan],
           ...,
           [-0.06446762, -0.06720979, -0.05627083, ..., -0.10106491,
            -0.13652472, -0.16418122],
           [-0.06440824, -0.06866269, -0.06032661, ..., -0.0962485 ,
            -0.11710767, -0.13457416],
           [-0.06403563, -0.06815775, -0.05836084, ...,         nan,
                    nan,         nan]], dtype=float32)

  • sea_water_temperature
    (trajectory, time)
    float32
    -1.725 -1.715 -1.716 ... nan nan

    minval :

    -1.7741855

    maxval :

    271.65103

    array([[ -1.7252067,  -1.7147746,  -1.7164727, ...,  -1.5144022,
             -1.5112064, 271.6371   ],
           [ -1.7234782,  -1.7166959,  -1.7193401, ...,         nan,
                    nan,         nan],
           [ -1.7241776,  -1.7221906,  -1.7287687, ...,         nan,
                    nan,         nan],
           ...,
           [ -1.7227112,  -1.7195339,  -1.7249702, ...,  -1.5164993,
             -1.5125393, 271.63388  ],
           [ -1.7214614,  -1.7106429,  -1.7149098, ...,  -1.654883 ,
             -1.6450124, 271.51227  ],
           [ -1.7228723,  -1.7201538,  -1.7206268, ...,         nan,
                    nan,         nan]], dtype=float32)

  • sea_water_salinity
    (trajectory, time)
    float32
    34.13 34.14 34.14 ... nan nan nan

    minval :

    34.09042

    maxval :

    34.148785

    array([[34.13499 , 34.13639 , 34.135418, ..., 34.141777, 34.140724,
            34.14069 ],
           [34.135326, 34.135906, 34.135315, ...,       nan,       nan,
                  nan],
           [34.135136, 34.135403, 34.13258 , ...,       nan,       nan,
                  nan],
           ...,
           [34.135445, 34.13582 , 34.13451 , ..., 34.141113, 34.14051 ,
            34.14032 ],
           [34.13578 , 34.137543, 34.1361  , ..., 34.116623, 34.117657,
            34.119728],
           [34.135452, 34.135513, 34.135124, ...,       nan,       nan,
                  nan]], dtype=float32)

  • sea_floor_depth_below_sea_level
    (trajectory, time)
    float32
    170.4 168.0 168.4 ... nan nan nan

    minval :

    42.480755

    maxval :

    212.24054

    array([[170.4251 , 168.00531, 168.41661, ...,  70.33888,  67.64344,
             66.17612],
           [170.93388, 166.38481, 167.92552, ...,       nan,       nan,
                  nan],
           [170.35571, 164.76051, 158.58131, ...,       nan,       nan,
                  nan],
           ...,
           [170.93738, 165.99576, 165.43837, ...,  69.60985,  68.93237,
             69.22015],
           [171.87785, 172.50331, 170.66501, ...,  65.52206,  65.69405,
             64.77807],
           [171.14177, 165.01735, 167.29327, ...,       nan,       nan,
                  nan]], dtype=float32)

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

    minval :

    8.212984e-05

    maxval :

    0.074306734

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

• Attributes: (162)

  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-04-03 15:18:21.343398

  date_created :

  2025-04-03T15:18:21.343414

  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-03-27 15:18:17.209134

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

  True

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

  seed_lon :

  19.1909

  seed_lat :

  79.5986

  seed_radius :

  50.0

  seed_time :

  2025-03-27 15:18:17.209134

  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-04-03 15:18:21.455845

  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-03-29 15:18:17.209134

  time_coverage_duration :

  P2DT0H0M0S

  time_coverage_resolution :

  P0DT0H30M0S

  performance :

  -------------------- Reader performance: -------------------- global_landmask 0:00:01.5 total 0:00:00.0 preparing 0:00:01.5 reading 0:00:00.0 masking -------------------- https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be 0:02:43.7 total 0:00:00.0 preparing 0:02:40.9 reading 0:00:01.6 interpolation 0:00:00.1 interpolation_time 0:00:02.7 rotating vectors 0:00:00.0 masking -------------------- Performance: 5:03.7 total time 4.2 configuration 0.1 preparing main loop 0.0 moving elements to ocean 4:59.3 main loop 5.4 updating elements 0.0 oil weathering 0.0 updating viscosities 0.0 updating densities 4.9 vertical mixing 0.0 cleaning up --------------------

  geospatial_bounds_crs :

  EPSG:4326

  geospatial_bounds_vertical_crs :

  EPSG:5831

  geospatial_lat_min :

  79.419716

  geospatial_lat_max :

  79.61108

  geospatial_lat_units :

  degrees_north

  geospatial_lat_resolution :

  point

  geospatial_lon_min :

  18.717579

  geospatial_lon_max :

  19.329117

  geospatial_lon_units :

  degrees_east

  geospatial_lon_resolution :

  point

  runtime :

  0:05:03.759016

  geospatial_vertical_min :

  -121.51813

  geospatial_vertical_max :

  0.0

  geospatial_vertical_positive :

  up

Print and plot results

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

===========================
--------------------
Reader performance:
--------------------
global_landmask
 0:00:01.5  total
 0:00:00.0  preparing
 0:00:01.5  reading
 0:00:00.0  masking
--------------------
https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
 0:02:43.7  total
 0:00:00.0  preparing
 0:02:40.9  reading
 0:00:01.6  interpolation
 0:00:00.1  interpolation_time
 0:00:02.7  rotating vectors
 0:00:00.0  masking
--------------------
Performance:
 5:03.7 total time
    4.2 configuration
    0.1 preparing main loop
      0.0 moving elements to ocean
 4:59.3 main loop
      5.4 updating elements
        0.0 oil weathering
          0.0 updating viscosities
          0.0 updating densities
        4.9 vertical mixing
    0.0 cleaning up
--------------------
===========================
Model:  OpenOil     (OpenDrift version 1.13.1)
        2135 active Oil particles  (865 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-03-27 15:18:17.209134 UTC
        Present: 2025-03-29 15:18:17.209134 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.11/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)
15:23:34 INFO    opendrift:4600: Saving animation to /root/project/docs/source/gallery/animations/example_oil_ice_0.gif...
15:24:09 INFO    opendrift:3036: Time to make animation: 0:00:48.910043

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-03-27 15:18 to 2025-03-29 15:18 UTC (49 steps)'}>, <Figure size 477.789x1100 with 2 Axes>)