Openoil sample output netCDF file

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

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

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

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

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

o.add_reader([rc, reader_norkyst, reader_arome])
14:01:22 INFO    opendrift.models.basemodel:539: OpenDriftSimulation initialised (version 1.10.7 / v1.10.6-119-g1da5bec)
14:01:22 INFO    opendrift.readers.reader_netCDF_CF_generic:92: Opening dataset: /root/project/tests/test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc
/opt/conda/envs/opendrift/lib/python3.11/site-packages/pyproj/crs/crs.py:1286: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems
  proj = self._crs.to_proj4(version=version)
14:01:22 INFO    opendrift.readers.reader_netCDF_CF_generic:92: Opening dataset: /root/project/tests/test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc
/opt/conda/envs/opendrift/lib/python3.11/site-packages/pyproj/crs/crs.py:1286: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems
  proj = self._crs.to_proj4(version=version)

Seed oil elements at defined position and time

time = [reader_arome.start_time,
        reader_arome.start_time + timedelta(hours=1)]
o.seed_elements(lon=5.05, lat=59.95, radius=1500, number=100,
                time=time, z=0, m3_per_hour=2, oil_type='TAU 1999')
14:01:22 INFO    opendrift.models.openoil.openoil:1462: setting oil_type to: TAU 1999
14:01:22 INFO    opendrift.models.openoil.adios.dirjs:71: Querying ADIOS database for oil: TAU 1999
14:01:22 INFO    opendrift.models.openoil.openoil:1652: Using density 851.003285 and viscosity 0.001249255043769228 of oiltype TAU 1999

Running model

