Model landmask

Comparing two simulation runs, with landmask from ocean model and GSHHG

from datetime import timedelta

from opendrift import test_data_folder as tdf
from opendrift.readers import reader_ROMS_native
from opendrift.models.oceandrift import OceanDrift

lon = 14.75; lat = 68.1

o = OceanDrift(loglevel=20)

reader_nordic = reader_ROMS_native.Reader(tdf +
    '2Feb2016_Nordic_sigma_3d/Nordic-4km_SLEVELS_avg_00_subset2Feb2016.nc')

07:41:32 INFO    opendrift:513: OpenDriftSimulation initialised (version 1.14.2 / v1.14.2-90-g1dd1995)
07:41:32 INFO    opendrift.readers:61: Opening file with xr.open_dataset
07:41:32 WARNING opendrift.readers.reader_ROMS_native:282: Duplicate variables for land_binary_mask, selecting mask_psi, and discarding mask_rho
07:41:32 WARNING opendrift.readers.reader_ROMS_native:282: Duplicate variables for latitude, selecting lat_psi, and discarding lat_rho
07:41:32 WARNING opendrift.readers.reader_ROMS_native:282: Duplicate variables for latitude, selecting lat_psi, and discarding lat_u
07:41:32 WARNING opendrift.readers.reader_ROMS_native:282: Duplicate variables for latitude, selecting lat_psi, and discarding lat_v
07:41:32 WARNING opendrift.readers.reader_ROMS_native:282: Duplicate variables for longitude, selecting lon_psi, and discarding lon_rho
07:41:32 WARNING opendrift.readers.reader_ROMS_native:282: Duplicate variables for longitude, selecting lon_psi, and discarding lon_u
07:41:32 WARNING opendrift.readers.reader_ROMS_native:282: Duplicate variables for longitude, selecting lon_psi, and discarding lon_v
07:41:32 WARNING opendrift.readers.reader_ROMS_native:282: Duplicate variables for ocean_s_coordinate_g2, selecting s_rho, and discarding s_w
07:41:32 INFO    opendrift.readers.reader_ROMS_native:289: The following variables without standard_name are discarded: ['AICEnudass', 'Akk_bak', 'Akp_bak', 'Akt_bak', 'Akv_bak', 'Charnok_alpha', 'CrgBan_cw', 'Cs_r', 'Cs_w', 'FSobc_in', 'FSobc_out', 'Falpha', 'Fbeta', 'Fgamma', 'LtracerSrc', 'M2nudg', 'M2obc_in', 'M2obc_out', 'M3nudg', 'M3obc_in', 'M3obc_out', 'Tcline', 'Tnudg', 'Tobc_in', 'Tobc_out', 'Vstretching', 'Vtransform', 'Znudg', 'Zob', 'Zos', 'Zos_hsig_alpha', 'angle', 'dstart', 'dt', 'dtfast', 'el', 'f', 'gamma2', 'gls_Kmin', 'gls_Pmin', 'gls_c1', 'gls_c2', 'gls_c3m', 'gls_c3p', 'gls_cmu0', 'gls_m', 'gls_n', 'gls_p', 'gls_sigk', 'gls_sigp', 'hc', 'mask_u', 'mask_v', 'nAVG', 'nHIS', 'nRST', 'nSTA', 'ndefAVG', 'ndefHIS', 'ndtfast', 'ntimes', 'ntsAVG', 'pm', 'pn', 'rdrg', 'rdrg2', 'rho0', 'spherical', 'swrad', 'sz_alpha', 'theta_b', 'theta_s', 'ubar', 'vbar', 'xl', 'zeta_detided', 'ocean_time']
07:41:32 WARNING opendrift.readers.basereader.structured:50: No proj string or projection could be derived, using 'fakeproj'. This assumes that the variables are structured and gridded approximately equidistantly on the surface (i.e. in meters). This must be guaranteed by the user. You can get rid of this warning by supplying a valid projection to the reader.
07:41:32 INFO    opendrift.readers.basereader.structured:90: Making interpolator for lon,lat to x,y conversion...

First run, with default GSHHG vector landmask

o.add_reader([reader_nordic])
time = reader_nordic.start_time
o.set_config('general:use_auto_landmask', True)
o.seed_elements(lon, lat, radius=3000, number=1000, time=time)
o.run(end_time=reader_nordic.end_time, time_step=1800, time_step_output=3*3600)

