Ship drift

from datetime import datetime
from opendrift.models.shipdrift import ShipDrift

o = ShipDrift(loglevel=20)

o.add_readers_from_list([
    'https://thredds.met.no/thredds/dodsC/cmems/topaz6/dataset-topaz6-arc-15min-3km-be.ncml',
    'https://thredds.met.no/thredds/dodsC/mepslatest/meps_lagged_6_h_latest_2_5km_latest.nc',
    'https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml'
    ])

21:56:07 INFO    opendrift:509: OpenDriftSimulation initialised (version 1.13.0 / v1.13.0-61-gf300067)

Seed ship elements at defined position and time Note: beam/length ratio is larger than allowed, but is then clipped internally

o.seed_elements(lon=5.0, lat=63.0, radius=1000, number=1000,
                time=datetime.utcnow(),
                length=80.0, beam=20.0, height=9.0, draft=4.0)

21:56:07 WARNING opendrift.models.shipdrift:181: Ratio of beam to length should be in range 0.12 to 0.18, given range is 0.25-0.25. Using border value.
21:56:07 INFO    opendrift.models.basemodel.environment:206: Adding a dynamical landmask with max. priority based on assumed maximum speed of 2.0 m/s. Adding a customised landmask may be faster...
21:56:11 INFO    opendrift.models.basemodel.environment:233: Fallback values will be used for the following variables which have no readers:
21:56:11 INFO    opendrift.models.basemodel.environment:236:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
21:56:11 INFO    opendrift.models.basemodel.environment:236:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
21:56:11 INFO    opendrift.models.basemodel.environment:236:    sea_surface_wave_significant_height: 0.000000
21:56:11 INFO    opendrift.models.basemodel.environment:236:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000

Running model

o.run(steps=24, stop_on_error=True)

