Seafloor oil spill

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

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

# Norkyst
reader_norkyst = reader_netCDF_CF_generic.Reader(o.test_data_folder() + '14Jan2016_NorKyst_z_3d/NorKyst-800m_ZDEPTHS_his_00_3Dsubset.nc')
#reader_norkyst = reader_netCDF_CF_generic.Reader('https://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')

o.add_reader([reader_norkyst])
o.set_config('environment:fallback:x_wind', 3)
o.set_config('environment:fallback:y_wind', 7)
o.set_config('drift:vertical_mixing', True)
20:22:42 INFO    opendrift.models.basemodel:515: OpenDriftSimulation initialised (version 1.12.0 / v1.12.0-26-g390e945)
20:22:42 INFO    opendrift.readers:58: Opening file with xr.open_dataset
20:22:42 INFO    opendrift.readers.reader_netCDF_CF_generic:337: Detected dimensions: {'x': 'X', 'y': 'Y', 'z': 'depth', 'time': 'time'}

Setting the range of droplet sizes for the seafloor release

o.set_config('seed:droplet_size_distribution','uniform')
o.set_config('seed:droplet_diameter_min_subsea', 0.0001)
o.set_config('seed:droplet_diameter_max_subsea', 0.0005)

# 'normal' and 'lognormal' distributions can also be specified
# o.set_config('seed:droplet_size_distribution','lognormal')
# o.set_config('seed:droplet_diameter_mu',0.001)  # 1 mm
# o.set_config('seed:droplet_diameter_sigma',0.0008) # 0.8 mm

Seeding some particles

time = [reader_norkyst.start_time,
        reader_norkyst.start_time + timedelta(hours=1)]
o.seed_elements(lon=4.5, lat=62.0, z='seafloor', radius=0, number=3000,
                time=time, oil_type='GENERIC DIESEL')
20:22:42 INFO    opendrift.models.openoil.openoil:1658: Using uniform droplet size distribution between 0.0001 and 0.0005 m for elements seeded below sea surface.
20:22:42 INFO    opendrift.models.openoil.adios.dirjs:86: Querying ADIOS database for oil: GENERIC DIESEL
20:22:42 INFO    opendrift.models.openoil.openoil:1715: Using density 855.84763 and viscosity 8.644570579524653e-06 of oiltype GENERIC DIESEL
20:22:42 INFO    opendrift.models.basemodel.environment:218: Adding a dynamical landmask with max. priority based on assumed maximum speed of 1.3 m/s. Adding a customised landmask may be faster...
20:22:46 INFO    opendrift.models.basemodel.environment:245: Fallback values will be used for the following variables which have no readers:
20:22:46 INFO    opendrift.models.basemodel.environment:248:    x_wind: 3.000000
20:22:46 INFO    opendrift.models.basemodel.environment:248:    y_wind: 7.000000
20:22:46 INFO    opendrift.models.basemodel.environment:248:    upward_sea_water_velocity: 0.000000
20:22:46 INFO    opendrift.models.basemodel.environment:248:    sea_surface_wave_significant_height: 0.000000
20:22:46 INFO    opendrift.models.basemodel.environment:248:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
20:22:46 INFO    opendrift.models.basemodel.environment:248:    sea_surface_wave_stokes_drift_y_velocity: 0.000000
20:22:46 INFO    opendrift.models.basemodel.environment:248:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
20:22:46 INFO    opendrift.models.basemodel.environment:248:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
20:22:46 INFO    opendrift.models.basemodel.environment:248:    sea_ice_area_fraction: 0.000000
20:22:46 INFO    opendrift.models.basemodel.environment:248:    sea_ice_x_velocity: 0.000000
20:22:46 INFO    opendrift.models.basemodel.environment:248:    sea_ice_y_velocity: 0.000000
20:22:46 INFO    opendrift.models.basemodel.environment:248:    ocean_vertical_diffusivity: 0.020000
20:22:46 INFO    opendrift.models.basemodel.environment:248:    ocean_mixed_layer_thickness: 50.000000
20:22:46 WARNING opendrift.models.basemodel.environment:458: Simulation has no simulation_extent, cannot check reader coverage