07:41:32 INFO    opendrift.models.basemodel.environment:206: Adding a global landmask from GSHHG
07:41:38 INFO    opendrift.models.basemodel.environment:229: Fallback values will be used for the following variables which have no readers:
07:41:38 INFO    opendrift.models.basemodel.environment:232:    x_wind: 0.000000
07:41:38 INFO    opendrift.models.basemodel.environment:232:    y_wind: 0.000000
07:41:38 INFO    opendrift.models.basemodel.environment:232:    upward_sea_water_velocity: 0.000000
07:41:38 INFO    opendrift.models.basemodel.environment:232:    ocean_vertical_diffusivity: 0.000000
07:41:38 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_significant_height: 0.000000
07:41:38 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
07:41:38 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
07:41:38 INFO    opendrift.models.basemodel.environment:232:    ocean_mixed_layer_thickness: 50.000000
07:41:38 INFO    opendrift:1732: Skipping environment variable ocean_vertical_diffusivity because of condition ['drift:vertical_mixing', 'is', False]
07:41:38 INFO    opendrift:1732: Skipping environment variable ocean_mixed_layer_thickness because of condition ['drift:vertical_mixing', 'is', False]
07:41:38 INFO    opendrift:1743: Storing previous values of element property lon because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
07:41:38 INFO    opendrift:1743: Storing previous values of element property lat because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
07:41:38 INFO    opendrift:1751: Storing previous values of environment variable sea_surface_height because of condition ['drift:vertical_advection', 'is', True]
07:41:38 INFO    opendrift:899: Using existing reader for land_binary_mask
07:41:40 INFO    opendrift:930: Moving 3 out of 1000 points from land to water
07:41:40 INFO    opendrift:2035: 2016-02-02 12:00:00 - step 1 of 96 - 1000 active elements (0 deactivated)
07:41:40 INFO    opendrift.readers.reader_ROMS_native:339: Using mask_rho for mask_rho
07:41:40 INFO    opendrift.readers.reader_ROMS_native:360: Using mask_u for mask_u
07:41:40 INFO    opendrift.readers.reader_ROMS_native:381: Using mask_v for mask_v
07:41:40 INFO    opendrift.readers.reader_ROMS_native:404: Using angle from Dataset.
07:41:40 INFO    opendrift:2035: 2016-02-02 12:30:00 - step 2 of 96 - 1000 active elements (0 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 13:00:00 - step 3 of 96 - 996 active elements (4 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 13:30:00 - step 4 of 96 - 995 active elements (5 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 14:00:00 - step 5 of 96 - 995 active elements (5 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 14:30:00 - step 6 of 96 - 994 active elements (6 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 15:00:00 - step 7 of 96 - 994 active elements (6 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 15:30:00 - step 8 of 96 - 993 active elements (7 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 16:00:00 - step 9 of 96 - 992 active elements (8 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 16:30:00 - step 10 of 96 - 992 active elements (8 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 17:00:00 - step 11 of 96 - 992 active elements (8 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 17:30:00 - step 12 of 96 - 989 active elements (11 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 18:00:00 - step 13 of 96 - 988 active elements (12 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 18:30:00 - step 14 of 96 - 988 active elements (12 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 19:00:00 - step 15 of 96 - 988 active elements (12 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 19:30:00 - step 16 of 96 - 987 active elements (13 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 20:00:00 - step 17 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 20:30:00 - step 18 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 21:00:00 - step 19 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 21:30:00 - step 20 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 22:00:00 - step 21 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 22:30:00 - step 22 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 23:00:00 - step 23 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-02 23:30:00 - step 24 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 00:00:00 - step 25 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 00:30:00 - step 26 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 01:00:00 - step 27 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 01:30:00 - step 28 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 02:00:00 - step 29 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 02:30:00 - step 30 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 03:00:00 - step 31 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 03:30:00 - step 32 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 04:00:00 - step 33 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 04:30:00 - step 34 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 05:00:00 - step 35 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 05:30:00 - step 36 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 06:00:00 - step 37 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 06:30:00 - step 38 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 07:00:00 - step 39 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 07:30:00 - step 40 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 08:00:00 - step 41 of 96 - 986 active elements (14 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 08:30:00 - step 42 of 96 - 985 active elements (15 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 09:00:00 - step 43 of 96 - 985 active elements (15 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 09:30:00 - step 44 of 96 - 984 active elements (16 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 10:00:00 - step 45 of 96 - 984 active elements (16 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 10:30:00 - step 46 of 96 - 982 active elements (18 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 11:00:00 - step 47 of 96 - 982 active elements (18 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 11:30:00 - step 48 of 96 - 980 active elements (20 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 12:00:00 - step 49 of 96 - 979 active elements (21 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 12:30:00 - step 50 of 96 - 977 active elements (23 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 13:00:00 - step 51 of 96 - 976 active elements (24 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 13:30:00 - step 52 of 96 - 975 active elements (25 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 14:00:00 - step 53 of 96 - 975 active elements (25 deactivated)
07:41:42 INFO    opendrift:2035: 2016-02-03 14:30:00 - step 54 of 96 - 975 active elements (25 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 15:00:00 - step 55 of 96 - 974 active elements (26 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 15:30:00 - step 56 of 96 - 971 active elements (29 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 16:00:00 - step 57 of 96 - 969 active elements (31 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 16:30:00 - step 58 of 96 - 964 active elements (36 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 17:00:00 - step 59 of 96 - 958 active elements (42 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 17:30:00 - step 60 of 96 - 956 active elements (44 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 18:00:00 - step 61 of 96 - 953 active elements (47 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 18:30:00 - step 62 of 96 - 945 active elements (55 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 19:00:00 - step 63 of 96 - 942 active elements (58 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 19:30:00 - step 64 of 96 - 939 active elements (61 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 20:00:00 - step 65 of 96 - 935 active elements (65 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 20:30:00 - step 66 of 96 - 931 active elements (69 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 21:00:00 - step 67 of 96 - 930 active elements (70 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 21:30:00 - step 68 of 96 - 928 active elements (72 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 22:00:00 - step 69 of 96 - 924 active elements (76 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 22:30:00 - step 70 of 96 - 918 active elements (82 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 23:00:00 - step 71 of 96 - 913 active elements (87 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-03 23:30:00 - step 72 of 96 - 908 active elements (92 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-04 00:00:00 - step 73 of 96 - 899 active elements (101 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-04 00:30:00 - step 74 of 96 - 890 active elements (110 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-04 01:00:00 - step 75 of 96 - 882 active elements (118 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-04 01:30:00 - step 76 of 96 - 877 active elements (123 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-04 02:00:00 - step 77 of 96 - 870 active elements (130 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-04 02:30:00 - step 78 of 96 - 863 active elements (137 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-04 03:00:00 - step 79 of 96 - 859 active elements (141 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-04 03:30:00 - step 80 of 96 - 847 active elements (153 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-04 04:00:00 - step 81 of 96 - 845 active elements (155 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-04 04:30:00 - step 82 of 96 - 840 active elements (160 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-04 05:00:00 - step 83 of 96 - 833 active elements (167 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-04 05:30:00 - step 84 of 96 - 823 active elements (177 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-04 06:00:00 - step 85 of 96 - 810 active elements (190 deactivated)
07:41:43 INFO    opendrift:2035: 2016-02-04 06:30:00 - step 86 of 96 - 799 active elements (201 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-04 07:00:00 - step 87 of 96 - 789 active elements (211 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-04 07:30:00 - step 88 of 96 - 777 active elements (223 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-04 08:00:00 - step 89 of 96 - 772 active elements (228 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-04 08:30:00 - step 90 of 96 - 757 active elements (243 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-04 09:00:00 - step 91 of 96 - 751 active elements (249 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-04 09:30:00 - step 92 of 96 - 747 active elements (253 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-04 10:00:00 - step 93 of 96 - 740 active elements (260 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-04 10:30:00 - step 94 of 96 - 738 active elements (262 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-04 11:00:00 - step 95 of 96 - 731 active elements (269 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-04 11:30:00 - step 96 of 96 - 723 active elements (277 deactivated)

<xarray.Dataset> Size: 1MB
Dimensions:                                   (trajectory: 1000, time: 17)
Coordinates:
  * trajectory                                (trajectory) int64 8kB 0 1 ... 999
  * time                                      (time) datetime64[ns] 136B 2016...
Data variables: (12/21)
    status                                    (trajectory, time) float32 68kB ...
    moving                                    (trajectory, time) float32 68kB ...
    age_seconds                               (trajectory, time) float32 68kB ...
    origin_marker                             (trajectory, time) float32 68kB ...
    lon                                       (trajectory, time) float32 68kB ...
    lat                                       (trajectory, time) float32 68kB ...
    ...                                        ...
    upward_sea_water_velocity                 (trajectory, time) float32 68kB ...
    sea_surface_wave_significant_height       (trajectory, time) float32 68kB ...
    sea_surface_wave_stokes_drift_x_velocity  (trajectory, time) float32 68kB ...
    sea_surface_wave_stokes_drift_y_velocity  (trajectory, time) float32 68kB ...
    sea_floor_depth_below_sea_level           (trajectory, time) float32 68kB ...
    land_binary_mask                          (trajectory, time) float32 68kB ...
Attributes: (12/119)
    Conventions:                                                           CF...
    standard_name_vocabulary:                                              CF...
    featureType:                                                           tr...
    title:                                                                 Op...
    summary:                                                               Ou...
    keywords:                                                              tr...
    ...                                                                    ...
    geospatial_lon_units:                                                  de...
    geospatial_lon_resolution:                                             point
    runtime:                                                               0:...
    geospatial_vertical_min:                                               0.0
    geospatial_vertical_max:                                               0.0
    geospatial_vertical_positive:                                          up

xarray.Dataset

• Dimensions:
  • trajectory: 1000
  • time: 17
• 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], shape=(1000,))

  • time
    (time)
    datetime64[ns]
    2016-02-02T12:00:00 ... 2016-02-...

    standard_name :

    time

    long_name :

    time

    dtype :

    <class 'numpy.float64'>

    array(['2016-02-02T12:00:00.000000000', '2016-02-02T15:00:00.000000000',
           '2016-02-02T18:00:00.000000000', '2016-02-02T21:00:00.000000000',
           '2016-02-03T00:00:00.000000000', '2016-02-03T03:00:00.000000000',
           '2016-02-03T06:00:00.000000000', '2016-02-03T09:00:00.000000000',
           '2016-02-03T12:00:00.000000000', '2016-02-03T15:00:00.000000000',
           '2016-02-03T18:00:00.000000000', '2016-02-03T21:00:00.000000000',
           '2016-02-04T00:00:00.000000000', '2016-02-04T03:00:00.000000000',
           '2016-02-04T06:00:00.000000000', '2016-02-04T09:00:00.000000000',
           '2016-02-04T12:00:00.000000000'], dtype='datetime64[ns]')

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

    dtype :

    <class 'numpy.int32'>

    valid_range :

    [0 1]

    flag_values :

    [0 1]

    flag_meanings :

    active stranded

    array([[ 0.,  0.,  0., ...,  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., ...,  1., nan, nan]],
          shape=(1000, 17), dtype=float32)

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

    dtype :

    <class 'numpy.int32'>

    minval :

    0

    maxval :

    1

    array([[ 1.,  1.,  1., ...,  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., ...,  0., nan, nan]],
          shape=(1000, 17), dtype=float32)

  • age_seconds
    (trajectory, time)
    float32
    0.0 1.08e+04 2.16e+04 ... nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    s

    minval :

    0.0

    maxval :

    172800.0

    array([[     0.,  10800.,  21600., ..., 151200., 162000., 172800.],
           [     0.,  10800.,  21600., ..., 151200., 162000., 172800.],
           [     0.,  10800.,  21600., ..., 151200., 162000., 172800.],
           ...,
           [     0.,  10800.,  21600., ..., 151200., 162000., 172800.],
           [     0.,  10800.,  21600., ..., 151200., 162000., 172800.],
           [     0.,  10800.,  21600., ..., 145800.,     nan,     nan]],
          shape=(1000, 17), dtype=float32)

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

    dtype :

    <class 'numpy.int32'>

    unit :

    description :

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

    flag_values :

    [0]

    flag_meanings :

    Seed_0

    minval :

    0

    maxval :

    0

    array([[ 0.,  0.,  0., ...,  0.,  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., nan, nan]],
          shape=(1000, 17), dtype=float32)

  • lon
    (trajectory, time)
    float32
    14.88 14.9 14.93 ... 15.03 nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    axis :

    X

    minval :

    14.301589

    maxval :

    15.369776

    array([[14.877172 , 14.897906 , 14.926504 , ..., 15.220199 , 15.225651 ,
            15.223183 ],
           [14.778859 , 14.79981  , 14.833956 , ..., 15.251485 , 15.281697 ,
            15.309673 ],
           [14.8204775, 14.823901 , 14.831969 , ..., 15.189053 , 15.219408 ,
            15.249455 ],
           ...,
           [14.756787 , 14.771563 , 14.788212 , ..., 15.161574 , 15.191191 ,
            15.223537 ],
           [14.667431 , 14.676323 , 14.689936 , ..., 14.991064 , 15.017084 ,
            15.037321 ],
           [14.724239 , 14.726044 , 14.735568 , ..., 15.029565 ,        nan,
                   nan]], shape=(1000, 17), dtype=float32)

  • lat
    (trajectory, time)
    float32
    68.11 68.12 68.13 ... 68.22 nan nan

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    axis :

    Y

    minval :

    68.01945

    maxval :

    68.29827

    array([[68.11491 , 68.12478 , 68.13388 , ..., 68.161194, 68.15143 ,
            68.139275],
           [68.124   , 68.1395  , 68.15173 , ..., 68.218414, 68.21322 ,
            68.20647 ],
           [68.08862 , 68.100494, 68.11341 , ..., 68.21265 , 68.21133 ,
            68.20674 ],
           ...,
           [68.10426 , 68.11746 , 68.13246 , ..., 68.23475 , 68.23589 ,
            68.23565 ],
           [68.06927 , 68.07746 , 68.085815, ..., 68.22157 , 68.231766,
            68.240906],
           [68.064735, 68.076256, 68.0875  , ..., 68.22296 ,       nan,
                  nan]], shape=(1000, 17), dtype=float32)

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

    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., nan, nan]],
          shape=(1000, 17), dtype=float32)

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

    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],
           ...,
           [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]],
          shape=(1000, 17), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    1

    description :

    Elements are moved with this fraction of the current vector, in addition to currents and Stokes drift

    minval :

    1.0

    maxval :

    1.0

    array([[ 1.,  1.,  1., ...,  1.,  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., nan, nan]],
          shape=(1000, 17), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    m/s

    description :

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

    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., nan, nan]],
          shape=(1000, 17), dtype=float32)

  • x_sea_water_velocity
    (trajectory, time)
    float32
    0.06472 0.09985 0.1156 ... nan nan

    minval :

    -0.1452677

    maxval :

    0.19952533

    array([[ 0.06472018,  0.09985437,  0.11561183, ...,  0.03219877,
             0.00421871, -0.02390574],
           [ 0.0662813 ,  0.11068583,  0.15806125, ...,  0.11883222,
             0.11283366,  0.10125312],
           [ 0.00541813,  0.02397641,  0.04624955, ...,  0.11740991,
             0.11551356,  0.11522176],
           ...,
           [ 0.05207663,  0.06197634,  0.07827664, ...,  0.10772207,
             0.12265687,  0.12121784],
           [ 0.0252429 ,  0.04573552,  0.06140871, ...,  0.11053084,
             0.08583724,  0.06984734],
           [-0.0058907 ,  0.02566006,  0.050713  , ...,  0.09346184,
                    nan,         nan]], shape=(1000, 17), dtype=float32)

  • y_sea_water_velocity
    (trajectory, time)
    float32
    0.1046 0.0975 0.09143 ... nan nan

    minval :

    -0.15142582

    maxval :

    0.17843722

    array([[ 0.10463738,  0.09750166,  0.09143377, ..., -0.09034176,
            -0.11543446, -0.13537003],
           [ 0.16189604,  0.14258766,  0.10423066, ..., -0.04678254,
            -0.06395643, -0.07522707],
           [ 0.11846498,  0.12889196,  0.13551562, ..., -0.00109887,
            -0.03224672, -0.06321396],
           ...,
           [ 0.13096647,  0.14587188,  0.15556085, ...,  0.01838147,
             0.00224334, -0.00822518],
           [ 0.08452237,  0.08524593,  0.08806531, ...,  0.11010524,
             0.0987203 ,  0.09017027],
           [ 0.11773563,  0.1173166 ,  0.11494715, ...,  0.09777828,
                    nan,         nan]], shape=(1000, 17), dtype=float32)

  • sea_surface_height
    (trajectory, time)
    float32
    0.4329 0.4166 0.4003 ... nan nan

    minval :

    0.15286748

    maxval :

    0.43597323

    array([[0.43292335, 0.4165746 , 0.40026605, ..., 0.20985389, 0.19413781,
            0.18113978],
           [0.42671713, 0.40887395, 0.3922047 , ..., 0.20401895, 0.18931794,
            0.1761517 ],
           [0.43128845, 0.414747  , 0.3976727 , ..., 0.20556459, 0.19034782,
            0.17798916],
           ...,
           [0.4267774 , 0.40993518, 0.39258167, ..., 0.20152742, 0.18653151,
            0.1737539 ],
           [0.42316896, 0.4068736 , 0.39046595, ..., 0.19807771, 0.1820827 ,
            0.16869767],
           [0.4257891 , 0.4094062 , 0.3928822 , ..., 0.20668049,        nan,
                   nan]], shape=(1000, 17), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  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., nan, nan]],
          shape=(1000, 17), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  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., nan, nan]],
          shape=(1000, 17), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  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., nan, nan]],
          shape=(1000, 17), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  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., nan, nan]],
          shape=(1000, 17), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  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., nan, nan]],
          shape=(1000, 17), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  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., nan, nan]],
          shape=(1000, 17), dtype=float32)

  • sea_floor_depth_below_sea_level
    (trajectory, time)
    float32
    300.2 295.4 293.5 ... 163.4 nan nan

    minval :

    82.55389

    maxval :

    374.22107

    array([[300.2193 , 295.36343, 293.53345, ..., 341.46396, 352.88043,
            363.92258],
           [244.83344, 227.90446, 221.28603, ..., 237.8977 , 252.32233,
            258.5296 ],
           [309.48007, 299.2717 , 285.0905 , ..., 240.24332, 249.39017,
            261.99194],
           ...,
           [265.4019 , 251.81877, 235.02663, ..., 185.85023, 189.51863,
            194.08647],
           [276.11667, 269.16086, 263.58356, ..., 165.85138, 151.94714,
            140.25618],
           [302.48672, 291.23718, 280.44785, ..., 163.43324,       nan,
                  nan]], shape=(1000, 17), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

    array([[ 0.,  0.,  0., ...,  0.,  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., nan, nan]],
          shape=(1000, 17), 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(['2016-02-02 12:00:00', '2016-02-02 15:00:00',
                   '2016-02-02 18:00:00', '2016-02-02 21:00:00',
                   '2016-02-03 00:00:00', '2016-02-03 03:00:00',
                   '2016-02-03 06:00:00', '2016-02-03 09:00:00',
                   '2016-02-03 12:00:00', '2016-02-03 15:00:00',
                   '2016-02-03 18:00:00', '2016-02-03 21:00:00',
                   '2016-02-04 00:00:00', '2016-02-04 03:00:00',
                   '2016-02-04 06:00:00', '2016-02-04 09:00:00',
                   '2016-02-04 12:00:00'],
                  dtype='datetime64[ns]', name='time', freq='3h'))

• Attributes: (119)

  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-07-17 07:41:38.307829

  date_created :

  2025-07-17T07:41:38.307839

  source :

  Output from simulation with OpenDrift

  model_url :

  https://github.com/OpenDrift/opendrift

  opendrift_class :

  OceanDrift

  opendrift_module :

  opendrift.models.oceandrift

  readers :

  odict_keys(['roms native', 'global_landmask'])

  time_coverage_start :

  2016-02-02 12:00:00

  time_step_calculation :

  0:30:00

  time_step_output :

  3:00: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:x_wind :

  None

  config_environment:fallback:x_wind :

  0

  config_environment:constant:y_wind :

  None

  config_environment:fallback:y_wind :

  0

  config_environment:constant:upward_sea_water_velocity :

  None

  config_environment:fallback:upward_sea_water_velocity :

  0

  config_environment:constant:ocean_vertical_diffusivity :

  None

  config_environment:fallback:ocean_vertical_diffusivity :

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

  10000

  config_environment:constant:land_binary_mask :

  None

  config_environment:fallback:land_binary_mask :

  None

  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 :

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

  0.02

  config_seed:current_drift_factor :

  1

  config_seed:terminal_velocity :

  0.0

  config_drift:vertical_advection :

  True

  config_drift:vertical_advection_at_surface :

  False

  config_drift:water_column_stretching :

  False

  config_drift:vertical_advection_correction :

  False

  config_drift:vertical_mixing :

  False

  config_drift:vertical_mixing_at_surface :

  False

  config_vertical_mixing:timestep :

  60

  config_vertical_mixing:diffusivitymodel :

  environment

  config_vertical_mixing:background_diffusivity :

  1.2e-05

  config_vertical_mixing:TSprofiles :

  False

  config_drift:wind_drift_depth :

  0.1

  config_drift:stokes_drift :

  True

  config_drift:stokes_drift_profile :

  Phillips

  config_drift:use_tabularised_stokes_drift :

  False

  config_drift:tabularised_stokes_drift_fetch :

  25000

  config_general:seafloor_action :

  lift_to_seafloor

  config_drift:truncate_ocean_model_below_m :

  None

  config_seed:seafloor :

  False

  opendrift_version :

  1.14.2

  seed_geojson :

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

  simulation_time :

  2025-07-17 07:41:40.121847

  reader_x_sea_water_velocity :

  ['roms native']

  reader_y_sea_water_velocity :

  ['roms native']

  reader_sea_surface_height :

  ['roms native']

  reader_x_wind :

  0

  reader_y_wind :

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

  ['roms native']

  reader_land_binary_mask :

  ['global_landmask']

  time_coverage_end :

  2016-02-04 12:00:00

  time_coverage_duration :

  P2DT0H0M0S

  time_coverage_resolution :

  P0DT0H30M0S

  performance :

  -------------------- Reader performance: -------------------- roms native 0:00:00.3 total 0:00:00.0 preparing 0:00:00.1 reading 0:00:00.0 interpolation 0:00:00.0 interpolation_time 0:00:00.1 rotating vectors 0:00:00.0 masking -------------------- global_landmask 0:00:01.8 total 0:00:00.0 preparing 0:00:01.8 reading 0:00:00.0 masking -------------------- Performance: 12.1 total time 6.0 configuration 1.8 preparing main loop 1.8 moving elements to ocean 4.2 main loop 0.0 updating elements 0.0 cleaning up --------------------

  geospatial_bounds_crs :

  EPSG:4326

  geospatial_bounds_vertical_crs :

  EPSG:5831

  geospatial_lat_min :

  68.01945

  geospatial_lat_max :

  68.29827

  geospatial_lat_units :

  degrees_north

  geospatial_lat_resolution :

  point

  geospatial_lon_min :

  14.301589

  geospatial_lon_max :

  15.369776

  geospatial_lon_units :

  degrees_east

  geospatial_lon_resolution :

  point

  runtime :

  0:00:12.126664

  geospatial_vertical_min :

  0.0

  geospatial_vertical_max :

  0.0

  geospatial_vertical_positive :

  up

Second run, with landmask from ocean model

o2 = OceanDrift(loglevel=20)
o2.add_reader([reader_nordic])
lon = 14.75; lat = 68.1
o2.set_config('general:use_auto_landmask', False)
o2.seed_elements(lon, lat, radius=3000, number=1000, time=time)
o2.run(end_time=reader_nordic.end_time, time_step=1800, time_step_output=3*3600)

07:41:44 INFO    opendrift:513: OpenDriftSimulation initialised (version 1.14.2 / v1.14.2-90-g1dd1995)
07:41:44 INFO    opendrift.models.basemodel.environment:229: Fallback values will be used for the following variables which have no readers:
07:41:44 INFO    opendrift.models.basemodel.environment:232:    x_wind: 0.000000
07:41:44 INFO    opendrift.models.basemodel.environment:232:    y_wind: 0.000000
07:41:44 INFO    opendrift.models.basemodel.environment:232:    upward_sea_water_velocity: 0.000000
07:41:44 INFO    opendrift.models.basemodel.environment:232:    ocean_vertical_diffusivity: 0.000000
07:41:44 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_significant_height: 0.000000
07:41:44 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
07:41:44 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
07:41:44 INFO    opendrift.models.basemodel.environment:232:    ocean_mixed_layer_thickness: 50.000000
07:41:44 INFO    opendrift:1732: Skipping environment variable ocean_vertical_diffusivity because of condition ['drift:vertical_mixing', 'is', False]
07:41:44 INFO    opendrift:1732: Skipping environment variable ocean_mixed_layer_thickness because of condition ['drift:vertical_mixing', 'is', False]
07:41:44 INFO    opendrift:1743: Storing previous values of element property lon because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
07:41:44 INFO    opendrift:1743: Storing previous values of element property lat because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
07:41:44 INFO    opendrift:1751: Storing previous values of environment variable sea_surface_height because of condition ['drift:vertical_advection', 'is', True]
07:41:44 INFO    opendrift:899: Using existing reader for land_binary_mask
07:41:44 INFO    opendrift:928: All points are in ocean
07:41:44 INFO    opendrift:2035: 2016-02-02 12:00:00 - step 1 of 96 - 1000 active elements (0 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 12:30:00 - step 2 of 96 - 1000 active elements (0 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 13:00:00 - step 3 of 96 - 1000 active elements (0 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 13:30:00 - step 4 of 96 - 999 active elements (1 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 14:00:00 - step 5 of 96 - 998 active elements (2 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 14:30:00 - step 6 of 96 - 998 active elements (2 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 15:00:00 - step 7 of 96 - 997 active elements (3 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 15:30:00 - step 8 of 96 - 995 active elements (5 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 16:00:00 - step 9 of 96 - 992 active elements (8 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 16:30:00 - step 10 of 96 - 991 active elements (9 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 17:00:00 - step 11 of 96 - 989 active elements (11 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 17:30:00 - step 12 of 96 - 986 active elements (14 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 18:00:00 - step 13 of 96 - 981 active elements (19 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 18:30:00 - step 14 of 96 - 975 active elements (25 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 19:00:00 - step 15 of 96 - 970 active elements (30 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 19:30:00 - step 16 of 96 - 964 active elements (36 deactivated)
07:41:44 INFO    opendrift:2035: 2016-02-02 20:00:00 - step 17 of 96 - 957 active elements (43 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-02 20:30:00 - step 18 of 96 - 953 active elements (47 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-02 21:00:00 - step 19 of 96 - 949 active elements (51 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-02 21:30:00 - step 20 of 96 - 947 active elements (53 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-02 22:00:00 - step 21 of 96 - 944 active elements (56 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-02 22:30:00 - step 22 of 96 - 933 active elements (67 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-02 23:00:00 - step 23 of 96 - 926 active elements (74 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-02 23:30:00 - step 24 of 96 - 920 active elements (80 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 00:00:00 - step 25 of 96 - 912 active elements (88 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 00:30:00 - step 26 of 96 - 906 active elements (94 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 01:00:00 - step 27 of 96 - 894 active elements (106 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 01:30:00 - step 28 of 96 - 889 active elements (111 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 02:00:00 - step 29 of 96 - 881 active elements (119 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 02:30:00 - step 30 of 96 - 875 active elements (125 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 03:00:00 - step 31 of 96 - 864 active elements (136 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 03:30:00 - step 32 of 96 - 856 active elements (144 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 04:00:00 - step 33 of 96 - 848 active elements (152 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 04:30:00 - step 34 of 96 - 843 active elements (157 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 05:00:00 - step 35 of 96 - 836 active elements (164 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 05:30:00 - step 36 of 96 - 829 active elements (171 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 06:00:00 - step 37 of 96 - 826 active elements (174 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 06:30:00 - step 38 of 96 - 818 active elements (182 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 07:00:00 - step 39 of 96 - 815 active elements (185 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 07:30:00 - step 40 of 96 - 814 active elements (186 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 08:00:00 - step 41 of 96 - 809 active elements (191 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 08:30:00 - step 42 of 96 - 806 active elements (194 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 09:00:00 - step 43 of 96 - 797 active elements (203 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 09:30:00 - step 44 of 96 - 792 active elements (208 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 10:00:00 - step 45 of 96 - 786 active elements (214 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 10:30:00 - step 46 of 96 - 780 active elements (220 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 11:00:00 - step 47 of 96 - 774 active elements (226 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 11:30:00 - step 48 of 96 - 771 active elements (229 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 12:00:00 - step 49 of 96 - 765 active elements (235 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 12:30:00 - step 50 of 96 - 762 active elements (238 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 13:00:00 - step 51 of 96 - 756 active elements (244 deactivated)
07:41:45 INFO    opendrift:2035: 2016-02-03 13:30:00 - step 52 of 96 - 752 active elements (248 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 14:00:00 - step 53 of 96 - 751 active elements (249 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 14:30:00 - step 54 of 96 - 749 active elements (251 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 15:00:00 - step 55 of 96 - 747 active elements (253 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 15:30:00 - step 56 of 96 - 743 active elements (257 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 16:00:00 - step 57 of 96 - 742 active elements (258 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 16:30:00 - step 58 of 96 - 741 active elements (259 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 17:00:00 - step 59 of 96 - 739 active elements (261 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 17:30:00 - step 60 of 96 - 738 active elements (262 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 18:00:00 - step 61 of 96 - 732 active elements (268 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 18:30:00 - step 62 of 96 - 732 active elements (268 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 19:00:00 - step 63 of 96 - 732 active elements (268 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 19:30:00 - step 64 of 96 - 730 active elements (270 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 20:00:00 - step 65 of 96 - 730 active elements (270 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 20:30:00 - step 66 of 96 - 729 active elements (271 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 21:00:00 - step 67 of 96 - 728 active elements (272 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 21:30:00 - step 68 of 96 - 726 active elements (274 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 22:00:00 - step 69 of 96 - 724 active elements (276 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 22:30:00 - step 70 of 96 - 724 active elements (276 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 23:00:00 - step 71 of 96 - 724 active elements (276 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-03 23:30:00 - step 72 of 96 - 724 active elements (276 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 00:00:00 - step 73 of 96 - 724 active elements (276 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 00:30:00 - step 74 of 96 - 723 active elements (277 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 01:00:00 - step 75 of 96 - 721 active elements (279 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 01:30:00 - step 76 of 96 - 721 active elements (279 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 02:00:00 - step 77 of 96 - 719 active elements (281 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 02:30:00 - step 78 of 96 - 719 active elements (281 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 03:00:00 - step 79 of 96 - 719 active elements (281 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 03:30:00 - step 80 of 96 - 718 active elements (282 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 04:00:00 - step 81 of 96 - 718 active elements (282 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 04:30:00 - step 82 of 96 - 718 active elements (282 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 05:00:00 - step 83 of 96 - 717 active elements (283 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 05:30:00 - step 84 of 96 - 717 active elements (283 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 06:00:00 - step 85 of 96 - 716 active elements (284 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 06:30:00 - step 86 of 96 - 716 active elements (284 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 07:00:00 - step 87 of 96 - 716 active elements (284 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 07:30:00 - step 88 of 96 - 715 active elements (285 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 08:00:00 - step 89 of 96 - 715 active elements (285 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 08:30:00 - step 90 of 96 - 715 active elements (285 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 09:00:00 - step 91 of 96 - 715 active elements (285 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 09:30:00 - step 92 of 96 - 715 active elements (285 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 10:00:00 - step 93 of 96 - 715 active elements (285 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 10:30:00 - step 94 of 96 - 714 active elements (286 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 11:00:00 - step 95 of 96 - 714 active elements (286 deactivated)
07:41:46 INFO    opendrift:2035: 2016-02-04 11:30:00 - step 96 of 96 - 714 active elements (286 deactivated)

<xarray.Dataset> Size: 1MB
Dimensions:                                   (trajectory: 1000, time: 17)
Coordinates:
  * trajectory                                (trajectory) int64 8kB 0 1 ... 999
  * time                                      (time) datetime64[ns] 136B 2016...
Data variables: (12/21)
    status                                    (trajectory, time) float32 68kB ...
    moving                                    (trajectory, time) float32 68kB ...
    age_seconds                               (trajectory, time) float32 68kB ...
    origin_marker                             (trajectory, time) float32 68kB ...
    lon                                       (trajectory, time) float32 68kB ...
    lat                                       (trajectory, time) float32 68kB ...
    ...                                        ...
    upward_sea_water_velocity                 (trajectory, time) float32 68kB ...
    sea_surface_wave_significant_height       (trajectory, time) float32 68kB ...
    sea_surface_wave_stokes_drift_x_velocity  (trajectory, time) float32 68kB ...
    sea_surface_wave_stokes_drift_y_velocity  (trajectory, time) float32 68kB ...
    sea_floor_depth_below_sea_level           (trajectory, time) float32 68kB ...
    land_binary_mask                          (trajectory, time) float32 68kB ...
Attributes: (12/119)
    Conventions:                                                           CF...
    standard_name_vocabulary:                                              CF...
    featureType:                                                           tr...
    title:                                                                 Op...
    summary:                                                               Ou...
    keywords:                                                              tr...
    ...                                                                    ...
    geospatial_lon_units:                                                  de...
    geospatial_lon_resolution:                                             point
    runtime:                                                               0:...
    geospatial_vertical_min:                                               0.0
    geospatial_vertical_max:                                               0.0
    geospatial_vertical_positive:                                          up

xarray.Dataset

• Dimensions:
  • trajectory: 1000
  • time: 17
• 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], shape=(1000,))

  • time
    (time)
    datetime64[ns]
    2016-02-02T12:00:00 ... 2016-02-...

    standard_name :

    time

    long_name :

    time

    dtype :

    <class 'numpy.float64'>

    array(['2016-02-02T12:00:00.000000000', '2016-02-02T15:00:00.000000000',
           '2016-02-02T18:00:00.000000000', '2016-02-02T21:00:00.000000000',
           '2016-02-03T00:00:00.000000000', '2016-02-03T03:00:00.000000000',
           '2016-02-03T06:00:00.000000000', '2016-02-03T09:00:00.000000000',
           '2016-02-03T12:00:00.000000000', '2016-02-03T15:00:00.000000000',
           '2016-02-03T18:00:00.000000000', '2016-02-03T21:00:00.000000000',
           '2016-02-04T00:00:00.000000000', '2016-02-04T03:00:00.000000000',
           '2016-02-04T06:00:00.000000000', '2016-02-04T09:00:00.000000000',
           '2016-02-04T12:00:00.000000000'], dtype='datetime64[ns]')

• Data variables: (21)
  • 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 1]

    flag_values :

    [0 1]

    flag_meanings :

    active stranded

    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.]], shape=(1000, 17), 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 :

    0

    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.]], shape=(1000, 17), dtype=float32)

  • age_seconds
    (trajectory, time)
    float32
    0.0 1.08e+04 ... 1.62e+05 1.728e+05

    dtype :

    <class 'numpy.float32'>

    units :

    s

    minval :

    0.0

    maxval :

    172800.0

    array([[     0.,  10800.,  21600., ..., 151200., 162000., 172800.],
           [     0.,  10800.,  21600., ..., 151200., 162000., 172800.],
           [     0.,  10800.,  21600., ..., 151200., 162000., 172800.],
           ...,
           [     0.,  10800.,  21600., ..., 151200., 162000., 172800.],
           [     0.,  10800.,  21600., ..., 151200., 162000., 172800.],
           [     0.,  10800.,  21600., ..., 151200., 162000., 172800.]],
          shape=(1000, 17), 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.]], shape=(1000, 17), dtype=float32)

  • lon
    (trajectory, time)
    float32
    14.88 14.9 14.93 ... 15.06 15.08

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_east

    standard_name :

    longitude

    long_name :

    longitude

    axis :

    X

    minval :

    14.301589

    maxval :

    15.369776

    array([[14.877172 , 14.897906 , 14.926504 , ..., 15.220199 , 15.225651 ,
            15.223183 ],
           [14.778859 , 14.79981  , 14.833956 , ..., 15.251485 , 15.281697 ,
            15.309673 ],
           [14.8204775, 14.823901 , 14.831969 , ..., 15.189053 , 15.219408 ,
            15.249455 ],
           ...,
           [14.756787 , 14.771563 , 14.788212 , ..., 15.161574 , 15.191191 ,
            15.223537 ],
           [14.667431 , 14.676323 , 14.689936 , ..., 14.991064 , 15.017084 ,
            15.037321 ],
           [14.724239 , 14.726044 , 14.735568 , ..., 15.044024 , 15.062516 ,
            15.078517 ]], shape=(1000, 17), dtype=float32)

  • lat
    (trajectory, time)
    float32
    68.11 68.12 68.13 ... 68.24 68.25

    dtype :

    <class 'numpy.float32'>

    units :

    degrees_north

    standard_name :

    latitude

    long_name :

    latitude

    axis :

    Y

    minval :

    68.01945

    maxval :

    68.2609

    array([[68.11491 , 68.12478 , 68.13388 , ..., 68.161194, 68.15143 ,
            68.139275],
           [68.124   , 68.1395  , 68.15173 , ..., 68.218414, 68.21322 ,
            68.20647 ],
           [68.08862 , 68.100494, 68.11341 , ..., 68.21265 , 68.21133 ,
            68.20674 ],
           ...,
           [68.10426 , 68.11746 , 68.13246 , ..., 68.23475 , 68.23589 ,
            68.23565 ],
           [68.06927 , 68.07746 , 68.085815, ..., 68.22157 , 68.231766,
            68.240906],
           [68.064735, 68.076256, 68.0875  , ..., 68.22921 , 68.23775 ,
            68.2451  ]], shape=(1000, 17), 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.]], shape=(1000, 17), dtype=float32)

  • wind_drift_factor
    (trajectory, time)
    float32
    0.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],
           ...,
           [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]],
          shape=(1000, 17), dtype=float32)

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

    dtype :

    <class 'numpy.float32'>

    units :

    1

    description :

    Elements are moved with this fraction of the current vector, in addition to currents and Stokes drift

    minval :

    1.0

    maxval :

    1.0

    array([[1., 1., 1., ..., 1., 1., 1.],
           [1., 1., 1., ..., 1., 1., 1.],
           [1., 1., 1., ..., 1., 1., 1.],
           ...,
           [1., 1., 1., ..., 1., 1., 1.],
           [1., 1., 1., ..., 1., 1., 1.],
           [1., 1., 1., ..., 1., 1., 1.]], shape=(1000, 17), dtype=float32)

  • terminal_velocity
    (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/s

    description :

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

    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.]], shape=(1000, 17), dtype=float32)

  • x_sea_water_velocity
    (trajectory, time)
    float32
    0.06472 0.09985 ... 0.06142 0.0615

    minval :

    -0.1452677

    maxval :

    0.20277517

    array([[ 0.06472018,  0.09985437,  0.11561183, ...,  0.03219877,
             0.00421871, -0.02390574],
           [ 0.0662813 ,  0.11068583,  0.15806125, ...,  0.11883222,
             0.11283366,  0.10125312],
           [ 0.00541813,  0.02397641,  0.04624955, ...,  0.11740991,
             0.11551356,  0.11522176],
           ...,
           [ 0.05207663,  0.06197634,  0.07827664, ...,  0.10772207,
             0.12265687,  0.12121784],
           [ 0.0252429 ,  0.04573552,  0.06140871, ...,  0.11053084,
             0.08583724,  0.06984734],
           [-0.0058907 ,  0.02566006,  0.050713  , ...,  0.08062977,
             0.06142236,  0.06150147]], shape=(1000, 17), dtype=float32)

  • y_sea_water_velocity
    (trajectory, time)
    float32
    0.1046 0.0975 ... 0.08137 0.07108

    minval :

    -0.15142582

    maxval :

    0.17843722

    array([[ 0.10463738,  0.09750166,  0.09143377, ..., -0.09034176,
            -0.11543446, -0.13537003],
           [ 0.16189604,  0.14258766,  0.10423066, ..., -0.04678254,
            -0.06395643, -0.07522707],
           [ 0.11846498,  0.12889196,  0.13551562, ..., -0.00109887,
            -0.03224672, -0.06321396],
           ...,
           [ 0.13096647,  0.14587188,  0.15556085, ...,  0.01838147,
             0.00224334, -0.00822518],
           [ 0.08452237,  0.08524593,  0.08806531, ...,  0.11010524,
             0.0987203 ,  0.09017027],
           [ 0.11773563,  0.1173166 ,  0.11494715, ...,  0.09242418,
             0.08137171,  0.07108025]], shape=(1000, 17), dtype=float32)

  • sea_surface_height
    (trajectory, time)
    float32
    0.4329 0.4166 ... 0.1829 0.1694

    minval :

    0.15286748

    maxval :

    0.43597323

    array([[0.43292335, 0.4165746 , 0.40026605, ..., 0.20985389, 0.19413781,
            0.18113978],
           [0.42671713, 0.40887395, 0.3922047 , ..., 0.20401895, 0.18931794,
            0.1761517 ],
           [0.43128845, 0.414747  , 0.3976727 , ..., 0.20556459, 0.19034782,
            0.17798916],
           ...,
           [0.4267774 , 0.40993518, 0.39258167, ..., 0.20152742, 0.18653151,
            0.1737539 ],
           [0.42316896, 0.4068736 , 0.39046595, ..., 0.19807771, 0.1820827 ,
            0.16869767],
           [0.4257891 , 0.4094062 , 0.3928822 , ..., 0.1988137 , 0.18287966,
            0.1694215 ]], shape=(1000, 17), dtype=float32)

  • x_wind
    (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.]], shape=(1000, 17), dtype=float32)

  • y_wind
    (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.]], shape=(1000, 17), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

    array([[0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.],
           ...,
           [0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.]], shape=(1000, 17), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

    array([[0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.],
           ...,
           [0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.]], shape=(1000, 17), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

    array([[0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.],
           ...,
           [0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.]], shape=(1000, 17), dtype=float32)

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

    minval :

    0.0

    maxval :

    0.0

    array([[0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.],
           ...,
           [0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.],
           [0., 0., 0., ..., 0., 0., 0.]], shape=(1000, 17), dtype=float32)

  • sea_floor_depth_below_sea_level
    (trajectory, time)
    float32
    300.2 295.4 293.5 ... 149.0 142.4

    minval :

    91.9484

    maxval :

    374.22107

    array([[300.2193 , 295.36343, 293.53345, ..., 341.46396, 352.88043,
            363.92258],
           [244.83344, 227.90446, 221.28603, ..., 237.8977 , 252.32233,
            258.5296 ],
           [309.48007, 299.2717 , 285.0905 , ..., 240.24332, 249.39017,
            261.99194],
           ...,
           [265.4019 , 251.81877, 235.02663, ..., 185.85023, 189.51863,
            194.08647],
           [276.11667, 269.16086, 263.58356, ..., 165.85138, 151.94714,
            140.25618],
           [302.48672, 291.23718, 280.44785, ..., 158.77815, 149.03651,
            142.36003]], shape=(1000, 17), 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.5

    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.]], shape=(1000, 17), 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(['2016-02-02 12:00:00', '2016-02-02 15:00:00',
                   '2016-02-02 18:00:00', '2016-02-02 21:00:00',
                   '2016-02-03 00:00:00', '2016-02-03 03:00:00',
                   '2016-02-03 06:00:00', '2016-02-03 09:00:00',
                   '2016-02-03 12:00:00', '2016-02-03 15:00:00',
                   '2016-02-03 18:00:00', '2016-02-03 21:00:00',
                   '2016-02-04 00:00:00', '2016-02-04 03:00:00',
                   '2016-02-04 06:00:00', '2016-02-04 09:00:00',
                   '2016-02-04 12:00:00'],
                  dtype='datetime64[ns]', name='time', freq='3h'))

• Attributes: (119)

  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-07-17 07:41:44.504618

  date_created :

  2025-07-17T07:41:44.504626

  source :

  Output from simulation with OpenDrift

  model_url :

  https://github.com/OpenDrift/opendrift

  opendrift_class :

  OceanDrift

  opendrift_module :

  opendrift.models.oceandrift

  readers :

  odict_keys(['roms native'])

  time_coverage_start :

  2016-02-02 12:00:00

  time_step_calculation :

  0:30:00

  time_step_output :

  3:00: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:x_wind :

  None

  config_environment:fallback:x_wind :

  0

  config_environment:constant:y_wind :

  None

  config_environment:fallback:y_wind :

  0

  config_environment:constant:upward_sea_water_velocity :

  None

  config_environment:fallback:upward_sea_water_velocity :

  0

  config_environment:constant:ocean_vertical_diffusivity :

  None

  config_environment:fallback:ocean_vertical_diffusivity :

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

  10000

  config_environment:constant:land_binary_mask :

  None

  config_environment:fallback:land_binary_mask :

  None

  config_general:use_auto_landmask :

  False

  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 :

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

  0.02

  config_seed:current_drift_factor :

  1

  config_seed:terminal_velocity :

  0.0

  config_drift:vertical_advection :

  True

  config_drift:vertical_advection_at_surface :

  False

  config_drift:water_column_stretching :

  False

  config_drift:vertical_advection_correction :

  False

  config_drift:vertical_mixing :

  False

  config_drift:vertical_mixing_at_surface :

  False

  config_vertical_mixing:timestep :

  60

  config_vertical_mixing:diffusivitymodel :

  environment

  config_vertical_mixing:background_diffusivity :

  1.2e-05

  config_vertical_mixing:TSprofiles :

  False

  config_drift:wind_drift_depth :

  0.1

  config_drift:stokes_drift :

  True

  config_drift:stokes_drift_profile :

  Phillips

  config_drift:use_tabularised_stokes_drift :

  False

  config_drift:tabularised_stokes_drift_fetch :

  25000

  config_general:seafloor_action :

  lift_to_seafloor

  config_drift:truncate_ocean_model_below_m :

  None

  config_seed:seafloor :

  False

  opendrift_version :

  1.14.2

  seed_geojson :

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

  simulation_time :

  2025-07-17 07:41:44.511159

  reader_x_sea_water_velocity :

  ['roms native']

  reader_y_sea_water_velocity :

  ['roms native']

  reader_sea_surface_height :

  ['roms native']

  reader_x_wind :

  0

  reader_y_wind :

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

  ['roms native']

  reader_land_binary_mask :

  ['roms native']

  time_coverage_end :

  2016-02-04 12:00:00

  time_coverage_duration :

  P2DT0H0M0S

  time_coverage_resolution :

  P0DT0H30M0S

  performance :

  -------------------- Reader performance: -------------------- roms native 0:00:00.6 total 0:00:00.0 preparing 0:00:00.2 reading 0:00:00.1 interpolation 0:00:00.0 interpolation_time 0:00:00.3 rotating vectors 0:00:00.0 masking -------------------- Performance: 2.3 total time 0.0 configuration 0.0 preparing main loop 0.0 moving elements to ocean 2.3 main loop 0.0 updating elements 0.0 cleaning up --------------------

  geospatial_bounds_crs :

  EPSG:4326

  geospatial_bounds_vertical_crs :

  EPSG:5831

  geospatial_lat_min :

  68.01945

  geospatial_lat_max :

  68.2609

  geospatial_lat_units :

  degrees_north

  geospatial_lat_resolution :

  point

  geospatial_lon_min :

  14.301589

  geospatial_lon_max :

  15.369776

  geospatial_lon_units :

  degrees_east

  geospatial_lon_resolution :

  point

  runtime :

  0:00:02.351429

  geospatial_vertical_min :

  0.0

  geospatial_vertical_max :

  0.0

  geospatial_vertical_positive :

  up

from matplotlib.colors import ListedColormap
import cartopy.feature as cfeature
cmap = ListedColormap([cfeature.COLORS['water'],
                       cfeature.COLORS['land']])

To only show the landmask from the model, hide the coastline landmask by doing:

o2.plot(background='land_binary_mask', hide_landmask=True, cmap=cmap)

(<GeoAxes: title={'center': 'OpenDrift - OceanDrift\n2016-02-02 12:00 to 2016-02-04 12:00 UTC (17 steps)'}>, <Figure size 1100x1014.37 with 2 Axes>)

Animation illustrating that red particles strand at ocean model land cells, and black particles strand at GSHHG land polygons

o.animation(compare=o2, background='land_binary_mask', cmap=cmap,
            legend=['Default GSHHG landmask', 'Ocean model landmask'])

07:42:12 INFO    opendrift:4553: Saving animation to /root/project/docs/source/gallery/animations/example_model_landmask_0.gif...
07:42:27 INFO    opendrift:2994: Time to make animation: 0:00:32.559208

o.plot(compare=o2, background='land_binary_mask', cmap=cmap,
       legend=['Default GSHHG landmask', 'Ocean model landmask'])

(<GeoAxes: title={'center': 'OpenDrift - OceanDrift\n2016-02-02 12:00 to 2016-02-04 12:00 UTC (17 steps)'}>, <Figure size 1100x1074.35 with 2 Axes>)