21:56:11 INFO    opendrift:919: Using existing reader for land_binary_mask
21:56:11 INFO    opendrift:930: All points are in ocean
21:56:11 INFO    opendrift:2056: 2025-03-05 21:56:07.344141 - step 1 of 24 - 1000 active elements (0 deactivated)
21:56:11 INFO    opendrift.readers:58: Opening file with xr.open_dataset
21:56:14 INFO    opendrift.readers.reader_netCDF_CF_generic:338: Detected dimensions: {'x': 'x', 'y': 'y', 'time': 'time'}
21:56:14 INFO    opendrift.readers:58: Opening file with xr.open_dataset
21:56:15 INFO    opendrift.readers.reader_netCDF_CF_generic:338: Detected dimensions: {'time': 'time', 'x': 'x', 'y': 'y'}
21:56:15 INFO    opendrift.readers:58: Opening file with xr.open_dataset
21:56:16 INFO    opendrift.readers.reader_netCDF_CF_generic:338: Detected dimensions: {'y': 'rlat', 'x': 'rlon', 'time': 'time'}
21:56:18 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:18 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:19 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:19 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:21 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:21 INFO    opendrift:2056: 2025-03-05 22:56:07.344141 - step 2 of 24 - 1000 active elements (0 deactivated)
21:56:22 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:22 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:23 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:23 INFO    opendrift:2056: 2025-03-05 23:56:07.344141 - step 3 of 24 - 1000 active elements (0 deactivated)
21:56:26 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:26 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:27 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:27 INFO    opendrift:2056: 2025-03-06 00:56:07.344141 - step 4 of 24 - 1000 active elements (0 deactivated)
21:56:29 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:29 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:29 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:29 INFO    opendrift:2056: 2025-03-06 01:56:07.344141 - step 5 of 24 - 1000 active elements (0 deactivated)
21:56:31 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:31 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:32 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:32 INFO    opendrift:2056: 2025-03-06 02:56:07.344141 - step 6 of 24 - 1000 active elements (0 deactivated)
21:56:33 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:33 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:34 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:34 INFO    opendrift:2056: 2025-03-06 03:56:07.344141 - step 7 of 24 - 1000 active elements (0 deactivated)
21:56:36 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:36 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:37 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:37 INFO    opendrift:2056: 2025-03-06 04:56:07.344141 - step 8 of 24 - 1000 active elements (0 deactivated)
21:56:38 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:38 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:39 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:39 INFO    opendrift:2056: 2025-03-06 05:56:07.344141 - step 9 of 24 - 1000 active elements (0 deactivated)
21:56:40 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:40 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:41 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:41 INFO    opendrift:2056: 2025-03-06 06:56:07.344141 - step 10 of 24 - 1000 active elements (0 deactivated)
21:56:42 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:42 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:43 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:43 INFO    opendrift:2056: 2025-03-06 07:56:07.344141 - step 11 of 24 - 1000 active elements (0 deactivated)
21:56:44 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:44 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:45 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:45 INFO    opendrift:2056: 2025-03-06 08:56:07.344141 - step 12 of 24 - 1000 active elements (0 deactivated)
21:56:46 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:46 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:47 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:47 INFO    opendrift:2056: 2025-03-06 09:56:07.344141 - step 13 of 24 - 1000 active elements (0 deactivated)
21:56:49 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:49 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:50 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:50 INFO    opendrift:2056: 2025-03-06 10:56:07.344141 - step 14 of 24 - 1000 active elements (0 deactivated)
21:56:51 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:51 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:52 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:52 INFO    opendrift:2056: 2025-03-06 11:56:07.344141 - step 15 of 24 - 1000 active elements (0 deactivated)
21:56:54 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:54 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:55 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:55 INFO    opendrift:2056: 2025-03-06 12:56:07.344141 - step 16 of 24 - 1000 active elements (0 deactivated)
21:56:56 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:56 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:57 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:56:57 INFO    opendrift:2056: 2025-03-06 13:56:07.344141 - step 17 of 24 - 1000 active elements (0 deactivated)
21:56:58 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:56:58 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:57:00 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:57:00 INFO    opendrift:2056: 2025-03-06 14:56:07.344141 - step 18 of 24 - 1000 active elements (0 deactivated)
21:57:01 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:57:01 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:57:02 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:57:02 INFO    opendrift:2056: 2025-03-06 15:56:07.344141 - step 19 of 24 - 1000 active elements (0 deactivated)
21:57:03 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:57:03 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:57:04 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:57:04 INFO    opendrift:2056: 2025-03-06 16:56:07.344141 - step 20 of 24 - 1000 active elements (0 deactivated)
21:57:06 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:57:06 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:57:07 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:57:07 INFO    opendrift:2056: 2025-03-06 17:56:07.344141 - step 21 of 24 - 1000 active elements (0 deactivated)
21:57:08 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:57:08 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:57:09 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:57:09 INFO    opendrift:2056: 2025-03-06 18:56:07.344141 - step 22 of 24 - 1000 active elements (0 deactivated)
21:57:10 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:57:10 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:57:11 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:57:11 INFO    opendrift:2056: 2025-03-06 19:56:07.344141 - step 23 of 24 - 1000 active elements (0 deactivated)
21:57:12 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:57:12 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:57:13 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction
21:57:13 INFO    opendrift:2056: 2025-03-06 20:56:07.344141 - step 24 of 24 - 1000 active elements (0 deactivated)
21:57:15 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:57:15 WARNING opendrift.readers.interpolation.structured:63: Ensemble data currently not extrapolated towards seafloor
21:57:16 INFO    opendrift.models.shipdrift:311: Using Stokes drift direction as wave direction

<xarray.Dataset> Size: 2MB
Dimensions:                                                                              (
                                                                                          trajectory: 1000,
                                                                                          time: 25)
Coordinates:
  * trajectory                                                                           (trajectory) int64 8kB ...
  * time                                                                                 (time) datetime64[ns] 200B ...
Data variables: (12/24)
    status                                                                               (trajectory, time) float32 100kB ...
    moving                                                                               (trajectory, time) float32 100kB ...
    age_seconds                                                                          (trajectory, time) float32 100kB ...
    origin_marker                                                                        (trajectory, time) float32 100kB ...
    lon                                                                                  (trajectory, time) float32 100kB ...
    lat                                                                                  (trajectory, time) float32 100kB ...
    ...                                                                                   ...
    x_sea_water_velocity                                                                 (trajectory, time) float32 100kB ...
    y_sea_water_velocity                                                                 (trajectory, time) float32 100kB ...
    sea_surface_wave_stokes_drift_x_velocity                                             (trajectory, time) float32 100kB ...
    sea_surface_wave_stokes_drift_y_velocity                                             (trajectory, time) float32 100kB ...
    sea_surface_wave_significant_height                                                  (trajectory, time) float32 100kB ...
    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment  (trajectory, time) float32 100kB ...
Attributes: (12/96)
    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: 1000
  • time: 25
