Openoil sample output netCDF file

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

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

rc = reader_constant.Reader({'x_wind': 5,
                             'y_wind': 4})

# Arome
reader_arome = reader_netCDF_CF_generic.Reader(o.test_data_folder() +
    '16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc')

# Norkyst
reader_norkyst = reader_netCDF_CF_generic.Reader(o.test_data_folder() +
    '16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc')

o.add_reader([rc, reader_norkyst, reader_arome])

15:28:09 INFO    opendrift:509: OpenDriftSimulation initialised (version 1.13.1 / v1.13.1-82-gbc545aa)
15:28:09 INFO    opendrift.readers:61: Opening file with xr.open_dataset
15:28:09 INFO    opendrift.readers.reader_netCDF_CF_generic:332: Detected dimensions: {'time': 'time', 'x': 'x', 'y': 'y'}
15:28:09 INFO    opendrift.readers:61: Opening file with xr.open_dataset
15:28:09 INFO    opendrift.readers.reader_netCDF_CF_generic:332: Detected dimensions: {'x': 'X', 'y': 'Y', 'z': 'depth', 'time': 'time'}

Seed oil elements at defined position and time

time = [reader_arome.start_time,
        reader_arome.start_time + timedelta(hours=1)]
o.seed_elements(lon=5.05, lat=59.95, radius=1500, number=100,
                time=time, z=0, m3_per_hour=2, oil_type='TAU 1999')

15:28:09 INFO    opendrift.models.openoil.adios.dirjs:86: Querying ADIOS database for oil: TAU 1999
15:28:09 INFO    opendrift.models.openoil.openoil:1710: Using density 851.002835 and viscosity 0.0012492557043601948 of oiltype TAU 1999
15:28:09 INFO    opendrift.models.basemodel.environment:206: Adding a global landmask from GSHHG
15:28:13 INFO    opendrift.models.basemodel.environment:229: Fallback values will be used for the following variables which have no readers:
15:28:13 INFO    opendrift.models.basemodel.environment:232:    sea_surface_height: 0.000000
15:28:13 INFO    opendrift.models.basemodel.environment:232:    upward_sea_water_velocity: 0.000000
15:28:13 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_significant_height: 0.000000
15:28:13 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
15:28:13 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
15:28:13 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
15:28:13 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
15:28:13 INFO    opendrift.models.basemodel.environment:232:    sea_ice_area_fraction: 0.000000
15:28:13 INFO    opendrift.models.basemodel.environment:232:    sea_ice_x_velocity: 0.000000
15:28:13 INFO    opendrift.models.basemodel.environment:232:    sea_ice_y_velocity: 0.000000
15:28:13 INFO    opendrift.models.basemodel.environment:232:    sea_water_temperature: 10.000000
15:28:13 INFO    opendrift.models.basemodel.environment:232:    sea_water_salinity: 34.000000
15:28:13 INFO    opendrift.models.basemodel.environment:232:    sea_floor_depth_below_sea_level: 10000.000000
15:28:13 INFO    opendrift.models.basemodel.environment:232:    ocean_vertical_diffusivity: 0.020000
15:28:13 INFO    opendrift.models.basemodel.environment:232:    ocean_mixed_layer_thickness: 50.000000

Running model

ncfile = 'openoil_sample_output.nc'
o.run(steps=4*8, time_step=900, time_step_output=3600, outfile=ncfile)