ncfile = 'openoil_sample_output.nc'
o.run(steps=4*8, time_step=900, time_step_output=3600, outfile=ncfile)
14:01:22 INFO    opendrift.models.basemodel:2591: Fallback values will be used for the following variables which have no readers:
14:01:22 INFO    opendrift.models.basemodel:2594:       upward_sea_water_velocity: 0.000000
14:01:22 INFO    opendrift.models.basemodel:2594:       sea_surface_wave_significant_height: 0.000000
14:01:22 INFO    opendrift.models.basemodel:2594:       sea_surface_wave_stokes_drift_x_velocity: 0.000000
14:01:22 INFO    opendrift.models.basemodel:2594:       sea_surface_wave_stokes_drift_y_velocity: 0.000000
14:01:22 INFO    opendrift.models.basemodel:2594:       sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
14:01:22 INFO    opendrift.models.basemodel:2594:       sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
14:01:22 INFO    opendrift.models.basemodel:2594:       sea_ice_area_fraction: 0.000000
14:01:22 INFO    opendrift.models.basemodel:2594:       sea_ice_x_velocity: 0.000000
14:01:22 INFO    opendrift.models.basemodel:2594:       sea_ice_y_velocity: 0.000000
14:01:22 INFO    opendrift.models.basemodel:2594:       sea_water_temperature: 10.000000
14:01:22 INFO    opendrift.models.basemodel:2594:       sea_water_salinity: 34.000000
14:01:22 INFO    opendrift.models.basemodel:2594:       sea_floor_depth_below_sea_level: 10000.000000
14:01:22 INFO    opendrift.models.basemodel:2594:       ocean_vertical_diffusivity: 0.020000
14:01:22 INFO    opendrift.models.basemodel:2594:       ocean_mixed_layer_thickness: 50.000000
14:01:22 INFO    opendrift.models.basemodel:2753: Adding a dynamical landmask with max. priority based on assumed maximum speed of 1.3 m/s. Adding a customised landmask may be faster...
14:01:29 INFO    opendrift.models.basemodel:1701: Using existing reader for land_binary_mask
14:01:29 INFO    opendrift.models.basemodel:1715: Moving 1 out of 100 points from land to water
14:01:29 INFO    opendrift.models.openoil.openoil:646: Oil-water surface tension is 0.030101 Nm
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 00:00:00 - step 1 of 32 - 25 active elements (0 deactivated)
/root/project/opendrift/readers/interpolation/interpolators.py:17: RuntimeWarning: overflow encountered in cast
  data[mask] = np.finfo(np.float64).min
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 00:15:00 - step 2 of 32 - 50 active elements (0 deactivated)
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 00:30:00 - step 3 of 32 - 75 active elements (0 deactivated)
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 00:45:00 - step 4 of 32 - 100 active elements (0 deactivated)
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 01:00:00 - step 5 of 32 - 100 active elements (0 deactivated)
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 01:15:00 - step 6 of 32 - 99 active elements (1 deactivated)
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 01:30:00 - step 7 of 32 - 99 active elements (1 deactivated)
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 01:45:00 - step 8 of 32 - 95 active elements (5 deactivated)
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 02:00:00 - step 9 of 32 - 92 active elements (8 deactivated)
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 02:15:00 - step 10 of 32 - 90 active elements (10 deactivated)
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 02:30:00 - step 11 of 32 - 87 active elements (13 deactivated)
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 02:45:00 - step 12 of 32 - 80 active elements (20 deactivated)
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 03:00:00 - step 13 of 32 - 75 active elements (25 deactivated)
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 03:15:00 - step 14 of 32 - 72 active elements (28 deactivated)
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 03:30:00 - step 15 of 32 - 63 active elements (37 deactivated)
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 03:45:00 - step 16 of 32 - 60 active elements (40 deactivated)
14:01:29 INFO    opendrift.models.basemodel:2882: 2015-11-16 04:00:00 - step 17 of 32 - 53 active elements (47 deactivated)
14:01:30 INFO    opendrift.models.basemodel:2882: 2015-11-16 04:15:00 - step 18 of 32 - 46 active elements (54 deactivated)
14:01:30 INFO    opendrift.models.basemodel:2882: 2015-11-16 04:30:00 - step 19 of 32 - 38 active elements (62 deactivated)
14:01:30 INFO    opendrift.models.basemodel:2882: 2015-11-16 04:45:00 - step 20 of 32 - 35 active elements (65 deactivated)
14:01:30 INFO    opendrift.models.basemodel:2882: 2015-11-16 05:00:00 - step 21 of 32 - 32 active elements (68 deactivated)
14:01:30 INFO    opendrift.models.basemodel:2882: 2015-11-16 05:15:00 - step 22 of 32 - 25 active elements (75 deactivated)
14:01:30 INFO    opendrift.models.basemodel:2882: 2015-11-16 05:30:00 - step 23 of 32 - 21 active elements (79 deactivated)
14:01:30 INFO    opendrift.models.basemodel:2882: 2015-11-16 05:45:00 - step 24 of 32 - 17 active elements (83 deactivated)
14:01:30 INFO    opendrift.models.basemodel:2882: 2015-11-16 06:00:00 - step 25 of 32 - 16 active elements (84 deactivated)
14:01:30 INFO    opendrift.models.basemodel:2882: 2015-11-16 06:15:00 - step 26 of 32 - 15 active elements (85 deactivated)
14:01:30 INFO    opendrift.models.basemodel:2882: 2015-11-16 06:30:00 - step 27 of 32 - 13 active elements (87 deactivated)
14:01:30 INFO    opendrift.models.basemodel:2882: 2015-11-16 06:45:00 - step 28 of 32 - 10 active elements (90 deactivated)
14:01:30 INFO    opendrift.models.basemodel:2882: 2015-11-16 07:00:00 - step 29 of 32 - 7 active elements (93 deactivated)
14:01:30 INFO    opendrift.models.basemodel:2882: 2015-11-16 07:15:00 - step 30 of 32 - 6 active elements (94 deactivated)
14:01:30 INFO    opendrift.models.basemodel:2882: 2015-11-16 07:30:00 - step 31 of 32 - 3 active elements (97 deactivated)
14:01:30 INFO    opendrift.models.basemodel:2882: 2015-11-16 07:45:00 - step 32 of 32 - 1 active elements (99 deactivated)
14:01:30 INFO    opendrift.export.io_netcdf:111: Wrote 9 steps to file openoil_sample_output.nc
/opt/conda/envs/opendrift/lib/python3.11/site-packages/numpy/ma/core.py:467: RuntimeWarning: invalid value encountered in cast
  fill_value = np.array(fill_value, copy=False, dtype=ndtype)

