ROMS native reader

import numpy as np
from opendrift.readers import reader_ROMS_native
from opendrift.models.oceandrift import OceanDrift

o = OceanDrift(loglevel=20)  # Set loglevel to 0 for debug information

17:12:30 INFO    opendrift:509: OpenDriftSimulation initialised (version 1.13.0 / v1.13.0-15-g93fccf7)

Creating and adding reader for Nordic 4km current dataset

nordic_native = reader_ROMS_native.Reader(o.test_data_folder() +
    '2Feb2016_Nordic_sigma_3d/Nordic-4km_SLEVELS_avg_00_subset2Feb2016.nc')
o.add_reader(nordic_native)

12:30 INFO    opendrift.readers:58: Opening file with xr.open_dataset
12:30 INFO    opendrift.readers.reader_ROMS_native:263: Read GLS parameters from file.
12:30 WARNING opendrift.readers.reader_ROMS_native:317: Duplicate variables for latitude, selecting lat_psi, and discarding lat_rho
12:30 WARNING opendrift.readers.reader_ROMS_native:317: Duplicate variables for latitude, selecting lat_psi, and discarding lat_u
12:30 WARNING opendrift.readers.reader_ROMS_native:317: Duplicate variables for latitude, selecting lat_psi, and discarding lat_v
12:30 WARNING opendrift.readers.reader_ROMS_native:317: Duplicate variables for longitude, selecting lon_psi, and discarding lon_rho
12:30 WARNING opendrift.readers.reader_ROMS_native:317: Duplicate variables for longitude, selecting lon_psi, and discarding lon_u
12:30 WARNING opendrift.readers.reader_ROMS_native:317: Duplicate variables for longitude, selecting lon_psi, and discarding lon_v
12:30 WARNING opendrift.readers.reader_ROMS_native:317: Duplicate variables for land_binary_mask, selecting mask_psi, and discarding mask_rho
12:30 WARNING opendrift.readers.reader_ROMS_native:317: Duplicate variables for ocean_s_coordinate_g2, selecting s_rho, and discarding s_w
12:30 INFO    opendrift.readers.reader_ROMS_native:324: 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', 'ocean_time', 'pm', 'pn', 'rdrg', 'rdrg2', 'rho0', 'spherical', 'swrad', 'sz_alpha', 'theta_b', 'theta_s', 'ubar', 'vbar', 'xl', 'zeta_detided']
12:30 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.
12:30 INFO    opendrift.readers.basereader.structured:90: Making interpolator for lon,lat to x,y conversion...

Seed elements at defined positions, depth and time

o.seed_elements(lon=12.0, lat=68.3, radius=0, number=10,
                z=np.linspace(0, -150, 10), time=nordic_native.start_time)

12:31 INFO    opendrift.models.basemodel.environment:206: Adding a dynamical landmask with max. priority based on assumed maximum speed of 2.0 m/s. Adding a customised landmask may be faster...
12:34 INFO    opendrift.models.basemodel.environment:233: Fallback values will be used for the following variables which have no readers:
12:34 INFO    opendrift.models.basemodel.environment:236:    x_wind: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    y_wind: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    upward_sea_water_velocity: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    ocean_vertical_diffusivity: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    sea_surface_wave_significant_height: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    sea_surface_swell_wave_to_direction: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    sea_surface_swell_wave_peak_period_from_variance_spectral_density: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    sea_surface_swell_wave_significant_height: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    sea_surface_wind_wave_to_direction: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    sea_surface_wind_wave_mean_period: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    sea_surface_wind_wave_significant_height: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    surface_downward_x_stress: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    surface_downward_y_stress: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    turbulent_kinetic_energy: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    turbulent_generic_length_scale: 0.000000
12:34 INFO    opendrift.models.basemodel.environment:236:    ocean_mixed_layer_thickness: 50.000000

Running model

o.run(time_step=3600)