15:28:13 INFO    opendrift:907: Using existing reader for land_binary_mask
15:28:14 INFO    opendrift:938: Moving 1 out of 100 points from land to water
15:28:14 INFO    opendrift.models.openoil.openoil:685: Oil-water surface tension is 0.030101 Nm
15:28:14 INFO    opendrift.models.openoil.openoil:694: Using max water fractions [0.5, 0.6] for temperatures [5.0, 15.0] for oiltype TAU 1999
15:28:14 INFO    opendrift.models.openoil.openoil:695: Corresponding max water fraction from GNOME is 0.6321836371465919
15:28:14 INFO    opendrift:2069: 2015-11-16 00:00:00 - step 1 of 32 - 25 active elements (0 deactivated)
15:28:15 INFO    opendrift:2069: 2015-11-16 00:15:00 - step 2 of 32 - 50 active elements (0 deactivated)
15:28:15 INFO    opendrift:2069: 2015-11-16 00:30:00 - step 3 of 32 - 75 active elements (0 deactivated)
15:28:15 INFO    opendrift:2069: 2015-11-16 00:45:00 - step 4 of 32 - 100 active elements (0 deactivated)
15:28:15 INFO    opendrift:2069: 2015-11-16 01:00:00 - step 5 of 32 - 100 active elements (0 deactivated)
15:28:16 INFO    opendrift:2069: 2015-11-16 01:15:00 - step 6 of 32 - 99 active elements (1 deactivated)
15:28:16 INFO    opendrift:2069: 2015-11-16 01:30:00 - step 7 of 32 - 99 active elements (1 deactivated)
15:28:16 INFO    opendrift:2069: 2015-11-16 01:45:00 - step 8 of 32 - 95 active elements (5 deactivated)
15:28:16 INFO    opendrift:2069: 2015-11-16 02:00:00 - step 9 of 32 - 92 active elements (8 deactivated)
15:28:16 INFO    opendrift:2069: 2015-11-16 02:15:00 - step 10 of 32 - 90 active elements (10 deactivated)
15:28:16 INFO    opendrift:2069: 2015-11-16 02:30:00 - step 11 of 32 - 87 active elements (13 deactivated)
15:28:16 INFO    opendrift:2069: 2015-11-16 02:45:00 - step 12 of 32 - 80 active elements (20 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 03:00:00 - step 13 of 32 - 75 active elements (25 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 03:15:00 - step 14 of 32 - 72 active elements (28 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 03:30:00 - step 15 of 32 - 63 active elements (37 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 03:45:00 - step 16 of 32 - 60 active elements (40 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 04:00:00 - step 17 of 32 - 53 active elements (47 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 04:15:00 - step 18 of 32 - 46 active elements (54 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 04:30:00 - step 19 of 32 - 38 active elements (62 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 04:45:00 - step 20 of 32 - 35 active elements (65 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 05:00:00 - step 21 of 32 - 32 active elements (68 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 05:15:00 - step 22 of 32 - 25 active elements (75 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 05:30:00 - step 23 of 32 - 21 active elements (79 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 05:45:00 - step 24 of 32 - 17 active elements (83 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 06:00:00 - step 25 of 32 - 16 active elements (84 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 06:15:00 - step 26 of 32 - 15 active elements (85 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 06:30:00 - step 27 of 32 - 13 active elements (87 deactivated)
15:28:17 INFO    opendrift:2069: 2015-11-16 06:45:00 - step 28 of 32 - 10 active elements (90 deactivated)
15:28:18 INFO    opendrift:2069: 2015-11-16 07:00:00 - step 29 of 32 - 7 active elements (93 deactivated)
15:28:18 INFO    opendrift:2069: 2015-11-16 07:15:00 - step 30 of 32 - 6 active elements (94 deactivated)
15:28:18 INFO    opendrift:2069: 2015-11-16 07:30:00 - step 31 of 32 - 3 active elements (97 deactivated)
15:28:18 INFO    opendrift:2069: 2015-11-16 07:45:00 - step 32 of 32 - 1 active elements (99 deactivated)

<xarray.Dataset> Size: 170kB
Dimensions:                                                                              (
                                                                                          trajectory: 100,
                                                                                          time: 9)
Coordinates:
  * trajectory                                                                           (trajectory) int32 400B ...
  * time                                                                                 (time) datetime64[ns] 72B ...
Data variables: (12/44)
    status                                                                               (trajectory, time) float64 7kB ...
    moving                                                                               (trajectory, time) float64 7kB ...
    age_seconds                                                                          (trajectory, time) float32 4kB ...
    origin_marker                                                                        (trajectory, time) float64 7kB ...
    lon                                                                                  (trajectory, time) float32 4kB ...
    lat                                                                                  (trajectory, time) float32 4kB ...
    ...                                                                                   ...
    sea_water_temperature                                                                (trajectory, time) float32 4kB ...
    sea_water_salinity                                                                   (trajectory, time) float32 4kB ...
    sea_floor_depth_below_sea_level                                                      (trajectory, time) float32 4kB ...
    ocean_vertical_diffusivity                                                           (trajectory, time) float32 4kB ...
    land_binary_mask                                                                     (trajectory, time) float32 4kB ...
    ocean_mixed_layer_thickness                                                          (trajectory, time) float32 4kB ...
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: 100
  • time: 9
• Coordinates: (2)
  • trajectory
    (trajectory)
    int32
    0 1 2 3 4 5 6 ... 94 95 96 97 98 99

    cf_role :

    trajectory_id

    array([ 0,  1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12, 13, 14, 15, 16, 17,
           18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,
           36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53,
           54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71,
           72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89,
           90, 91, 92, 93, 94, 95, 96, 97, 98, 99], dtype=int32)

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

    standard_name :

    time

    long_name :

    time

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

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

    valid_range :

    [0 1]

    flag_values :

    [0 1]

    flag_meanings :

    active stranded

    [900 values with dtype=float64]

  • moving
    (trajectory, time)
    float64
    ...

    minval :

    0

    maxval :

    1

    [900 values with dtype=float64]

  • age_seconds
    (trajectory, time)
    float32
    ...

    units :

    s

    minval :

    0.0

    maxval :

    27000.0

    [900 values with dtype=float32]

  • origin_marker
    (trajectory, time)
    float64
    ...

    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

    [900 values with dtype=float64]

  • lon
    (trajectory, time)
    float32
    ...

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    axis :

    X

    minval :

    4.9814644

    maxval :

    5.143353

    [900 values with dtype=float32]

  • lat
    (trajectory, time)
    float32
    ...

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    axis :

    Y

    minval :

    59.921787

    maxval :

    60.01548

    [900 values with dtype=float32]

  • z
    (trajectory, time)
    float32
    ...

    units :

    m

    standard_name :

    z

    long_name :

    vertical position

    axis :

    Z

    positive :

    up

    minval :

    0.0

    maxval :

    0.0

    [900 values with dtype=float32]

  • wind_drift_factor
    (trajectory, time)
    float32
    ...

    units :

    %

    description :

    Elements at the ocean surface are moved by this fraction of the wind vector, in addition to currents and Stokes drift

    minval :

    0.03

    maxval :

    0.03

    [900 values with dtype=float32]

  • current_drift_factor
    (trajectory, time)
    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

    [900 values with dtype=float32]

  • terminal_velocity
    (trajectory, time)
    float32
    ...

    units :

    m/s

    description :

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

    minval :

    0.0

    maxval :

    0.10783295

    [900 values with dtype=float32]

  • mass_oil
    (trajectory, time)
    float32
    ...

    units :

    kg

    minval :

    8.155731

    maxval :

    17.020056

    [900 values with dtype=float32]

  • viscosity
    (trajectory, time)
    float32
    ...

    units :

    m2/s

    description :

    Kinematic viscosity of oil (emulsion)

    minval :

    0.0012492557

    maxval :

    1.7628639

    [900 values with dtype=float32]

  • density
    (trajectory, time)
    float32
    ...

    units :

    kg/m^3

    minval :

    851.0028

    maxval :

    947.7759

    [900 values with dtype=float32]

  • bulltime
    (trajectory, time)
    float32
    ...

    units :

    s

    minval :

    0.0

    maxval :

    0.0

    [900 values with dtype=float32]

  • interfacial_area
    (trajectory, time)
    float32
    ...

    units :

    m2

    minval :

    0.0

    maxval :

    1141499.9

    [900 values with dtype=float32]

  • mass_dispersed
    (trajectory, time)
    float32
    ...

    units :

    kg

    minval :

    0.0

    maxval :

    0.3321804

    [900 values with dtype=float32]

  • mass_evaporated
    (trajectory, time)
    float32
    ...

    units :

    kg

    minval :

    0.0

    maxval :

    8.568908

    [900 values with dtype=float32]

  • mass_biodegraded
    (trajectory, time)
    float32
    ...

    units :

    kg

    minval :

    0.0

    maxval :

    0.0

    [900 values with dtype=float32]

  • biodegradation_half_time_droplet
    (trajectory, time)
    float32
    ...

    units :

    Days

    description :

    Biodegradation half time in days for submerged oil droplets

    minval :

    1.0

    maxval :

    1.0

    [900 values with dtype=float32]

  • biodegradation_half_time_slick
    (trajectory, time)
    float32
    ...

    units :

    Days

    description :

    Biodegradation half time in days for surface oil slick

    minval :

    3.0

    maxval :

    3.0

    [900 values with dtype=float32]

  • fraction_evaporated
    (trajectory, time)
    float32
    ...

    units :

    %

    minval :

    0.0

    maxval :

    0.5109257

    [900 values with dtype=float32]

  • water_fraction
    (trajectory, time)
    float32
    ...

    units :

    %

    minval :

    0.0

    maxval :

    0.55

    [900 values with dtype=float32]

  • oil_film_thickness
    (trajectory, time)
    float32
    ...

    units :

    m

    minval :

    0.001

    maxval :

    0.001

    [900 values with dtype=float32]

  • diameter
    (trajectory, time)
    float32
    ...

    units :

    m

    minval :

    0.0

    maxval :

    0.0026039516

    [900 values with dtype=float32]

  • x_sea_water_velocity
    (trajectory, time)
    float32
    ...

    minval :

    -0.41676733

    maxval :

    0.19736072

    [900 values with dtype=float32]

  • y_sea_water_velocity
    (trajectory, time)
    float32
    ...

    minval :

    -0.1364768

    maxval :

    0.6068177

    [900 values with dtype=float32]

  • x_wind
    (trajectory, time)
    float32
    ...

    minval :

    3.6137037

    maxval :

    6.4470487

    [900 values with dtype=float32]

  • y_wind
    (trajectory, time)
    float32
    ...

    minval :

    2.5616503

    maxval :

    5.5025315

    [900 values with dtype=float32]

  • sea_surface_height
    (trajectory, time)
    float32
    ...

    minval :

    0.0

    maxval :

    0.0

    [900 values with dtype=float32]

  • upward_sea_water_velocity
    (trajectory, time)
    float32
    ...

    minval :

    0.0

    maxval :

    0.0

    [900 values with dtype=float32]

  • sea_surface_wave_significant_height
    (trajectory, time)
    float32
    ...

    minval :

    0.54447395

    maxval :

    1.5387571

    [900 values with dtype=float32]

  • sea_surface_wave_stokes_drift_x_velocity
    (trajectory, time)
    float32
    ...

    minval :

    0.0

    maxval :

    0.0

    [900 values with dtype=float32]

  • sea_surface_wave_stokes_drift_y_velocity
    (trajectory, time)
    float32
    ...

    minval :

    0.0

    maxval :

    0.0

    [900 values with dtype=float32]

  • sea_surface_wave_period_at_variance_spectral_density_maximum
    (trajectory, time)
    float32
    ...

    minval :

    0.0

    maxval :

    0.0

    [900 values with dtype=float32]

  • sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment
    (trajectory, time)
    float32
    ...

    minval :

    4.0199265

    maxval :

    6.7579446

    [900 values with dtype=float32]

  • sea_ice_area_fraction
    (trajectory, time)
    float32
    ...

    minval :

    0.0

    maxval :

    0.0

    [900 values with dtype=float32]

  • sea_ice_x_velocity
    (trajectory, time)
    float32
    ...

    minval :

    0.0

    maxval :

    0.0

    [900 values with dtype=float32]

  • sea_ice_y_velocity
    (trajectory, time)
    float32
    ...

    minval :

    0.0

    maxval :

    0.0

    [900 values with dtype=float32]

  • sea_water_temperature
    (trajectory, time)
    float32
    ...

    minval :

    10.0

    maxval :

    283.15

    [900 values with dtype=float32]

  • sea_water_salinity
    (trajectory, time)
    float32
    ...

    minval :

    34.0

    maxval :

    34.0

    [900 values with dtype=float32]

  • sea_floor_depth_below_sea_level
    (trajectory, time)
    float32
    ...

    minval :

    10000.0

    maxval :

    10000.0

    [900 values with dtype=float32]

  • ocean_vertical_diffusivity
    (trajectory, time)
    float32
    ...

    minval :

    0.02

    maxval :

    0.02

    [900 values with dtype=float32]

  • land_binary_mask
    (trajectory, time)
    float32
    ...

    minval :

    0.0

    maxval :

    0.0

    [900 values with dtype=float32]

  • ocean_mixed_layer_thickness
    (trajectory, time)
    float32
    ...

    minval :

    50.0

    maxval :

    50.0

    [900 values with dtype=float32]

• Indexes: (2)
  • trajectory
    PandasIndex

    PandasIndex(Index([ 0,  1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12, 13, 14, 15, 16, 17,
           18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,
           36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53,
           54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71,
           72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89,
           90, 91, 92, 93, 94, 95, 96, 97, 98, 99],
          dtype='int32', name='trajectory'))

  • 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'],
                  dtype='datetime64[ns]', name='time', freq=None))

• 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-05-09 15:28:13.532386

  date_created :

  2025-05-09T15:28:13.532392

  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(['constant_reader', '/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: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: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 :

  TAU 1999

  opendrift_version :

  1.13.1

  seed_lon :

  5.05

  seed_lat :

  59.95

  seed_radius :

  1500.0

  seed_time :

  2015-11-16 00:00:00

  seed_number :

  100.0

  seed_z :

  0.0

  seed_m3_per_hour :

  1.0

  seed_oiltype :

  TAU 1999

  seed_geojson :

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

  simulation_time :

  2025-05-09 15:28:14.957762

  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 :

  5

  reader_y_wind :

  4

  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-16 08:00:00

  time_coverage_duration :

  P0DT8H0M0S

  time_coverage_resolution :

  P0DT0H15M0S

  performance :

  -------------------- Reader performance: -------------------- constant_reader 0:00:00.0 total 0:00:00.0 preparing 0:00:00.0 reading 0:00:00.0 masking -------------------- /root/project/tests/test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc 0:00:00.0 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.0 rotating vectors 0:00:00.0 masking -------------------- /root/project/tests/test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc -------------------- global_landmask 0:00:01.4 total 0:00:00.0 preparing 0:00:01.4 reading 0:00:00.0 masking -------------------- Performance: 9.1 total time 4.4 configuration 1.4 preparing main loop 1.4 moving elements to ocean 3.2 main loop 0.9 updating elements 0.0 oil weathering 0.0 updating viscosities 0.0 updating densities 0.0 evaporation 0.0 emulsification 0.0 dispersion 0.8 vertical mixing 0.0 cleaning up --------------------

  geospatial_bounds_crs :

  EPSG:4326

  geospatial_bounds_vertical_crs :

  EPSG:5831

  geospatial_lat_min :

  59.921787

  geospatial_lat_max :

  60.01548

  geospatial_lat_units :

  degrees_north

  geospatial_lat_resolution :

  point

  geospatial_lon_min :

  4.9814644

  geospatial_lon_max :

  5.143353

  geospatial_lon_units :

  degrees_east

  geospatial_lon_resolution :

  point

  runtime :

  0:00:09.195359

  geospatial_vertical_min :

  0.0

  geospatial_vertical_max :

  0.0

  geospatial_vertical_positive :

  up

Print and plot results

o.plot_oil_budget()
o.plot()

• 
• 

/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)
/opt/conda/envs/opendrift/lib/python3.13/site-packages/shapely/creation.py:218: RuntimeWarning: invalid value encountered in linestrings
  return lib.linestrings(coords, np.intc(handle_nan), out=out, **kwargs)

(<GeoAxes: title={'center': 'OpenDrift - OpenOil (TAU 1999)\n2015-11-16 00:00 to 2015-11-16 08:00 UTC (9 steps)'}>, <Figure size 1072.61x1100 with 1 Axes>)

ncdump of the output file

from subprocess import check_output
args = [ "ncdump", "-h", ncfile ]
ncdump = check_output(args).decode().strip()
print(ncdump)

# Cleaning up
os.remove(ncfile)

netcdf openoil_sample_output {
dimensions:
        trajectory = 100 ;
        time = 9 ;
variables:
        int status(trajectory, time) ;
                status:_FillValue = 2147483647 ;
                status:valid_range = 0, 1 ;
                status:flag_values = 0, 1 ;
                status:flag_meanings = "active stranded" ;
        int moving(trajectory, time) ;
                moving:_FillValue = 2147483647 ;
                moving:minval = 0 ;
                moving:maxval = 1 ;
        float age_seconds(trajectory, time) ;
                age_seconds:_FillValue = NaNf ;
                age_seconds:units = "s" ;
                age_seconds:minval = 0.f ;
                age_seconds:maxval = 27000.f ;
        int origin_marker(trajectory, time) ;
                origin_marker:_FillValue = 2147483647 ;
                origin_marker:unit = "" ;
                origin_marker: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." ;
                origin_marker:flag_values = 0 ;
                origin_marker:flag_meanings = "Seed_0" ;
                origin_marker:minval = 0 ;
                origin_marker:maxval = 0 ;
        float lon(trajectory, time) ;
                lon:_FillValue = NaNf ;
                lon:units = "degrees_east" ;
                lon:standard_name = "longitude" ;
                lon:long_name = "longitude" ;
                lon:axis = "X" ;
                lon:minval = 4.981464f ;
                lon:maxval = 5.143353f ;
        float lat(trajectory, time) ;
                lat:_FillValue = NaNf ;
                lat:units = "degrees_north" ;
                lat:standard_name = "latitude" ;
                lat:long_name = "latitude" ;
                lat:axis = "Y" ;
                lat:minval = 59.92179f ;
                lat:maxval = 60.01548f ;
        float z(trajectory, time) ;
                z:_FillValue = NaNf ;
                z:units = "m" ;
                z:standard_name = "z" ;
                z:long_name = "vertical position" ;
                z:axis = "Z" ;
                z:positive = "up" ;
                z:minval = 0.f ;
                z:maxval = 0.f ;
        float wind_drift_factor(trajectory, time) ;
                wind_drift_factor:_FillValue = NaNf ;
                wind_drift_factor:units = "%" ;
                wind_drift_factor:description = "Elements at the ocean surface are moved by this fraction of the wind vector, in addition to currents and Stokes drift" ;
                wind_drift_factor:minval = 0.03f ;
                wind_drift_factor:maxval = 0.03f ;
        float current_drift_factor(trajectory, time) ;
                current_drift_factor:_FillValue = NaNf ;
                current_drift_factor:units = "1" ;
                current_drift_factor:description = "Elements are moved with this fraction of the current vector, in addition to currents and Stokes drift" ;
                current_drift_factor:minval = 1.f ;
                current_drift_factor:maxval = 1.f ;
        float terminal_velocity(trajectory, time) ;
                terminal_velocity:_FillValue = NaNf ;
                terminal_velocity:units = "m/s" ;
                terminal_velocity:description = "Terminal rise/sinking velocity (buoyancy) in the ocean column" ;
                terminal_velocity:minval = 0.f ;
                terminal_velocity:maxval = 0.107833f ;
        float mass_oil(trajectory, time) ;
                mass_oil:_FillValue = NaNf ;
                mass_oil:units = "kg" ;
                mass_oil:minval = 8.155731f ;
                mass_oil:maxval = 17.02006f ;
        float viscosity(trajectory, time) ;
                viscosity:_FillValue = NaNf ;
                viscosity:units = "m2/s" ;
                viscosity:description = "Kinematic viscosity of oil (emulsion)" ;
                viscosity:minval = 0.001249256f ;
                viscosity:maxval = 1.762864f ;
        float density(trajectory, time) ;
                density:_FillValue = NaNf ;
                density:units = "kg/m^3" ;
                density:minval = 851.0028f ;
                density:maxval = 947.7759f ;
        float bulltime(trajectory, time) ;
                bulltime:_FillValue = NaNf ;
                bulltime:units = "s" ;
                bulltime:minval = 0.f ;
                bulltime:maxval = 0.f ;
        float interfacial_area(trajectory, time) ;
                interfacial_area:_FillValue = NaNf ;
                interfacial_area:units = "m2" ;
                interfacial_area:minval = 0.f ;
                interfacial_area:maxval = 1141500.f ;
        float mass_dispersed(trajectory, time) ;
                mass_dispersed:_FillValue = NaNf ;
                mass_dispersed:units = "kg" ;
                mass_dispersed:minval = 0.f ;
                mass_dispersed:maxval = 0.3321804f ;
        float mass_evaporated(trajectory, time) ;
                mass_evaporated:_FillValue = NaNf ;
                mass_evaporated:units = "kg" ;
                mass_evaporated:minval = 0.f ;
                mass_evaporated:maxval = 8.568908f ;
        float mass_biodegraded(trajectory, time) ;
                mass_biodegraded:_FillValue = NaNf ;
                mass_biodegraded:units = "kg" ;
                mass_biodegraded:minval = 0.f ;
                mass_biodegraded:maxval = 0.f ;
        float biodegradation_half_time_droplet(trajectory, time) ;
                biodegradation_half_time_droplet:_FillValue = NaNf ;
                biodegradation_half_time_droplet:units = "Days" ;
                biodegradation_half_time_droplet:description = "Biodegradation half time in days for submerged oil droplets" ;
                biodegradation_half_time_droplet:minval = 1.f ;
                biodegradation_half_time_droplet:maxval = 1.f ;
        float biodegradation_half_time_slick(trajectory, time) ;
                biodegradation_half_time_slick:_FillValue = NaNf ;
                biodegradation_half_time_slick:units = "Days" ;
                biodegradation_half_time_slick:description = "Biodegradation half time in days for surface oil slick" ;
                biodegradation_half_time_slick:minval = 3.f ;
                biodegradation_half_time_slick:maxval = 3.f ;
        float fraction_evaporated(trajectory, time) ;
                fraction_evaporated:_FillValue = NaNf ;
                fraction_evaporated:units = "%" ;
                fraction_evaporated:minval = 0.f ;
                fraction_evaporated:maxval = 0.5109257f ;
        float water_fraction(trajectory, time) ;
                water_fraction:_FillValue = NaNf ;
                water_fraction:units = "%" ;
                water_fraction:minval = 0.f ;
                water_fraction:maxval = 0.55f ;
        float oil_film_thickness(trajectory, time) ;
                oil_film_thickness:_FillValue = NaNf ;
                oil_film_thickness:units = "m" ;
                oil_film_thickness:minval = 0.001f ;
                oil_film_thickness:maxval = 0.001f ;
        float diameter(trajectory, time) ;
                diameter:_FillValue = NaNf ;
                diameter:units = "m" ;
                diameter:minval = 0.f ;
                diameter:maxval = 0.002603952f ;
        float x_sea_water_velocity(trajectory, time) ;
                x_sea_water_velocity:_FillValue = NaNf ;
                x_sea_water_velocity:minval = -0.4167673f ;
                x_sea_water_velocity:maxval = 0.1973607f ;
        float y_sea_water_velocity(trajectory, time) ;
                y_sea_water_velocity:_FillValue = NaNf ;
                y_sea_water_velocity:minval = -0.1364768f ;
                y_sea_water_velocity:maxval = 0.6068177f ;
        float x_wind(trajectory, time) ;
                x_wind:_FillValue = NaNf ;
                x_wind:minval = 3.613704f ;
                x_wind:maxval = 6.447049f ;
        float y_wind(trajectory, time) ;
                y_wind:_FillValue = NaNf ;
                y_wind:minval = 2.56165f ;
                y_wind:maxval = 5.502532f ;
        float sea_surface_height(trajectory, time) ;
                sea_surface_height:_FillValue = NaNf ;
                sea_surface_height:minval = 0.f ;
                sea_surface_height:maxval = 0.f ;
        float upward_sea_water_velocity(trajectory, time) ;
                upward_sea_water_velocity:_FillValue = NaNf ;
                upward_sea_water_velocity:minval = 0.f ;
                upward_sea_water_velocity:maxval = 0.f ;
        float sea_surface_wave_significant_height(trajectory, time) ;
                sea_surface_wave_significant_height:_FillValue = NaNf ;
                sea_surface_wave_significant_height:minval = 0.5444739f ;
                sea_surface_wave_significant_height:maxval = 1.538757f ;
        float sea_surface_wave_stokes_drift_x_velocity(trajectory, time) ;
                sea_surface_wave_stokes_drift_x_velocity:_FillValue = NaNf ;
                sea_surface_wave_stokes_drift_x_velocity:minval = 0.f ;
                sea_surface_wave_stokes_drift_x_velocity:maxval = 0.f ;
        float sea_surface_wave_stokes_drift_y_velocity(trajectory, time) ;
                sea_surface_wave_stokes_drift_y_velocity:_FillValue = NaNf ;
                sea_surface_wave_stokes_drift_y_velocity:minval = 0.f ;
                sea_surface_wave_stokes_drift_y_velocity:maxval = 0.f ;
        float sea_surface_wave_period_at_variance_spectral_density_maximum(trajectory, time) ;
                sea_surface_wave_period_at_variance_spectral_density_maximum:_FillValue = NaNf ;
                sea_surface_wave_period_at_variance_spectral_density_maximum:minval = 0.f ;
                sea_surface_wave_period_at_variance_spectral_density_maximum:maxval = 0.f ;
        float sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment(trajectory, time) ;
                sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment:_FillValue = NaNf ;
                sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment:minval = 4.019927f ;
                sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment:maxval = 6.757945f ;
        float sea_ice_area_fraction(trajectory, time) ;
                sea_ice_area_fraction:_FillValue = NaNf ;
                sea_ice_area_fraction:minval = 0.f ;
                sea_ice_area_fraction:maxval = 0.f ;
        float sea_ice_x_velocity(trajectory, time) ;
                sea_ice_x_velocity:_FillValue = NaNf ;
                sea_ice_x_velocity:minval = 0.f ;
                sea_ice_x_velocity:maxval = 0.f ;
        float sea_ice_y_velocity(trajectory, time) ;
                sea_ice_y_velocity:_FillValue = NaNf ;
                sea_ice_y_velocity:minval = 0.f ;
                sea_ice_y_velocity:maxval = 0.f ;
        float sea_water_temperature(trajectory, time) ;
                sea_water_temperature:_FillValue = NaNf ;
                sea_water_temperature:minval = 10.f ;
                sea_water_temperature:maxval = 283.15f ;
        float sea_water_salinity(trajectory, time) ;
                sea_water_salinity:_FillValue = NaNf ;
                sea_water_salinity:minval = 34.f ;
                sea_water_salinity:maxval = 34.f ;
        float sea_floor_depth_below_sea_level(trajectory, time) ;
                sea_floor_depth_below_sea_level:_FillValue = NaNf ;
                sea_floor_depth_below_sea_level:minval = 10000.f ;
                sea_floor_depth_below_sea_level:maxval = 10000.f ;
        float ocean_vertical_diffusivity(trajectory, time) ;
                ocean_vertical_diffusivity:_FillValue = NaNf ;
                ocean_vertical_diffusivity:minval = 0.02f ;
                ocean_vertical_diffusivity:maxval = 0.02f ;
        float land_binary_mask(trajectory, time) ;
                land_binary_mask:_FillValue = NaNf ;
                land_binary_mask:minval = 0.f ;
                land_binary_mask:maxval = 0.f ;
        float ocean_mixed_layer_thickness(trajectory, time) ;
                ocean_mixed_layer_thickness:_FillValue = NaNf ;
                ocean_mixed_layer_thickness:minval = 50.f ;
                ocean_mixed_layer_thickness:maxval = 50.f ;
        int trajectory(trajectory) ;
                trajectory:cf_role = "trajectory_id" ;
        double time(time) ;
                time:standard_name = "time" ;
                time:long_name = "time" ;
                time:units = "seconds since 1970-01-01" ;
                time:calendar = "proleptic_gregorian" ;

// global attributes:
                :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-05-09 15:28:13.532386" ;
                :date_created = "2025-05-09T15:28:13.532392" ;
                :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([\'constant_reader\', \'/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: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 = 0LL ;
                :config_environment\:constant\:upward_sea_water_velocity = "None" ;
                :config_environment\:fallback\:upward_sea_water_velocity = 0LL ;
                :config_environment\:constant\:sea_surface_wave_significant_height = "None" ;
                :config_environment\:fallback\:sea_surface_wave_significant_height = 0LL ;
                :config_environment\:constant\:sea_surface_wave_stokes_drift_x_velocity = "None" ;
                :config_environment\:fallback\:sea_surface_wave_stokes_drift_x_velocity = 0LL ;
                :config_environment\:constant\:sea_surface_wave_stokes_drift_y_velocity = "None" ;
                :config_environment\:fallback\:sea_surface_wave_stokes_drift_y_velocity = 0LL ;
                :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 = 0LL ;
                :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 = 0LL ;
                :config_environment\:constant\:sea_ice_area_fraction = "None" ;
                :config_environment\:fallback\:sea_ice_area_fraction = 0LL ;
                :config_environment\:constant\:sea_ice_x_velocity = "None" ;
                :config_environment\:fallback\:sea_ice_x_velocity = 0LL ;
                :config_environment\:constant\:sea_ice_y_velocity = "None" ;
                :config_environment\:fallback\:sea_ice_y_velocity = 0LL ;
                :config_environment\:constant\:sea_water_temperature = "None" ;
                :config_environment\:fallback\:sea_water_temperature = 10LL ;
                :config_environment\:constant\:sea_water_salinity = "None" ;
                :config_environment\:fallback\:sea_water_salinity = 34LL ;
                :config_environment\:constant\:sea_floor_depth_below_sea_level = "None" ;
                :config_environment\:fallback\:sea_floor_depth_below_sea_level = 10000LL ;
                :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 = 50LL ;
                :config_general\:use_auto_landmask = "True" ;
                :config_drift\:current_uncertainty = 0.05 ;
                :config_drift\:current_uncertainty_uniform = 0LL ;
                :config_drift\:max_speed = 1.3 ;
                :config_readers\:max_number_of_fails = 1LL ;
                :config_general\:simulation_name = "" ;
                :config_general\:coastline_action = "stranding" ;
                :config_general\:coastline_approximation_precision = 0.001 ;
                :config_general\:time_step_minutes = 60LL ;
                :config_general\:time_step_output_minutes = "None" ;
                :config_seed\:ocean_only = "True" ;
                :config_seed\:number = 1LL ;
                :config_drift\:max_age_seconds = "None" ;
                :config_drift\:advection_scheme = "euler" ;
                :config_drift\:horizontal_diffusivity = 0LL ;
                :config_drift\:profiles_depth = 50LL ;
                :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 = 0LL ;
                :config_seed\:z = 0LL ;
                :config_seed\:wind_drift_factor = 0.03 ;
                :config_seed\:current_drift_factor = 1LL ;
                :config_seed\:terminal_velocity = 0. ;
                :config_seed\:oil_film_thickness = 0.001 ;
                :config_drift\:vertical_advection = "False" ;
                :config_drift\:vertical_mixing = "True" ;
                :config_vertical_mixing\:timestep = 60LL ;
                :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 = 1LL ;
                :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 = "TAU 1999" ;
                :opendrift_version = "1.13.1" ;
                :seed_lon = "5.05" ;
                :seed_lat = "59.95" ;
                :seed_radius = "1500.0" ;
                :seed_time = "2015-11-16 00:00:00" ;
                :seed_number = "100.0" ;
                :seed_z = "0.0" ;
                :seed_m3_per_hour = "1.0" ;
                :seed_oiltype = "TAU 1999" ;
                :seed_geojson = "{\"features\": [], \"type\": \"FeatureCollection\"}" ;
                :simulation_time = "2025-05-09 15:28:14.957762" ;
                :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 = "5" ;
                :reader_y_wind = "4" ;
                :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-16 08:00:00" ;
                :time_coverage_duration = "P0DT8H0M0S" ;
                :time_coverage_resolution = "P0DT0H15M0S" ;
                :performance = "--------------------\nReader performance:\n--------------------\nconstant_reader\n 0:00:00.0  total\n 0:00:00.0  preparing\n 0:00:00.0  reading\n 0:00:00.0  masking\n--------------------\n/root/project/tests/test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc\n 0:00:00.0  total\n 0:00:00.0  preparing\n 0:00:00.0  reading\n 0:00:00.0  interpolation\n 0:00:00.0  interpolation_time\n 0:00:00.0  rotating vectors\n 0:00:00.0  masking\n--------------------\n/root/project/tests/test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc\n--------------------\nglobal_landmask\n 0:00:01.4  total\n 0:00:00.0  preparing\n 0:00:01.4  reading\n 0:00:00.0  masking\n--------------------\nPerformance:\n    9.1 total time\n    4.4 configuration\n    1.4 preparing main loop\n      1.4 moving elements to ocean\n    3.2 main loop\n      0.9 updating elements\n        0.0 oil weathering\n          0.0 updating viscosities\n          0.0 updating densities\n          0.0 evaporation\n          0.0 emulsification\n          0.0 dispersion\n        0.8 vertical mixing\n    0.0 cleaning up\n--------------------\n" ;
                :geospatial_bounds_crs = "EPSG:4326" ;
                :geospatial_bounds_vertical_crs = "EPSG:5831" ;
                :geospatial_lat_min = 59.92179f ;
                :geospatial_lat_max = 60.01548f ;
                :geospatial_lat_units = "degrees_north" ;
                :geospatial_lat_resolution = "point" ;
                :geospatial_lon_min = 4.981464f ;
                :geospatial_lon_max = 5.143353f ;
                :geospatial_lon_units = "degrees_east" ;
                :geospatial_lon_resolution = "point" ;
                :runtime = "0:00:09.195359" ;
                :geospatial_vertical_min = 0.f ;
                :geospatial_vertical_max = 0.f ;
                :geospatial_vertical_positive = "up" ;
}