Running model with a small timestep to resolve the boyant rising

o.run(duration=timedelta(hours=6), time_step=60, time_step_output=60)
20:22:48 INFO    opendrift.models.basemodel:935: Using existing reader for land_binary_mask
20:22:48 INFO    opendrift.models.basemodel:946: All points are in ocean
20:22:48 INFO    opendrift.models.openoil.openoil:684: Oil-water surface tension is 0.030314 Nm
20:22:48 INFO    opendrift.models.openoil.openoil:697: Max water fraction not available for GENERIC DIESEL, using default
20:22:48 WARNING opendrift.models.basemodel:729: 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.
20:22:48 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:00:00 - step 1 of 360 - 50 active elements (0 deactivated)
20:22:48 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:01:00 - step 2 of 360 - 100 active elements (0 deactivated)
20:22:48 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:02:00 - step 3 of 360 - 150 active elements (0 deactivated)
20:22:48 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:03:00 - step 4 of 360 - 200 active elements (0 deactivated)
20:22:48 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:04:00 - step 5 of 360 - 250 active elements (0 deactivated)
20:22:48 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:05:00 - step 6 of 360 - 300 active elements (0 deactivated)
20:22:48 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:06:00 - step 7 of 360 - 350 active elements (0 deactivated)
20:22:48 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:07:00 - step 8 of 360 - 400 active elements (0 deactivated)
20:22:48 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:08:00 - step 9 of 360 - 450 active elements (0 deactivated)
20:22:49 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:09:00 - step 10 of 360 - 500 active elements (0 deactivated)
20:22:49 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:10:00 - step 11 of 360 - 550 active elements (0 deactivated)
20:22:49 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:11:00 - step 12 of 360 - 600 active elements (0 deactivated)
20:22:49 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:12:00 - step 13 of 360 - 650 active elements (0 deactivated)
20:22:49 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:13:00 - step 14 of 360 - 700 active elements (0 deactivated)
20:22:49 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:14:00 - step 15 of 360 - 750 active elements (0 deactivated)
20:22:49 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:15:00 - step 16 of 360 - 800 active elements (0 deactivated)
20:22:49 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:16:00 - step 17 of 360 - 850 active elements (0 deactivated)
20:22:49 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:17:00 - step 18 of 360 - 900 active elements (0 deactivated)
20:22:49 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:18:00 - step 19 of 360 - 950 active elements (0 deactivated)
20:22:49 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:19:00 - step 20 of 360 - 1000 active elements (0 deactivated)
20:22:50 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:20:00 - step 21 of 360 - 1050 active elements (0 deactivated)
20:22:50 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:21:00 - step 22 of 360 - 1100 active elements (0 deactivated)
20:22:50 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:22:00 - step 23 of 360 - 1150 active elements (0 deactivated)
20:22:50 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:23:00 - step 24 of 360 - 1200 active elements (0 deactivated)
20:22:50 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:24:00 - step 25 of 360 - 1250 active elements (0 deactivated)
20:22:50 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:25:00 - step 26 of 360 - 1300 active elements (0 deactivated)
20:22:50 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:26:00 - step 27 of 360 - 1350 active elements (0 deactivated)
20:22:50 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:27:00 - step 28 of 360 - 1400 active elements (0 deactivated)
20:22:50 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:28:00 - step 29 of 360 - 1450 active elements (0 deactivated)
20:22:50 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:29:00 - step 30 of 360 - 1500 active elements (0 deactivated)
20:22:51 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:30:00 - step 31 of 360 - 1550 active elements (0 deactivated)
20:22:51 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:31:00 - step 32 of 360 - 1600 active elements (0 deactivated)
20:22:51 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:32:00 - step 33 of 360 - 1650 active elements (0 deactivated)
20:22:51 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:33:00 - step 34 of 360 - 1700 active elements (0 deactivated)
20:22:51 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:34:00 - step 35 of 360 - 1750 active elements (0 deactivated)
20:22:51 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:35:00 - step 36 of 360 - 1800 active elements (0 deactivated)
20:22:51 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:36:00 - step 37 of 360 - 1850 active elements (0 deactivated)
20:22:51 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:37:00 - step 38 of 360 - 1900 active elements (0 deactivated)
20:22:51 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:38:00 - step 39 of 360 - 1950 active elements (0 deactivated)
20:22:52 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:39:00 - step 40 of 360 - 2000 active elements (0 deactivated)
20:22:52 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:40:00 - step 41 of 360 - 2050 active elements (0 deactivated)
20:22:52 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:41:00 - step 42 of 360 - 2100 active elements (0 deactivated)
20:22:52 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:42:00 - step 43 of 360 - 2150 active elements (0 deactivated)
20:22:52 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:43:00 - step 44 of 360 - 2200 active elements (0 deactivated)
20:22:52 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:44:00 - step 45 of 360 - 2250 active elements (0 deactivated)
20:22:52 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:45:00 - step 46 of 360 - 2300 active elements (0 deactivated)
20:22:52 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:46:00 - step 47 of 360 - 2350 active elements (0 deactivated)
20:22:53 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:47:00 - step 48 of 360 - 2400 active elements (0 deactivated)
20:22:53 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:48:00 - step 49 of 360 - 2450 active elements (0 deactivated)
20:22:53 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:49:00 - step 50 of 360 - 2500 active elements (0 deactivated)
20:22:53 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:50:00 - step 51 of 360 - 2550 active elements (0 deactivated)
20:22:53 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:51:00 - step 52 of 360 - 2600 active elements (0 deactivated)
20:22:53 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:52:00 - step 53 of 360 - 2650 active elements (0 deactivated)
20:22:53 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:53:00 - step 54 of 360 - 2700 active elements (0 deactivated)
20:22:53 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:54:00 - step 55 of 360 - 2750 active elements (0 deactivated)
20:22:54 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:55:00 - step 56 of 360 - 2800 active elements (0 deactivated)
20:22:54 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:56:00 - step 57 of 360 - 2850 active elements (0 deactivated)
20:22:54 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:57:00 - step 58 of 360 - 2900 active elements (0 deactivated)
20:22:54 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:58:00 - step 59 of 360 - 2950 active elements (0 deactivated)
20:22:54 INFO    opendrift.models.basemodel:2035: 2016-01-14 00:59:00 - step 60 of 360 - 3000 active elements (0 deactivated)
20:22:54 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:00:00 - step 61 of 360 - 3000 active elements (0 deactivated)
20:22:54 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:01:00 - step 62 of 360 - 3000 active elements (0 deactivated)
20:22:55 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:02:00 - step 63 of 360 - 3000 active elements (0 deactivated)
20:22:55 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:03:00 - step 64 of 360 - 3000 active elements (0 deactivated)
20:22:55 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:04:00 - step 65 of 360 - 3000 active elements (0 deactivated)
20:22:55 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:05:00 - step 66 of 360 - 3000 active elements (0 deactivated)
20:22:55 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:06:00 - step 67 of 360 - 3000 active elements (0 deactivated)
20:22:55 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:07:00 - step 68 of 360 - 3000 active elements (0 deactivated)
20:22:55 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:08:00 - step 69 of 360 - 3000 active elements (0 deactivated)
20:22:55 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:09:00 - step 70 of 360 - 3000 active elements (0 deactivated)
20:22:56 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:10:00 - step 71 of 360 - 3000 active elements (0 deactivated)
20:22:56 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:11:00 - step 72 of 360 - 3000 active elements (0 deactivated)
20:22:56 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:12:00 - step 73 of 360 - 3000 active elements (0 deactivated)
20:22:56 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:13:00 - step 74 of 360 - 3000 active elements (0 deactivated)
20:22:56 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:14:00 - step 75 of 360 - 3000 active elements (0 deactivated)
20:22:56 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:15:00 - step 76 of 360 - 3000 active elements (0 deactivated)
20:22:56 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:16:00 - step 77 of 360 - 3000 active elements (0 deactivated)
20:22:56 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:17:00 - step 78 of 360 - 3000 active elements (0 deactivated)
20:22:57 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:18:00 - step 79 of 360 - 3000 active elements (0 deactivated)
20:22:57 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:19:00 - step 80 of 360 - 3000 active elements (0 deactivated)
20:22:57 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:20:00 - step 81 of 360 - 3000 active elements (0 deactivated)
20:22:57 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:21:00 - step 82 of 360 - 3000 active elements (0 deactivated)
20:22:57 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:22:00 - step 83 of 360 - 3000 active elements (0 deactivated)
20:22:57 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:23:00 - step 84 of 360 - 3000 active elements (0 deactivated)
20:22:57 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:24:00 - step 85 of 360 - 3000 active elements (0 deactivated)
20:22:57 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:25:00 - step 86 of 360 - 3000 active elements (0 deactivated)
20:22:58 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:26:00 - step 87 of 360 - 3000 active elements (0 deactivated)
20:22:58 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:27:00 - step 88 of 360 - 3000 active elements (0 deactivated)
20:22:58 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:28:00 - step 89 of 360 - 3000 active elements (0 deactivated)
20:22:58 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:29:00 - step 90 of 360 - 3000 active elements (0 deactivated)
20:22:58 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:30:00 - step 91 of 360 - 3000 active elements (0 deactivated)
20:22:58 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:31:00 - step 92 of 360 - 3000 active elements (0 deactivated)
20:22:58 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:32:00 - step 93 of 360 - 3000 active elements (0 deactivated)
20:22:58 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:33:00 - step 94 of 360 - 3000 active elements (0 deactivated)
20:22:58 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:34:00 - step 95 of 360 - 3000 active elements (0 deactivated)
20:22:59 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:35:00 - step 96 of 360 - 3000 active elements (0 deactivated)
20:22:59 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:36:00 - step 97 of 360 - 3000 active elements (0 deactivated)
20:22:59 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:37:00 - step 98 of 360 - 3000 active elements (0 deactivated)
20:22:59 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:38:00 - step 99 of 360 - 3000 active elements (0 deactivated)
20:22:59 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:39:00 - step 100 of 360 - 3000 active elements (0 deactivated)
20:22:59 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:40:00 - step 101 of 360 - 3000 active elements (0 deactivated)
20:22:59 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:41:00 - step 102 of 360 - 3000 active elements (0 deactivated)
20:22:59 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:42:00 - step 103 of 360 - 3000 active elements (0 deactivated)
20:22:59 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:43:00 - step 104 of 360 - 3000 active elements (0 deactivated)
20:23:00 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:44:00 - step 105 of 360 - 3000 active elements (0 deactivated)
20:23:00 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:45:00 - step 106 of 360 - 3000 active elements (0 deactivated)
20:23:00 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:46:00 - step 107 of 360 - 3000 active elements (0 deactivated)
20:23:00 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:47:00 - step 108 of 360 - 3000 active elements (0 deactivated)
20:23:00 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:48:00 - step 109 of 360 - 3000 active elements (0 deactivated)
20:23:00 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:49:00 - step 110 of 360 - 3000 active elements (0 deactivated)
20:23:00 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:50:00 - step 111 of 360 - 3000 active elements (0 deactivated)
20:23:00 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:51:00 - step 112 of 360 - 3000 active elements (0 deactivated)
20:23:00 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:52:00 - step 113 of 360 - 3000 active elements (0 deactivated)
20:23:00 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:53:00 - step 114 of 360 - 3000 active elements (0 deactivated)
20:23:01 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:54:00 - step 115 of 360 - 3000 active elements (0 deactivated)
20:23:01 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:55:00 - step 116 of 360 - 3000 active elements (0 deactivated)
20:23:01 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:56:00 - step 117 of 360 - 3000 active elements (0 deactivated)
20:23:01 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:57:00 - step 118 of 360 - 3000 active elements (0 deactivated)
20:23:01 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:58:00 - step 119 of 360 - 3000 active elements (0 deactivated)
20:23:01 INFO    opendrift.models.basemodel:2035: 2016-01-14 01:59:00 - step 120 of 360 - 3000 active elements (0 deactivated)
20:23:01 INFO    opendrift.models.basemodel:2035: 2016-01-14 02:00:00 - step 121 of 360 - 3000 active elements (0 deactivated)
20:23:02 INFO    opendrift.models.basemodel:2035: 2016-01-14 02:01:00 - step 122 of 360 - 3000 active elements (0 deactivated)
20:23:02 WARNING opendrift.models.basemodel:2327: Missing variables: ['x_sea_water_velocity', 'y_sea_water_velocity']
20:23:02 WARNING opendrift.models.basemodel:2133: The simulation stopped before requested end time was reached.
20:23:02 INFO    opendrift.models.basemodel:2135: ========================
20:23:02 INFO    opendrift.models.basemodel:2136: End of simulation:
20:23:02 INFO    opendrift.models.basemodel:2137: No more active or scheduled elements, quitting.
20:23:02 INFO    opendrift.models.basemodel:2138: Traceback (most recent call last):
  File "/root/project/opendrift/models/basemodel/__init__.py", line 2100, in run
    raise ValueError(
ValueError: No more active or scheduled elements, quitting.

20:23:02 INFO    opendrift.models.basemodel:2139: "Missing variables: ['x_sea_water_velocity', 'y_sea_water_velocity']", 'The simulation stopped before requested end time was reached.'

20:23:02 INFO    opendrift.models.basemodel:2142: ========================

Print and plot results

print(o)

o.animation_profile(markersize='z', color='z')
===========================
--------------------
Reader performance:
--------------------
global_landmask
 0:00:00.0  total
 0:00:00.0  preparing
 0:00:00.0  reading
 0:00:00.0  masking
--------------------
Performance:
   19.6 total time
    4.4 configuration
    1.3 preparing main loop
      0.0 moving elements to ocean
   13.8 main loop
      5.4 updating elements
        0.0 oil weathering
          0.0 updating viscosities
          0.0 updating densities
          0.0 evaporation
          0.0 emulsification
          0.0 dispersion
        4.9 vertical mixing
    0.0 cleaning up
--------------------
===========================
Model:  OpenOil     (OpenDrift version 1.12.0)
        0 active Oil particles  (3000 deactivated, 0 scheduled)
-------------------
Environment variables:
  -----
  land_binary_mask
     1) global_landmask
  -----
