Oil budget (NOAA)

from datetime import datetime
from opendrift import test_data_folder as tdf
from opendrift.readers import reader_netCDF_CF_generic
from opendrift.models.openoil import OpenOil

o = OpenOil(loglevel=20, weathering_model='noaa')

# Using constand wind and current
#o.set_config('environment:fallback:x_wind', 7)
#o.set_config('environment:fallback:x_sea_water_velocity', .7)
#o.set_config('environment:fallback:y_sea_water_velocity', .3)
#o.set_config('environment:fallback:land_binary_mask', 0)

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

o.add_reader([reader_arome, reader_norkyst])

17:03:22 INFO    opendrift:569: OpenDriftSimulation initialised (version 1.14.11 / v1.14.12-66-g07283eb)
17:03:22 INFO    opendrift.readers:67: Opening file with xr.open_dataset
17:03:23 INFO    opendrift.readers.reader_netCDF_CF_generic:340: Detected dimensions: {'time': 'time', 'x': 'x', 'y': 'y'}
17:03:23 INFO    opendrift.readers:67: Opening file with xr.open_dataset
17:03:23 INFO    opendrift.readers.reader_netCDF_CF_generic:340: Detected dimensions: {'x': 'X', 'y': 'Y', 'z': 'depth', 'time': 'time'}

Adjusting some configuration

o.set_config('processes:dispersion', False)
o.set_config('processes:evaporation', True)
o.set_config('processes:emulsification', True)
o.set_config('drift:vertical_mixing', True)

Seeding some oil particles

oil_type='MARTIN LINGE CRUDE 2016'
o.seed_elements(lon=4.88, lat=60.1, z=0, radius=3000, number=500,
                time=reader_norkyst.start_time, oil_type=oil_type)

17:03:23 INFO    opendrift.models.openoil.adios.dirjs:86: Querying ADIOS database for oil: MARTIN LINGE CRUDE 2016
17:03:23 INFO    opendrift.models.openoil.openoil:1731: Using density 930.00252 and viscosity 0.0004141770017829454 of oiltype MARTIN LINGE CRUDE 2016
17:03:23 INFO    opendrift.models.basemodel.environment:203: Adding a global landmask from GSHHG
17:03:28 INFO    opendrift.models.basemodel.environment:227: Fallback values will be used for the following variables which have no readers:
17:03:28 INFO    opendrift.models.basemodel.environment:230:    sea_surface_height: 0.000000
17:03:28 INFO    opendrift.models.basemodel.environment:230:    upward_sea_water_velocity: 0.000000
17:03:28 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_significant_height: 0.000000
17:03:28 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
17:03:28 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
17:03:28 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
17:03:28 INFO    opendrift.models.basemodel.environment:230:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
17:03:28 INFO    opendrift.models.basemodel.environment:230:    sea_ice_area_fraction: 0.000000
17:03:28 INFO    opendrift.models.basemodel.environment:230:    sea_ice_x_velocity: 0.000000
17:03:28 INFO    opendrift.models.basemodel.environment:230:    sea_ice_y_velocity: 0.000000
17:03:28 INFO    opendrift.models.basemodel.environment:230:    sea_water_temperature: 10.000000
17:03:28 INFO    opendrift.models.basemodel.environment:230:    sea_water_salinity: 34.000000
17:03:28 INFO    opendrift.models.basemodel.environment:230:    sea_floor_depth_below_sea_level: 10000.000000
17:03:28 INFO    opendrift.models.basemodel.environment:230:    horizontal_diffusivity: 0.000000
17:03:28 INFO    opendrift.models.basemodel.environment:230:    ocean_vertical_diffusivity: 0.020000
17:03:28 INFO    opendrift.models.basemodel.environment:230:    ocean_mixed_layer_thickness: 50.000000

Running model

o.run(steps=4*24, time_step=900, time_step_output=3600)