Print and plot results

o.plot_oil_budget()
o.plot()
  • TAU 1999 (851.0 kg/m3) - 2015-11-16 00:00 to 2015-11-16 08:00
  • OpenDrift - OpenOil (TAU 1999) 2015-11-16 00:00 to 2015-11-16 08:00 UTC (9 steps)
/opt/conda/envs/opendrift/lib/python3.11/site-packages/numpy/ma/core.py:1020: RuntimeWarning: overflow encountered in multiply
  result = self.f(da, db, *args, **kwargs)
/opt/conda/envs/opendrift/lib/python3.11/site-packages/numpy/ma/core.py:6975: RuntimeWarning: overflow encountered in power
  result = np.where(m, fa, umath.power(fa, fb)).view(basetype)
0.9128345689654023 ALPHA

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

ncdump of the output file

from subprocess import check_output
args = [ "ncdump", "-h", ncfile ]
ncdump = check_output(args).decode().strip()
print(ncdump)
netcdf openoil_sample_output {
dimensions:
        trajectory = 100 ;
        time = 9 ;
variables:
        int trajectory(trajectory) ;
                trajectory:cf_role = "trajectory_id" ;
                trajectory:units = "1" ;
        double time(time) ;
                time:units = "seconds since 1970-01-01 00:00:00" ;
                time:standard_name = "time" ;
                time:long_name = "time" ;
        int status(trajectory, time) ;
                status:valid_range = 0, 1 ;
                status:flag_values = 0, 1 ;
                status:flag_meanings = "active stranded" ;
                status:minval = 0 ;
                status:maxval = 1 ;
        int moving(trajectory, time) ;
                moving:minval = 0 ;
                moving:maxval = 1 ;
        float age_seconds(trajectory, time) ;
                age_seconds:units = "s" ;
                age_seconds:minval = 900.f ;
                age_seconds:maxval = 27000.f ;
        int origin_marker(trajectory, time) ;
                origin_marker:unit = "" ;
                origin_marker:description = "An integer kept constant during the simulation. Different values may be used for different seedings, to separate elements during analysis. With GUI, only a single seeding is possible." ;
                origin_marker:flag_values = 0LL ;
                origin_marker:flag_meanings = "Seed_0" ;
                origin_marker:minval = 0 ;
                origin_marker:maxval = 0 ;
        float lon(trajectory, time) ;
                lon:units = "degrees_east" ;
                lon:standard_name = "longitude" ;
                lon:long_name = "longitude" ;
                lon:minval = 4.981464f ;
                lon:maxval = 5.143353f ;
        float lat(trajectory, time) ;
                lat:units = "degrees_north" ;
                lat:standard_name = "latitude" ;
                lat:long_name = "latitude" ;
                lat:minval = 59.92179f ;
                lat:maxval = 60.01746f ;
        float z(trajectory, time) ;
                z:units = "m" ;
                z:standard_name = "z" ;
                z:long_name = "vertical position" ;
                z:axis = "Z" ;
                z:positive = "up" ;
                z:minval = 0.f ;
                z:maxval = 0.f ;
        float wind_drift_factor(trajectory, time) ;
                wind_drift_factor:units = "%" ;
                wind_drift_factor:description = "Elements at the ocean surface are moved by this fraction of the wind vector, in addition to currents and Stokes drift" ;
                wind_drift_factor:minval = 0.03f ;
                wind_drift_factor:maxval = 0.03f ;
        float current_drift_factor(trajectory, time) ;
                current_drift_factor:units = "1" ;
                current_drift_factor:description = "Elements are moved with this fraction of the current vector, in addition to currents and Stokes drift" ;
                current_drift_factor:minval = 1.f ;
                current_drift_factor:maxval = 1.f ;
        float terminal_velocity(trajectory, time) ;
                terminal_velocity:units = "m/s" ;
                terminal_velocity:description = "Terminal rise/sinking velocity (buoyancy) in the ocean column" ;
                terminal_velocity:minval = 0.f ;
                terminal_velocity:maxval = 0.1068134f ;
        float mass_oil(trajectory, time) ;
                mass_oil:units = "kg" ;
                mass_oil:minval = 8.185321f ;
                mass_oil:maxval = 17.02007f ;
        float viscosity(trajectory, time) ;
                viscosity:units = "m2/s" ;
                viscosity:description = "Kinematic viscosity of oil" ;
                viscosity:minval = 0.001249255f ;
                viscosity:maxval = 4.003683f ;
        float density(trajectory, time) ;
                density:units = "kg/m^3" ;
                density:minval = 851.0033f ;
                density:maxval = 970.0526f ;
        float bulltime(trajectory, time) ;
                bulltime:units = "s" ;
                bulltime:minval = 0.f ;
                bulltime:maxval = 0.f ;
        float interfacial_area(trajectory, time) ;
                interfacial_area:units = "m2" ;
                interfacial_area:minval = 0.f ;
                interfacial_area:maxval = 1296232.f ;
        float mass_dispersed(trajectory, time) ;
                mass_dispersed:units = "kg" ;
                mass_dispersed:minval = 0.f ;
                mass_dispersed:maxval = 0.3069384f ;
        float mass_evaporated(trajectory, time) ;
                mass_evaporated:units = "kg" ;
                mass_evaporated:minval = 0.f ;
                mass_evaporated:maxval = 8.568092f ;
        float mass_biodegraded(trajectory, time) ;
                mass_biodegraded:units = "kg" ;
                mass_biodegraded:minval = 0.f ;
                mass_biodegraded:maxval = 0.f ;
        float fraction_evaporated(trajectory, time) ;
                fraction_evaporated:units = "%" ;
                fraction_evaporated:minval = 0.f ;
                fraction_evaporated:maxval = 0.5100498f ;
        float water_fraction(trajectory, time) ;
                water_fraction:units = "%" ;
                water_fraction:minval = 0.f ;
                water_fraction:maxval = 0.6835831f ;
        float oil_film_thickness(trajectory, time) ;
                oil_film_thickness:units = "m" ;
                oil_film_thickness:minval = 0.001f ;
                oil_film_thickness:maxval = 0.001f ;
        float diameter(trajectory, time) ;
                diameter:units = "m" ;
                diameter:minval = 0.f ;
                diameter:maxval = 0.002621652f ;
        float x_sea_water_velocity(trajectory, time) ;
                x_sea_water_velocity:minval = -0.4167673f ;
                x_sea_water_velocity:maxval = 0.1973607f ;
        float y_sea_water_velocity(trajectory, time) ;
                y_sea_water_velocity:minval = -0.1364768f ;
                y_sea_water_velocity:maxval = 0.6068177f ;
        float x_wind(trajectory, time) ;
                x_wind:minval = 3.613704f ;
                x_wind:maxval = 6.447049f ;
        float y_wind(trajectory, time) ;
                y_wind:minval = 2.56165f ;
                y_wind:maxval = 5.502532f ;
        float upward_sea_water_velocity(trajectory, time) ;
                upward_sea_water_velocity:minval = 0.f ;
                upward_sea_water_velocity:maxval = 0.f ;
        float sea_surface_wave_significant_height(trajectory, time) ;
                sea_surface_wave_significant_height:minval = 0.5444739f ;
                sea_surface_wave_significant_height:maxval = 1.538757f ;
        float sea_surface_wave_stokes_drift_x_velocity(trajectory, time) ;
                sea_surface_wave_stokes_drift_x_velocity:minval = 0.f ;
                sea_surface_wave_stokes_drift_x_velocity:maxval = 0.f ;
        float sea_surface_wave_stokes_drift_y_velocity(trajectory, time) ;
                sea_surface_wave_stokes_drift_y_velocity:minval = 0.f ;
                sea_surface_wave_stokes_drift_y_velocity:maxval = 0.f ;
        float sea_surface_wave_period_at_variance_spectral_density_maximum(trajectory, time) ;
                sea_surface_wave_period_at_variance_spectral_density_maximum:minval = 0.f ;
                sea_surface_wave_period_at_variance_spectral_density_maximum:maxval = 0.f ;
        float sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment(trajectory, time) ;
                sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment:minval = 4.019927f ;
                sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment:maxval = 6.757945f ;
        float sea_ice_area_fraction(trajectory, time) ;
                sea_ice_area_fraction:minval = 0.f ;
                sea_ice_area_fraction:maxval = 0.f ;
        float sea_ice_x_velocity(trajectory, time) ;
                sea_ice_x_velocity:minval = 0.f ;
                sea_ice_x_velocity:maxval = 0.f ;
        float sea_ice_y_velocity(trajectory, time) ;
                sea_ice_y_velocity:minval = 0.f ;
                sea_ice_y_velocity:maxval = 0.f ;
        float sea_water_temperature(trajectory, time) ;
                sea_water_temperature:minval = 10.f ;
                sea_water_temperature:maxval = 283.15f ;
        float sea_water_salinity(trajectory, time) ;
                sea_water_salinity:minval = 34.f ;
                sea_water_salinity:maxval = 34.f ;
        float sea_floor_depth_below_sea_level(trajectory, time) ;
                sea_floor_depth_below_sea_level:minval = 10000.f ;
                sea_floor_depth_below_sea_level:maxval = 10000.f ;
        float ocean_vertical_diffusivity(trajectory, time) ;
                ocean_vertical_diffusivity:minval = 0.02f ;
                ocean_vertical_diffusivity:maxval = 0.02f ;
        float land_binary_mask(trajectory, time) ;
                land_binary_mask:minval = 0.f ;
                land_binary_mask:maxval = 1.f ;
        float ocean_mixed_layer_thickness(trajectory, time) ;
                ocean_mixed_layer_thickness:minval = 50.f ;
                ocean_mixed_layer_thickness:maxval = 50.f ;

// global attributes:
                :Conventions = "CF-1.6" ;
                :standard_name_vocabulary = "CF-1.6" ;
                :featureType = "trajectory" ;
                :history = "Created 2023-09-01 14:01:29.049958" ;
                :source = "Output from simulation with OpenDrift" ;
                :model_url = "https://github.com/OpenDrift/opendrift" ;
                :opendrift_class = "OpenOil" ;
                :opendrift_module = "opendrift.models.openoil.openoil" ;
                :readers = "odict_keys([\'constant_reader\', \'/root/project/tests/test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc\', \'/root/project/tests/test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc\', \'global_landmask\'])" ;
                :time_coverage_start = "2015-11-16 00:00:00" ;
                :time_step_calculation = "0:15:00" ;
                :time_step_output = "1:00:00" ;
                :config_general\:use_auto_landmask = "True" ;
                :config_general\:coastline_action = "stranding" ;
                :config_general\:time_step_minutes = 60LL ;
                :config_general\:time_step_output_minutes = "None" ;
                :config_seed\:ocean_only = "True" ;
                :config_seed\:number = 1LL ;
                :config_drift\:max_age_seconds = "None" ;
                :config_drift\:advection_scheme = "euler" ;
                :config_drift\:current_uncertainty = 0.05 ;
                :config_drift\:current_uncertainty_uniform = 0LL ;
                :config_drift\:horizontal_diffusivity = 0LL ;
                :config_drift\:wind_uncertainty = 0.5 ;
                :config_drift\:relative_wind = "False" ;
                :config_drift\:deactivate_north_of = "None" ;
                :config_drift\:deactivate_south_of = "None" ;
                :config_drift\:deactivate_east_of = "None" ;
                :config_drift\:deactivate_west_of = "None" ;
                :config_readers\:max_number_of_fails = 1LL ;
                :config_seed\:origin_marker = 0LL ;
                :config_seed\:z = 0LL ;
                :config_seed\:wind_drift_factor = 0.03 ;
                :config_seed\:current_drift_factor = 1LL ;
                :config_seed\:terminal_velocity = 0. ;
                :config_seed\:oil_film_thickness = 0.001 ;
                :config_environment\:constant\:x_sea_water_velocity = "None" ;
                :config_environment\:fallback\:x_sea_water_velocity = "None" ;
                :config_environment\:constant\:y_sea_water_velocity = "None" ;
                :config_environment\:fallback\:y_sea_water_velocity = "None" ;
                :config_environment\:constant\:x_wind = "None" ;
                :config_environment\:fallback\:x_wind = "None" ;
                :config_environment\:constant\:y_wind = "None" ;
                :config_environment\:fallback\:y_wind = "None" ;
                :config_environment\:constant\:upward_sea_water_velocity = "None" ;
                :config_environment\:fallback\:upward_sea_water_velocity = 0LL ;
                :config_environment\:constant\:sea_surface_wave_significant_height = "None" ;
                :config_environment\:fallback\:sea_surface_wave_significant_height = 0LL ;
                :config_environment\:constant\:sea_surface_wave_stokes_drift_x_velocity = "None" ;
                :config_environment\:fallback\:sea_surface_wave_stokes_drift_x_velocity = 0LL ;
                :config_environment\:constant\:sea_surface_wave_stokes_drift_y_velocity = "None" ;
                :config_environment\:fallback\:sea_surface_wave_stokes_drift_y_velocity = 0LL ;
                :config_environment\:constant\:sea_surface_wave_period_at_variance_spectral_density_maximum = "None" ;
                :config_environment\:fallback\:sea_surface_wave_period_at_variance_spectral_density_maximum = 0LL ;
                :config_environment\:constant\:sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment = "None" ;
                :config_environment\:fallback\:sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment = 0LL ;
                :config_environment\:constant\:sea_ice_area_fraction = "None" ;
                :config_environment\:fallback\:sea_ice_area_fraction = 0LL ;
                :config_environment\:constant\:sea_ice_x_velocity = "None" ;
                :config_environment\:fallback\:sea_ice_x_velocity = 0LL ;
                :config_environment\:constant\:sea_ice_y_velocity = "None" ;
                :config_environment\:fallback\:sea_ice_y_velocity = 0LL ;
                :config_environment\:constant\:sea_water_temperature = "None" ;
                :config_environment\:fallback\:sea_water_temperature = 10LL ;
                :config_environment\:constant\:sea_water_salinity = "None" ;
                :config_environment\:fallback\:sea_water_salinity = 34LL ;
                :config_environment\:constant\:sea_floor_depth_below_sea_level = "None" ;
                :config_environment\:fallback\:sea_floor_depth_below_sea_level = 10000LL ;
                :config_environment\:constant\:ocean_vertical_diffusivity = "None" ;
                :config_environment\:fallback\:ocean_vertical_diffusivity = 0.02 ;
                :config_environment\:constant\:land_binary_mask = "None" ;
                :config_environment\:fallback\:land_binary_mask = "None" ;
                :config_environment\:constant\:ocean_mixed_layer_thickness = "None" ;
                :config_environment\:fallback\:ocean_mixed_layer_thickness = 50LL ;
                :config_drift\:vertical_advection = "False" ;
                :config_drift\:vertical_mixing = "True" ;
                :config_vertical_mixing\:timestep = 60LL ;
                :config_vertical_mixing\:diffusivitymodel = "environment" ;
                :config_vertical_mixing\:background_diffusivity = 1.2e-05 ;
                :config_vertical_mixing\:TSprofiles = "False" ;
                :config_drift\:wind_drift_depth = 0.1 ;
                :config_drift\:stokes_drift = "True" ;
                :config_drift\:stokes_drift_profile = "Phillips" ;
                :config_drift\:use_tabularised_stokes_drift = "False" ;
                :config_drift\:tabularised_stokes_drift_fetch = "25000" ;
                :config_general\:seafloor_action = "lift_to_seafloor" ;
                :config_drift\:truncate_ocean_model_below_m = "None" ;
                :config_seed\:seafloor = "False" ;
                :config_seed\:m3_per_hour = 1LL ;
                :config_seed\:droplet_size_distribution = "uniform" ;
                :config_seed\:droplet_diameter_mu = 0.001 ;
                :config_seed\:droplet_diameter_sigma = 0.0005 ;
                :config_seed\:droplet_diameter_min_subsea = 0.0005 ;
                :config_seed\:droplet_diameter_max_subsea = 0.005 ;
                :config_processes\:dispersion = "True" ;
                :config_processes\:evaporation = "True" ;
                :config_processes\:emulsification = "True" ;
                :config_processes\:biodegradation = "False" ;
                :config_processes\:update_oilfilm_thickness = "False" ;
                :config_wave_entrainment\:droplet_size_distribution = "Johansen et al. (2015)" ;
                :config_wave_entrainment\:entrainment_rate = "Li et al. (2017)" ;
                :config_seed\:oil_type = "TAU 1999" ;
                :opendrift_version = "1.10.7" ;
                :seed_lon = "5.05" ;
                :seed_lat = "59.95" ;
                :seed_radius = "1500.0" ;
                :seed_time = "2015-11-16 00:00:00" ;
                :seed_number = "100.0" ;
                :seed_z = "0.0" ;
                :seed_m3_per_hour = "2.0" ;
                :seed_oiltype = "TAU 1999" ;
                :seed_geojson = "{\"features\": [], \"type\": \"FeatureCollection\"}" ;
                :steps_exported = 9LL ;
                :time_coverage_end = "2015-11-16 08:00:00" ;
                :performance = "--------------------\nReader performance:\n--------------------\nconstant_reader\n 0:00:00.0  total\n 0:00:00.0  preparing\n 0:00:00.0  reading\n 0:00:00.0  masking\n--------------------\n/root/project/tests/test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc\n 0:00:00.3  total\n 0:00:00.0  preparing\n 0:00:00.0  reading\n 0:00:00.0  interpolation\n 0:00:00.0  interpolation_time\n 0:00:00.2  rotating vectors\n 0:00:00.0  masking\n--------------------\n/root/project/tests/test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc\n--------------------\nglobal_landmask\n 0:00:00.0  total\n 0:00:00.0  preparing\n 0:00:00.0  reading\n 0:00:00.0  masking\n--------------------\nPerformance:\n      0 total time\n    0.4 configuration\n    6.3 preparing main loop\n      6.3 making dynamical landmask\n      0.0 moving elements to ocean\n      0.6 readers\n        0.0 global_landmask\n        0.1 postprocessing\n    1.8 main loop\n        0.0 constant_reader\n        0.3 /root/project/tests/test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc\n      1.0 updating elements\n        0.0 oil weathering\n          0.0 updating viscosities\n          0.0 updating densities\n          0.0 evaporation\n          0.0 emulsification\n          0.0 dispersion\n        0.9 vertical mixing\n      0 cleaning up\n--------------------\n" ;
                :simulation_time = "2023-09-01 14:01:29.078373" ;
                :reader_x_sea_water_velocity = "[\'/root/project/tests/test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc\']" ;
                :reader_y_sea_water_velocity = "[\'/root/project/tests/test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc\']" ;
                :reader_x_wind = "5" ;
                :reader_y_wind = "4" ;
                :reader_upward_sea_water_velocity = "0" ;
                :reader_sea_surface_wave_significant_height = "0" ;
                :reader_sea_surface_wave_stokes_drift_x_velocity = "0" ;
                :reader_sea_surface_wave_stokes_drift_y_velocity = "0" ;
                :reader_sea_surface_wave_period_at_variance_spectral_density_maximum = "0" ;
                :reader_sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment = "0" ;
                :reader_sea_ice_area_fraction = "0" ;
                :reader_sea_ice_x_velocity = "0" ;
                :reader_sea_ice_y_velocity = "0" ;
                :reader_sea_water_temperature = "10" ;
                :reader_sea_water_salinity = "34" ;
                :reader_sea_floor_depth_below_sea_level = "10000" ;
                :reader_ocean_vertical_diffusivity = "0.02" ;
                :reader_land_binary_mask = "[\'global_landmask\']" ;
                :reader_ocean_mixed_layer_thickness = "50" ;
                :geospatial_lat_min = 0. ;
                :geospatial_lat_max = 0. ;
                :geospatial_lat_units = "degrees_north" ;
                :geospatial_lat_resolution = "point" ;
                :geospatial_lon_min = 0. ;
                :geospatial_lon_max = 0. ;
                :geospatial_lon_units = "degrees_east" ;
                :geospatial_lon_resolution = "point" ;
                :runtime = "0:00:08.689466" ;
}

Download sample file

Total running time of the script: (1 minutes 51.811 seconds)

Gallery generated by Sphinx-Gallery