Readers not added for the following variables:
  ocean_mixed_layer_thickness
  ocean_vertical_diffusivity
  sea_floor_depth_below_sea_level
  sea_ice_area_fraction
  sea_ice_x_velocity
  sea_ice_y_velocity
  sea_surface_height
  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_water_salinity
  sea_water_temperature
  upward_sea_water_velocity
  x_sea_water_velocity
  x_wind
  y_sea_water_velocity
  y_wind

Discarded readers:
  /root/project/tests/test_data/14Jan2016_NorKyst_z_3d/NorKyst-800m_ZDEPTHS_his_00_3Dsubset.nc (ends before simuation is finished)

Time:
        Start: 2016-01-14 00:00:00 UTC
        Present: 2016-01-14 02:01:00 UTC
        Calculation steps: 122 * 0:01:00 - total time: 2:01:00
        Output steps: 123 * 0:01:00
-------------------
"Missing variables: ['x_sea_water_velocity', 'y_sea_water_velocity']", 'The simulation stopped before requested end time was reached.'
===========================

20:23:02 INFO    opendrift.models.basemodel:4613: Saving animation to /root/project/docs/source/gallery/animations/example_oilspill_seafloor_0.gif...
20:23:18 INFO    opendrift.models.basemodel:3251: Time to make animation: 0:00:16.166482
../_images/example_oilspill_seafloor_0.gif
o.animate_vertical_distribution(bins=30, subsamplingstep=5)
20:23:20 INFO    opendrift.models.basemodel:4613: Saving animation to /root/project/docs/source/gallery/animations/example_oilspill_seafloor_1.gif...
20:23:42 INFO    opendrift.models.oceandrift:713: Time to make animation: 0:00:24.311407
../_images/example_oilspill_seafloor_1.gif
o.plot_oil_budget()
GENERIC DIESEL (855.8 kg/m3) - 2016-01-14 00:00 to 2016-01-14 02:01

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

Gallery generated by Sphinx-Gallery