17:03:28 INFO    opendrift:1911: Skipping environment variable upward_sea_water_velocity because of condition ['drift:vertical_advection', 'is', False]
17:03:28 INFO    opendrift:1922: Storing previous values of element property lon because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
17:03:28 INFO    opendrift:1922: Storing previous values of element property lat because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
17:03:28 INFO    opendrift:958: Using existing reader for land_binary_mask to move elements to ocean
17:03:28 INFO    opendrift:990: All points are in ocean
17:03:28 INFO    opendrift.models.openoil.openoil:697: Oil-water surface tension is 0.033729 Nm
17:03:28 INFO    opendrift.models.openoil.openoil:706: Using max water fractions [0.72, 0.78] for temperatures [5.0, 15.0] for oiltype MARTIN LINGE CRUDE 2016
17:03:28 INFO    opendrift.models.openoil.openoil:707: Corresponding max water fraction from GNOME is 0.731347796266117
17:03:28 INFO    opendrift:2238: 2015-11-16 00:00:00 - step 1 of 96 - 500 active elements (0 deactivated)
17:03:28 INFO    opendrift:2238: 2015-11-16 00:15:00 - step 2 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 00:30:00 - step 3 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 00:45:00 - step 4 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 01:00:00 - step 5 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 01:15:00 - step 6 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 01:30:00 - step 7 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 01:45:00 - step 8 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 02:00:00 - step 9 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 02:15:00 - step 10 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 02:30:00 - step 11 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 02:45:00 - step 12 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 03:00:00 - step 13 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 03:15:00 - step 14 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 03:30:00 - step 15 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 03:45:00 - step 16 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 04:00:00 - step 17 of 96 - 500 active elements (0 deactivated)
17:03:29 INFO    opendrift:2238: 2015-11-16 04:15:00 - step 18 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 04:30:00 - step 19 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 04:45:00 - step 20 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 05:00:00 - step 21 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 05:15:00 - step 22 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 05:30:00 - step 23 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 05:45:00 - step 24 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 06:00:00 - step 25 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 06:15:00 - step 26 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 06:30:00 - step 27 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 06:45:00 - step 28 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 07:00:00 - step 29 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 07:15:00 - step 30 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 07:30:00 - step 31 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 07:45:00 - step 32 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 08:00:00 - step 33 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 08:15:00 - step 34 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 08:30:00 - step 35 of 96 - 500 active elements (0 deactivated)
17:03:30 INFO    opendrift:2238: 2015-11-16 08:45:00 - step 36 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 09:00:00 - step 37 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 09:15:00 - step 38 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 09:30:00 - step 39 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 09:45:00 - step 40 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 10:00:00 - step 41 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 10:15:00 - step 42 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 10:30:00 - step 43 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 10:45:00 - step 44 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 11:00:00 - step 45 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 11:15:00 - step 46 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 11:30:00 - step 47 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 11:45:00 - step 48 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 12:00:00 - step 49 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 12:15:00 - step 50 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 12:30:00 - step 51 of 96 - 500 active elements (0 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 12:45:00 - step 52 of 96 - 499 active elements (1 deactivated)
17:03:31 INFO    opendrift:2238: 2015-11-16 13:00:00 - step 53 of 96 - 499 active elements (1 deactivated)
17:03:32 INFO    opendrift:2238: 2015-11-16 13:15:00 - step 54 of 96 - 499 active elements (1 deactivated)
17:03:32 INFO    opendrift:2238: 2015-11-16 13:30:00 - step 55 of 96 - 499 active elements (1 deactivated)
17:03:32 INFO    opendrift:2238: 2015-11-16 13:45:00 - step 56 of 96 - 499 active elements (1 deactivated)
17:03:32 INFO    opendrift:2238: 2015-11-16 14:00:00 - step 57 of 96 - 499 active elements (1 deactivated)
17:03:32 INFO    opendrift:2238: 2015-11-16 14:15:00 - step 58 of 96 - 497 active elements (3 deactivated)
17:03:32 INFO    opendrift:2238: 2015-11-16 14:30:00 - step 59 of 96 - 495 active elements (5 deactivated)
17:03:32 INFO    opendrift:2238: 2015-11-16 14:45:00 - step 60 of 96 - 494 active elements (6 deactivated)
17:03:32 INFO    opendrift:2238: 2015-11-16 15:00:00 - step 61 of 96 - 494 active elements (6 deactivated)
17:03:32 INFO    opendrift:2238: 2015-11-16 15:15:00 - step 62 of 96 - 494 active elements (6 deactivated)
17:03:32 INFO    opendrift:2238: 2015-11-16 15:30:00 - step 63 of 96 - 490 active elements (10 deactivated)
17:03:32 INFO    opendrift:2238: 2015-11-16 15:45:00 - step 64 of 96 - 489 active elements (11 deactivated)
17:03:32 INFO    opendrift:2238: 2015-11-16 16:00:00 - step 65 of 96 - 488 active elements (12 deactivated)
17:03:32 INFO    opendrift:2238: 2015-11-16 16:15:00 - step 66 of 96 - 485 active elements (15 deactivated)
17:03:32 INFO    opendrift:2238: 2015-11-16 16:30:00 - step 67 of 96 - 485 active elements (15 deactivated)
17:03:33 INFO    opendrift:2238: 2015-11-16 16:45:00 - step 68 of 96 - 483 active elements (17 deactivated)
17:03:33 INFO    opendrift:2238: 2015-11-16 17:00:00 - step 69 of 96 - 480 active elements (20 deactivated)
17:03:33 INFO    opendrift:2238: 2015-11-16 17:15:00 - step 70 of 96 - 477 active elements (23 deactivated)
17:03:33 INFO    opendrift:2238: 2015-11-16 17:30:00 - step 71 of 96 - 471 active elements (29 deactivated)
17:03:33 INFO    opendrift:2238: 2015-11-16 17:45:00 - step 72 of 96 - 463 active elements (37 deactivated)
17:03:33 INFO    opendrift:2238: 2015-11-16 18:00:00 - step 73 of 96 - 457 active elements (43 deactivated)
17:03:33 INFO    opendrift:2238: 2015-11-16 18:15:00 - step 74 of 96 - 451 active elements (49 deactivated)
17:03:33 INFO    opendrift:2238: 2015-11-16 18:30:00 - step 75 of 96 - 446 active elements (54 deactivated)
17:03:33 INFO    opendrift:2238: 2015-11-16 18:45:00 - step 76 of 96 - 437 active elements (63 deactivated)
17:03:33 INFO    opendrift:2238: 2015-11-16 19:00:00 - step 77 of 96 - 431 active elements (69 deactivated)
17:03:33 INFO    opendrift:2238: 2015-11-16 19:15:00 - step 78 of 96 - 425 active elements (75 deactivated)
17:03:33 INFO    opendrift:2238: 2015-11-16 19:30:00 - step 79 of 96 - 416 active elements (84 deactivated)
17:03:33 INFO    opendrift:2238: 2015-11-16 19:45:00 - step 80 of 96 - 414 active elements (86 deactivated)
17:03:33 WARNING opendrift.readers.basereader.structured:324: Data block from /root/project/tests/test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc not large enough to cover element positions within timestep. Buffer size (4) must be increased. See `Variables.set_buffer_size`.
17:03:34 INFO    opendrift:2238: 2015-11-16 20:00:00 - step 81 of 96 - 411 active elements (89 deactivated)
17:03:34 INFO    opendrift:2238: 2015-11-16 20:15:00 - step 82 of 96 - 405 active elements (95 deactivated)
17:03:34 INFO    opendrift:2238: 2015-11-16 20:30:00 - step 83 of 96 - 402 active elements (98 deactivated)
17:03:34 INFO    opendrift:2238: 2015-11-16 20:45:00 - step 84 of 96 - 398 active elements (102 deactivated)
17:03:34 INFO    opendrift:2238: 2015-11-16 21:00:00 - step 85 of 96 - 394 active elements (106 deactivated)
17:03:34 INFO    opendrift:2238: 2015-11-16 21:15:00 - step 86 of 96 - 392 active elements (108 deactivated)
17:03:34 INFO    opendrift:2238: 2015-11-16 21:30:00 - step 87 of 96 - 385 active elements (115 deactivated)
17:03:34 INFO    opendrift:2238: 2015-11-16 21:45:00 - step 88 of 96 - 378 active elements (122 deactivated)
17:03:34 INFO    opendrift:2238: 2015-11-16 22:00:00 - step 89 of 96 - 372 active elements (128 deactivated)
17:03:34 INFO    opendrift:2238: 2015-11-16 22:15:00 - step 90 of 96 - 366 active elements (134 deactivated)
17:03:34 INFO    opendrift:2238: 2015-11-16 22:30:00 - step 91 of 96 - 360 active elements (140 deactivated)
17:03:34 INFO    opendrift:2238: 2015-11-16 22:45:00 - step 92 of 96 - 356 active elements (144 deactivated)
17:03:34 INFO    opendrift:2238: 2015-11-16 23:00:00 - step 93 of 96 - 351 active elements (149 deactivated)
17:03:35 INFO    opendrift:2238: 2015-11-16 23:15:00 - step 94 of 96 - 344 active elements (156 deactivated)
17:03:35 INFO    opendrift:2238: 2015-11-16 23:30:00 - step 95 of 96 - 333 active elements (167 deactivated)
17:03:35 INFO    opendrift:2238: 2015-11-16 23:45:00 - step 96 of 96 - 329 active elements (171 deactivated)

<xarray.Dataset> Size: 2MB
Dimensions:                                                                              (
                                                                                          trajectory: 500,
                                                                                          time: 25)
Coordinates:
  * trajectory                                                                           (trajectory) int64 4kB ...
  * time                                                                                 (time) datetime64[ns] 200B ...
Data variables: (12/44)
    status                                                                               (trajectory, time) float32 50kB ...
    moving                                                                               (trajectory, time) float32 50kB ...
    age_seconds                                                                          (trajectory, time) float32 50kB ...
    origin_marker                                                                        (trajectory, time) float32 50kB ...
    lon                                                                                  (trajectory, time) float32 50kB ...
    lat                                                                                  (trajectory, time) float32 50kB ...
    ...                                                                                   ...
    sea_water_salinity                                                                   (trajectory, time) float32 50kB ...
    sea_floor_depth_below_sea_level                                                      (trajectory, time) float32 50kB ...
    horizontal_diffusivity                                                               (trajectory, time) float32 50kB ...
    ocean_vertical_diffusivity                                                           (trajectory, time) float32 50kB ...
    land_binary_mask                                                                     (trajectory, time) float32 50kB ...
    ocean_mixed_layer_thickness                                                          (trajectory, time) float32 50kB ...
Attributes: (12/167)
    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: 500
  • time: 25
• Coordinates: (2)
  • trajectory
    (trajectory)
    int64
    0 1 2 3 4 5 ... 495 496 497 498 499

    cf_role :

    trajectory_id

    dtype :

    <class 'numpy.int32'>

    array([  0,   1,   2, ..., 497, 498, 499], shape=(500,))

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

    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'], dtype='datetime64[ns]')

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

    dtype :

    <class 'numpy.int32'>

    valid_range :

    [0 1]

    flag_values :

    [0 1]

    flag_meanings :

    active stranded

    array([[ 0.,  0.,  0., ...,  1., nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  1., nan],
           [ 0.,  0.,  0., ...,  1., nan, nan]],
          shape=(500, 25), dtype=float32)

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

    dtype :

    <class 'numpy.int32'>

    minval :

    0

    maxval :

    1

    array([[ 1.,  1.,  1., ...,  0., nan, nan],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ..., nan, nan, nan],
           ...,
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ...,  1.,  0., nan],
           [ 1.,  1.,  1., ...,  0., nan, nan]],
          shape=(500, 25), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    s

    minval :

    0.0

    maxval :

    86400.0

    array([[    0.,  3600.,  7200., ..., 78300.,    nan,    nan],
           [    0.,  3600.,  7200., ..., 79200., 82800., 86400.],
           [    0.,  3600.,  7200., ...,    nan,    nan,    nan],
           ...,
           [    0.,  3600.,  7200., ..., 79200., 82800., 86400.],
           [    0.,  3600.,  7200., ..., 79200., 80100.,    nan],
           [    0.,  3600.,  7200., ..., 78300.,    nan,    nan]],
          shape=(500, 25), dtype=float32)

  • origin_marker
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 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., nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0., nan, nan]],
          shape=(500, 25), dtype=float32)

  • lon
    (trajectory, time)
    float32
    4.975 4.956 4.95 ... 4.83 nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    axis :

    X

    minval :

    4.4456744

    maxval :

    5.039443

    array([[4.975158 , 4.955796 , 4.94978  , ..., 4.9056845,       nan,
                  nan],
           [4.9015784, 4.8976707, 4.9055676, ..., 4.6928587, 4.6793685,
            4.663991 ],
           [4.932827 , 4.927009 , 4.9286203, ...,       nan,       nan,
                  nan],
           ...,
           [4.8522673, 4.8636813, 4.872975 , ..., 4.848733 , 4.858781 ,
            4.8733835],
           [4.8251524, 4.8337426, 4.840465 , ..., 4.9533596, 4.9553313,
                  nan],
           [4.875803 , 4.87272  , 4.869493 , ..., 4.8297772,       nan,
                  nan]], shape=(500, 25), dtype=float32)

  • lat
    (trajectory, time)
    float32
    60.11 60.12 60.13 ... 60.49 nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    axis :

    Y

    minval :

    60.017963

    maxval :

    60.857864

    array([[60.11027 , 60.118885, 60.12617 , ..., 60.407158,       nan,
                  nan],
           [60.099075, 60.111034, 60.120705, ..., 60.52879 , 60.55287 ,
            60.576794],
           [60.12951 , 60.137405, 60.143715, ...,       nan,       nan,
                  nan],
           ...,
           [60.102535, 60.118565, 60.131996, ..., 60.47669 , 60.48363 ,
            60.49152 ],
           [60.069088, 60.086052, 60.104355, ..., 60.417583, 60.419556,
                  nan],
           [60.09036 , 60.104614, 60.122284, ..., 60.489468,       nan,
                  nan]], shape=(500, 25), dtype=float32)

  • z
    (trajectory, time)
    float32
    0.0 0.0 -3.178 ... -26.31 nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    m

    standard_name :

    z

    long_name :

    vertical position

    axis :

    Z

    positive :

    up

    minval :

    -49.75206

    maxval :

    0.0

    array([[  0.       ,   0.       ,  -3.1781898, ..., -13.103731 ,
                    nan,         nan],
           [  0.       ,   0.       , -21.838715 , ..., -44.324642 ,
            -17.489655 , -13.550843 ],
           [  0.       ,  -3.9620168, -11.072356 , ...,         nan,
                    nan,         nan],
           ...,
           [  0.       , -26.904764 , -14.055898 , ...,  -5.0461526,
            -10.819094 ,  -2.7984245],
           [  0.       , -23.107323 ,   0.       , ..., -48.159054 ,
            -45.757393 ,         nan],
           [  0.       ,   0.       ,  -2.9442632, ..., -26.307129 ,
                    nan,         nan]], shape=(500, 25), dtype=float32)

  • wind_drift_factor
    (trajectory, time)
    float32
    0.03 0.03 0.03 ... 0.03 nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    1

    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,  nan,  nan],
           [0.03, 0.03, 0.03, ..., 0.03, 0.03, 0.03],
           [0.03, 0.03, 0.03, ...,  nan,  nan,  nan],
           ...,
           [0.03, 0.03, 0.03, ..., 0.03, 0.03, 0.03],
           [0.03, 0.03, 0.03, ..., 0.03, 0.03,  nan],
           [0.03, 0.03, 0.03, ..., 0.03,  nan,  nan]],
          shape=(500, 25), dtype=float32)

  • current_drift_factor
    (trajectory, time)
    float32
    1.0 1.0 1.0 1.0 ... 1.0 1.0 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., nan, nan],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ..., nan, nan, nan],
           ...,
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ...,  1.,  1., nan],
           [ 1.,  1.,  1., ...,  1., nan, nan]],
          shape=(500, 25), dtype=float32)

  • terminal_velocity
    (trajectory, time)
    float32
    0.0 0.0 0.006431 ... nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    m/s

    description :

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

    minval :

    0.0

    maxval :

    0.08517955

    array([[0.        , 0.        , 0.00643122, ..., 0.00756295,        nan,
                   nan],
           [0.        , 0.06005226, 0.00242527, ..., 0.00029005, 0.00029005,
            0.00029005],
           [0.        , 0.07118084, 0.00134592, ...,        nan,        nan,
                   nan],
           ...,
           [0.        , 0.00060199, 0.00060199, ..., 0.00060199, 0.00060199,
            0.00060199],
           [0.        , 0.00803695, 0.00803695, ..., 0.00141008, 0.00141008,
                   nan],
           [0.        , 0.0355003 , 0.01889891, ..., 0.00337315,        nan,
                   nan]], shape=(500, 25), dtype=float32)

  • mass_oil
    (trajectory, time)
    float32
    1.86 1.823 1.817 ... 1.815 nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    1.7781155

    maxval :

    1.860005

    array([[1.860005 , 1.8231219, 1.8170761, ..., 1.7934159,       nan,
                  nan],
           [1.860005 , 1.8224559, 1.81851  , ..., 1.8139673, 1.8139673,
            1.8139673],
           [1.860005 , 1.8228952, 1.8206681, ...,       nan,       nan,
                  nan],
           ...,
           [1.860005 , 1.8298123, 1.8298123, ..., 1.8298123, 1.8298123,
            1.8298123],
           [1.860005 , 1.8288217, 1.8288217, ..., 1.802007 , 1.802007 ,
                  nan],
           [1.860005 , 1.8249747, 1.8178513, ..., 1.8153247,       nan,
                  nan]], shape=(500, 25), dtype=float32)

  • viscosity
    (trajectory, time)
    float32
    0.0004142 0.0005839 ... nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    m2/s

    description :

    Kinematic viscosity of oil (emulsion)

    minval :

    0.000414177

    maxval :

    0.00075598055

    array([[0.00041418, 0.00058386, 0.00061311, ..., 0.00069628,        nan,
                   nan],
           [0.00041418, 0.00058687, 0.00060841, ..., 0.00062345, 0.00062345,
            0.00062345],
           [0.00041418, 0.00058403, 0.00060139, ...,        nan,        nan,
                   nan],
           ...,
           [0.00041418, 0.00057254, 0.00057254, ..., 0.00057254, 0.00057254,
            0.00057254],
           [0.00041418, 0.00057559, 0.00057559, ..., 0.00066485, 0.00066485,
                   nan],
           [0.00041418, 0.00058762, 0.00060431, ..., 0.00061891,        nan,
                   nan]], shape=(500, 25), dtype=float32)

  • density
    (trajectory, time)
    float32
    930.0 930.0 930.0 ... 930.0 nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    kg/m^3

    minval :

    930.0025

    maxval :

    930.004

    array([[930.0025, 930.004 , 930.004 , ..., 930.004 ,      nan,      nan],
           [930.0025, 930.004 , 930.004 , ..., 930.004 , 930.004 , 930.004 ],
           [930.0025, 930.004 , 930.004 , ...,      nan,      nan,      nan],
           ...,
           [930.0025, 930.004 , 930.004 , ..., 930.004 , 930.004 , 930.004 ],
           [930.0025, 930.004 , 930.004 , ..., 930.004 , 930.004 ,      nan],
           [930.0025, 930.004 , 930.004 , ..., 930.004 ,      nan,      nan]],
          shape=(500, 25), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    s

    minval :

    0.0

    maxval :

    0.0

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

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

    dtype :

    <class 'numpy.float32'>

    units :

    m2

    minval :

    0.0

    maxval :

    0.0

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

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

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.0

    maxval :

    0.0

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

  • mass_evaporated
    (trajectory, time)
    float32
    0.0 0.03688 0.04293 ... nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.0

    maxval :

    0.08188945

    array([[0.        , 0.03688309, 0.04292899, ..., 0.06658912,        nan,
                   nan],
           [0.        , 0.03754914, 0.04149498, ..., 0.04603761, 0.04603761,
            0.04603761],
           [0.        , 0.0371099 , 0.03933695, ...,        nan,        nan,
                   nan],
           ...,
           [0.        , 0.03019271, 0.03019271, ..., 0.03019271, 0.03019271,
            0.03019271],
           [0.        , 0.03118336, 0.03118336, ..., 0.05799801, 0.05799801,
                   nan],
           [0.        , 0.03503038, 0.04215366, ..., 0.04468033,        nan,
                   nan]], shape=(500, 25), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    kg

    minval :

    0.0

    maxval :

    0.0

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

  • biodegradation_half_time_droplet
    (trajectory, time)
    float32
    1.0 1.0 1.0 1.0 ... 1.0 1.0 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., nan, nan],
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ..., nan, nan, nan],
           ...,
           [ 1.,  1.,  1., ...,  1.,  1.,  1.],
           [ 1.,  1.,  1., ...,  1.,  1., nan],
           [ 1.,  1.,  1., ...,  1., nan, nan]],
          shape=(500, 25), dtype=float32)

  • biodegradation_half_time_slick
    (trajectory, time)
    float32
    3.0 3.0 3.0 3.0 ... 3.0 3.0 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., nan, nan],
           [ 3.,  3.,  3., ...,  3.,  3.,  3.],
           [ 3.,  3.,  3., ..., nan, nan, nan],
           ...,
           [ 3.,  3.,  3., ...,  3.,  3.,  3.],
           [ 3.,  3.,  3., ...,  3.,  3., nan],
           [ 3.,  3.,  3., ...,  3., nan, nan]],
          shape=(500, 25), dtype=float32)

  • fraction_evaporated
    (trajectory, time)
    float32
    0.0 0.01819 0.02308 ... nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    1

    minval :

    0.0

    maxval :

    0.044026468

    array([[0.        , 0.01819093, 0.02308004, ..., 0.0358005 ,        nan,
                   nan],
           [0.        , 0.01870579, 0.02230907, ..., 0.02475134, 0.02475134,
            0.02475134],
           [0.        , 0.01822083, 0.02114884, ...,        nan,        nan,
                   nan],
           ...,
           [0.        , 0.0162326 , 0.0162326 , ..., 0.0162326 , 0.0162326 ,
            0.0162326 ],
           [0.        , 0.0167652 , 0.0167652 , ..., 0.03118164, 0.03118164,
                   nan],
           [0.        , 0.01883349, 0.02163387, ..., 0.02402162,        nan,
                   nan]], shape=(500, 25), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    1

    minval :

    0.0

    maxval :

    0.0

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

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

    dtype :

    <class 'numpy.float32'>

    units :

    m

    minval :

    0.001

    maxval :

    0.001

    array([[0.001, 0.001, 0.001, ..., 0.001,   nan,   nan],
           [0.001, 0.001, 0.001, ..., 0.001, 0.001, 0.001],
           [0.001, 0.001, 0.001, ...,   nan,   nan,   nan],
           ...,
           [0.001, 0.001, 0.001, ..., 0.001, 0.001, 0.001],
           [0.001, 0.001, 0.001, ..., 0.001, 0.001,   nan],
           [0.001, 0.001, 0.001, ..., 0.001,   nan,   nan]],
          shape=(500, 25), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    m

    minval :

    0.0

    maxval :

    0.0029595434

    array([[0.0000000e+00, 0.0000000e+00, 4.1356386e-04, ..., 4.4847824e-04,
                      nan,           nan],
           [0.0000000e+00, 1.4710003e-03, 2.5396590e-04, ..., 8.7827131e-05,
            8.7827131e-05, 8.7827131e-05],
           [0.0000000e+00, 2.0673671e-03, 1.8919344e-04, ...,           nan,
                      nan,           nan],
           ...,
           [0.0000000e+00, 1.2652927e-04, 1.2652927e-04, ..., 1.2652927e-04,
            1.2652927e-04, 1.2652927e-04],
           [0.0000000e+00, 4.6231863e-04, 4.6231863e-04, ..., 1.9364996e-04,
            1.9364996e-04,           nan],
           [0.0000000e+00, 9.7165431e-04, 7.5132353e-04, ..., 2.9951177e-04,
                      nan,           nan]], shape=(500, 25), dtype=float32)

  • x_sea_water_velocity
    (trajectory, time)
    float32
    -0.01505 -0.01725 ... nan nan

    minval :

    -0.5844203

    maxval :

    0.51184624

    array([[-1.50461467e-02, -1.72463357e-02, -2.34047431e-04, ...,
             1.29799470e-01,             nan,             nan],
           [ 2.53742516e-01,  2.46543452e-01,  1.00513153e-01, ...,
            -2.11540937e-01, -2.42315277e-01, -2.54002839e-01],
           [ 1.00974575e-01,  9.35190096e-02,  6.90075532e-02, ...,
                        nan,             nan,             nan],
           ...,
           [ 2.38440365e-01,  1.08330570e-01,  8.79487693e-02, ...,
             1.69723049e-01,  2.09516898e-01,  2.77311057e-01],
           [ 1.48586914e-01,  2.31852010e-01,  9.68209505e-02, ...,
             1.20701529e-01,  9.72386450e-02,             nan],
           [ 1.64890945e-01,  9.59495381e-02,  1.18995897e-01, ...,
             1.09491639e-01,             nan,             nan]],
          shape=(500, 25), dtype=float32)

  • y_sea_water_velocity
    (trajectory, time)
    float32
    0.09734 0.1887 0.08978 ... nan nan

    minval :

    -0.26684105

    maxval :

    1.3990502

    array([[0.09733942, 0.18868227, 0.08977842, ..., 0.23964447,        nan,
                   nan],
           [0.14077923, 0.1713915 , 0.17812297, ..., 0.63268757, 0.711916  ,
            0.713946  ],
           [0.15881118, 0.18909352, 0.08885948, ...,        nan,        nan,
                   nan],
           ...,
           [0.4853535 , 0.37503147, 0.39718142, ..., 0.21499637, 0.25574386,
            0.18410066],
           [0.42777485, 0.48688143, 0.49083236, ..., 0.24432613, 0.2261381 ,
                   nan],
           [0.21463431, 0.34166428, 0.3611865 , ..., 0.22694081,        nan,
                   nan]], shape=(500, 25), dtype=float32)

  • x_wind
    (trajectory, time)
    float32
    -8.657 -10.1 -8.91 ... nan nan

    minval :

    -15.133531

    maxval :

    20.457739

    array([[ -8.657231 , -10.104213 ,  -8.909777 , ...,   4.0527873,
                    nan,         nan],
           [ -8.152298 ,  -9.309067 ,  -7.783069 , ...,   4.47048  ,
              4.752893 ,   4.2170067],
           [ -7.6702156, -10.078344 ,  -8.270208 , ...,         nan,
                    nan,         nan],
           ...,
           [ -7.1102133,  -9.170782 ,  -8.548274 , ...,   4.7771564,
              4.7016654,   4.607498 ],
           [ -7.665609 ,  -8.580118 ,  -7.9489083, ...,   2.5318959,
              2.7571797,         nan],
           [ -8.061482 ,  -7.630372 ,  -8.101029 , ...,   3.4641743,
                    nan,         nan]], shape=(500, 25), dtype=float32)

  • y_wind
    (trajectory, time)
    float32
    4.204 7.347 9.057 ... 15.09 nan nan

    minval :

    -2.3921654

    maxval :

    18.849596

    array([[ 4.2043343,  7.3468723,  9.056636 , ..., 14.656202 ,        nan,
                   nan],
           [ 6.5876226,  7.7234125,  8.995858 , ..., 14.879236 , 16.516947 ,
            16.411995 ],
           [ 6.1272945,  6.357647 ,  9.321331 , ...,        nan,        nan,
                   nan],
           ...,
           [ 7.2361493,  7.4399505, 10.434131 , ..., 14.867546 , 15.590568 ,
            15.917381 ],
           [ 7.2700877,  9.282365 , 10.121726 , ..., 11.372083 , 11.941132 ,
                   nan],
           [ 7.1522975,  8.999196 , 11.018624 , ..., 15.087062 ,        nan,
                   nan]], shape=(500, 25), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

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

  • sea_surface_wave_significant_height
    (trajectory, time)
    float32
    2.279 3.839 3.971 ... 5.895 nan nan

    minval :

    1.2112613

    maxval :

    10.894341

    array([[2.2785523, 3.8393629, 3.9706063, ..., 5.688242 ,       nan,
                  nan],
           [2.7024755, 3.5992212, 3.4809394, ..., 5.937871 , 7.2668285,
            7.0635624],
           [2.3708482, 3.4930205, 3.8199756, ...,       nan,       nan,
                  nan],
           ...,
           [2.5317578, 3.43062  , 4.4758244, ..., 5.9990826, 6.523218 ,
            6.754965 ],
           [2.7457469, 3.9306054, 4.0746083, ..., 3.3390749, 3.6947389,
                  nan],
           [2.857114 , 3.424519 , 4.601104 , ..., 5.89465  ,       nan,
                  nan]], shape=(500, 25), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

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

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

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0., nan, nan],
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ..., nan, nan, nan],
           ...,
           [ 0.,  0.,  0., ...,  0.,  0.,  0.],
           [ 0.,  0.,  0., ...,  0.,  0., nan],
           [ 0.,  0.,  0., ...,  0., nan, nan]],
          shape=(500, 25), 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 nan nan

    minval :

    0.0

    maxval :

    0.0

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

  • sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment
    (trajectory, time)
    float32
    8.224 10.67 10.86 ... 13.23 nan nan

    minval :

    5.9958186

    maxval :

    17.981672

    array([[ 8.223545 , 10.674785 , 10.855704 , ..., 12.993274 ,        nan,
                   nan],
           [ 8.955926 , 10.335557 , 10.164309 , ..., 13.275319 , 14.685949 ,
            14.479098 ],
           [ 8.388445 , 10.1819315, 10.6477995, ...,        nan,        nan,
                   nan],
           ...,
           [ 8.668435 , 10.090574 , 11.525668 , ..., 13.343569 , 13.914273 ,
            14.159279 ],
           [ 9.02734  , 10.800885 , 10.996957 , ...,  9.955031 , 10.471804 ,
                   nan],
           [ 9.208594 , 10.081598 , 11.685859 , ..., 13.226916 ,        nan,
                   nan]], shape=(500, 25), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

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

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

    minval :

    0.0

    maxval :

    0.0

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

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

    minval :

    0.0

    maxval :

    0.0

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

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

    minval :

    10.0

    maxval :

    283.15

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

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

    minval :

    34.0

    maxval :

    34.0

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

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

    minval :

    10000.0

    maxval :

    10000.0

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

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

    minval :

    0.0

    maxval :

    0.0

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

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

    minval :

    0.02

    maxval :

    0.02

    array([[0.02, 0.02, 0.02, ..., 0.02,  nan,  nan],
           [0.02, 0.02, 0.02, ..., 0.02, 0.02, 0.02],
           [0.02, 0.02, 0.02, ...,  nan,  nan,  nan],
           ...,
           [0.02, 0.02, 0.02, ..., 0.02, 0.02, 0.02],
           [0.02, 0.02, 0.02, ..., 0.02, 0.02,  nan],
           [0.02, 0.02, 0.02, ..., 0.02,  nan,  nan]],
          shape=(500, 25), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

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

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

    minval :

    50.0

    maxval :

    50.0

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

• Attributes: (167)

  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 2026-10-26 17:03:28.890327

  date_created :

  2026-10-26T17:03:28.890336

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

  time_coverage_start :

  2015-11-16 00:00:00

  time_step_calculation :

  0:15: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:horizontal_diffusivity :

  None

  config_environment:fallback:horizontal_diffusivity :

  0

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

  50

  config_drift:wind_uncertainty :

  0.5

  config_drift:relative_wind :

  False

  config_drift:deactivate_north_of :

  None

  config_drift:deactivate_south_of :

  None

  config_drift:deactivate_east_of :

  None

  config_drift:deactivate_west_of :

  None

  config_seed:origin_marker :

  0

  config_seed:z :

  0

  config_seed:wind_drift_factor :

  0.03

  config_seed:current_drift_factor :

  1

  config_seed:terminal_velocity :

  0.0

  config_seed:oil_film_thickness :

  0.001

  config_drift:vertical_advection :

  False

  config_drift:vertical_advection_at_surface :

  False

  config_drift:water_column_stretching :

  False

  config_drift:vertical_advection_correction :

  False

  config_drift:vertical_mixing :

  True

  config_drift:vertical_mixing_at_surface :

  False

  config_vertical_mixing:timestep :

  60

  config_vertical_mixing:diffusivitymodel :

  environment

  config_vertical_mixing:background_diffusivity :

  1.2e-05

  config_vertical_mixing:TSprofiles :

  False

  config_drift:wind_drift_depth :

  0.1

  config_drift:stokes_drift :

  True

  config_drift:stokes_drift_profile :

  Phillips

  config_drift:use_tabularised_stokes_drift :

  False

  config_drift:tabularised_stokes_drift_fetch :

  25000

  config_general:seafloor_action :

  lift_to_seafloor

  config_drift:truncate_ocean_model_below_m :

  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 :

  False

  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 :

  MARTIN LINGE CRUDE 2016

  opendrift_version :

  1.14.11

  seed_lon :

  4.88

  seed_lat :

  60.1

  seed_radius :

  3000.0

  seed_time :

  2015-11-16 00:00:00

  seed_number :

  500.0

  seed_z :

  0.0

  seed_m3_per_hour :

  1.0

  seed_oiltype :

  MARTIN LINGE CRUDE 2016

  seed_geojson :

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

  simulation_time :

  2026-10-26 17:03:28.897943

  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_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_horizontal_diffusivity :

  0

  reader_ocean_vertical_diffusivity :

  0.02

  reader_land_binary_mask :

  ['global_landmask']

  reader_ocean_mixed_layer_thickness :

  50

  time_coverage_end :

  2015-11-17 00:00:00

  time_coverage_duration :

  P1DT0H0M0S

  time_coverage_resolution :

  P0DT0H15M0S

  performance :

  -------------------- Reader performance: -------------------- /root/project/tests/test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc 0:00:00.2 total 0:00:00.0 preparing 0:00:00.0 reading 0:00:00.0 interpolation 0:00:00.0 interpolation_time 0:00:00.1 rotating vectors 0:00:00.0 masking -------------------- /root/project/tests/test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc 0:00:00.2 total 0:00:00.0 preparing 0:00:00.1 reading 0:00:00.0 interpolation 0:00:00.0 interpolation_time 0:00:00.1 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.8 total time 6.3 configuration 0.0 preparing main loop 0.0 moving elements to ocean 6.3 main loop 2.9 updating elements 0.1 oil weathering 0.0 updating viscosities 0.0 updating densities 0.0 evaporation 0.0 emulsification 2.7 vertical mixing 0.0 cleaning up --------------------

  geospatial_bounds_crs :

  EPSG:4326

  geospatial_bounds_vertical_crs :

  EPSG:5831

  geospatial_lat_min :

  60.017963

  geospatial_lat_max :

  60.857864

  geospatial_lat_units :

  degrees_north

  geospatial_lat_resolution :

  point

  geospatial_lon_min :

  4.4456744

  geospatial_lon_max :

  5.039443

  geospatial_lon_units :

  degrees_east

  geospatial_lon_resolution :

  point

  runtime :

  0:00:12.801694

  geospatial_vertical_min :

  -49.75206

  geospatial_vertical_max :

  0.0

  geospatial_vertical_positive :

  up

Print and plot results

o.plot_oil_budget()
o.plot(fast=True)
o.animation(fast=True)
#o.plot_oil_budget('oil_budget_MartinLingeCrude.png')

• 
• 

17:03:35 WARNING opendrift:2614: Plotting fast. This will make your plots less accurate.
17:03:36 WARNING opendrift:2614: Plotting fast. This will make your plots less accurate.
17:03:37 INFO    opendrift:4889: Saving animation to /root/project/docs/source/gallery/animations/example_oil_budget_noaa_0.gif...
17:03:45 INFO    opendrift:3297: Time to make animation: 0:00:08.405387

o.plot_property('fraction_evaporated')
o.plot_property('density')
o.plot_property('water_fraction')
o.plot_property('viscosity')
o.plot_property('interfacial_area')
o.plot_property('z')

• 
• 
• 
• 
• 
•