12:34 INFO    opendrift:1841: Duration, steps or end time not specified, running until end of first reader: 2016-02-04 12:00:00
12:34 INFO    opendrift:1841: Duration, steps or end time not specified, running until end of first reader: 2016-02-04 12:00:00
12:34 INFO    opendrift:919: Using existing reader for land_binary_mask
12:34 INFO    opendrift:930: All points are in ocean
12:34 WARNING opendrift:713: Seafloor check not being run because environment is missing. This will happen the first time the function is run but if it happens subsequently there is probably a problem.
12:34 INFO    opendrift:2056: 2016-02-02 12:00:00 - step 1 of 48 - 10 active elements (0 deactivated)
12:34 INFO    opendrift.readers.reader_ROMS_native:374: Using mask_rho for mask_rho
12:34 INFO    opendrift.readers.reader_ROMS_native:425: Using zeta for sea surface height
12:35 INFO    opendrift.readers.reader_ROMS_native:395: Using mask_u for mask_u
12:35 INFO    opendrift.readers.reader_ROMS_native:646: Time: 0:00:00.049197
12:35 INFO    opendrift.readers.reader_ROMS_native:416: Using mask_v for mask_v
12:35 INFO    opendrift.readers.reader_ROMS_native:439: Using angle from Dataset.
12:35 INFO    opendrift:2056: 2016-02-02 13:00:00 - step 2 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-02 14:00:00 - step 3 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-02 15:00:00 - step 4 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-02 16:00:00 - step 5 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-02 17:00:00 - step 6 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-02 18:00:00 - step 7 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-02 19:00:00 - step 8 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-02 20:00:00 - step 9 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-02 21:00:00 - step 10 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-02 22:00:00 - step 11 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-02 23:00:00 - step 12 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-03 00:00:00 - step 13 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-03 01:00:00 - step 14 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-03 02:00:00 - step 15 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-03 03:00:00 - step 16 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-03 04:00:00 - step 17 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-03 05:00:00 - step 18 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-03 06:00:00 - step 19 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-03 07:00:00 - step 20 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-03 08:00:00 - step 21 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-03 09:00:00 - step 22 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-03 10:00:00 - step 23 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-03 11:00:00 - step 24 of 48 - 10 active elements (0 deactivated)
12:35 INFO    opendrift:2056: 2016-02-03 12:00:00 - step 25 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-03 13:00:00 - step 26 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-03 14:00:00 - step 27 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-03 15:00:00 - step 28 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-03 16:00:00 - step 29 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-03 17:00:00 - step 30 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-03 18:00:00 - step 31 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-03 19:00:00 - step 32 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-03 20:00:00 - step 33 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-03 21:00:00 - step 34 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-03 22:00:00 - step 35 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-03 23:00:00 - step 36 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-04 00:00:00 - step 37 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-04 01:00:00 - step 38 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-04 02:00:00 - step 39 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-04 03:00:00 - step 40 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-04 04:00:00 - step 41 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-04 05:00:00 - step 42 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-04 06:00:00 - step 43 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-04 07:00:00 - step 44 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-04 08:00:00 - step 45 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-04 09:00:00 - step 46 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-04 10:00:00 - step 47 of 48 - 10 active elements (0 deactivated)
12:36 INFO    opendrift:2056: 2016-02-04 11:00:00 - step 48 of 48 - 10 active elements (0 deactivated)

<xarray.Dataset> Size: 69kB
Dimensions:                                                                              (
                                                                                          trajectory: 10,
                                                                                          time: 49)
Coordinates:
  * trajectory                                                                           (trajectory) int64 80B ...
  * time                                                                                 (time) datetime64[ns] 392B ...
Data variables: (12/35)
    status                                                                               (trajectory, time) float32 2kB ...
    moving                                                                               (trajectory, time) float32 2kB ...
    age_seconds                                                                          (trajectory, time) float32 2kB ...
    origin_marker                                                                        (trajectory, time) float32 2kB ...
    lon                                                                                  (trajectory, time) float32 2kB ...
    lat                                                                                  (trajectory, time) float32 2kB ...
    ...                                                                                   ...
    surface_downward_y_stress                                                            (trajectory, time) float32 2kB ...
    turbulent_kinetic_energy                                                             (trajectory, time) float32 2kB ...
    turbulent_generic_length_scale                                                       (trajectory, time) float32 2kB ...
    ocean_mixed_layer_thickness                                                          (trajectory, time) float32 2kB ...
    sea_floor_depth_below_sea_level                                                      (trajectory, time) float32 2kB ...
    land_binary_mask                                                                     (trajectory, time) float32 2kB ...
Attributes: (12/153)
    Conventions:                                                             ...
    standard_name_vocabulary:                                                ...
    featureType:                                                             ...
    title:                                                                   ...
    summary:                                                                 ...
    keywords:                                                                ...
    ...                                                                                                               ...
    geospatial_lon_units:                                                    ...
    geospatial_lon_resolution:                                               ...
    runtime:                                                                 ...
    geospatial_vertical_min:                                                 ...
    geospatial_vertical_max:                                                 ...
    geospatial_vertical_positive:                                            ...

xarray.Dataset

Dimensions:
- trajectory: 10
- time: 49

Coordinates: (2)

