Note
Go to the end to download the full example code
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()
/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)