Compare

Comparing two oil drift simulation runs, with and without wind

from datetime import timedelta
from opendrift.readers import reader_netCDF_CF_generic
from opendrift.models.openoil import OpenOil

o = OpenOil(loglevel=20)  # Set loglevel to 0 for debug information

# Arome atmospheric model
reader_arome = reader_netCDF_CF_generic.Reader(o.test_data_folder() +
    '16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc')
# Norkyst ocean model
reader_norkyst = reader_netCDF_CF_generic.Reader(o.test_data_folder() +
    '16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc')

o.add_reader([reader_norkyst, reader_arome])

# Seeding some particles
lon = 4.5; lat = 60.0; # Outside Bergen
time = [reader_arome.start_time,
        reader_arome.start_time + timedelta(hours=30)]
o.seed_elements(lon, lat, radius=50, number=5000, time=time,
                oil_type='GENERIC HEAVY CRUDE',
                wind_drift_factor=0.03) # 3% wind drift

# Running model
o.run(steps=66, time_step=1800, time_step_output=3600)

# Second run, for comparison
o2 = OpenOil(loglevel=20)  # Set loglevel to 0 for debug information
o2.add_reader([reader_norkyst, reader_arome])
o2.seed_elements(lon, lat, radius=50, number=5000, time=time,
                 oil_type='GENERIC HEAVY CRUDE',
                 wind_drift_factor=0.0) # No wind drift
o2.run(steps=66, time_step=1800, time_step_output=3600)