trajectory
(trajectory)
int64
0 1 2 3 4 5 6 7 8 9
cf_role :
trajectory_id
dtype :
<class 'numpy.int32'>
```
array([0, 1, 2, 3, 4, 5, 6, 7, 8, 9])
```

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-02T13:00:00.000000000',
       '2016-02-02T14:00:00.000000000', '2016-02-02T15:00:00.000000000',
       '2016-02-02T16:00:00.000000000', '2016-02-02T17:00:00.000000000',
       '2016-02-02T18:00:00.000000000', '2016-02-02T19:00:00.000000000',
       '2016-02-02T20:00:00.000000000', '2016-02-02T21:00:00.000000000',
       '2016-02-02T22:00:00.000000000', '2016-02-02T23:00:00.000000000',
       '2016-02-03T00:00:00.000000000', '2016-02-03T01:00:00.000000000',
       '2016-02-03T02:00:00.000000000', '2016-02-03T03:00:00.000000000',
       '2016-02-03T04:00:00.000000000', '2016-02-03T05:00:00.000000000',
       '2016-02-03T06:00:00.000000000', '2016-02-03T07:00:00.000000000',
       '2016-02-03T08:00:00.000000000', '2016-02-03T09:00:00.000000000',
       '2016-02-03T10:00:00.000000000', '2016-02-03T11:00:00.000000000',
       '2016-02-03T12:00:00.000000000', '2016-02-03T13:00:00.000000000',
       '2016-02-03T14:00:00.000000000', '2016-02-03T15:00:00.000000000',
       '2016-02-03T16:00:00.000000000', '2016-02-03T17:00:00.000000000',
       '2016-02-03T18:00:00.000000000', '2016-02-03T19:00:00.000000000',
       '2016-02-03T20:00:00.000000000', '2016-02-03T21:00:00.000000000',
       '2016-02-03T22:00:00.000000000', '2016-02-03T23:00:00.000000000',
       '2016-02-04T00:00:00.000000000', '2016-02-04T01:00:00.000000000',
       '2016-02-04T02:00:00.000000000', '2016-02-04T03:00:00.000000000',
       '2016-02-04T04:00:00.000000000', '2016-02-04T05:00:00.000000000',
       '2016-02-04T06:00:00.000000000', '2016-02-04T07:00:00.000000000',
       '2016-02-04T08:00:00.000000000', '2016-02-04T09:00:00.000000000',
       '2016-02-04T10:00:00.000000000', '2016-02-04T11:00:00.000000000',
       '2016-02-04T12:00:00.000000000'], dtype='datetime64[ns]')

Data variables: (35)

Indexes: (2)

trajectory

PandasIndex

PandasIndex(Index([0, 1, 2, 3, 4, 5, 6, 7, 8, 9], dtype='int64', name='trajectory'))

time

PandasIndex

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

Attributes: (153)
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:12:34.958072
date_created :
2025-02-18T17:12:34.958079
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 :
1:00:00
time_step_output :
1: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:sea_surface_wave_period_at_variance_spectral_density_maximum :
None
config_environment:fallback:sea_surface_wave_period_at_variance_spectral_density_maximum :
0
config_environment:constant:sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment :
None
config_environment:fallback:sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment :
0
config_environment:constant:sea_surface_swell_wave_to_direction :
None
config_environment:fallback:sea_surface_swell_wave_to_direction :
0
config_environment:constant:sea_surface_swell_wave_peak_period_from_variance_spectral_density :
None
config_environment:fallback:sea_surface_swell_wave_peak_period_from_variance_spectral_density :
0
config_environment:constant:sea_surface_swell_wave_significant_height :
None
config_environment:fallback:sea_surface_swell_wave_significant_height :
0
config_environment:constant:sea_surface_wind_wave_to_direction :
None
config_environment:fallback:sea_surface_wind_wave_to_direction :
0
config_environment:constant:sea_surface_wind_wave_mean_period :
None
config_environment:fallback:sea_surface_wind_wave_mean_period :
0
config_environment:constant:sea_surface_wind_wave_significant_height :
None
config_environment:fallback:sea_surface_wind_wave_significant_height :
0
config_environment:constant:surface_downward_x_stress :
None
config_environment:fallback:surface_downward_x_stress :
0
config_environment:constant:surface_downward_y_stress :
None
config_environment:fallback:surface_downward_y_stress :
0
config_environment:constant:turbulent_kinetic_energy :
None
config_environment:fallback:turbulent_kinetic_energy :
0
config_environment:constant:turbulent_generic_length_scale :
None
config_environment:fallback:turbulent_generic_length_scale :
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_mixing :
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.13.0
seed_geojson :
{"features": [], "type": "FeatureCollection"}
simulation_time :
2025-02-18 17:12:34.965015
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_ocean_vertical_diffusivity :
0
reader_sea_surface_wave_significant_height :
0
reader_sea_surface_wave_stokes_drift_x_velocity :
0
reader_sea_surface_wave_stokes_drift_y_velocity :
0
reader_sea_surface_wave_period_at_variance_spectral_density_maximum :
0
reader_sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment :
0
reader_sea_surface_swell_wave_to_direction :
0
reader_sea_surface_swell_wave_peak_period_from_variance_spectral_density :
0
reader_sea_surface_swell_wave_significant_height :
0
reader_sea_surface_wind_wave_to_direction :
0
reader_sea_surface_wind_wave_mean_period :
0
reader_sea_surface_wind_wave_significant_height :
0
reader_surface_downward_x_stress :
0
reader_surface_downward_y_stress :
0
reader_turbulent_kinetic_energy :
0
reader_turbulent_generic_length_scale :
0
reader_ocean_mixed_layer_thickness :
50
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 :
P0DT1H0M0S
performance :
-------------------- Reader performance: -------------------- roms native 0:00:00.4 total 0:00:00.0 preparing 0:00:00.3 reading 0:00:00.1 interpolation 0:00:00.0 interpolation_time 0:00:00.1 rotating vectors 0:00:00.0 masking -------------------- global_landmask 0:00:00.0 total 0:00:00.0 preparing 0:00:00.0 reading 0:00:00.0 masking -------------------- Performance: 6.1 total time 4.3 configuration 0.0 preparing main loop 0.0 moving elements to ocean 1.8 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.3
geospatial_lat_max :
68.59404
geospatial_lat_units :
degrees_north
geospatial_lat_resolution :
point
geospatial_lon_min :
12.0
geospatial_lon_max :
12.947481
geospatial_lon_units :
degrees_east
geospatial_lon_resolution :
point
runtime :
0:00:06.164972
geospatial_vertical_min :
-144.72884
geospatial_vertical_max :
0.0
geospatial_vertical_positive :
up