• Coordinates: (2)
  • trajectory
    (trajectory)
    int64
    0 1 2 3 4 5 ... 995 996 997 998 999

    cf_role :

    trajectory_id

    dtype :

    <class 'numpy.int32'>

    array([  0,   1,   2, ..., 997, 998, 999])

  • time
    (time)
    datetime64[ns]
    2025-03-05T21:56:07.344141 ... 2...

    standard_name :

    time

    long_name :

    time

    dtype :

    <class 'numpy.float64'>

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

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

    dtype :

    <class 'numpy.int32'>

    valid_range :

    [0 0]

    flag_values :

    [0]

    flag_meanings :

    active

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

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

    dtype :

    <class 'numpy.int32'>

    minval :

    1

    maxval :

    1

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

  • age_seconds
    (trajectory, time)
    float32
    0.0 3.6e+03 ... 8.28e+04 8.64e+04

    dtype :

    <class 'numpy.float32'>

    units :

    s

    minval :

    0.0

    maxval :

    86400.0

    array([[    0.,  3600.,  7200., ..., 79200., 82800., 86400.],
           [    0.,  3600.,  7200., ..., 79200., 82800., 86400.],
           [    0.,  3600.,  7200., ..., 79200., 82800., 86400.],
           ...,
           [    0.,  3600.,  7200., ..., 79200., 82800., 86400.],
           [    0.,  3600.,  7200., ..., 79200., 82800., 86400.],
           [    0.,  3600.,  7200., ..., 79200., 82800., 86400.]],
          dtype=float32)

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

    dtype :

    <class 'numpy.int32'>

    unit :

    description :

    An integer kept constant during the simulation. Different values may be used for different seedings, to separate elements during analysis. With GUI, only a single seeding is possible.

    flag_values :

    [0]

    flag_meanings :

    Seed_0

    minval :

    0

    maxval :

    0

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

  • lon
    (trajectory, time)
    float32
    5.035 5.106 5.204 ... 6.658 6.679

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    axis :

    X

    minval :

    4.9398756

    maxval :

    7.024521

    array([[5.0348186, 5.1055193, 5.20421  , ..., 6.5028   , 6.527642 ,
            6.54208  ],
           [5.007899 , 5.088261 , 5.2176037, ..., 6.6896753, 6.692416 ,
            6.6563616],
           [5.0193124, 5.1226597, 5.221775 , ..., 6.768137 , 6.790417 ,
            6.8115716],
           ...,
           [5.0018587, 5.113972 , 5.230477 , ..., 6.8485565, 6.8897295,
            6.9019556],
           [4.9773602, 5.055765 , 5.1563244, ..., 6.5270514, 6.56716  ,
            6.620329 ],
           [4.9929357, 5.097072 , 5.2047887, ..., 6.666166 , 6.6579285,
            6.678587 ]], dtype=float32)

  • lat
    (trajectory, time)
    float32
    63.0 63.04 63.08 ... 63.11 63.1

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    axis :

    Y

    minval :

    62.973133

    maxval :

    63.456875

    array([[63.004982, 63.042458, 63.076187, ..., 63.286846, 63.297382,
            63.29263 ],
           [63.008007, 63.033394, 63.055378, ..., 63.135845, 63.131954,
            63.118702],
           [62.99621 , 63.019745, 63.037422, ..., 63.22511 , 63.23228 ,
            63.23584 ],
           ...,
           [63.001423, 63.018272, 63.037785, ..., 63.17078 , 63.149956,
            63.14799 ],
           [62.989754, 63.02732 , 63.058994, ..., 63.340073, 63.334885,
            63.344246],
           [62.98824 , 63.00404 , 63.03024 , ..., 63.114162, 63.11162 ,
            63.100002]], dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    m

    standard_name :

    z

    long_name :

    vertical position

    axis :

    Z

    positive :

    up

    minval :

    0.0

    maxval :

    0.0

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

  • orientation
    (trajectory, time)
    float32
    0.0 0.0 0.0 0.0 ... 1.0 1.0 1.0 1.0

    dtype :

    <class 'numpy.uint8'>

    units :

    1

    minval :

    0

    maxval :

    1

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

  • length
    (trajectory, time)
    float32
    80.0 80.0 80.0 ... 80.0 80.0 80.0

    dtype :

    <class 'numpy.float32'>

    units :

    m

    min :

    1

    max :

    500

    description :

    Length of ship

    level :

    1

    minval :

    80.0

    maxval :

    80.0

    array([[80., 80., 80., ..., 80., 80., 80.],
           [80., 80., 80., ..., 80., 80., 80.],
           [80., 80., 80., ..., 80., 80., 80.],
           ...,
           [80., 80., 80., ..., 80., 80., 80.],
           [80., 80., 80., ..., 80., 80., 80.],
           [80., 80., 80., ..., 80., 80., 80.]], dtype=float32)

  • height
    (trajectory, time)
    float32
    9.0 9.0 9.0 9.0 ... 9.0 9.0 9.0 9.0

    dtype :

    <class 'numpy.float32'>

    units :

    m

    min :

    1

    max :

    100

    description :

    Total height of ship (above and below waterline)

    level :

    1

    minval :

    9.0

    maxval :

    9.0

    array([[9., 9., 9., ..., 9., 9., 9.],
           [9., 9., 9., ..., 9., 9., 9.],
           [9., 9., 9., ..., 9., 9., 9.],
           ...,
           [9., 9., 9., ..., 9., 9., 9.],
           [9., 9., 9., ..., 9., 9., 9.],
           [9., 9., 9., ..., 9., 9., 9.]], dtype=float32)

  • draft
    (trajectory, time)
    float32
    4.0 4.0 4.0 4.0 ... 4.0 4.0 4.0 4.0

    dtype :

    <class 'numpy.float32'>

    units :

    m

    min :

    1

    max :

    30

    description :

    Draft of ship (depth below water)

    level :

    1

    minval :

    4.0

    maxval :

    4.0

    array([[4., 4., 4., ..., 4., 4., 4.],
           [4., 4., 4., ..., 4., 4., 4.],
           [4., 4., 4., ..., 4., 4., 4.],
           ...,
           [4., 4., 4., ..., 4., 4., 4.],
           [4., 4., 4., ..., 4., 4., 4.],
           [4., 4., 4., ..., 4., 4., 4.]], dtype=float32)

  • beam
    (trajectory, time)
    float32
    20.0 20.0 20.0 ... 20.0 20.0 20.0

    dtype :

    <class 'numpy.float32'>

    min :

    1

    max :

    70

    units :

    m

    description :

    Beam (width) of ship

    level :

    1

    minval :

    20.0

    maxval :

    20.0

    array([[20., 20., 20., ..., 20., 20., 20.],
           [20., 20., 20., ..., 20., 20., 20.],
           [20., 20., 20., ..., 20., 20., 20.],
           ...,
           [20., 20., 20., ..., 20., 20., 20.],
           [20., 20., 20., ..., 20., 20., 20.],
           [20., 20., 20., ..., 20., 20., 20.]], dtype=float32)

  • wind_drag_coeff
    (trajectory, time)
    float32
    0.815 0.815 0.815 ... 0.815 0.815

    dtype :

    <class 'numpy.float32'>

    units :

    1

    minval :

    0.815

    maxval :

    0.815

    array([[0.815, 0.815, 0.815, ..., 0.815, 0.815, 0.815],
           [0.815, 0.815, 0.815, ..., 0.815, 0.815, 0.815],
           [0.815, 0.815, 0.815, ..., 0.815, 0.815, 0.815],
           ...,
           [0.815, 0.815, 0.815, ..., 0.815, 0.815, 0.815],
           [0.815, 0.815, 0.815, ..., 0.815, 0.815, 0.815],
           [0.815, 0.815, 0.815, ..., 0.815, 0.815, 0.815]], dtype=float32)

  • water_drag_coeff
    (trajectory, time)
    float32
    1.32 1.32 1.32 ... 1.32 1.32 1.32

    dtype :

    <class 'numpy.float32'>

    units :

    1

    minval :

    1.32

    maxval :

    1.32

    array([[1.32, 1.32, 1.32, ..., 1.32, 1.32, 1.32],
           [1.32, 1.32, 1.32, ..., 1.32, 1.32, 1.32],
           [1.32, 1.32, 1.32, ..., 1.32, 1.32, 1.32],
           ...,
           [1.32, 1.32, 1.32, ..., 1.32, 1.32, 1.32],
           [1.32, 1.32, 1.32, ..., 1.32, 1.32, 1.32],
           [1.32, 1.32, 1.32, ..., 1.32, 1.32, 1.32]], dtype=float32)

  • jibeProbability
    (trajectory, time)
    float32
    0.04 0.04 0.04 ... 0.04 0.04 0.04

    dtype :

    <class 'numpy.float32'>

    units :

    1/h

    minval :

    0.04

    maxval :

    0.04

    array([[0.04, 0.04, 0.04, ..., 0.04, 0.04, 0.04],
           [0.04, 0.04, 0.04, ..., 0.04, 0.04, 0.04],
           [0.04, 0.04, 0.04, ..., 0.04, 0.04, 0.04],
           ...,
           [0.04, 0.04, 0.04, ..., 0.04, 0.04, 0.04],
           [0.04, 0.04, 0.04, ..., 0.04, 0.04, 0.04],
           [0.04, 0.04, 0.04, ..., 0.04, 0.04, 0.04]], dtype=float32)

  • x_wind
    (trajectory, time)
    float32
    11.27 12.76 ... -0.3356 -0.3356

    minval :

    -5.515845

    maxval :

    21.608671

    array([[11.26649   , 12.755262  , 15.606082  , ..., -2.9460158 ,
            -0.9176974 , -0.9176974 ],
           [10.235095  , 11.768533  , 14.177288  , ..., -1.2511168 ,
            -0.45435467, -0.45435467],
           [14.037307  , 15.443307  , 15.768857  , ...,  1.5024804 ,
             1.0739654 ,  1.0739654 ],
           ...,
           [14.404433  , 13.189688  , 12.900795  , ..., -2.261806  ,
            -4.384785  , -4.384785  ],
           [13.214402  , 13.858331  , 15.556949  , ...,  2.955427  ,
             4.0727596 ,  4.0727596 ],
           [11.801314  , 11.758478  , 12.249408  , ..., -0.1957619 ,
            -0.33560616, -0.33560616]], dtype=float32)

  • y_wind
    (trajectory, time)
    float32
    14.27 13.75 15.49 ... 3.286 3.286

    minval :

    -2.7679462

    maxval :

    18.456747

    array([[14.270587 , 13.751572 , 15.48721  , ...,  2.9490972,  6.723205 ,
             6.723205 ],
           [15.209942 , 16.431871 , 17.038383 , ...,  2.6208057,  3.8412652,
             3.8412652],
           [14.03058  , 14.437188 , 16.507475 , ...,  4.0992794,  4.096858 ,
             4.096858 ],
           ...,
           [11.279642 , 13.149542 , 15.666657 , ...,  3.5364757,  2.0792866,
             2.0792866],
           [11.7544775, 13.515211 , 15.396202 , ...,  3.9081109,  4.446479 ,
             4.446479 ],
           [14.796185 , 15.431989 , 15.122896 , ...,  3.376833 ,  3.2858713,
             3.2858713]], dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

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

  • x_sea_water_velocity
    (trajectory, time)
    float32
    0.9062 0.7758 ... 0.04043 0.04043

    minval :

    -0.094790354

    maxval :

    1.1157072

    array([[0.9061539 , 0.77578676, 0.6478728 , ..., 0.27829927, 0.2952455 ,
            0.2952455 ],
           [0.96216726, 0.80344546, 0.7093493 , ..., 0.08362284, 0.09854357,
            0.09854357],
           [0.9268899 , 0.80017096, 0.7634205 , ..., 0.22249196, 0.23150867,
            0.23150867],
           ...,
           [0.97517794, 0.80621254, 0.75034004, ..., 0.2618705 , 0.27915207,
            0.27915207],
           [0.88906956, 0.8626184 , 0.75992954, ..., 0.13480474, 0.17101021,
            0.17101021],
           [0.9002457 , 0.8402197 , 0.7916693 , ..., 0.01799555, 0.04043007,
            0.04043007]], dtype=float32)

  • y_sea_water_velocity
    (trajectory, time)
    float32
    0.2186 0.3887 ... -0.05879 -0.05879

    minval :

    -0.39625064

    maxval :

    0.6322862

    array([[ 0.21857361,  0.3886664 ,  0.4881873 , ..., -0.04878   ,
            -0.02964123, -0.02964123],
           [ 0.21840107,  0.35094404,  0.44764873, ..., -0.10092552,
            -0.04677178, -0.04677178],
           [ 0.2092679 ,  0.33228442,  0.41195402, ...,  0.01087107,
             0.04693509,  0.04693509],
           ...,
           [ 0.20876144,  0.32372177,  0.4076674 , ..., -0.01874658,
             0.0082694 ,  0.0082694 ],
           [ 0.28875318,  0.3147557 ,  0.46824053, ..., -0.01643222,
             0.00759739,  0.00759739],
           [ 0.25490025,  0.270845  ,  0.39732778, ..., -0.11277741,
            -0.05878534, -0.05878534]], dtype=float32)

  • sea_surface_wave_stokes_drift_x_velocity
    (trajectory, time)
    float32
    0.3997 0.3944 ... 0.07573 0.07573

    minval :

    0.071776874

    maxval :

    0.47799468

    array([[0.39965108, 0.39444688, 0.3858561 , ..., 0.08754208, 0.0809886 ,
            0.0809886 ],
           [0.41238713, 0.4008502 , 0.39690506, ..., 0.08446582, 0.07643831,
            0.07643831],
           [0.40476054, 0.3883635 , 0.39784232, ..., 0.08890618, 0.08160526,
            0.08160526],
           ...,
           [0.41221154, 0.3895004 , 0.39891446, ..., 0.08709908, 0.07502812,
            0.07502812],
           [0.39489928, 0.40320444, 0.3889383 , ..., 0.08870544, 0.08160088,
            0.08160088],
           [0.39535415, 0.3880487 , 0.39454132, ..., 0.08403654, 0.0757299 ,
            0.0757299 ]], dtype=float32)

  • sea_surface_wave_stokes_drift_y_velocity
    (trajectory, time)
    float32
    0.2433 0.3034 ... -0.01989 -0.01989

    minval :

    -0.028467521

    maxval :

    0.38437542

    array([[ 2.4329168e-01,  3.0337444e-01,  3.8049132e-01, ...,
            -4.2348015e-03,  1.3689036e-04,  1.3689036e-04],
           [ 2.3597899e-01,  2.9953647e-01,  3.7602741e-01, ...,
            -2.1273935e-02, -1.9808084e-02, -1.9808084e-02],
           [ 2.3177755e-01,  2.9422623e-01,  3.7262398e-01, ...,
            -9.7126691e-03, -6.0040383e-03, -6.0040383e-03],
           ...,
           [ 2.2981362e-01,  2.9383296e-01,  3.7445980e-01, ...,
            -1.4666517e-02, -1.6440056e-02, -1.6440056e-02],
           [ 2.2889642e-01,  3.0091119e-01,  3.7058756e-01, ...,
             8.8897906e-04,  3.3322501e-03,  3.3322501e-03],
           [ 2.2957169e-01,  2.9196367e-01,  3.6657324e-01, ...,
            -2.2290871e-02, -1.9892445e-02, -1.9892445e-02]], dtype=float32)

  • sea_surface_wave_significant_height
    (trajectory, time)
    float32
    6.308 6.69 7.188 ... 3.361 3.361

    minval :

    3.2360716

    maxval :

    9.69218

    array([[6.3078666, 6.689787 , 7.187976 , ..., 3.533606 , 3.3298557,
            3.3298557],
           [6.3194246, 6.7229156, 7.1681747, ..., 3.5899444, 3.3630445,
            3.3630445],
           [6.3257775, 6.679742 , 7.1414833, ..., 3.5146956, 3.3214138,
            3.3214138],
           ...,
           [6.3265295, 6.685296 , 7.146053 , ..., 3.500289 , 3.3090627,
            3.3090627],
           [6.356126 , 6.727033 , 7.1268277, ..., 3.671039 , 3.4527273,
            3.4527273],
           [6.3513007, 6.682521 , 7.126814 , ..., 3.6049707, 3.361217 ,
            3.361217 ]], dtype=float32)

  • sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment
    (trajectory, time)
    float32
    7.221 7.326 7.428 ... 7.995 7.995

    minval :

    7.151821

    maxval :

    9.294474

    array([[7.221003 , 7.3264785, 7.4281898, ..., 7.846495 , 7.7007422,
            7.7007422],
           [7.2020483, 7.33431  , 7.3966174, ..., 8.096302 , 8.005932 ,
            8.005932 ],
           [7.240784 , 7.365852 , 7.39515  , ..., 7.828516 , 7.7119536,
            7.7119536],
           ...,
           [7.2204494, 7.3662806, 7.388879 , ..., 7.907409 , 8.00433  ,
            8.00433  ],
           [7.29451  , 7.3260627, 7.4240465, ..., 7.852135 , 7.713879 ,
            7.713879 ],
           [7.2906227, 7.3745756, 7.4152393, ..., 8.110569 , 7.994925 ,
            7.994925 ]], dtype=float32)

• Indexes: (2)
  • trajectory
    PandasIndex

    PandasIndex(Index([  0,   1,   2,   3,   4,   5,   6,   7,   8,   9,
           ...
           990, 991, 992, 993, 994, 995, 996, 997, 998, 999],
          dtype='int64', name='trajectory', length=1000))

  • time
    PandasIndex

    PandasIndex(DatetimeIndex(['2025-03-05 21:56:07.344141', '2025-03-05 22:56:07.344141',
                   '2025-03-05 23:56:07.344141', '2025-03-06 00:56:07.344141',
                   '2025-03-06 01:56:07.344141', '2025-03-06 02:56:07.344141',
                   '2025-03-06 03:56:07.344141', '2025-03-06 04:56:07.344141',
                   '2025-03-06 05:56:07.344141', '2025-03-06 06:56:07.344141',
                   '2025-03-06 07:56:07.344141', '2025-03-06 08:56:07.344141',
                   '2025-03-06 09:56:07.344141', '2025-03-06 10:56:07.344141',
                   '2025-03-06 11:56:07.344141', '2025-03-06 12:56:07.344141',
                   '2025-03-06 13:56:07.344141', '2025-03-06 14:56:07.344141',
                   '2025-03-06 15:56:07.344141', '2025-03-06 16:56:07.344141',
                   '2025-03-06 17:56:07.344141', '2025-03-06 18:56:07.344141',
                   '2025-03-06 19:56:07.344141', '2025-03-06 20:56:07.344141',
                   '2025-03-06 21:56:07.344141'],
                  dtype='datetime64[ns]', name='time', freq='h'))

• Attributes: (96)

  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-03-05 21:56:11.695381

  date_created :

  2025-03-05T21:56:11.695395

  source :

  Output from simulation with OpenDrift

  model_url :

  https://github.com/OpenDrift/opendrift

  opendrift_class :

  ShipDrift

  opendrift_module :

  opendrift.models.shipdrift

  readers :

  odict_keys(['LazyReader: https://thredds.met.no/thredds/dodsC/cmems/topaz6/dataset-topaz6-arc-15min-3km-be.ncml', 'LazyReader: https://thredds.met.no/thredds/dodsC/mepslatest/meps_lagged_6_h_latest_2_5km_latest.nc', 'LazyReader: https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml', 'global_landmask'])

  time_coverage_start :

  2025-03-05 21:56:07.344141

  time_step_calculation :

  1:00:00

  time_step_output :

  1:00:00

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

  None

  config_environment:fallback:land_binary_mask :

  None

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

  None

  config_environment:fallback:sea_surface_wave_significant_height :

  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_general:use_auto_landmask :

  True

  config_drift:current_uncertainty :

  0

  config_drift:current_uncertainty_uniform :

  0

  config_drift:max_speed :

  2.0

  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 :

  100.0

  config_drift:profiles_depth :

  50

  config_drift:wind_uncertainty :

  0

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

  random

  config_seed:length :

  80

  config_seed:height :

  8

  config_seed:draft :

  4.0

  config_seed:beam :

  10

  config_seed:wind_drag_coeff :

  1

  config_seed:water_drag_coeff :

  1

  config_seed:jibeProbability :

  0.04

  opendrift_version :

  1.13.0

  seed_geojson :

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

  simulation_time :

  2025-03-05 21:56:11.700631

  reader_x_wind :

  ['https://thredds.met.no/thredds/dodsC/mepslatest/meps_lagged_6_h_latest_2_5km_latest.nc']

  reader_y_wind :

  ['https://thredds.met.no/thredds/dodsC/mepslatest/meps_lagged_6_h_latest_2_5km_latest.nc']

  reader_land_binary_mask :

  ['global_landmask']

  reader_x_sea_water_velocity :

  ['https://thredds.met.no/thredds/dodsC/cmems/topaz6/dataset-topaz6-arc-15min-3km-be.ncml', 'https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml']

  reader_y_sea_water_velocity :

  ['https://thredds.met.no/thredds/dodsC/cmems/topaz6/dataset-topaz6-arc-15min-3km-be.ncml', 'https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml']

  reader_sea_surface_wave_stokes_drift_x_velocity :

  ['https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml']

  reader_sea_surface_wave_stokes_drift_y_velocity :

  ['https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml']

  reader_sea_surface_wave_significant_height :

  ['https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml']

  reader_sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment :

  ['https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml']

  time_coverage_end :

  2025-03-06 21:56:07.344141

  time_coverage_duration :

  P1DT0H0M0S

  time_coverage_resolution :

  P0DT1H0M0S

  performance :

  -------------------- Reader performance: -------------------- global_landmask 0:00:00.0 total 0:00:00.0 preparing 0:00:00.0 reading 0:00:00.0 masking -------------------- https://thredds.met.no/thredds/dodsC/cmems/topaz6/dataset-topaz6-arc-15min-3km-be.ncml 0:00:16.3 total 0:00:00.0 preparing 0:00:16.1 reading 0:00:00.0 interpolation 0:00:00.0 interpolation_time 0:00:00.2 rotating vectors 0:00:00.0 masking -------------------- https://thredds.met.no/thredds/dodsC/mepslatest/meps_lagged_6_h_latest_2_5km_latest.nc 0:00:22.8 total 0:00:00.0 preparing 0:00:22.6 reading 0:00:00.3 interpolation 0:00:00.0 interpolation_time 0:00:00.2 rotating vectors 0:00:00.0 masking -------------------- https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml 0:00:18.7 total 0:00:00.0 preparing 0:00:18.4 reading 0:00:00.0 interpolation 0:00:00.0 interpolation_time 0:00:00.2 rotating vectors 0:00:00.0 masking -------------------- Performance: 1:08.7 total time 4.3 configuration 0.0 preparing main loop 0.0 moving elements to ocean 1:04.3 main loop 0.6 updating elements 0.0 cleaning up --------------------

  geospatial_bounds_crs :

  EPSG:4326

  geospatial_bounds_vertical_crs :

  EPSG:5831

  geospatial_lat_min :

  62.973133

  geospatial_lat_max :

  63.456875

  geospatial_lat_units :

  degrees_north

  geospatial_lat_resolution :

  point

  geospatial_lon_min :

  4.9398756

  geospatial_lon_max :

  7.024521

  geospatial_lon_units :

  degrees_east

  geospatial_lon_resolution :

  point

  runtime :

  0:01:08.712187

  geospatial_vertical_min :

  0.0

  geospatial_vertical_max :

  0.0

  geospatial_vertical_positive :

  up

Print and plot results

print(o)
o.plot(linecolor='orientation')
#o.animation(color='orientation', markersize=20, filename="orientation.gif", legend=['left','right'],colorbar=False,cmap='bwr')
o.animation(color='orientation', legend=['left','right'], markersize=20, colorbar=False, cmap='bwr')

===========================
--------------------
Reader performance:
--------------------
global_landmask
 0:00:00.0  total
 0:00:00.0  preparing
 0:00:00.0  reading
 0:00:00.0  masking
--------------------
https://thredds.met.no/thredds/dodsC/cmems/topaz6/dataset-topaz6-arc-15min-3km-be.ncml
 0:00:16.3  total
 0:00:00.0  preparing
 0:00:16.1  reading
 0:00:00.0  interpolation
 0:00:00.0  interpolation_time
 0:00:00.2  rotating vectors
 0:00:00.0  masking
--------------------
https://thredds.met.no/thredds/dodsC/mepslatest/meps_lagged_6_h_latest_2_5km_latest.nc
 0:00:22.8  total
 0:00:00.0  preparing
 0:00:22.6  reading
 0:00:00.3  interpolation
 0:00:00.0  interpolation_time
 0:00:00.2  rotating vectors
 0:00:00.0  masking
--------------------
https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml
 0:00:18.7  total
 0:00:00.0  preparing
 0:00:18.4  reading
 0:00:00.0  interpolation
 0:00:00.0  interpolation_time
 0:00:00.2  rotating vectors
 0:00:00.0  masking
--------------------
Performance:
 1:08.7 total time
    4.3 configuration
    0.0 preparing main loop
      0.0 moving elements to ocean
 1:04.3 main loop
      0.6 updating elements
    0.0 cleaning up
--------------------
===========================
Model:  ShipDrift     (OpenDrift version 1.13.0)
        1000 active ShipObject particles  (0 deactivated, 0 scheduled)
-------------------
Environment variables:
  -----
  land_binary_mask
     1) global_landmask
  -----
  x_sea_water_velocity
  y_sea_water_velocity
     1) https://thredds.met.no/thredds/dodsC/cmems/topaz6/dataset-topaz6-arc-15min-3km-be.ncml
     2) https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml
  -----
  x_wind
  y_wind
     1) https://thredds.met.no/thredds/dodsC/mepslatest/meps_lagged_6_h_latest_2_5km_latest.nc
  -----
  sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment
  sea_surface_wave_significant_height
  sea_surface_wave_stokes_drift_x_velocity
  sea_surface_wave_stokes_drift_y_velocity
     1) https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml

Discarded readers:

Time:
        Start: 2025-03-05 21:56:07.344141 UTC
        Present: 2025-03-06 21:56:07.344141 UTC
        Calculation steps: 24 * 1:00:00 - total time: 1 day, 0:00:00
        Output steps: 25 * 1:00:00
===========================

21:57:49 INFO    opendrift:4600: Saving animation to /root/project/docs/source/gallery/animations/example_shipdrift_0.gif...
21:58:04 INFO    opendrift:3036: Time to make animation: 0:00:15.837047