15:22:37 INFO    opendrift:509: OpenDriftSimulation initialised (version 1.13.1 / v1.13.1-12-gca7f755)
15:22:37 INFO    opendrift.readers:61: Opening file with xr.open_dataset
15:22:38 INFO    opendrift.readers.reader_netCDF_CF_generic:332: Detected dimensions: {'time': 'time', 'x': 'x', 'y': 'y'}
15:22:38 INFO    opendrift.readers:61: Opening file with xr.open_dataset
15:22:38 INFO    opendrift.readers.reader_netCDF_CF_generic:332: Detected dimensions: {'x': 'X', 'y': 'Y', 'z': 'depth', 'time': 'time'}
15:22:38 INFO    opendrift.models.openoil.adios.dirjs:86: Querying ADIOS database for oil: GENERIC HEAVY CRUDE
15:22:38 INFO    opendrift.models.openoil.openoil:1710: Using density 963.50937 and viscosity 0.007984481780964631 of oiltype GENERIC HEAVY CRUDE
15:22:38 INFO    opendrift.models.basemodel.environment:206: Adding a global landmask from GSHHG
15:22:42 INFO    opendrift.models.basemodel.environment:229: Fallback values will be used for the following variables which have no readers:
15:22:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_height: 0.000000
15:22:42 INFO    opendrift.models.basemodel.environment:232:    upward_sea_water_velocity: 0.000000
15:22:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_significant_height: 0.000000
15:22:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
15:22:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
15:22:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
15:22:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
15:22:42 INFO    opendrift.models.basemodel.environment:232:    sea_ice_area_fraction: 0.000000
15:22:42 INFO    opendrift.models.basemodel.environment:232:    sea_ice_x_velocity: 0.000000
15:22:42 INFO    opendrift.models.basemodel.environment:232:    sea_ice_y_velocity: 0.000000
15:22:42 INFO    opendrift.models.basemodel.environment:232:    sea_water_temperature: 10.000000
15:22:42 INFO    opendrift.models.basemodel.environment:232:    sea_water_salinity: 34.000000
15:22:42 INFO    opendrift.models.basemodel.environment:232:    sea_floor_depth_below_sea_level: 10000.000000
15:22:42 INFO    opendrift.models.basemodel.environment:232:    ocean_vertical_diffusivity: 0.020000
15:22:42 INFO    opendrift.models.basemodel.environment:232:    ocean_mixed_layer_thickness: 50.000000
15:22:42 INFO    opendrift:919: Using existing reader for land_binary_mask
15:22:42 INFO    opendrift:930: All points are in ocean
15:22:42 INFO    opendrift.models.openoil.openoil:685: Oil-water surface tension is 0.035048 Nm
15:22:42 INFO    opendrift.models.openoil.openoil:698: Max water fraction not available for GENERIC HEAVY CRUDE, using default
15:22:42 INFO    opendrift:2056: 2015-11-16 00:00:00 - step 1 of 66 - 84 active elements (0 deactivated)
15:22:42 INFO    opendrift:2056: 2015-11-16 00:30:00 - step 2 of 66 - 167 active elements (0 deactivated)
15:22:42 INFO    opendrift:2056: 2015-11-16 01:00:00 - step 3 of 66 - 250 active elements (0 deactivated)
15:22:42 INFO    opendrift:2056: 2015-11-16 01:30:00 - step 4 of 66 - 334 active elements (0 deactivated)
15:22:43 INFO    opendrift:2056: 2015-11-16 02:00:00 - step 5 of 66 - 417 active elements (0 deactivated)
15:22:43 INFO    opendrift:2056: 2015-11-16 02:30:00 - step 6 of 66 - 500 active elements (0 deactivated)
15:22:43 INFO    opendrift:2056: 2015-11-16 03:00:00 - step 7 of 66 - 584 active elements (0 deactivated)
15:22:43 INFO    opendrift:2056: 2015-11-16 03:30:00 - step 8 of 66 - 667 active elements (0 deactivated)
15:22:43 INFO    opendrift:2056: 2015-11-16 04:00:00 - step 9 of 66 - 750 active elements (0 deactivated)
15:22:43 INFO    opendrift:2056: 2015-11-16 04:30:00 - step 10 of 66 - 834 active elements (0 deactivated)
15:22:43 INFO    opendrift:2056: 2015-11-16 05:00:00 - step 11 of 66 - 917 active elements (0 deactivated)
15:22:43 INFO    opendrift:2056: 2015-11-16 05:30:00 - step 12 of 66 - 1000 active elements (0 deactivated)
15:22:43 INFO    opendrift:2056: 2015-11-16 06:00:00 - step 13 of 66 - 1084 active elements (0 deactivated)
15:22:44 INFO    opendrift:2056: 2015-11-16 06:30:00 - step 14 of 66 - 1167 active elements (0 deactivated)
15:22:44 INFO    opendrift:2056: 2015-11-16 07:00:00 - step 15 of 66 - 1250 active elements (0 deactivated)
15:22:44 INFO    opendrift:2056: 2015-11-16 07:30:00 - step 16 of 66 - 1334 active elements (0 deactivated)
15:22:44 INFO    opendrift:2056: 2015-11-16 08:00:00 - step 17 of 66 - 1417 active elements (0 deactivated)
15:22:44 INFO    opendrift:2056: 2015-11-16 08:30:00 - step 18 of 66 - 1500 active elements (0 deactivated)
15:22:44 INFO    opendrift:2056: 2015-11-16 09:00:00 - step 19 of 66 - 1584 active elements (0 deactivated)
15:22:45 INFO    opendrift:2056: 2015-11-16 09:30:00 - step 20 of 66 - 1667 active elements (0 deactivated)
15:22:45 INFO    opendrift:2056: 2015-11-16 10:00:00 - step 21 of 66 - 1750 active elements (0 deactivated)
15:22:45 INFO    opendrift:2056: 2015-11-16 10:30:00 - step 22 of 66 - 1833 active elements (0 deactivated)
15:22:45 INFO    opendrift:2056: 2015-11-16 11:00:00 - step 23 of 66 - 1917 active elements (0 deactivated)
15:22:45 INFO    opendrift:2056: 2015-11-16 11:30:00 - step 24 of 66 - 2000 active elements (0 deactivated)
15:22:45 INFO    opendrift:2056: 2015-11-16 12:00:00 - step 25 of 66 - 2083 active elements (0 deactivated)
15:22:45 INFO    opendrift:2056: 2015-11-16 12:30:00 - step 26 of 66 - 2167 active elements (0 deactivated)
15:22:46 INFO    opendrift:2056: 2015-11-16 13:00:00 - step 27 of 66 - 2250 active elements (0 deactivated)
15:22:46 INFO    opendrift:2056: 2015-11-16 13:30:00 - step 28 of 66 - 2333 active elements (0 deactivated)
15:22:46 INFO    opendrift:2056: 2015-11-16 14:00:00 - step 29 of 66 - 2417 active elements (0 deactivated)
15:22:46 INFO    opendrift:2056: 2015-11-16 14:30:00 - step 30 of 66 - 2500 active elements (0 deactivated)
15:22:46 INFO    opendrift:2056: 2015-11-16 15:00:00 - step 31 of 66 - 2583 active elements (0 deactivated)
15:22:47 INFO    opendrift:2056: 2015-11-16 15:30:00 - step 32 of 66 - 2667 active elements (0 deactivated)
15:22:47 INFO    opendrift:2056: 2015-11-16 16:00:00 - step 33 of 66 - 2750 active elements (0 deactivated)
15:22:47 INFO    opendrift:2056: 2015-11-16 16:30:00 - step 34 of 66 - 2833 active elements (0 deactivated)
15:22:47 INFO    opendrift:2056: 2015-11-16 17:00:00 - step 35 of 66 - 2917 active elements (0 deactivated)
15:22:47 INFO    opendrift:2056: 2015-11-16 17:30:00 - step 36 of 66 - 3000 active elements (0 deactivated)
15:22:48 INFO    opendrift:2056: 2015-11-16 18:00:00 - step 37 of 66 - 3083 active elements (0 deactivated)
15:22:48 INFO    opendrift:2056: 2015-11-16 18:30:00 - step 38 of 66 - 3167 active elements (0 deactivated)
15:22:48 INFO    opendrift:2056: 2015-11-16 19:00:00 - step 39 of 66 - 3250 active elements (0 deactivated)
15:22:48 INFO    opendrift:2056: 2015-11-16 19:30:00 - step 40 of 66 - 3333 active elements (0 deactivated)
15:22:48 INFO    opendrift:2056: 2015-11-16 20:00:00 - step 41 of 66 - 3416 active elements (0 deactivated)
15:22:49 INFO    opendrift:2056: 2015-11-16 20:30:00 - step 42 of 66 - 3500 active elements (0 deactivated)
15:22:49 INFO    opendrift:2056: 2015-11-16 21:00:00 - step 43 of 66 - 3583 active elements (0 deactivated)
15:22:49 INFO    opendrift:2056: 2015-11-16 21:30:00 - step 44 of 66 - 3666 active elements (0 deactivated)
15:22:49 INFO    opendrift:2056: 2015-11-16 22:00:00 - step 45 of 66 - 3750 active elements (0 deactivated)
15:22:50 INFO    opendrift:2056: 2015-11-16 22:30:00 - step 46 of 66 - 3833 active elements (0 deactivated)
15:22:53 INFO    opendrift:2056: 2015-11-16 23:00:00 - step 47 of 66 - 3915 active elements (1 deactivated)
15:22:53 INFO    opendrift:2056: 2015-11-16 23:30:00 - step 48 of 66 - 3996 active elements (4 deactivated)
15:22:53 INFO    opendrift:2056: 2015-11-17 00:00:00 - step 49 of 66 - 4073 active elements (10 deactivated)
15:22:54 INFO    opendrift:2056: 2015-11-17 00:30:00 - step 50 of 66 - 4147 active elements (19 deactivated)
15:22:54 INFO    opendrift:2056: 2015-11-17 01:00:00 - step 51 of 66 - 4184 active elements (66 deactivated)
15:22:54 INFO    opendrift:2056: 2015-11-17 01:30:00 - step 52 of 66 - 4171 active elements (162 deactivated)
15:22:54 INFO    opendrift:2056: 2015-11-17 02:00:00 - step 53 of 66 - 4116 active elements (300 deactivated)
15:22:55 INFO    opendrift:2056: 2015-11-17 02:30:00 - step 54 of 66 - 4024 active elements (476 deactivated)
15:22:55 INFO    opendrift:2056: 2015-11-17 03:00:00 - step 55 of 66 - 3887 active elements (696 deactivated)
15:22:56 INFO    opendrift:2056: 2015-11-17 03:30:00 - step 56 of 66 - 3792 active elements (874 deactivated)
15:22:56 INFO    opendrift:2056: 2015-11-17 04:00:00 - step 57 of 66 - 3692 active elements (1058 deactivated)
15:22:56 INFO    opendrift:2056: 2015-11-17 04:30:00 - step 58 of 66 - 3596 active elements (1237 deactivated)
15:22:57 INFO    opendrift:2056: 2015-11-17 05:00:00 - step 59 of 66 - 3563 active elements (1353 deactivated)
15:22:57 INFO    opendrift:2056: 2015-11-17 05:30:00 - step 60 of 66 - 3563 active elements (1436 deactivated)
15:22:57 INFO    opendrift:2056: 2015-11-17 06:00:00 - step 61 of 66 - 3478 active elements (1522 deactivated)
15:22:57 INFO    opendrift:2056: 2015-11-17 06:30:00 - step 62 of 66 - 3392 active elements (1608 deactivated)
15:22:58 INFO    opendrift:2056: 2015-11-17 07:00:00 - step 63 of 66 - 3302 active elements (1698 deactivated)
15:22:58 INFO    opendrift:2056: 2015-11-17 07:30:00 - step 64 of 66 - 3173 active elements (1827 deactivated)
15:22:58 INFO    opendrift:2056: 2015-11-17 08:00:00 - step 65 of 66 - 2933 active elements (2067 deactivated)
15:22:59 INFO    opendrift:2056: 2015-11-17 08:30:00 - step 66 of 66 - 2713 active elements (2287 deactivated)
15:22:59 INFO    opendrift:509: OpenDriftSimulation initialised (version 1.13.1 / v1.13.1-12-gca7f755)
15:22:59 INFO    opendrift.models.openoil.adios.dirjs:86: Querying ADIOS database for oil: GENERIC HEAVY CRUDE
15:22:59 INFO    opendrift.models.openoil.openoil:1710: Using density 963.50937 and viscosity 0.007984481780964631 of oiltype GENERIC HEAVY CRUDE
15:22:59 INFO    opendrift.models.basemodel.environment:206: Adding a global landmask from GSHHG
15:22:59 INFO    opendrift.models.basemodel.environment:229: Fallback values will be used for the following variables which have no readers:
15:22:59 INFO    opendrift.models.basemodel.environment:232:    sea_surface_height: 0.000000
15:22:59 INFO    opendrift.models.basemodel.environment:232:    upward_sea_water_velocity: 0.000000
15:22:59 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_significant_height: 0.000000
15:22:59 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
15:22:59 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
15:22:59 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
15:22:59 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
15:22:59 INFO    opendrift.models.basemodel.environment:232:    sea_ice_area_fraction: 0.000000
15:22:59 INFO    opendrift.models.basemodel.environment:232:    sea_ice_x_velocity: 0.000000
15:22:59 INFO    opendrift.models.basemodel.environment:232:    sea_ice_y_velocity: 0.000000
15:22:59 INFO    opendrift.models.basemodel.environment:232:    sea_water_temperature: 10.000000
15:22:59 INFO    opendrift.models.basemodel.environment:232:    sea_water_salinity: 34.000000
15:22:59 INFO    opendrift.models.basemodel.environment:232:    sea_floor_depth_below_sea_level: 10000.000000
15:22:59 INFO    opendrift.models.basemodel.environment:232:    ocean_vertical_diffusivity: 0.020000
15:22:59 INFO    opendrift.models.basemodel.environment:232:    ocean_mixed_layer_thickness: 50.000000
15:22:59 INFO    opendrift:919: Using existing reader for land_binary_mask
15:22:59 INFO    opendrift:930: All points are in ocean
15:22:59 INFO    opendrift.models.openoil.openoil:685: Oil-water surface tension is 0.035048 Nm
15:22:59 INFO    opendrift.models.openoil.openoil:698: Max water fraction not available for GENERIC HEAVY CRUDE, using default
15:22:59 INFO    opendrift:2056: 2015-11-16 00:00:00 - step 1 of 66 - 84 active elements (0 deactivated)
15:22:59 INFO    opendrift:2056: 2015-11-16 00:30:00 - step 2 of 66 - 167 active elements (0 deactivated)
15:22:59 INFO    opendrift:2056: 2015-11-16 01:00:00 - step 3 of 66 - 250 active elements (0 deactivated)
15:22:59 INFO    opendrift:2056: 2015-11-16 01:30:00 - step 4 of 66 - 334 active elements (0 deactivated)
15:22:59 INFO    opendrift:2056: 2015-11-16 02:00:00 - step 5 of 66 - 417 active elements (0 deactivated)
15:23:00 INFO    opendrift:2056: 2015-11-16 02:30:00 - step 6 of 66 - 500 active elements (0 deactivated)
15:23:00 INFO    opendrift:2056: 2015-11-16 03:00:00 - step 7 of 66 - 584 active elements (0 deactivated)
15:23:00 INFO    opendrift:2056: 2015-11-16 03:30:00 - step 8 of 66 - 667 active elements (0 deactivated)
15:23:00 INFO    opendrift:2056: 2015-11-16 04:00:00 - step 9 of 66 - 750 active elements (0 deactivated)
15:23:00 INFO    opendrift:2056: 2015-11-16 04:30:00 - step 10 of 66 - 834 active elements (0 deactivated)
15:23:00 INFO    opendrift:2056: 2015-11-16 05:00:00 - step 11 of 66 - 917 active elements (0 deactivated)
15:23:00 INFO    opendrift:2056: 2015-11-16 05:30:00 - step 12 of 66 - 1000 active elements (0 deactivated)
15:23:00 INFO    opendrift:2056: 2015-11-16 06:00:00 - step 13 of 66 - 1084 active elements (0 deactivated)
15:23:00 INFO    opendrift:2056: 2015-11-16 06:30:00 - step 14 of 66 - 1167 active elements (0 deactivated)
15:23:00 INFO    opendrift:2056: 2015-11-16 07:00:00 - step 15 of 66 - 1250 active elements (0 deactivated)
15:23:01 INFO    opendrift:2056: 2015-11-16 07:30:00 - step 16 of 66 - 1334 active elements (0 deactivated)
15:23:01 INFO    opendrift:2056: 2015-11-16 08:00:00 - step 17 of 66 - 1417 active elements (0 deactivated)
15:23:01 INFO    opendrift:2056: 2015-11-16 08:30:00 - step 18 of 66 - 1500 active elements (0 deactivated)
15:23:01 INFO    opendrift:2056: 2015-11-16 09:00:00 - step 19 of 66 - 1584 active elements (0 deactivated)
15:23:01 INFO    opendrift:2056: 2015-11-16 09:30:00 - step 20 of 66 - 1667 active elements (0 deactivated)
15:23:01 INFO    opendrift:2056: 2015-11-16 10:00:00 - step 21 of 66 - 1750 active elements (0 deactivated)
15:23:01 INFO    opendrift:2056: 2015-11-16 10:30:00 - step 22 of 66 - 1833 active elements (0 deactivated)
15:23:02 INFO    opendrift:2056: 2015-11-16 11:00:00 - step 23 of 66 - 1917 active elements (0 deactivated)
15:23:02 INFO    opendrift:2056: 2015-11-16 11:30:00 - step 24 of 66 - 2000 active elements (0 deactivated)
15:23:02 INFO    opendrift:2056: 2015-11-16 12:00:00 - step 25 of 66 - 2083 active elements (0 deactivated)
15:23:02 INFO    opendrift:2056: 2015-11-16 12:30:00 - step 26 of 66 - 2167 active elements (0 deactivated)
15:23:02 INFO    opendrift:2056: 2015-11-16 13:00:00 - step 27 of 66 - 2250 active elements (0 deactivated)
15:23:02 INFO    opendrift:2056: 2015-11-16 13:30:00 - step 28 of 66 - 2333 active elements (0 deactivated)
15:23:03 INFO    opendrift:2056: 2015-11-16 14:00:00 - step 29 of 66 - 2417 active elements (0 deactivated)
15:23:03 INFO    opendrift:2056: 2015-11-16 14:30:00 - step 30 of 66 - 2500 active elements (0 deactivated)
15:23:03 INFO    opendrift:2056: 2015-11-16 15:00:00 - step 31 of 66 - 2583 active elements (0 deactivated)
15:23:03 INFO    opendrift:2056: 2015-11-16 15:30:00 - step 32 of 66 - 2667 active elements (0 deactivated)
15:23:03 INFO    opendrift:2056: 2015-11-16 16:00:00 - step 33 of 66 - 2750 active elements (0 deactivated)
15:23:03 INFO    opendrift:2056: 2015-11-16 16:30:00 - step 34 of 66 - 2833 active elements (0 deactivated)
15:23:04 INFO    opendrift:2056: 2015-11-16 17:00:00 - step 35 of 66 - 2917 active elements (0 deactivated)
15:23:04 INFO    opendrift:2056: 2015-11-16 17:30:00 - step 36 of 66 - 3000 active elements (0 deactivated)
15:23:04 INFO    opendrift:2056: 2015-11-16 18:00:00 - step 37 of 66 - 3083 active elements (0 deactivated)
15:23:04 INFO    opendrift:2056: 2015-11-16 18:30:00 - step 38 of 66 - 3167 active elements (0 deactivated)
15:23:04 INFO    opendrift:2056: 2015-11-16 19:00:00 - step 39 of 66 - 3250 active elements (0 deactivated)
15:23:04 INFO    opendrift:2056: 2015-11-16 19:30:00 - step 40 of 66 - 3333 active elements (0 deactivated)
15:23:05 INFO    opendrift:2056: 2015-11-16 20:00:00 - step 41 of 66 - 3416 active elements (0 deactivated)
15:23:05 INFO    opendrift:2056: 2015-11-16 20:30:00 - step 42 of 66 - 3500 active elements (0 deactivated)
15:23:05 INFO    opendrift:2056: 2015-11-16 21:00:00 - step 43 of 66 - 3583 active elements (0 deactivated)
15:23:06 INFO    opendrift:2056: 2015-11-16 21:30:00 - step 44 of 66 - 3666 active elements (0 deactivated)
15:23:06 INFO    opendrift:2056: 2015-11-16 22:00:00 - step 45 of 66 - 3750 active elements (0 deactivated)
15:23:06 INFO    opendrift:2056: 2015-11-16 22:30:00 - step 46 of 66 - 3833 active elements (0 deactivated)
15:23:06 INFO    opendrift:2056: 2015-11-16 23:00:00 - step 47 of 66 - 3916 active elements (0 deactivated)
15:23:07 INFO    opendrift:2056: 2015-11-16 23:30:00 - step 48 of 66 - 4000 active elements (0 deactivated)
15:23:07 INFO    opendrift:2056: 2015-11-17 00:00:00 - step 49 of 66 - 4083 active elements (0 deactivated)
15:23:07 INFO    opendrift:2056: 2015-11-17 00:30:00 - step 50 of 66 - 4166 active elements (0 deactivated)
15:23:08 INFO    opendrift:2056: 2015-11-17 01:00:00 - step 51 of 66 - 4250 active elements (0 deactivated)
15:23:08 INFO    opendrift:2056: 2015-11-17 01:30:00 - step 52 of 66 - 4333 active elements (0 deactivated)
15:23:08 INFO    opendrift:2056: 2015-11-17 02:00:00 - step 53 of 66 - 4416 active elements (0 deactivated)
15:23:08 INFO    opendrift:2056: 2015-11-17 02:30:00 - step 54 of 66 - 4500 active elements (0 deactivated)
15:23:09 INFO    opendrift:2056: 2015-11-17 03:00:00 - step 55 of 66 - 4583 active elements (0 deactivated)
15:23:09 INFO    opendrift:2056: 2015-11-17 03:30:00 - step 56 of 66 - 4666 active elements (0 deactivated)
15:23:09 INFO    opendrift:2056: 2015-11-17 04:00:00 - step 57 of 66 - 4750 active elements (0 deactivated)
15:23:09 INFO    opendrift:2056: 2015-11-17 04:30:00 - step 58 of 66 - 4833 active elements (0 deactivated)
15:23:10 INFO    opendrift:2056: 2015-11-17 05:00:00 - step 59 of 66 - 4916 active elements (0 deactivated)
15:23:10 INFO    opendrift:2056: 2015-11-17 05:30:00 - step 60 of 66 - 4999 active elements (0 deactivated)
15:23:10 INFO    opendrift:2056: 2015-11-17 06:00:00 - step 61 of 66 - 5000 active elements (0 deactivated)
15:23:10 INFO    opendrift:2056: 2015-11-17 06:30:00 - step 62 of 66 - 5000 active elements (0 deactivated)
15:23:11 INFO    opendrift:2056: 2015-11-17 07:00:00 - step 63 of 66 - 5000 active elements (0 deactivated)
15:23:11 INFO    opendrift:2056: 2015-11-17 07:30:00 - step 64 of 66 - 5000 active elements (0 deactivated)
15:23:11 INFO    opendrift:2056: 2015-11-17 08:00:00 - step 65 of 66 - 5000 active elements (0 deactivated)
15:23:12 INFO    opendrift:2056: 2015-11-17 08:30:00 - step 66 of 66 - 5000 active elements (0 deactivated)