Print and plot results, with lines colored by particle depth

print(o)
o.plot(linecolor='z', fast=True)
#o.animation()

OpenDrift - OceanDrift 2016-02-02 12:00 to 2016-02-04 12:00 UTC (49 steps)

===========================
--------------------
Reader performance:
--------------------
roms native
 0:00:00.4  total
 0:00:00.0  preparing
 0:00:00.3  reading
 0:00:00.1  interpolation
 0:00:00.0  interpolation_time
 0:00:00.1  rotating vectors
 0:00:00.0  masking
--------------------
global_landmask
 0:00:00.0  total
 0:00:00.0  preparing
 0:00:00.0  reading
 0:00:00.0  masking
--------------------
Performance:
    6.1 total time
    4.3 configuration
    0.0 preparing main loop
      0.0 moving elements to ocean
    1.8 main loop
      0.0 updating elements
    0.0 cleaning up
--------------------
===========================
Model:  OceanDrift     (OpenDrift version 1.13.0)
        10 active Lagrangian3DArray particles  (0 deactivated, 0 scheduled)
-------------------
Environment variables:
  -----
  sea_floor_depth_below_sea_level
  sea_surface_height
  x_sea_water_velocity
  y_sea_water_velocity
     1) roms native
  -----
  land_binary_mask
     1) global_landmask
  -----
Readers not added for the following variables:
  ocean_mixed_layer_thickness
  ocean_vertical_diffusivity
  sea_surface_swell_wave_peak_period_from_variance_spectral_density
  sea_surface_swell_wave_significant_height
  sea_surface_swell_wave_to_direction
  sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment
  sea_surface_wave_period_at_variance_spectral_density_maximum
  sea_surface_wave_significant_height
  sea_surface_wave_stokes_drift_x_velocity
  sea_surface_wave_stokes_drift_y_velocity
  sea_surface_wind_wave_mean_period
  sea_surface_wind_wave_significant_height
  sea_surface_wind_wave_to_direction
  surface_downward_x_stress
  surface_downward_y_stress
  turbulent_generic_length_scale
  turbulent_kinetic_energy
  upward_sea_water_velocity
  x_wind
  y_wind

Discarded readers:

Time:
        Start: 2016-02-02 12:00:00 UTC
        Present: 2016-02-04 12:00:00 UTC
        Calculation steps: 48 * 1:00:00 - total time: 2 days, 0:00:00
        Output steps: 49 * 1:00:00
===========================

17:12:37 WARNING opendrift:2429: Plotting fast. This will make your plots less accurate.

(<GeoAxes: title={'center': 'OpenDrift - OceanDrift\n2016-02-02 12:00 to 2016-02-04 12:00 UTC (49 steps)'}>, <Figure size 1017.46x1100 with 2 Axes>)

Total running time of the script: (0 minutes 12.875 seconds)

Gallery generated by Sphinx-Gallery