Note
Go to the end to download the full example code.
Sediment drift with resuspension
from datetime import timedelta, datetime
from opendrift.readers import reader_oscillating
from opendrift.models.sedimentdrift import SedimentDrift
Constructing an artificial current field where x- and y-components are oscilating with different amplitude and period
reader_oscx = reader_oscillating.Reader('x_sea_water_velocity',
amplitude=0.6, zero_time=datetime.utcnow())
reader_oscy = reader_oscillating.Reader('y_sea_water_velocity',
amplitude=.3, period_seconds=3600*5, zero_time=datetime.utcnow())
o = SedimentDrift(loglevel=50) # 0 for debug output
if True:
o.add_reader([reader_oscx, reader_oscy])
o.set_config('environment:fallback:y_wind', -6)
o.set_config('environment:fallback:x_wind', -3)
o.set_config('environment:fallback:sea_floor_depth_below_sea_level', 30) # 100m depth
else: # Using live data from Thredds instead of oscillating currents
o.add_readers_from_list([
'https://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be'])
Set threshold for bottom resuspension
o.set_config('vertical_mixing:resuspension_threshold', .5)
# Adding some horizontal and vertical diffusion
o.set_config('drift:current_uncertainty', 0.1)
o.set_config('drift:wind_uncertainty', 1)
o.set_config('vertical_mixing:diffusivitymodel', 'windspeed_Large1994')
#o.set_config('vertical_mixing:diffusivitymodel', 'environment')
Seeding sediments
o.seed_elements(lon=4.65, lat=60, number=10000,
time=[datetime.utcnow(), datetime.utcnow()+timedelta(hours=6)],
terminal_velocity=-.01) # 1 cm/s settling speed
o.run(time_step=1800, time_step_output=1800, duration=timedelta(hours=72))
<xarray.Dataset> Size: 145MB
Dimensions: (
trajectory: 10000,
time: 145)
Coordinates:
* trajectory (trajectory) int64 80kB ...
* time (time) datetime64[ns] 1kB ...
Data variables: (12/25)
status (trajectory, time) float32 6MB ...
moving (trajectory, time) float32 6MB ...
age_seconds (trajectory, time) float32 6MB ...
origin_marker (trajectory, time) float32 6MB ...
lon (trajectory, time) float32 6MB ...
lat (trajectory, time) float32 6MB ...
... ...
sea_surface_wave_period_at_variance_spectral_density_maximum (trajectory, time) float32 6MB ...
sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment (trajectory, time) float32 6MB ...
land_binary_mask (trajectory, time) float32 6MB ...
ocean_vertical_diffusivity (trajectory, time) float32 6MB ...
ocean_mixed_layer_thickness (trajectory, time) float32 6MB ...
sea_floor_depth_below_sea_level (trajectory, time) float32 6MB ...
Attributes: (12/122)
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
- trajectory: 10000
- time: 145
- trajectory(trajectory)int640 1 2 3 4 ... 9996 9997 9998 9999
- cf_role :
- trajectory_id
- dtype :
- <class 'numpy.int32'>
array([ 0, 1, 2, ..., 9997, 9998, 9999])
- time(time)datetime64[ns]2025-02-18T17:08:32.688998 ... 2...
- standard_name :
- time
- long_name :
- time
- dtype :
- <class 'numpy.float64'>
array(['2025-02-18T17:08:32.688998000', '2025-02-18T17:38:32.688998000', '2025-02-18T18:08:32.688998000', '2025-02-18T18:38:32.688998000', '2025-02-18T19:08:32.688998000', '2025-02-18T19:38:32.688998000', '2025-02-18T20:08:32.688998000', '2025-02-18T20:38:32.688998000', '2025-02-18T21:08:32.688998000', '2025-02-18T21:38:32.688998000', '2025-02-18T22:08:32.688998000', '2025-02-18T22:38:32.688998000', '2025-02-18T23:08:32.688998000', '2025-02-18T23:38:32.688998000', '2025-02-19T00:08:32.688998000', '2025-02-19T00:38:32.688998000', '2025-02-19T01:08:32.688998000', '2025-02-19T01:38:32.688998000', '2025-02-19T02:08:32.688998000', '2025-02-19T02:38:32.688998000', '2025-02-19T03:08:32.688998000', '2025-02-19T03:38:32.688998000', '2025-02-19T04:08:32.688998000', '2025-02-19T04:38:32.688998000', '2025-02-19T05:08:32.688998000', '2025-02-19T05:38:32.688998000', '2025-02-19T06:08:32.688998000', '2025-02-19T06:38:32.688998000', '2025-02-19T07:08:32.688998000', '2025-02-19T07:38:32.688998000', '2025-02-19T08:08:32.688998000', '2025-02-19T08:38:32.688998000', '2025-02-19T09:08:32.688998000', '2025-02-19T09:38:32.688998000', '2025-02-19T10:08:32.688998000', '2025-02-19T10:38:32.688998000', '2025-02-19T11:08:32.688998000', '2025-02-19T11:38:32.688998000', '2025-02-19T12:08:32.688998000', '2025-02-19T12:38:32.688998000', '2025-02-19T13:08:32.688998000', '2025-02-19T13:38:32.688998000', '2025-02-19T14:08:32.688998000', '2025-02-19T14:38:32.688998000', '2025-02-19T15:08:32.688998000', '2025-02-19T15:38:32.688998000', '2025-02-19T16:08:32.688998000', '2025-02-19T16:38:32.688998000', '2025-02-19T17:08:32.688998000', '2025-02-19T17:38:32.688998000', '2025-02-19T18:08:32.688998000', '2025-02-19T18:38:32.688998000', '2025-02-19T19:08:32.688998000', '2025-02-19T19:38:32.688998000', '2025-02-19T20:08:32.688998000', '2025-02-19T20:38:32.688998000', '2025-02-19T21:08:32.688998000', '2025-02-19T21:38:32.688998000', '2025-02-19T22:08:32.688998000', '2025-02-19T22:38:32.688998000', '2025-02-19T23:08:32.688998000', '2025-02-19T23:38:32.688998000', '2025-02-20T00:08:32.688998000', '2025-02-20T00:38:32.688998000', '2025-02-20T01:08:32.688998000', '2025-02-20T01:38:32.688998000', '2025-02-20T02:08:32.688998000', '2025-02-20T02:38:32.688998000', '2025-02-20T03:08:32.688998000', '2025-02-20T03:38:32.688998000', '2025-02-20T04:08:32.688998000', '2025-02-20T04:38:32.688998000', '2025-02-20T05:08:32.688998000', '2025-02-20T05:38:32.688998000', '2025-02-20T06:08:32.688998000', '2025-02-20T06:38:32.688998000', '2025-02-20T07:08:32.688998000', '2025-02-20T07:38:32.688998000', '2025-02-20T08:08:32.688998000', '2025-02-20T08:38:32.688998000', '2025-02-20T09:08:32.688998000', '2025-02-20T09:38:32.688998000', '2025-02-20T10:08:32.688998000', '2025-02-20T10:38:32.688998000', '2025-02-20T11:08:32.688998000', '2025-02-20T11:38:32.688998000', '2025-02-20T12:08:32.688998000', '2025-02-20T12:38:32.688998000', '2025-02-20T13:08:32.688998000', '2025-02-20T13:38:32.688998000', '2025-02-20T14:08:32.688998000', '2025-02-20T14:38:32.688998000', '2025-02-20T15:08:32.688998000', '2025-02-20T15:38:32.688998000', '2025-02-20T16:08:32.688998000', '2025-02-20T16:38:32.688998000', '2025-02-20T17:08:32.688998000', '2025-02-20T17:38:32.688998000', '2025-02-20T18:08:32.688998000', '2025-02-20T18:38:32.688998000', '2025-02-20T19:08:32.688998000', '2025-02-20T19:38:32.688998000', '2025-02-20T20:08:32.688998000', '2025-02-20T20:38:32.688998000', '2025-02-20T21:08:32.688998000', '2025-02-20T21:38:32.688998000', '2025-02-20T22:08:32.688998000', '2025-02-20T22:38:32.688998000', '2025-02-20T23:08:32.688998000', '2025-02-20T23:38:32.688998000', '2025-02-21T00:08:32.688998000', '2025-02-21T00:38:32.688998000', '2025-02-21T01:08:32.688998000', '2025-02-21T01:38:32.688998000', '2025-02-21T02:08:32.688998000', '2025-02-21T02:38:32.688998000', '2025-02-21T03:08:32.688998000', '2025-02-21T03:38:32.688998000', '2025-02-21T04:08:32.688998000', '2025-02-21T04:38:32.688998000', '2025-02-21T05:08:32.688998000', '2025-02-21T05:38:32.688998000', '2025-02-21T06:08:32.688998000', '2025-02-21T06:38:32.688998000', '2025-02-21T07:08:32.688998000', '2025-02-21T07:38:32.688998000', '2025-02-21T08:08:32.688998000', '2025-02-21T08:38:32.688998000', '2025-02-21T09:08:32.688998000', '2025-02-21T09:38:32.688998000', '2025-02-21T10:08:32.688998000', '2025-02-21T10:38:32.688998000', '2025-02-21T11:08:32.688998000', '2025-02-21T11:38:32.688998000', '2025-02-21T12:08:32.688998000', '2025-02-21T12:38:32.688998000', '2025-02-21T13:08:32.688998000', '2025-02-21T13:38:32.688998000', '2025-02-21T14:08:32.688998000', '2025-02-21T14:38:32.688998000', '2025-02-21T15:08:32.688998000', '2025-02-21T15:38:32.688998000', '2025-02-21T16:08:32.688998000', '2025-02-21T16:38:32.688998000', '2025-02-21T17:08:32.688998000'], dtype='datetime64[ns]')
- status(trajectory, time)float320.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- dtype :
- <class 'numpy.int32'>
- valid_range :
- [0 0]
- flag_values :
- [0]
- flag_meanings :
- active
array([[ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], ..., [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.]], dtype=float32) - moving(trajectory, time)float321.0 1.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- dtype :
- <class 'numpy.int32'>
- minval :
- 0
- maxval :
- 1
array([[ 1., 1., 0., ..., 0., 0., 1.], [ 1., 0., 0., ..., 0., 0., 0.], [ 1., 1., 1., ..., 0., 0., 0.], ..., [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.]], dtype=float32) - age_seconds(trajectory, time)float320.0 1.8e+03 ... 2.376e+05 2.394e+05
- dtype :
- <class 'numpy.float32'>
- units :
- s
- minval :
- 0.0
- maxval :
- 259200.0
array([[ 0., 1800., 3600., ..., 255600., 257400., 259200.], [ 0., 1800., 3600., ..., 255600., 257400., 259200.], [ 0., 1800., 3600., ..., 255600., 257400., 259200.], ..., [ nan, nan, nan, ..., 235800., 237600., 239400.], [ nan, nan, nan, ..., 235800., 237600., 239400.], [ nan, nan, nan, ..., 235800., 237600., 239400.]], dtype=float32) - origin_marker(trajectory, time)float320.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- dtype :
- <class 'numpy.int32'>
- unit :
- description :
- An integer kept constant during the simulation. Different values may be used for different seedings, to separate elements during analysis. With GUI, only a single seeding is possible.
- flag_values :
- [0]
- flag_meanings :
- Seed_0
- minval :
- 0
- maxval :
- 0
array([[ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], ..., [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.]], dtype=float32) - lon(trajectory, time)float324.65 4.654 4.655 ... 4.98 4.98 4.98
- dtype :
- <class 'numpy.float32'>
- units :
- degrees_east
- standard_name :
- longitude
- long_name :
- longitude
- axis :
- X
- minval :
- 4.402975
- maxval :
- 5.258207
array([[4.65 , 4.6537676, 4.654725 , ..., 5.06053 , 5.06053 , 5.06053 ], [4.65 , 4.649303 , 4.649303 , ..., 5.0203714, 5.0203714, 5.0203714], [4.65 , 4.6506047, 4.648504 , ..., 5.1139092, 5.1139092, 5.1139092], ..., [ nan, nan, nan, ..., 4.9113955, 4.9113955, 4.9113955], [ nan, nan, nan, ..., 4.99289 , 4.99289 , 4.99289 ], [ nan, nan, nan, ..., 4.979738 , 4.979738 , 4.979738 ]], dtype=float32) - lat(trajectory, time)float3260.0 60.0 60.0 ... 59.98 59.98
- dtype :
- <class 'numpy.float32'>
- units :
- degrees_north
- standard_name :
- latitude
- long_name :
- latitude
- axis :
- Y
- minval :
- 59.90722
- maxval :
- 60.098274
array([[60. , 59.99876 , 59.999546, ..., 60.00605 , 60.00605 , 60.00605 ], [60. , 59.99921 , 59.99921 , ..., 59.979485, 59.979485, 59.979485], [60. , 59.99663 , 59.99691 , ..., 59.96593 , 59.96593 , 59.96593 ], ..., [ nan, nan, nan, ..., 59.983097, 59.983097, 59.983097], [ nan, nan, nan, ..., 59.971367, 59.971367, 59.971367], [ nan, nan, nan, ..., 59.983715, 59.983715, 59.983715]], dtype=float32) - z(trajectory, time)float320.0 -17.91 -30.14 ... -30.07 -30.07
- dtype :
- <class 'numpy.float32'>
- units :
- m
- standard_name :
- z
- long_name :
- vertical position
- axis :
- Z
- positive :
- up
- minval :
- -30.599997
- maxval :
- 0.0
array([[ 0. , -17.914232, -30.14472 , ..., -30.320045, -30.320045, -30.310045], [ 0. , -30.098269, -30.098269, ..., -30.24326 , -30.24326 , -30.24326 ], [ 0. , -17.514284, -26.29749 , ..., -30.179518, -30.179518, -30.179518], ..., [ nan, nan, nan, ..., -30.108414, -30.108414, -30.108414], [ nan, nan, nan, ..., -30.276293, -30.276293, -30.276293], [ nan, nan, nan, ..., -30.07484 , -30.07484 , -30.07484 ]], dtype=float32) - wind_drift_factor(trajectory, time)float320.02 0.02 0.02 ... 0.02 0.02 0.02
- dtype :
- <class 'numpy.float32'>
- units :
- 1
- description :
- Elements at surface are moved with this fraction of the vind vector, in addition to currents and Stokes drift
- minval :
- 0.02
- maxval :
- 0.02
array([[0.02, 0.02, 0.02, ..., 0.02, 0.02, 0.02], [0.02, 0.02, 0.02, ..., 0.02, 0.02, 0.02], [0.02, 0.02, 0.02, ..., 0.02, 0.02, 0.02], ..., [ nan, nan, nan, ..., 0.02, 0.02, 0.02], [ nan, nan, nan, ..., 0.02, 0.02, 0.02], [ nan, nan, nan, ..., 0.02, 0.02, 0.02]], dtype=float32) - current_drift_factor(trajectory, time)float321.0 1.0 1.0 1.0 ... 1.0 1.0 1.0 1.0
- dtype :
- <class 'numpy.float32'>
- units :
- 1
- description :
- Elements are moved with this fraction of the current vector, in addition to currents and Stokes drift
- minval :
- 1.0
- maxval :
- 1.0
array([[ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], [ 1., 1., 1., ..., 1., 1., 1.], ..., [nan, nan, nan, ..., 1., 1., 1.], [nan, nan, nan, ..., 1., 1., 1.], [nan, nan, nan, ..., 1., 1., 1.]], dtype=float32) - terminal_velocity(trajectory, time)float32-0.01 -0.01 -0.01 ... -0.01 -0.01
- dtype :
- <class 'numpy.float32'>
- units :
- m/s
- minval :
- -0.01
- maxval :
- -0.01
array([[-0.01, -0.01, -0.01, ..., -0.01, -0.01, -0.01], [-0.01, -0.01, -0.01, ..., -0.01, -0.01, -0.01], [-0.01, -0.01, -0.01, ..., -0.01, -0.01, -0.01], ..., [ nan, nan, nan, ..., -0.01, -0.01, -0.01], [ nan, nan, nan, ..., -0.01, -0.01, -0.01], [ nan, nan, nan, ..., -0.01, -0.01, -0.01]], dtype=float32) - settled(trajectory, time)float320.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- dtype :
- <class 'numpy.uint8'>
- units :
- 1
- minval :
- 0
- maxval :
- 0
array([[ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], ..., [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.]], dtype=float32) - x_sea_water_velocity(trajectory, time)float320.1764 0.02968 ... -0.08856
- minval :
- -1.0181199
- maxval :
- 1.0368029
array([[ 0.17641148, 0.02968115, 0.10041095, ..., 0.31032035, 0.09746848, 0.09746848], [ 0.04002196, 0.08297363, 0.11504264, ..., 0.21695481, -0.02193601, -0.02193601], [ 0.09788004, -0.06513204, 0.00947064, ..., 0.1016138 , 0.08584205, 0.08584205], ..., [ nan, nan, nan, ..., -0.0161072 , 0.07293599, 0.07293599], [ nan, nan, nan, ..., 0.14972244, 0.13996296, 0.13996296], [ nan, nan, nan, ..., 0.05969569, -0.0885559 , -0.0885559 ]], dtype=float32) - y_sea_water_velocity(trajectory, time)float320.005172 0.04867 ... 0.3579 0.3579
- minval :
- -0.7497912
- maxval :
- 0.71832526
array([[ 0.00517245, 0.04866524, 0.1517709 , ..., 0.2162844 , 0.5807213 , 0.5807213 ], [ 0.08674262, -0.07284362, 0.23437195, ..., 0.24177946, 0.26389802, 0.26389802], [-0.0848171 , 0.01733063, 0.15258947, ..., 0.19646555, 0.32843137, 0.32843137], ..., [ nan, nan, nan, ..., 0.18403319, 0.3729422 , 0.3729422 ], [ nan, nan, nan, ..., 0.23300299, 0.15925106, 0.15925106], [ nan, nan, nan, ..., 0.21413097, 0.35787895, 0.35787895]], dtype=float32) - sea_surface_height(trajectory, time)float320.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- minval :
- 0.0
- maxval :
- 0.0
array([[ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], ..., [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.]], dtype=float32) - upward_sea_water_velocity(trajectory, time)float320.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- minval :
- 0.0
- maxval :
- 0.0
array([[ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], ..., [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.]], dtype=float32) - x_wind(trajectory, time)float32-2.981 -4.244 ... -3.338 -3.338
- minval :
- -8.002298
- maxval :
- 1.8089983
array([[-2.981078 , -4.243658 , -2.8865137, ..., -4.9380593, -3.34425 , -3.34425 ], [-3.082034 , -2.827815 , -2.0777311, ..., -3.5350225, -2.32015 , -2.32015 ], [-3.957158 , -2.3757708, -1.5997821, ..., -2.5414686, -4.458367 , -4.458367 ], ..., [ nan, nan, nan, ..., -2.7222812, -1.0365078, -1.0365078], [ nan, nan, nan, ..., -2.0151415, -3.3912978, -3.3912978], [ nan, nan, nan, ..., -2.4119408, -3.3377228, -3.3377228]], dtype=float32) - y_wind(trajectory, time)float32-4.097 -5.883 ... -5.671 -5.671
- minval :
- -10.789874
- maxval :
- -1.0925639
array([[-4.0973816, -5.883371 , -6.234717 , ..., -7.2351575, -5.630373 , -5.630373 ], [-6.7756643, -6.090957 , -4.563066 , ..., -6.344755 , -7.9961433, -7.9961433], [-6.18809 , -5.2804637, -5.106999 , ..., -4.6859593, -7.3998036, -7.3998036], ..., [ nan, nan, nan, ..., -5.7076945, -6.708225 , -6.708225 ], [ nan, nan, nan, ..., -5.610828 , -5.963494 , -5.963494 ], [ nan, nan, nan, ..., -6.3788004, -5.67101 , -5.67101 ]], dtype=float32) - sea_surface_wave_stokes_drift_x_velocity(trajectory, time)float320.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- minval :
- 0.0
- maxval :
- 0.0
array([[ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], ..., [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.]], dtype=float32) - sea_surface_wave_stokes_drift_y_velocity(trajectory, time)float320.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- minval :
- 0.0
- maxval :
- 0.0
array([[ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], ..., [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.]], dtype=float32) - sea_surface_wave_period_at_variance_spectral_density_maximum(trajectory, time)float320.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- minval :
- 0.0
- maxval :
- 0.0
array([[ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], ..., [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.]], dtype=float32) - sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment(trajectory, time)float324.33 6.198 5.871 ... 5.623 5.623
- minval :
- 1.5306671
- maxval :
- 9.854761
array([[4.329676 , 6.1984577, 5.870633 , ..., 7.484886 , 5.5956483, 5.5956483], [6.3604164, 5.738093 , 4.284174 , ..., 6.2060866, 7.114278 , 7.114278 ], [6.276234 , 4.947642 , 4.5728745, ..., 4.555 , 7.3818607, 7.3818607], ..., [ nan, nan, nan, ..., 5.403376 , 5.8 , 5.8 ], [ nan, nan, nan, ..., 5.09412 , 5.861952 , 5.861952 ], [ nan, nan, nan, ..., 5.8271236, 5.622699 , 5.622699 ]], dtype=float32) - land_binary_mask(trajectory, time)float320.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0
- minval :
- 0.0
- maxval :
- 0.0
array([[ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], [ 0., 0., 0., ..., 0., 0., 0.], ..., [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.], [nan, nan, nan, ..., 0., 0., 0.]], dtype=float32) - ocean_vertical_diffusivity(trajectory, time)float320.02 0.02 0.02 ... 0.02 0.02 0.02
- minval :
- 0.02
- maxval :
- 0.02
array([[0.02, 0.02, 0.02, ..., 0.02, 0.02, 0.02], [0.02, 0.02, 0.02, ..., 0.02, 0.02, 0.02], [0.02, 0.02, 0.02, ..., 0.02, 0.02, 0.02], ..., [ nan, nan, nan, ..., 0.02, 0.02, 0.02], [ nan, nan, nan, ..., 0.02, 0.02, 0.02], [ nan, nan, nan, ..., 0.02, 0.02, 0.02]], dtype=float32) - ocean_mixed_layer_thickness(trajectory, time)float3250.0 50.0 50.0 ... 50.0 50.0 50.0
- minval :
- 50.0
- maxval :
- 50.0
array([[50., 50., 50., ..., 50., 50., 50.], [50., 50., 50., ..., 50., 50., 50.], [50., 50., 50., ..., 50., 50., 50.], ..., [nan, nan, nan, ..., 50., 50., 50.], [nan, nan, nan, ..., 50., 50., 50.], [nan, nan, nan, ..., 50., 50., 50.]], dtype=float32) - sea_floor_depth_below_sea_level(trajectory, time)float3230.0 30.0 30.0 ... 30.0 30.0 30.0
- minval :
- 30.0
- maxval :
- 30.0
array([[30., 30., 30., ..., 30., 30., 30.], [30., 30., 30., ..., 30., 30., 30.], [30., 30., 30., ..., 30., 30., 30.], ..., [nan, nan, nan, ..., 30., 30., 30.], [nan, nan, nan, ..., 30., 30., 30.], [nan, nan, nan, ..., 30., 30., 30.]], dtype=float32)
- trajectoryPandasIndex
PandasIndex(Index([ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, ... 9990, 9991, 9992, 9993, 9994, 9995, 9996, 9997, 9998, 9999], dtype='int64', name='trajectory', length=10000)) - timePandasIndex
PandasIndex(DatetimeIndex(['2025-02-18 17:08:32.688998', '2025-02-18 17:38:32.688998', '2025-02-18 18:08:32.688998', '2025-02-18 18:38:32.688998', '2025-02-18 19:08:32.688998', '2025-02-18 19:38:32.688998', '2025-02-18 20:08:32.688998', '2025-02-18 20:38:32.688998', '2025-02-18 21:08:32.688998', '2025-02-18 21:38:32.688998', ... '2025-02-21 12:38:32.688998', '2025-02-21 13:08:32.688998', '2025-02-21 13:38:32.688998', '2025-02-21 14:08:32.688998', '2025-02-21 14:38:32.688998', '2025-02-21 15:08:32.688998', '2025-02-21 15:38:32.688998', '2025-02-21 16:08:32.688998', '2025-02-21 16:38:32.688998', '2025-02-21 17:08:32.688998'], dtype='datetime64[ns]', name='time', length=145, freq='30min'))
- 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-02-18 17:08:38.205032
- date_created :
- 2025-02-18T17:08:38.205060
- source :
- Output from simulation with OpenDrift
- model_url :
- https://github.com/OpenDrift/opendrift
- opendrift_class :
- SedimentDrift
- opendrift_module :
- opendrift.models.sedimentdrift
- readers :
- odict_keys(['oscillating_reader', 'oscillating_reader_0', 'global_landmask'])
- time_coverage_start :
- 2025-02-18 17:08:32.688998
- time_step_calculation :
- 0:30:00
- time_step_output :
- 0:30:00
- config_environment:constant:x_sea_water_velocity :
- None
- config_environment:fallback:x_sea_water_velocity :
- 0
- config_environment:constant:y_sea_water_velocity :
- None
- config_environment:fallback:y_sea_water_velocity :
- 0
- 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:x_wind :
- None
- config_environment:fallback:x_wind :
- -3.0
- config_environment:constant:y_wind :
- None
- config_environment:fallback:y_wind :
- -6.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:land_binary_mask :
- None
- config_environment:fallback:land_binary_mask :
- None
- config_environment:constant:ocean_vertical_diffusivity :
- None
- config_environment:fallback:ocean_vertical_diffusivity :
- 0.02
- config_environment:constant:ocean_mixed_layer_thickness :
- None
- config_environment:fallback:ocean_mixed_layer_thickness :
- 50
- config_environment:constant:sea_floor_depth_below_sea_level :
- None
- config_environment:fallback:sea_floor_depth_below_sea_level :
- 30.0
- config_general:use_auto_landmask :
- True
- config_drift:current_uncertainty :
- 0.1
- 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 :
- previous
- 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 :
- 1.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:wind_drift_factor :
- 0.02
- config_seed:current_drift_factor :
- 1
- config_seed:terminal_velocity :
- -0.001
- config_seed:settled :
- 0
- config_drift:vertical_advection :
- True
- config_drift:vertical_mixing :
- True
- config_vertical_mixing:timestep :
- 60
- config_vertical_mixing:diffusivitymodel :
- windspeed_Large1994
- 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_vertical_mixing:resuspension_threshold :
- 0.5
- opendrift_version :
- 1.13.0
- seed_geojson :
- {"features": [], "type": "FeatureCollection"}
- simulation_time :
- 2025-02-18 17:08:38.261175
- reader_x_sea_water_velocity :
- ['oscillating_reader']
- reader_y_sea_water_velocity :
- ['oscillating_reader_0']
- reader_sea_surface_height :
- 0
- reader_upward_sea_water_velocity :
- 0
- reader_x_wind :
- -3.0
- reader_y_wind :
- -6.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_land_binary_mask :
- ['global_landmask']
- reader_ocean_vertical_diffusivity :
- 0.02
- reader_ocean_mixed_layer_thickness :
- 50
- reader_sea_floor_depth_below_sea_level :
- 30.0
- time_coverage_end :
- 2025-02-21 17:08:32.688998
- time_coverage_duration :
- P3DT0H0M0S
- time_coverage_resolution :
- P0DT0H30M0S
- performance :
- -------------------- Reader performance: -------------------- oscillating_reader 0:00:00.0 total 0:00:00.0 preparing 0:00:00.0 reading 0:00:00.0 masking -------------------- oscillating_reader_0 0:00:00.0 total 0:00:00.0 preparing 0:00:00.0 reading 0:00:00.0 masking -------------------- global_landmask 0:00:03.1 total 0:00:00.0 preparing 0:00:03.1 reading 0:00:00.0 masking -------------------- Performance: 43.5 total time 5.1 configuration 0.5 preparing main loop 0.0 moving elements to ocean 37.8 main loop 12.7 updating elements 11.8 vertical mixing 0.0 cleaning up --------------------
- geospatial_bounds_crs :
- EPSG:4326
- geospatial_bounds_vertical_crs :
- EPSG:5831
- geospatial_lat_min :
- 59.90722
- geospatial_lat_max :
- 60.098274
- geospatial_lat_units :
- degrees_north
- geospatial_lat_resolution :
- point
- geospatial_lon_min :
- 4.402975
- geospatial_lon_max :
- 5.258207
- geospatial_lon_units :
- degrees_east
- geospatial_lon_resolution :
- point
- runtime :
- 0:00:43.598158
- geospatial_vertical_min :
- -30.599997
- geospatial_vertical_max :
- 0.0
- geospatial_vertical_positive :
- up
Plotting the depth vs time
o.plot_property('z')
Animate sediment particles
o.animation(color='moving', colorbar=False, legend=['Sedimented', 'Moving'], fast=True, buffer=.01)
#o.animation_profile()
Total running time of the script: (2 minutes 26.850 seconds)