<xarray.Dataset> Size: 30MB
Dimensions:                                                                              (
                                                                                          trajectory: 5000,
                                                                                          time: 34)
Coordinates:
  * trajectory                                                                           (trajectory) int64 40kB ...
  * time                                                                                 (time) datetime64[ns] 272B ...
Data variables: (12/44)
    status                                                                               (trajectory, time) float32 680kB ...
    moving                                                                               (trajectory, time) float32 680kB ...
    age_seconds                                                                          (trajectory, time) float32 680kB ...
    origin_marker                                                                        (trajectory, time) float32 680kB ...
    lon                                                                                  (trajectory, time) float32 680kB ...
    lat                                                                                  (trajectory, time) float32 680kB ...
    ...                                                                                   ...
    sea_water_temperature                                                                (trajectory, time) float32 680kB ...
    sea_water_salinity                                                                   (trajectory, time) float32 680kB ...
    sea_floor_depth_below_sea_level                                                      (trajectory, time) float32 680kB ...
    ocean_vertical_diffusivity                                                           (trajectory, time) float32 680kB ...
    land_binary_mask                                                                     (trajectory, time) float32 680kB ...
    ocean_mixed_layer_thickness                                                          (trajectory, time) float32 680kB ...
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: 5000
  • time: 34
• Coordinates: (2)
  • trajectory
    (trajectory)
    int64
    0 1 2 3 4 ... 4996 4997 4998 4999

    cf_role :

    trajectory_id

    dtype :

    <class 'numpy.int32'>

    array([   0,    1,    2, ..., 4997, 4998, 4999])

  • time
    (time)
    datetime64[ns]
    2015-11-16 ... 2015-11-17T09:00:00

    standard_name :

    time

    long_name :

    time

    dtype :

    <class 'numpy.float64'>

    array(['2015-11-16T00:00:00.000000000', '2015-11-16T01:00:00.000000000',
           '2015-11-16T02:00:00.000000000', '2015-11-16T03:00:00.000000000',
           '2015-11-16T04:00:00.000000000', '2015-11-16T05:00:00.000000000',
           '2015-11-16T06:00:00.000000000', '2015-11-16T07:00:00.000000000',
           '2015-11-16T08:00:00.000000000', '2015-11-16T09:00:00.000000000',
           '2015-11-16T10:00:00.000000000', '2015-11-16T11:00:00.000000000',
           '2015-11-16T12:00:00.000000000', '2015-11-16T13:00:00.000000000',
           '2015-11-16T14:00:00.000000000', '2015-11-16T15:00:00.000000000',
           '2015-11-16T16:00:00.000000000', '2015-11-16T17:00:00.000000000',
           '2015-11-16T18:00:00.000000000', '2015-11-16T19:00:00.000000000',
           '2015-11-16T20:00:00.000000000', '2015-11-16T21:00:00.000000000',
           '2015-11-16T22:00:00.000000000', '2015-11-16T23:00:00.000000000',
           '2015-11-17T00:00:00.000000000', '2015-11-17T01:00:00.000000000',
           '2015-11-17T02:00:00.000000000', '2015-11-17T03:00:00.000000000',
           '2015-11-17T04:00:00.000000000', '2015-11-17T05:00:00.000000000',
           '2015-11-17T06:00:00.000000000', '2015-11-17T07:00:00.000000000',
           '2015-11-17T08:00:00.000000000', '2015-11-17T09:00:00.000000000'],
          dtype='datetime64[ns]')

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

    dtype :

    <class 'numpy.int32'>

    valid_range :

    [0 0]

    flag_values :

    [0]

    flag_meanings :

    active

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

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

    dtype :

    <class 'numpy.int32'>

    minval :

    1

    maxval :

    1

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

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

    dtype :

    <class 'numpy.float32'>

    units :

    s

    minval :

    0.0

    maxval :

    118800.0

    array([[     0.,   3600.,   7200., ..., 111600., 115200., 118800.],
           [     0.,   3600.,   7200., ..., 111600., 115200., 118800.],
           [     0.,   3600.,   7200., ..., 111600., 115200., 118800.],
           ...,
           [    nan,     nan,     nan, ...,   5400.,   9000.,  12600.],
           [    nan,     nan,     nan, ...,   5400.,   9000.,  12600.],
           [    nan,     nan,     nan, ...,   3600.,   7200.,  10800.]],
          dtype=float32)

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

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

  • lon
    (trajectory, time)
    float32
    4.502 4.519 4.533 ... 4.478 4.461

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    axis :

    X

    minval :

    4.4580245

    maxval :

    4.965604

    array([[4.5015807, 4.5188684, 4.532512 , ..., 4.6674204, 4.636593 ,
            4.604879 ],
           [4.5003586, 4.518657 , 4.5365477, ..., 4.638658 , 4.6016707,
            4.5720296],
           [4.500877 , 4.523524 , 4.5372343, ..., 4.724823 , 4.701701 ,
            4.6744027],
           ...,
           [      nan,       nan,       nan, ..., 4.493745 , 4.4800034,
            4.4613028],
           [      nan,       nan,       nan, ..., 4.4925485, 4.4793754,
            4.458472 ],
           [      nan,       nan,       nan, ..., 4.4918027, 4.478256 ,
            4.4612594]], dtype=float32)

  • lat
    (trajectory, time)
    float32
    60.0 60.01 60.03 ... 60.04 60.06

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    axis :

    Y

    minval :

    59.998623

    maxval :

    60.72471

    array([[60.000137, 60.01362 , 60.030224, ..., 60.570435, 60.603756,
            60.64003 ],
           [59.99967 , 60.014584, 60.030117, ..., 60.589813, 60.623653,
            60.661713],
           [59.99931 , 60.014126, 60.031284, ..., 60.53417 , 60.56159 ,
            60.598244],
           ...,
           [      nan,       nan,       nan, ..., 60.022236, 60.04377 ,
            60.069237],
           [      nan,       nan,       nan, ..., 60.021767, 60.04164 ,
            60.069458],
           [      nan,       nan,       nan, ..., 60.016193, 60.037376,
            60.064835]], dtype=float32)

  • z
    (trajectory, time)
    float32
    0.0 0.0 -5.794 0.0 ... 0.0 0.0 0.0

    dtype :

    <class 'numpy.float32'>

    units :

    m

    standard_name :

    z

    long_name :

    vertical position

    axis :

    Z

    positive :

    up

    minval :

    -48.601967

    maxval :

    0.0

    array([[ 0.      ,  0.      , -5.794138, ...,  0.      ,  0.      ,
             0.      ],
           [ 0.      ,  0.      ,  0.      , ...,  0.      ,  0.      ,
             0.      ],
           [ 0.      ,  0.      ,  0.      , ...,  0.      ,  0.      ,
             0.      ],
           ...,
           [      nan,       nan,       nan, ...,  0.      ,  0.      ,
             0.      ],
           [      nan,       nan,       nan, ...,  0.      ,  0.      ,
             0.      ],
           [      nan,       nan,       nan, ...,  0.      ,  0.      ,
             0.      ]], dtype=float32)

  • wind_drift_factor
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

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

    maxval :

    0.0

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

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

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

  • terminal_velocity
    (trajectory, time)
    float32
    0.0 0.05632 0.05344 ... 0.0 0.0 0.0

    dtype :

    <class 'numpy.float32'>

    units :

    m/s

    description :

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

    minval :

    0.0

    maxval :

    0.0689397

    array([[0.        , 0.0563175 , 0.05343587, ..., 0.02486652, 0.02486652,
            0.02486652],
           [0.        , 0.        , 0.        , ..., 0.00542277, 0.00542277,
            0.00542277],
           [0.        , 0.05577845, 0.05577845, ..., 0.00515751, 0.00515751,
            0.00515751],
           ...,
           [       nan,        nan,        nan, ..., 0.05601068, 0.03853785,
            0.03853785],
           [       nan,        nan,        nan, ..., 0.06047002, 0.04460654,
            0.04460654],
           [       nan,        nan,        nan, ..., 0.        , 0.        ,
            0.        ]], dtype=float32)

  • mass_oil
    (trajectory, time)
    float32
    5.781 4.029 3.108 ... 2.365 1.959

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.6992431

    maxval :

    5.7810564

    array([[5.7810564 , 4.028608  , 3.107647  , ..., 0.8742196 , 0.8635639 ,
            0.8557872 ],
           [5.7810564 , 3.9906292 , 2.9068484 , ..., 0.8510412 , 0.8422486 ,
            0.8344977 ],
           [5.7810564 , 3.910879  , 2.8984222 , ..., 0.8595499 , 0.85010624,
            0.8403391 ],
           ...,
           [       nan,        nan,        nan, ..., 2.7398813 , 2.1641777 ,
            1.8454107 ],
           [       nan,        nan,        nan, ..., 2.468495  , 2.0290995 ,
            1.7843633 ],
           [       nan,        nan,        nan, ..., 3.2246544 , 2.36527   ,
            1.9593251 ]], dtype=float32)

  • viscosity
    (trajectory, time)
    float32
    0.007984 0.02602 ... 0.2224 0.3903

    dtype :

    <class 'numpy.float32'>

    units :

    m2/s

    description :

    Kinematic viscosity of oil (emulsion)

    minval :

    0.007984482

    maxval :

    4.173632

    array([[0.00798448, 0.02602329, 0.03963313, ..., 3.668112  , 3.8097868 ,
            3.945326  ],
           [0.00798448, 0.02698573, 0.04249505, ..., 3.7260685 , 3.8688478 ,
            3.9934309 ],
           [0.00798448, 0.02698792, 0.04155284, ..., 3.7669709 , 3.9080415 ,
            4.049785  ],
           ...,
           [       nan,        nan,        nan, ..., 0.04408459, 0.32331142,
            0.4588824 ],
           [       nan,        nan,        nan, ..., 0.0491449 , 0.36879322,
            0.4913542 ],
           [       nan,        nan,        nan, ..., 0.03300022, 0.2224153 ,
            0.39032677]], dtype=float32)

  • density
    (trajectory, time)
    float32
    963.5 964.0 964.0 ... 989.7 992.4

    dtype :

    <class 'numpy.float32'>

    units :

    kg/m^3

    minval :

    963.5094

    maxval :

    992.3557

    array([[963.5094, 964.027 , 964.027 , ..., 992.3557, 992.3557, 992.3557],
           [963.5094, 964.027 , 964.027 , ..., 992.3557, 992.3557, 992.3557],
           [963.5094, 964.027 , 964.027 , ..., 992.3557, 992.3557, 992.3557],
           ...,
           [     nan,      nan,      nan, ..., 964.027 , 992.3557, 992.3557],
           [     nan,      nan,      nan, ..., 964.027 , 992.3557, 992.3557],
           [     nan,      nan,      nan, ..., 964.027 , 989.7101, 992.3557]],
          dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    s

    minval :

    0.0

    maxval :

    0.0

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

  • interfacial_area
    (trajectory, time)
    float32
    0.0 0.0 0.0 ... 7.957e+05 1.433e+06

    dtype :

    <class 'numpy.float32'>

    units :

    m2

    minval :

    0.0

    maxval :

    4912964.5

    array([[      0. ,       0. ,       0. , ..., 3124967.2, 3125008.8,
            3125173.8],
           [      0. ,       0. ,       0. , ..., 3160851.8, 3160980.5,
            3161113. ],
           [      0. ,       0. ,       0. , ..., 3079762. , 3079828.5,
            3079849.8],
           ...,
           [      nan,       nan,       nan, ...,  393601.5, 1130572.4,
            1684908. ],
           [      nan,       nan,       nan, ...,  424400.6, 1148134.5,
            1689387. ],
           [      nan,       nan,       nan, ...,       0. ,  795656.6,
            1433152.1]], dtype=float32)

  • mass_dispersed
    (trajectory, time)
    float32
    0.0 1.178 2.082 ... 2.813 3.206

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.0

    maxval :

    4.520344

    array([[0.       , 1.1775606, 2.082319 , ..., 4.201028 , 4.2092056,
            4.215042 ],
           [0.       , 1.2146484, 2.280606 , ..., 4.238018 , 4.2445116,
            4.250171 ],
           [0.       , 1.2932009, 2.2889338, ..., 4.219775 , 4.226879 ,
            4.234229 ],
           ...,
           [      nan,       nan,       nan, ..., 2.4471393, 3.0077236,
            3.3150659],
           [      nan,       nan,       nan, ..., 2.715673 , 3.141966 ,
            3.3769321],
           [      nan,       nan,       nan, ..., 1.9700981, 2.8134215,
            3.2062848]], dtype=float32)

  • mass_evaporated
    (trajectory, time)
    float32
    0.0 0.5749 0.5911 ... 0.6024 0.6154

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.0

    maxval :

    0.72040606

    array([[0.        , 0.57488805, 0.59109056, ..., 0.70580906, 0.70828664,
            0.7102271 ],
           [0.        , 0.57577884, 0.593602  , ..., 0.6919971 , 0.6942961 ,
            0.6963875 ],
           [0.        , 0.57697666, 0.59370035, ..., 0.7017314 , 0.7040712 ,
            0.7064882 ],
           ...,
           [       nan,        nan,        nan, ..., 0.5940359 , 0.6091552 ,
            0.62057984],
           [       nan,        nan,        nan, ..., 0.5968885 , 0.6099908 ,
            0.619761  ],
           [       nan,        nan,        nan, ..., 0.5863037 , 0.60236496,
            0.6154466 ]], dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.0

    maxval :

    0.0

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

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

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

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

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

  • fraction_evaporated
    (trajectory, time)
    float32
    0.0 0.0996 0.1417 ... 0.1831 0.2209

    dtype :

    <class 'numpy.float32'>

    units :

    %

    minval :

    0.0

    maxval :

    0.45781043

    array([[0.        , 0.09960305, 0.1416704 , ..., 0.44489947, 0.4486891 ,
            0.45218495],
           [0.        , 0.1032347 , 0.14864263, ..., 0.44646713, 0.45022744,
            0.45339686],
           [0.        , 0.10324281, 0.14640045, ..., 0.44755888, 0.4512354 ,
            0.45479816],
           ...,
           [       nan,        nan,        nan, ..., 0.1523149 , 0.20201784,
            0.23703563],
           [       nan,        nan,        nan, ..., 0.16318117, 0.21517986,
            0.24387276],
           [       nan,        nan,        nan, ..., 0.12335528, 0.18305375,
            0.22085467]], dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    %

    minval :

    0.0

    maxval :

    0.45019522

    array([[0.        , 0.        , 0.        , ..., 0.45019522, 0.45019522,
            0.45019522],
           [0.        , 0.        , 0.        , ..., 0.45019522, 0.45019522,
            0.45019522],
           [0.        , 0.        , 0.        , ..., 0.45019522, 0.45019522,
            0.45019522],
           ...,
           [       nan,        nan,        nan, ..., 0.39613616, 0.45019522,
            0.45019522],
           [       nan,        nan,        nan, ..., 0.41429162, 0.45019522,
            0.45019522],
           [       nan,        nan,        nan, ..., 0.        , 0.45019522,
            0.45019522]], dtype=float32)

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

    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, ..., 0.001, 0.001, 0.001],
           [0.001, 0.001, 0.001, ..., 0.001, 0.001, 0.001],
           ...,
           [  nan,   nan,   nan, ..., 0.001, 0.001, 0.001],
           [  nan,   nan,   nan, ..., 0.001, 0.001, 0.001],
           [  nan,   nan,   nan, ..., 0.001, 0.001, 0.001]], dtype=float32)

  • diameter
    (trajectory, time)
    float32
    0.0 0.002002 0.001802 ... 0.0 0.0

    dtype :

    <class 'numpy.float32'>

    units :

    m

    minval :

    0.0

    maxval :

    0.002999979

    array([[0.        , 0.00200201, 0.00180238, ..., 0.00159069, 0.00159069,
            0.00159069],
           [0.        , 0.        , 0.        , ..., 0.00074283, 0.00074283,
            0.00074283],
           [0.        , 0.00196387, 0.00196387, ..., 0.00072443, 0.00072443,
            0.00072443],
           ...,
           [       nan,        nan,        nan, ..., 0.00198025, 0.00198025,
            0.00198025],
           [       nan,        nan,        nan, ..., 0.00230813, 0.00230813,
            0.00230813],
           [       nan,        nan,        nan, ..., 0.        , 0.        ,
            0.        ]], dtype=float32)

  • x_sea_water_velocity
    (trajectory, time)
    float32
    0.2924 0.2011 ... -0.2531 -0.2729

    minval :

    -0.77124804

    maxval :

    1.1388335

    array([[ 0.29241988,  0.20110632,  0.28901815, ..., -0.47774726,
            -0.517898  , -0.4464849 ],
           [ 0.25904477,  0.26544204,  0.20585248, ..., -0.6569439 ,
            -0.5237557 , -0.3770282 ],
           [ 0.3933881 ,  0.20141809,  0.19791125, ..., -0.380811  ,
            -0.4921796 , -0.33902702],
           ...,
           [        nan,         nan,         nan, ..., -0.20507982,
            -0.21505156, -0.36355355],
           [        nan,         nan,         nan, ..., -0.12386315,
            -0.36160484, -0.2852007 ],
           [        nan,         nan,         nan, ..., -0.21303938,
            -0.25312763, -0.27285883]], dtype=float32)

  • y_sea_water_velocity
    (trajectory, time)
    float32
    0.3641 0.518 ... 0.8258 0.8737

    minval :

    -0.38523096

    maxval :

    1.3946526

    array([[0.36409444, 0.51796234, 0.5094887 , ..., 1.0736189 , 1.0909232 ,
            1.1545241 ],
           [0.50797606, 0.4592814 , 0.54768944, ..., 1.0296686 , 1.1706446 ,
            1.1853416 ],
           [0.5073473 , 0.5106243 , 0.47898614, ..., 0.7873976 , 1.1541189 ,
            1.1151277 ],
           ...,
           [       nan,        nan,        nan, ..., 0.6251363 , 0.7509022 ,
            0.8254842 ],
           [       nan,        nan,        nan, ..., 0.54717773, 0.7915698 ,
            0.93019843],
           [       nan,        nan,        nan, ..., 0.58333266, 0.8258305 ,
            0.87366754]], dtype=float32)

  • x_wind
    (trajectory, time)
    float32
    -8.748 -9.447 ... 6.493 7.163

    minval :

    -15.208915

    maxval :

    17.42325

    array([[ -8.748272 ,  -9.446504 ,  -8.983033 , ...,   4.8908305,
              3.7031386,   2.8885705],
           [ -8.874136 , -10.024847 ,  -9.447379 , ...,   4.574147 ,
              3.8485045,   3.1367726],
           [ -8.287212 ,  -9.56559  ,  -8.56232  , ...,   3.9473453,
              4.003879 ,   5.4945593],
           ...,
           [        nan,         nan,         nan, ...,   7.971263 ,
              6.9061637,   6.4345965],
           [        nan,         nan,         nan, ...,   7.128233 ,
              6.0720625,   6.757368 ],
           [        nan,         nan,         nan, ...,   6.501858 ,
              6.4932537,   7.1634073]], dtype=float32)

  • y_wind
    (trajectory, time)
    float32
    9.138 10.41 14.09 ... 16.27 16.6

    minval :

    -2.8727612

    maxval :

    21.256886

    array([[ 9.137849, 10.412806, 14.09369 , ..., 18.708199, 17.759605,
            16.924578],
           [ 9.804504, 11.794177, 14.172353, ..., 17.7678  , 17.40429 ,
            17.550827],
           [10.307045, 10.645775, 14.072082, ..., 18.322939, 18.912813,
            19.594948],
           ...,
           [      nan,       nan,       nan, ..., 16.229025, 16.359663,
            15.92441 ],
           [      nan,       nan,       nan, ..., 16.103004, 14.965876,
            15.693896],
           [      nan,       nan,       nan, ..., 15.16429 , 16.26992 ,
            16.59994 ]], dtype=float32)

  • sea_surface_height
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    minval :

    0.0

    maxval :

    0.0

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

  • upward_sea_water_velocity
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    minval :

    0.0

    maxval :

    0.0

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

  • sea_surface_wave_significant_height
    (trajectory, time)
    float32
    3.937 4.863 6.871 ... 7.549 8.041

    minval :

    1.667064

    maxval :

    12.234488

    array([[ 3.9368005,  4.862508 ,  6.8714433, ...,  9.198356 ,  8.096272 ,
             7.251715 ],
           [ 4.302013 ,  5.8941646,  7.136671 , ...,  8.28079  ,  7.815918 ,
             7.8196235],
           [ 4.3028607,  5.0388927,  6.6748843, ...,  8.642266 ,  9.193646 ,
            10.188143 ],
           ...,
           [       nan,        nan,        nan, ...,  8.042289 ,  7.7572083,
             7.256774 ],
           [       nan,        nan,        nan, ...,  7.6289134,  6.4168453,
             7.182225 ],
           [       nan,        nan,        nan, ...,  6.696854 ,  7.549066 ,
             8.04106  ]], dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

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

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

    minval :

    0.0

    maxval :

    0.0

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

  • sea_surface_wave_period_at_variance_spectral_density_maximum
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0

    minval :

    0.0

    maxval :

    0.0

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

  • sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment
    (trajectory, time)
    float32
    10.81 12.01 14.28 ... 14.97 15.45

    minval :

    7.0340548

    maxval :

    19.055593

    array([[10.809393, 12.013227, 14.280833, ..., 16.522839, 15.501444,
            14.67067 ],
           [11.299663, 13.226371, 14.553833, ..., 15.677091, 15.230689,
            15.234301],
           [11.300776, 12.229175, 14.075099, ..., 16.015606, 16.51861 ,
            17.3891  ],
           ...,
           [      nan,       nan,       nan, ..., 15.449677, 15.173379,
            14.675787],
           [      nan,       nan,       nan, ..., 15.047381, 13.800359,
            14.600211],
           [      nan,       nan,       nan, ..., 14.098244, 14.968429,
            15.448498]], dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

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

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

    minval :

    0.0

    maxval :

    0.0

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

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

    minval :

    0.0

    maxval :

    0.0

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

  • sea_water_temperature
    (trajectory, time)
    float32
    10.0 10.0 10.0 ... 10.0 10.0 283.1

    minval :

    10.0

    maxval :

    283.15

    array([[ 10.  ,  10.  ,  10.  , ...,  10.  ,  10.  , 283.15],
           [ 10.  ,  10.  ,  10.  , ...,  10.  ,  10.  , 283.15],
           [ 10.  ,  10.  ,  10.  , ...,  10.  ,  10.  , 283.15],
           ...,
           [   nan,    nan,    nan, ...,  10.  ,  10.  , 283.15],
           [   nan,    nan,    nan, ...,  10.  ,  10.  , 283.15],
           [   nan,    nan,    nan, ...,  10.  ,  10.  , 283.15]],
          dtype=float32)

  • sea_water_salinity
    (trajectory, time)
    float32
    34.0 34.0 34.0 ... 34.0 34.0 34.0

    minval :

    34.0

    maxval :

    34.0

    array([[34., 34., 34., ..., 34., 34., 34.],
           [34., 34., 34., ..., 34., 34., 34.],
           [34., 34., 34., ..., 34., 34., 34.],
           ...,
           [nan, nan, nan, ..., 34., 34., 34.],
           [nan, nan, nan, ..., 34., 34., 34.],
           [nan, nan, nan, ..., 34., 34., 34.]], dtype=float32)

  • sea_floor_depth_below_sea_level
    (trajectory, time)
    float32
    1e+04 1e+04 1e+04 ... 1e+04 1e+04

    minval :

    10000.0

    maxval :

    10000.0

    array([[10000., 10000., 10000., ..., 10000., 10000., 10000.],
           [10000., 10000., 10000., ..., 10000., 10000., 10000.],
           [10000., 10000., 10000., ..., 10000., 10000., 10000.],
           ...,
           [   nan,    nan,    nan, ..., 10000., 10000., 10000.],
           [   nan,    nan,    nan, ..., 10000., 10000., 10000.],
           [   nan,    nan,    nan, ..., 10000., 10000., 10000.]],
          dtype=float32)

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

    minval :

    0.02

    maxval :

    0.02

    array([[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, ..., 0.02, 0.02, 0.02],
           ...,
           [ nan,  nan,  nan, ..., 0.02, 0.02, 0.02],
           [ nan,  nan,  nan, ..., 0.02, 0.02, 0.02],
           [ nan,  nan,  nan, ..., 0.02, 0.02, 0.02]], dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

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

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

    minval :

    50.0

    maxval :

    50.0

    array([[50., 50., 50., ..., 50., 50., 50.],
           [50., 50., 50., ..., 50., 50., 50.],
           [50., 50., 50., ..., 50., 50., 50.],
           ...,
           [nan, nan, nan, ..., 50., 50., 50.],
           [nan, nan, nan, ..., 50., 50., 50.],
           [nan, nan, nan, ..., 50., 50., 50.]], dtype=float32)

• Indexes: (2)
  • trajectory
    PandasIndex

    PandasIndex(Index([   0,    1,    2,    3,    4,    5,    6,    7,    8,    9,
           ...
           4990, 4991, 4992, 4993, 4994, 4995, 4996, 4997, 4998, 4999],
          dtype='int64', name='trajectory', length=5000))

  • time
    PandasIndex

    PandasIndex(DatetimeIndex(['2015-11-16 00:00:00', '2015-11-16 01:00:00',
                   '2015-11-16 02:00:00', '2015-11-16 03:00:00',
                   '2015-11-16 04:00:00', '2015-11-16 05:00:00',
                   '2015-11-16 06:00:00', '2015-11-16 07:00:00',
                   '2015-11-16 08:00:00', '2015-11-16 09:00:00',
                   '2015-11-16 10:00:00', '2015-11-16 11:00:00',
                   '2015-11-16 12:00:00', '2015-11-16 13:00:00',
                   '2015-11-16 14:00:00', '2015-11-16 15:00:00',
                   '2015-11-16 16:00:00', '2015-11-16 17:00:00',
                   '2015-11-16 18:00:00', '2015-11-16 19:00:00',
                   '2015-11-16 20:00:00', '2015-11-16 21:00:00',
                   '2015-11-16 22:00:00', '2015-11-16 23:00:00',
                   '2015-11-17 00:00:00', '2015-11-17 01:00:00',
                   '2015-11-17 02:00:00', '2015-11-17 03:00:00',
                   '2015-11-17 04:00:00', '2015-11-17 05:00:00',
                   '2015-11-17 06:00:00', '2015-11-17 07:00:00',
                   '2015-11-17 08:00:00', '2015-11-17 09:00:00'],
                  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:22:59.655113

  date_created :

  2025-04-03T15:22:59.655131

  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(['/root/project/tests/test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc', '/root/project/tests/test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc', 'global_landmask'])

  time_coverage_start :

  2015-11-16 00:00:00

  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 :

  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 :

  None

  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 :

  True

  config_processes:evaporation :

  True

  config_processes:emulsification :

  True

  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 HEAVY CRUDE

  opendrift_version :

  1.13.1

  seed_lon :

  4.5

  seed_lat :

  60.0

  seed_radius :

  50.0

  seed_time :

  2015-11-16 00:00:00

  seed_number :

  5000.0

  seed_z :

  0.0

  seed_m3_per_hour :

  1.0

  seed_oiltype :

  GENERIC HEAVY CRUDE

  seed_geojson :

  {"features": [], "type": "FeatureCollection"}

  simulation_time :

  2025-04-03 15:22:59.660663

  reader_x_sea_water_velocity :

  ['/root/project/tests/test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc']

  reader_y_sea_water_velocity :

  ['/root/project/tests/test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc']

  reader_x_wind :

  ['/root/project/tests/test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc']

  reader_y_wind :

  ['/root/project/tests/test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc']

  reader_sea_surface_height :

  0

  reader_upward_sea_water_velocity :

  0

  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 :

  0

  reader_sea_ice_x_velocity :

  0

  reader_sea_ice_y_velocity :

  0

  reader_sea_water_temperature :

  10

  reader_sea_water_salinity :

  34

  reader_sea_floor_depth_below_sea_level :

  10000

  reader_ocean_vertical_diffusivity :

  0.02

  reader_land_binary_mask :

  ['global_landmask']

  reader_ocean_mixed_layer_thickness :

  50

  time_coverage_end :

  2015-11-17 09:00:00

  time_coverage_duration :

  P1DT9H0M0S

  time_coverage_resolution :

  P0DT0H30M0S

  performance :

  -------------------- Reader performance: -------------------- /root/project/tests/test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc 0:00:01.7 total 0:00:00.0 preparing 0:00:00.3 reading 0:00:00.1 interpolation 0:00:00.0 interpolation_time 0:00:01.3 rotating vectors 0:00:00.0 masking -------------------- /root/project/tests/test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc 0:00:01.6 total 0:00:00.0 preparing 0:00:00.2 reading 0:00:00.0 interpolation 0:00:00.0 interpolation_time 0:00:01.3 rotating vectors 0:00:00.0 masking -------------------- global_landmask 0:00:00.0 total 0:00:00.0 preparing 0:00:00.0 reading 0:00:00.0 masking -------------------- Performance: 12.9 total time 0.1 configuration 0.0 preparing main loop 0.0 moving elements to ocean 12.7 main loop 8.5 updating elements 3.9 oil weathering 0.0 updating viscosities 0.0 updating densities 3.8 evaporation 0.0 emulsification 0.0 dispersion 4.3 vertical mixing 0.0 cleaning up --------------------

  geospatial_bounds_crs :

  EPSG:4326

  geospatial_bounds_vertical_crs :

  EPSG:5831

  geospatial_lat_min :

  59.998623

  geospatial_lat_max :

  60.72471

  geospatial_lat_units :

  degrees_north

  geospatial_lat_resolution :

  point

  geospatial_lon_min :

  4.4580245

  geospatial_lon_max :

  4.965604

  geospatial_lon_units :

  degrees_east

  geospatial_lon_resolution :

  point

  runtime :

  0:00:12.914209

  geospatial_vertical_min :

  -48.601967

  geospatial_vertical_max :

  0.0

  geospatial_vertical_positive :

  up

Animate and compare the two runs. We see that there is much more stranding of oil when wind is considered.

o.animation(fast=True, compare=o2,
            legend=['Current + 3 % wind drift', 'Current only'])

15:23:12 WARNING opendrift:2429: Plotting fast. This will make your plots less accurate.
15:23:13 INFO    opendrift:4600: Saving animation to /root/project/docs/source/gallery/animations/example_compare_0.gif...
15:23:39 INFO    opendrift:3036: Time to make animation: 0:00:26.708819

o.plot(fast=True, compare=o2,
       legend=['Current + 3 % wind drift', 'Current only'])

15:23:39 WARNING opendrift:2429: Plotting fast. This will make your plots less accurate.

(<GeoAxes: title={'center': 'OpenDrift - OpenOil (GENERIC HEAVY CRUDE)\n2015-11-16 00:00 to 2015-11-17 09:00 UTC (34 steps)'}>, <Figure size 623.045x1100 with 1 Axes>)