Oil film thickness

from datetime import datetime, timedelta
import matplotlib.pyplot as plt
import numpy as np
from opendrift.models.openoil import OpenOil


number = 10000
timestep = timedelta(minutes=10)
timestep_output = timedelta(minutes=60)
duration = timedelta(hours=20)
mass_oil = 2000  # mass oil per particle
oil_type = 'GENERIC DIESEL'
#oil_type = 'GENERIC BUNKER C'

First run, where surface oil thickness is updated

o1 = OpenOil(loglevel=20, weathering_model='noaa')
15:17:37 INFO    opendrift:568: OpenDriftSimulation initialised (version 1.14.3 / v1.14.3-4-ga7b402c)

Northwards wind, eastwards current

o1.set_config('environment:fallback:land_binary_mask', 0)
o1.set_config('environment:fallback:x_wind', 0)
o1.set_config('environment:fallback:y_wind', 7)
o1.set_config('environment:fallback:sea_surface_wave_stokes_drift_x_velocity', 0)
o1.set_config('environment:fallback:sea_surface_wave_stokes_drift_y_velocity', .3)
o1.set_config('environment:fallback:x_sea_water_velocity', .1)
o1.set_config('environment:fallback:y_sea_water_velocity', 0)

Using Johansen droplet spectrum, which depends on oil film thickness

o1.set_config('wave_entrainment:droplet_size_distribution',
             'Johansen et al. (2015)')
o1.set_config('drift:wind_uncertainty', 2)
o1.set_config('drift:current_uncertainty', .1)
o1.set_config('processes:dispersion', False)
o1.set_config('processes:update_oilfilm_thickness', True)

o1.seed_elements(lon=4.5, lat=60, number=number,
                mass_oil=mass_oil, radius=1000,
                oil_type=oil_type,
                time=datetime.now())
o1.run(time_step=timestep, time_step_output=timestep_output,
       duration=duration)
15:17:37 INFO    opendrift.models.openoil.adios.dirjs:86: Querying ADIOS database for oil: GENERIC DIESEL
15:17:37 INFO    opendrift.models.openoil.openoil:1717: Using density 855.84763 and viscosity 8.644570579524623e-06 of oiltype GENERIC DIESEL
15:17:37 INFO    opendrift.models.basemodel.environment:206: Adding a global landmask from GSHHG
15:17:42 INFO    opendrift.models.basemodel.environment:229: Fallback values will be used for the following variables which have no readers:
15:17:42 INFO    opendrift.models.basemodel.environment:232:    x_sea_water_velocity: 0.100000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    y_sea_water_velocity: 0.000000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    x_wind: 0.000000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    y_wind: 7.000000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_height: 0.000000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    upward_sea_water_velocity: 0.000000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_significant_height: 0.000000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_y_velocity: 0.300000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    sea_ice_area_fraction: 0.000000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    sea_ice_x_velocity: 0.000000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    sea_ice_y_velocity: 0.000000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    sea_water_temperature: 10.000000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    sea_water_salinity: 34.000000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    sea_floor_depth_below_sea_level: 10000.000000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    ocean_vertical_diffusivity: 0.020000
15:17:42 INFO    opendrift.models.basemodel.environment:232:    ocean_mixed_layer_thickness: 50.000000
15:17:42 INFO    opendrift:1778: Skipping environment variable upward_sea_water_velocity because of condition ['drift:vertical_advection', 'is', False]
15:17:42 INFO    opendrift:1789: Storing previous values of element property lon because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
15:17:42 INFO    opendrift:1789: Storing previous values of element property lat because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
15:17:42 INFO    opendrift:944: Using existing reader for land_binary_mask to move elements to ocean
15:17:42 INFO    opendrift:974: All points are in ocean
15:17:42 INFO    opendrift.models.openoil.openoil:692: Oil-water surface tension is 0.030314 Nm
15:17:42 INFO    opendrift.models.openoil.openoil:705: Max water fraction not available for GENERIC DIESEL, using default
15:17:42 INFO    opendrift:2085: 2025-09-12 15:17:37.943161 - step 1 of 120 - 10000 active elements (0 deactivated)
15:17:42 INFO    opendrift:2085: 2025-09-12 15:27:37.943161 - step 2 of 120 - 10000 active elements (0 deactivated)
15:17:43 INFO    opendrift:2085: 2025-09-12 15:37:37.943161 - step 3 of 120 - 10000 active elements (0 deactivated)
15:17:43 INFO    opendrift:2085: 2025-09-12 15:47:37.943161 - step 4 of 120 - 10000 active elements (0 deactivated)
15:17:43 INFO    opendrift:2085: 2025-09-12 15:57:37.943161 - step 5 of 120 - 10000 active elements (0 deactivated)
15:17:43 INFO    opendrift:2085: 2025-09-12 16:07:37.943161 - step 6 of 120 - 10000 active elements (0 deactivated)
15:17:43 INFO    opendrift:2085: 2025-09-12 16:17:37.943161 - step 7 of 120 - 10000 active elements (0 deactivated)
15:17:44 INFO    opendrift:2085: 2025-09-12 16:27:37.943161 - step 8 of 120 - 10000 active elements (0 deactivated)
15:17:44 INFO    opendrift:2085: 2025-09-12 16:37:37.943161 - step 9 of 120 - 10000 active elements (0 deactivated)
15:17:44 INFO    opendrift:2085: 2025-09-12 16:47:37.943161 - step 10 of 120 - 10000 active elements (0 deactivated)
15:17:44 INFO    opendrift:2085: 2025-09-12 16:57:37.943161 - step 11 of 120 - 10000 active elements (0 deactivated)
15:17:44 INFO    opendrift:2085: 2025-09-12 17:07:37.943161 - step 12 of 120 - 10000 active elements (0 deactivated)
15:17:44 INFO    opendrift:2085: 2025-09-12 17:17:37.943161 - step 13 of 120 - 10000 active elements (0 deactivated)
15:17:45 INFO    opendrift:2085: 2025-09-12 17:27:37.943161 - step 14 of 120 - 10000 active elements (0 deactivated)
15:17:45 INFO    opendrift:2085: 2025-09-12 17:37:37.943161 - step 15 of 120 - 10000 active elements (0 deactivated)
15:17:45 INFO    opendrift:2085: 2025-09-12 17:47:37.943161 - step 16 of 120 - 10000 active elements (0 deactivated)
15:17:45 INFO    opendrift:2085: 2025-09-12 17:57:37.943161 - step 17 of 120 - 10000 active elements (0 deactivated)
15:17:45 INFO    opendrift:2085: 2025-09-12 18:07:37.943161 - step 18 of 120 - 10000 active elements (0 deactivated)
15:17:45 INFO    opendrift:2085: 2025-09-12 18:17:37.943161 - step 19 of 120 - 10000 active elements (0 deactivated)
15:17:46 INFO    opendrift:2085: 2025-09-12 18:27:37.943161 - step 20 of 120 - 10000 active elements (0 deactivated)
15:17:46 INFO    opendrift:2085: 2025-09-12 18:37:37.943161 - step 21 of 120 - 10000 active elements (0 deactivated)
15:17:46 INFO    opendrift:2085: 2025-09-12 18:47:37.943161 - step 22 of 120 - 10000 active elements (0 deactivated)
15:17:46 INFO    opendrift:2085: 2025-09-12 18:57:37.943161 - step 23 of 120 - 10000 active elements (0 deactivated)
15:17:46 INFO    opendrift:2085: 2025-09-12 19:07:37.943161 - step 24 of 120 - 10000 active elements (0 deactivated)
15:17:46 INFO    opendrift:2085: 2025-09-12 19:17:37.943161 - step 25 of 120 - 10000 active elements (0 deactivated)
15:17:47 INFO    opendrift:2085: 2025-09-12 19:27:37.943161 - step 26 of 120 - 10000 active elements (0 deactivated)
15:17:47 INFO    opendrift:2085: 2025-09-12 19:37:37.943161 - step 27 of 120 - 10000 active elements (0 deactivated)
15:17:47 INFO    opendrift:2085: 2025-09-12 19:47:37.943161 - step 28 of 120 - 10000 active elements (0 deactivated)
15:17:47 INFO    opendrift:2085: 2025-09-12 19:57:37.943161 - step 29 of 120 - 10000 active elements (0 deactivated)
15:17:47 INFO    opendrift:2085: 2025-09-12 20:07:37.943161 - step 30 of 120 - 10000 active elements (0 deactivated)
15:17:47 INFO    opendrift:2085: 2025-09-12 20:17:37.943161 - step 31 of 120 - 10000 active elements (0 deactivated)
15:17:47 INFO    opendrift:2085: 2025-09-12 20:27:37.943161 - step 32 of 120 - 10000 active elements (0 deactivated)
15:17:48 INFO    opendrift:2085: 2025-09-12 20:37:37.943161 - step 33 of 120 - 10000 active elements (0 deactivated)
15:17:48 INFO    opendrift:2085: 2025-09-12 20:47:37.943161 - step 34 of 120 - 10000 active elements (0 deactivated)
15:17:48 INFO    opendrift:2085: 2025-09-12 20:57:37.943161 - step 35 of 120 - 10000 active elements (0 deactivated)
15:17:48 INFO    opendrift:2085: 2025-09-12 21:07:37.943161 - step 36 of 120 - 10000 active elements (0 deactivated)
15:17:48 INFO    opendrift:2085: 2025-09-12 21:17:37.943161 - step 37 of 120 - 10000 active elements (0 deactivated)
15:17:48 INFO    opendrift:2085: 2025-09-12 21:27:37.943161 - step 38 of 120 - 10000 active elements (0 deactivated)
15:17:49 INFO    opendrift:2085: 2025-09-12 21:37:37.943161 - step 39 of 120 - 10000 active elements (0 deactivated)
15:17:49 INFO    opendrift:2085: 2025-09-12 21:47:37.943161 - step 40 of 120 - 10000 active elements (0 deactivated)
15:17:49 INFO    opendrift:2085: 2025-09-12 21:57:37.943161 - step 41 of 120 - 10000 active elements (0 deactivated)
15:17:49 INFO    opendrift:2085: 2025-09-12 22:07:37.943161 - step 42 of 120 - 10000 active elements (0 deactivated)
15:17:49 INFO    opendrift:2085: 2025-09-12 22:17:37.943161 - step 43 of 120 - 10000 active elements (0 deactivated)
15:17:49 INFO    opendrift:2085: 2025-09-12 22:27:37.943161 - step 44 of 120 - 10000 active elements (0 deactivated)
15:17:49 INFO    opendrift:2085: 2025-09-12 22:37:37.943161 - step 45 of 120 - 10000 active elements (0 deactivated)
15:17:50 INFO    opendrift:2085: 2025-09-12 22:47:37.943161 - step 46 of 120 - 10000 active elements (0 deactivated)
15:17:50 INFO    opendrift:2085: 2025-09-12 22:57:37.943161 - step 47 of 120 - 10000 active elements (0 deactivated)
15:17:50 INFO    opendrift:2085: 2025-09-12 23:07:37.943161 - step 48 of 120 - 10000 active elements (0 deactivated)
15:17:50 INFO    opendrift:2085: 2025-09-12 23:17:37.943161 - step 49 of 120 - 10000 active elements (0 deactivated)
15:17:50 INFO    opendrift:2085: 2025-09-12 23:27:37.943161 - step 50 of 120 - 10000 active elements (0 deactivated)
15:17:50 INFO    opendrift:2085: 2025-09-12 23:37:37.943161 - step 51 of 120 - 10000 active elements (0 deactivated)
15:17:51 INFO    opendrift:2085: 2025-09-12 23:47:37.943161 - step 52 of 120 - 10000 active elements (0 deactivated)
15:17:51 INFO    opendrift:2085: 2025-09-12 23:57:37.943161 - step 53 of 120 - 10000 active elements (0 deactivated)
15:17:51 INFO    opendrift:2085: 2025-09-13 00:07:37.943161 - step 54 of 120 - 10000 active elements (0 deactivated)
15:17:51 INFO    opendrift:2085: 2025-09-13 00:17:37.943161 - step 55 of 120 - 10000 active elements (0 deactivated)
15:17:51 INFO    opendrift:2085: 2025-09-13 00:27:37.943161 - step 56 of 120 - 10000 active elements (0 deactivated)
15:17:51 INFO    opendrift:2085: 2025-09-13 00:37:37.943161 - step 57 of 120 - 10000 active elements (0 deactivated)
15:17:51 INFO    opendrift:2085: 2025-09-13 00:47:37.943161 - step 58 of 120 - 10000 active elements (0 deactivated)
15:17:52 INFO    opendrift:2085: 2025-09-13 00:57:37.943161 - step 59 of 120 - 10000 active elements (0 deactivated)
15:17:52 INFO    opendrift:2085: 2025-09-13 01:07:37.943161 - step 60 of 120 - 10000 active elements (0 deactivated)
15:17:52 INFO    opendrift:2085: 2025-09-13 01:17:37.943161 - step 61 of 120 - 10000 active elements (0 deactivated)
15:17:52 INFO    opendrift:2085: 2025-09-13 01:27:37.943161 - step 62 of 120 - 10000 active elements (0 deactivated)
15:17:52 INFO    opendrift:2085: 2025-09-13 01:37:37.943161 - step 63 of 120 - 10000 active elements (0 deactivated)
15:17:52 INFO    opendrift:2085: 2025-09-13 01:47:37.943161 - step 64 of 120 - 10000 active elements (0 deactivated)
15:17:52 INFO    opendrift:2085: 2025-09-13 01:57:37.943161 - step 65 of 120 - 10000 active elements (0 deactivated)
15:17:53 INFO    opendrift:2085: 2025-09-13 02:07:37.943161 - step 66 of 120 - 10000 active elements (0 deactivated)
15:17:53 INFO    opendrift:2085: 2025-09-13 02:17:37.943161 - step 67 of 120 - 10000 active elements (0 deactivated)
15:17:53 INFO    opendrift:2085: 2025-09-13 02:27:37.943161 - step 68 of 120 - 10000 active elements (0 deactivated)
15:17:53 INFO    opendrift:2085: 2025-09-13 02:37:37.943161 - step 69 of 120 - 10000 active elements (0 deactivated)
15:17:53 INFO    opendrift:2085: 2025-09-13 02:47:37.943161 - step 70 of 120 - 10000 active elements (0 deactivated)
15:17:53 INFO    opendrift:2085: 2025-09-13 02:57:37.943161 - step 71 of 120 - 10000 active elements (0 deactivated)
15:17:53 INFO    opendrift:2085: 2025-09-13 03:07:37.943161 - step 72 of 120 - 10000 active elements (0 deactivated)
15:17:54 INFO    opendrift:2085: 2025-09-13 03:17:37.943161 - step 73 of 120 - 10000 active elements (0 deactivated)
15:17:54 INFO    opendrift:2085: 2025-09-13 03:27:37.943161 - step 74 of 120 - 10000 active elements (0 deactivated)
15:17:54 INFO    opendrift:2085: 2025-09-13 03:37:37.943161 - step 75 of 120 - 10000 active elements (0 deactivated)
15:17:54 INFO    opendrift:2085: 2025-09-13 03:47:37.943161 - step 76 of 120 - 10000 active elements (0 deactivated)
15:17:54 INFO    opendrift:2085: 2025-09-13 03:57:37.943161 - step 77 of 120 - 10000 active elements (0 deactivated)
15:17:54 INFO    opendrift:2085: 2025-09-13 04:07:37.943161 - step 78 of 120 - 10000 active elements (0 deactivated)
15:17:54 INFO    opendrift:2085: 2025-09-13 04:17:37.943161 - step 79 of 120 - 10000 active elements (0 deactivated)
15:17:55 INFO    opendrift:2085: 2025-09-13 04:27:37.943161 - step 80 of 120 - 10000 active elements (0 deactivated)
15:17:55 INFO    opendrift:2085: 2025-09-13 04:37:37.943161 - step 81 of 120 - 10000 active elements (0 deactivated)
15:17:55 INFO    opendrift:2085: 2025-09-13 04:47:37.943161 - step 82 of 120 - 10000 active elements (0 deactivated)
15:17:55 INFO    opendrift:2085: 2025-09-13 04:57:37.943161 - step 83 of 120 - 10000 active elements (0 deactivated)
15:17:55 INFO    opendrift:2085: 2025-09-13 05:07:37.943161 - step 84 of 120 - 10000 active elements (0 deactivated)
15:17:55 INFO    opendrift:2085: 2025-09-13 05:17:37.943161 - step 85 of 120 - 10000 active elements (0 deactivated)
15:17:56 INFO    opendrift:2085: 2025-09-13 05:27:37.943161 - step 86 of 120 - 10000 active elements (0 deactivated)
15:17:56 INFO    opendrift:2085: 2025-09-13 05:37:37.943161 - step 87 of 120 - 10000 active elements (0 deactivated)
15:17:56 INFO    opendrift:2085: 2025-09-13 05:47:37.943161 - step 88 of 120 - 10000 active elements (0 deactivated)
15:17:56 INFO    opendrift:2085: 2025-09-13 05:57:37.943161 - step 89 of 120 - 10000 active elements (0 deactivated)
15:17:56 INFO    opendrift:2085: 2025-09-13 06:07:37.943161 - step 90 of 120 - 10000 active elements (0 deactivated)
15:17:56 INFO    opendrift:2085: 2025-09-13 06:17:37.943161 - step 91 of 120 - 10000 active elements (0 deactivated)
15:17:56 INFO    opendrift:2085: 2025-09-13 06:27:37.943161 - step 92 of 120 - 10000 active elements (0 deactivated)
15:17:57 INFO    opendrift:2085: 2025-09-13 06:37:37.943161 - step 93 of 120 - 10000 active elements (0 deactivated)
15:17:57 INFO    opendrift:2085: 2025-09-13 06:47:37.943161 - step 94 of 120 - 10000 active elements (0 deactivated)
15:17:57 INFO    opendrift:2085: 2025-09-13 06:57:37.943161 - step 95 of 120 - 10000 active elements (0 deactivated)
15:17:57 INFO    opendrift:2085: 2025-09-13 07:07:37.943161 - step 96 of 120 - 10000 active elements (0 deactivated)
15:17:57 INFO    opendrift:2085: 2025-09-13 07:17:37.943161 - step 97 of 120 - 10000 active elements (0 deactivated)
15:17:57 INFO    opendrift:2085: 2025-09-13 07:27:37.943161 - step 98 of 120 - 10000 active elements (0 deactivated)
15:17:58 INFO    opendrift:2085: 2025-09-13 07:37:37.943161 - step 99 of 120 - 10000 active elements (0 deactivated)
15:17:58 INFO    opendrift:2085: 2025-09-13 07:47:37.943161 - step 100 of 120 - 10000 active elements (0 deactivated)
15:17:58 INFO    opendrift:2085: 2025-09-13 07:57:37.943161 - step 101 of 120 - 10000 active elements (0 deactivated)
15:17:58 INFO    opendrift:2085: 2025-09-13 08:07:37.943161 - step 102 of 120 - 10000 active elements (0 deactivated)
15:17:58 INFO    opendrift:2085: 2025-09-13 08:17:37.943161 - step 103 of 120 - 10000 active elements (0 deactivated)
15:17:58 INFO    opendrift:2085: 2025-09-13 08:27:37.943161 - step 104 of 120 - 10000 active elements (0 deactivated)
15:17:58 INFO    opendrift:2085: 2025-09-13 08:37:37.943161 - step 105 of 120 - 10000 active elements (0 deactivated)
15:17:59 INFO    opendrift:2085: 2025-09-13 08:47:37.943161 - step 106 of 120 - 10000 active elements (0 deactivated)
15:17:59 INFO    opendrift:2085: 2025-09-13 08:57:37.943161 - step 107 of 120 - 10000 active elements (0 deactivated)
15:17:59 INFO    opendrift:2085: 2025-09-13 09:07:37.943161 - step 108 of 120 - 10000 active elements (0 deactivated)
15:17:59 INFO    opendrift:2085: 2025-09-13 09:17:37.943161 - step 109 of 120 - 10000 active elements (0 deactivated)
15:17:59 INFO    opendrift:2085: 2025-09-13 09:27:37.943161 - step 110 of 120 - 10000 active elements (0 deactivated)
15:17:59 INFO    opendrift:2085: 2025-09-13 09:37:37.943161 - step 111 of 120 - 10000 active elements (0 deactivated)
15:18:00 INFO    opendrift:2085: 2025-09-13 09:47:37.943161 - step 112 of 120 - 10000 active elements (0 deactivated)
15:18:00 INFO    opendrift:2085: 2025-09-13 09:57:37.943161 - step 113 of 120 - 10000 active elements (0 deactivated)
15:18:00 INFO    opendrift:2085: 2025-09-13 10:07:37.943161 - step 114 of 120 - 10000 active elements (0 deactivated)
15:18:00 INFO    opendrift:2085: 2025-09-13 10:17:37.943161 - step 115 of 120 - 10000 active elements (0 deactivated)
15:18:00 INFO    opendrift:2085: 2025-09-13 10:27:37.943161 - step 116 of 120 - 10000 active elements (0 deactivated)
15:18:00 INFO    opendrift:2085: 2025-09-13 10:37:37.943161 - step 117 of 120 - 10000 active elements (0 deactivated)
15:18:01 INFO    opendrift:2085: 2025-09-13 10:47:37.943161 - step 118 of 120 - 10000 active elements (0 deactivated)
15:18:01 INFO    opendrift:2085: 2025-09-13 10:57:37.943161 - step 119 of 120 - 10000 active elements (0 deactivated)
15:18:01 INFO    opendrift:2085: 2025-09-13 11:07:37.943161 - step 120 of 120 - 10000 active elements (0 deactivated)
<xarray.Dataset> Size: 36MB
Dimensions:                                                                              (
                                                                                          trajectory: 10000,
                                                                                          time: 21)
Coordinates:
  * trajectory                                                                           (trajectory) int64 80kB ...
  * time                                                                                 (time) datetime64[ns] 168B ...
Data variables: (12/43)
    status                                                                               (trajectory, time) float32 840kB ...
    moving                                                                               (trajectory, time) float32 840kB ...
    age_seconds                                                                          (trajectory, time) float32 840kB ...
    origin_marker                                                                        (trajectory, time) float32 840kB ...
    lon                                                                                  (trajectory, time) float32 840kB ...
    lat                                                                                  (trajectory, time) float32 840kB ...
    ...                                                                                   ...
    sea_water_temperature                                                                (trajectory, time) float32 840kB ...
    sea_water_salinity                                                                   (trajectory, time) float32 840kB ...
    sea_floor_depth_below_sea_level                                                      (trajectory, time) float32 840kB ...
    ocean_vertical_diffusivity                                                           (trajectory, time) float32 840kB ...
    land_binary_mask                                                                     (trajectory, time) float32 840kB ...
    ocean_mixed_layer_thickness                                                          (trajectory, time) float32 840kB ...
Attributes: (12/165)
    Conventions:                                                             ...
    standard_name_vocabulary:                                                ...
    featureType:                                                             ...
    title:                                                                   ...
    summary:                                                                 ...
    keywords:                                                                ...
    ...                                                                                                               ...
    geospatial_lon_units:                                                    ...
    geospatial_lon_resolution:                                               ...
    runtime:                                                                 ...
    geospatial_vertical_min:                                                 ...
    geospatial_vertical_max:                                                 ...
    geospatial_vertical_positive:                                            ...


Animation shows how oil thickness evolves, and decreases due to evaporation and spreading

unitfactor=1e6  # show film thickness in micrometers
o1.animation(color='oil_film_thickness', fast=True,
             vmin=1e-7*unitfactor, vmax=1e-4*unitfactor,
             unitfactor=unitfactor, surface_only=True)
15:18:01 WARNING opendrift:2444: Plotting fast. This will make your plots less accurate.
15:18:03 INFO    opendrift:4612: Saving animation to /root/project/docs/source/gallery/animations/example_oil_thickness_0.gif...
15:18:19 INFO    opendrift:3053: Time to make animation: 0:00:17.961968
../_images/example_oil_thickness_0.gif

Second run, identical but without updating surface oil thickness

o2 = OpenOil(loglevel=20, weathering_model='noaa')
o2.set_config('environment:fallback:land_binary_mask', 0)
o2.set_config('environment:fallback:x_wind', 0)
o2.set_config('environment:fallback:y_wind', 7)
o2.set_config('environment:fallback:sea_surface_wave_stokes_drift_x_velocity', 0)
o2.set_config('environment:fallback:sea_surface_wave_stokes_drift_y_velocity', .3)
o2.set_config('environment:fallback:x_sea_water_velocity', .1)
o2.set_config('environment:fallback:y_sea_water_velocity', 0)

o2.set_config('wave_entrainment:droplet_size_distribution',
             'Johansen et al. (2015)')
o2.set_config('drift:wind_uncertainty', 2)
o2.set_config('drift:current_uncertainty', .1)
o2.set_config('processes:dispersion', False)
o2.set_config('processes:update_oilfilm_thickness', False)

o2.seed_elements(lon=4.5, lat=60, number=number,
                mass_oil=mass_oil, radius=1000,
                oil_type=oil_type,
                time=datetime.now())
o2.run(time_step=timestep, time_step_output=timestep_output,
       duration=duration)
15:18:19 INFO    opendrift:568: OpenDriftSimulation initialised (version 1.14.3 / v1.14.3-4-ga7b402c)
15:18:19 INFO    opendrift.models.openoil.adios.dirjs:86: Querying ADIOS database for oil: GENERIC DIESEL
15:18:19 INFO    opendrift.models.openoil.openoil:1717: Using density 855.84763 and viscosity 8.644570579524623e-06 of oiltype GENERIC DIESEL
15:18:19 INFO    opendrift.models.basemodel.environment:206: Adding a global landmask from GSHHG
15:18:19 INFO    opendrift.models.basemodel.environment:229: Fallback values will be used for the following variables which have no readers:
15:18:19 INFO    opendrift.models.basemodel.environment:232:    x_sea_water_velocity: 0.100000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    y_sea_water_velocity: 0.000000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    x_wind: 0.000000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    y_wind: 7.000000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    sea_surface_height: 0.000000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    upward_sea_water_velocity: 0.000000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_significant_height: 0.000000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_stokes_drift_y_velocity: 0.300000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    sea_ice_area_fraction: 0.000000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    sea_ice_x_velocity: 0.000000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    sea_ice_y_velocity: 0.000000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    sea_water_temperature: 10.000000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    sea_water_salinity: 34.000000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    sea_floor_depth_below_sea_level: 10000.000000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    ocean_vertical_diffusivity: 0.020000
15:18:19 INFO    opendrift.models.basemodel.environment:232:    ocean_mixed_layer_thickness: 50.000000
15:18:19 INFO    opendrift:1778: Skipping environment variable upward_sea_water_velocity because of condition ['drift:vertical_advection', 'is', False]
15:18:19 INFO    opendrift:1789: Storing previous values of element property lon because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
15:18:19 INFO    opendrift:1789: Storing previous values of element property lat because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
15:18:19 INFO    opendrift:944: Using existing reader for land_binary_mask to move elements to ocean
15:18:19 INFO    opendrift:974: All points are in ocean
15:18:19 INFO    opendrift.models.openoil.openoil:692: Oil-water surface tension is 0.030314 Nm
15:18:19 INFO    opendrift.models.openoil.openoil:705: Max water fraction not available for GENERIC DIESEL, using default
15:18:19 INFO    opendrift:2085: 2025-09-12 15:18:19.812677 - step 1 of 120 - 10000 active elements (0 deactivated)
15:18:20 INFO    opendrift:2085: 2025-09-12 15:28:19.812677 - step 2 of 120 - 10000 active elements (0 deactivated)
15:18:20 INFO    opendrift:2085: 2025-09-12 15:38:19.812677 - step 3 of 120 - 10000 active elements (0 deactivated)
15:18:20 INFO    opendrift:2085: 2025-09-12 15:48:19.812677 - step 4 of 120 - 10000 active elements (0 deactivated)
15:18:20 INFO    opendrift:2085: 2025-09-12 15:58:19.812677 - step 5 of 120 - 10000 active elements (0 deactivated)
15:18:21 INFO    opendrift:2085: 2025-09-12 16:08:19.812677 - step 6 of 120 - 10000 active elements (0 deactivated)
15:18:21 INFO    opendrift:2085: 2025-09-12 16:18:19.812677 - step 7 of 120 - 10000 active elements (0 deactivated)
15:18:21 INFO    opendrift:2085: 2025-09-12 16:28:19.812677 - step 8 of 120 - 10000 active elements (0 deactivated)
15:18:21 INFO    opendrift:2085: 2025-09-12 16:38:19.812677 - step 9 of 120 - 10000 active elements (0 deactivated)
15:18:21 INFO    opendrift:2085: 2025-09-12 16:48:19.812677 - step 10 of 120 - 10000 active elements (0 deactivated)
15:18:22 INFO    opendrift:2085: 2025-09-12 16:58:19.812677 - step 11 of 120 - 10000 active elements (0 deactivated)
15:18:22 INFO    opendrift:2085: 2025-09-12 17:08:19.812677 - step 12 of 120 - 10000 active elements (0 deactivated)
15:18:22 INFO    opendrift:2085: 2025-09-12 17:18:19.812677 - step 13 of 120 - 10000 active elements (0 deactivated)
15:18:22 INFO    opendrift:2085: 2025-09-12 17:28:19.812677 - step 14 of 120 - 10000 active elements (0 deactivated)
15:18:23 INFO    opendrift:2085: 2025-09-12 17:38:19.812677 - step 15 of 120 - 10000 active elements (0 deactivated)
15:18:23 INFO    opendrift:2085: 2025-09-12 17:48:19.812677 - step 16 of 120 - 10000 active elements (0 deactivated)
15:18:23 INFO    opendrift:2085: 2025-09-12 17:58:19.812677 - step 17 of 120 - 10000 active elements (0 deactivated)
15:18:23 INFO    opendrift:2085: 2025-09-12 18:08:19.812677 - step 18 of 120 - 10000 active elements (0 deactivated)
15:18:23 INFO    opendrift:2085: 2025-09-12 18:18:19.812677 - step 19 of 120 - 10000 active elements (0 deactivated)
15:18:24 INFO    opendrift:2085: 2025-09-12 18:28:19.812677 - step 20 of 120 - 10000 active elements (0 deactivated)
15:18:24 INFO    opendrift:2085: 2025-09-12 18:38:19.812677 - step 21 of 120 - 10000 active elements (0 deactivated)
15:18:24 INFO    opendrift:2085: 2025-09-12 18:48:19.812677 - step 22 of 120 - 10000 active elements (0 deactivated)
15:18:24 INFO    opendrift:2085: 2025-09-12 18:58:19.812677 - step 23 of 120 - 10000 active elements (0 deactivated)
15:18:24 INFO    opendrift:2085: 2025-09-12 19:08:19.812677 - step 24 of 120 - 10000 active elements (0 deactivated)
15:18:24 INFO    opendrift:2085: 2025-09-12 19:18:19.812677 - step 25 of 120 - 10000 active elements (0 deactivated)
15:18:25 INFO    opendrift:2085: 2025-09-12 19:28:19.812677 - step 26 of 120 - 10000 active elements (0 deactivated)
15:18:25 INFO    opendrift:2085: 2025-09-12 19:38:19.812677 - step 27 of 120 - 10000 active elements (0 deactivated)
15:18:25 INFO    opendrift:2085: 2025-09-12 19:48:19.812677 - step 28 of 120 - 10000 active elements (0 deactivated)
15:18:25 INFO    opendrift:2085: 2025-09-12 19:58:19.812677 - step 29 of 120 - 10000 active elements (0 deactivated)
15:18:25 INFO    opendrift:2085: 2025-09-12 20:08:19.812677 - step 30 of 120 - 10000 active elements (0 deactivated)
15:18:26 INFO    opendrift:2085: 2025-09-12 20:18:19.812677 - step 31 of 120 - 10000 active elements (0 deactivated)
15:18:26 INFO    opendrift:2085: 2025-09-12 20:28:19.812677 - step 32 of 120 - 10000 active elements (0 deactivated)
15:18:26 INFO    opendrift:2085: 2025-09-12 20:38:19.812677 - step 33 of 120 - 10000 active elements (0 deactivated)
15:18:26 INFO    opendrift:2085: 2025-09-12 20:48:19.812677 - step 34 of 120 - 10000 active elements (0 deactivated)
15:18:26 INFO    opendrift:2085: 2025-09-12 20:58:19.812677 - step 35 of 120 - 10000 active elements (0 deactivated)
15:18:26 INFO    opendrift:2085: 2025-09-12 21:08:19.812677 - step 36 of 120 - 10000 active elements (0 deactivated)
15:18:27 INFO    opendrift:2085: 2025-09-12 21:18:19.812677 - step 37 of 120 - 10000 active elements (0 deactivated)
15:18:27 INFO    opendrift:2085: 2025-09-12 21:28:19.812677 - step 38 of 120 - 10000 active elements (0 deactivated)
15:18:27 INFO    opendrift:2085: 2025-09-12 21:38:19.812677 - step 39 of 120 - 10000 active elements (0 deactivated)
15:18:27 INFO    opendrift:2085: 2025-09-12 21:48:19.812677 - step 40 of 120 - 10000 active elements (0 deactivated)
15:18:27 INFO    opendrift:2085: 2025-09-12 21:58:19.812677 - step 41 of 120 - 10000 active elements (0 deactivated)
15:18:27 INFO    opendrift:2085: 2025-09-12 22:08:19.812677 - step 42 of 120 - 10000 active elements (0 deactivated)
15:18:28 INFO    opendrift:2085: 2025-09-12 22:18:19.812677 - step 43 of 120 - 10000 active elements (0 deactivated)
15:18:28 INFO    opendrift:2085: 2025-09-12 22:28:19.812677 - step 44 of 120 - 10000 active elements (0 deactivated)
15:18:28 INFO    opendrift:2085: 2025-09-12 22:38:19.812677 - step 45 of 120 - 10000 active elements (0 deactivated)
15:18:28 INFO    opendrift:2085: 2025-09-12 22:48:19.812677 - step 46 of 120 - 10000 active elements (0 deactivated)
15:18:28 INFO    opendrift:2085: 2025-09-12 22:58:19.812677 - step 47 of 120 - 10000 active elements (0 deactivated)
15:18:29 INFO    opendrift:2085: 2025-09-12 23:08:19.812677 - step 48 of 120 - 10000 active elements (0 deactivated)
15:18:29 INFO    opendrift:2085: 2025-09-12 23:18:19.812677 - step 49 of 120 - 10000 active elements (0 deactivated)
15:18:29 INFO    opendrift:2085: 2025-09-12 23:28:19.812677 - step 50 of 120 - 10000 active elements (0 deactivated)
15:18:29 INFO    opendrift:2085: 2025-09-12 23:38:19.812677 - step 51 of 120 - 10000 active elements (0 deactivated)
15:18:29 INFO    opendrift:2085: 2025-09-12 23:48:19.812677 - step 52 of 120 - 10000 active elements (0 deactivated)
15:18:29 INFO    opendrift:2085: 2025-09-12 23:58:19.812677 - step 53 of 120 - 10000 active elements (0 deactivated)
15:18:30 INFO    opendrift:2085: 2025-09-13 00:08:19.812677 - step 54 of 120 - 10000 active elements (0 deactivated)
15:18:30 INFO    opendrift:2085: 2025-09-13 00:18:19.812677 - step 55 of 120 - 10000 active elements (0 deactivated)
15:18:30 INFO    opendrift:2085: 2025-09-13 00:28:19.812677 - step 56 of 120 - 10000 active elements (0 deactivated)
15:18:30 INFO    opendrift:2085: 2025-09-13 00:38:19.812677 - step 57 of 120 - 10000 active elements (0 deactivated)
15:18:30 INFO    opendrift:2085: 2025-09-13 00:48:19.812677 - step 58 of 120 - 10000 active elements (0 deactivated)
15:18:30 INFO    opendrift:2085: 2025-09-13 00:58:19.812677 - step 59 of 120 - 10000 active elements (0 deactivated)
15:18:31 INFO    opendrift:2085: 2025-09-13 01:08:19.812677 - step 60 of 120 - 10000 active elements (0 deactivated)
15:18:31 INFO    opendrift:2085: 2025-09-13 01:18:19.812677 - step 61 of 120 - 10000 active elements (0 deactivated)
15:18:31 INFO    opendrift:2085: 2025-09-13 01:28:19.812677 - step 62 of 120 - 10000 active elements (0 deactivated)
15:18:31 INFO    opendrift:2085: 2025-09-13 01:38:19.812677 - step 63 of 120 - 10000 active elements (0 deactivated)
15:18:31 INFO    opendrift:2085: 2025-09-13 01:48:19.812677 - step 64 of 120 - 10000 active elements (0 deactivated)
15:18:32 INFO    opendrift:2085: 2025-09-13 01:58:19.812677 - step 65 of 120 - 10000 active elements (0 deactivated)
15:18:32 INFO    opendrift:2085: 2025-09-13 02:08:19.812677 - step 66 of 120 - 10000 active elements (0 deactivated)
15:18:32 INFO    opendrift:2085: 2025-09-13 02:18:19.812677 - step 67 of 120 - 10000 active elements (0 deactivated)
15:18:32 INFO    opendrift:2085: 2025-09-13 02:28:19.812677 - step 68 of 120 - 10000 active elements (0 deactivated)
15:18:33 INFO    opendrift:2085: 2025-09-13 02:38:19.812677 - step 69 of 120 - 10000 active elements (0 deactivated)
15:18:33 INFO    opendrift:2085: 2025-09-13 02:48:19.812677 - step 70 of 120 - 10000 active elements (0 deactivated)
15:18:33 INFO    opendrift:2085: 2025-09-13 02:58:19.812677 - step 71 of 120 - 10000 active elements (0 deactivated)
15:18:33 INFO    opendrift:2085: 2025-09-13 03:08:19.812677 - step 72 of 120 - 10000 active elements (0 deactivated)
15:18:33 INFO    opendrift:2085: 2025-09-13 03:18:19.812677 - step 73 of 120 - 10000 active elements (0 deactivated)
15:18:34 INFO    opendrift:2085: 2025-09-13 03:28:19.812677 - step 74 of 120 - 10000 active elements (0 deactivated)
15:18:34 INFO    opendrift:2085: 2025-09-13 03:38:19.812677 - step 75 of 120 - 10000 active elements (0 deactivated)
15:18:34 INFO    opendrift:2085: 2025-09-13 03:48:19.812677 - step 76 of 120 - 10000 active elements (0 deactivated)
15:18:34 INFO    opendrift:2085: 2025-09-13 03:58:19.812677 - step 77 of 120 - 10000 active elements (0 deactivated)
15:18:34 INFO    opendrift:2085: 2025-09-13 04:08:19.812677 - step 78 of 120 - 10000 active elements (0 deactivated)
15:18:34 INFO    opendrift:2085: 2025-09-13 04:18:19.812677 - step 79 of 120 - 10000 active elements (0 deactivated)
15:18:35 INFO    opendrift:2085: 2025-09-13 04:28:19.812677 - step 80 of 120 - 10000 active elements (0 deactivated)
15:18:35 INFO    opendrift:2085: 2025-09-13 04:38:19.812677 - step 81 of 120 - 10000 active elements (0 deactivated)
15:18:35 INFO    opendrift:2085: 2025-09-13 04:48:19.812677 - step 82 of 120 - 10000 active elements (0 deactivated)
15:18:35 INFO    opendrift:2085: 2025-09-13 04:58:19.812677 - step 83 of 120 - 10000 active elements (0 deactivated)
15:18:35 INFO    opendrift:2085: 2025-09-13 05:08:19.812677 - step 84 of 120 - 10000 active elements (0 deactivated)
15:18:36 INFO    opendrift:2085: 2025-09-13 05:18:19.812677 - step 85 of 120 - 10000 active elements (0 deactivated)
15:18:36 INFO    opendrift:2085: 2025-09-13 05:28:19.812677 - step 86 of 120 - 10000 active elements (0 deactivated)
15:18:36 INFO    opendrift:2085: 2025-09-13 05:38:19.812677 - step 87 of 120 - 10000 active elements (0 deactivated)
15:18:36 INFO    opendrift:2085: 2025-09-13 05:48:19.812677 - step 88 of 120 - 10000 active elements (0 deactivated)
15:18:36 INFO    opendrift:2085: 2025-09-13 05:58:19.812677 - step 89 of 120 - 10000 active elements (0 deactivated)
15:18:37 INFO    opendrift:2085: 2025-09-13 06:08:19.812677 - step 90 of 120 - 10000 active elements (0 deactivated)
15:18:37 INFO    opendrift:2085: 2025-09-13 06:18:19.812677 - step 91 of 120 - 10000 active elements (0 deactivated)
15:18:37 INFO    opendrift:2085: 2025-09-13 06:28:19.812677 - step 92 of 120 - 10000 active elements (0 deactivated)
15:18:37 INFO    opendrift:2085: 2025-09-13 06:38:19.812677 - step 93 of 120 - 10000 active elements (0 deactivated)
15:18:37 INFO    opendrift:2085: 2025-09-13 06:48:19.812677 - step 94 of 120 - 10000 active elements (0 deactivated)
15:18:38 INFO    opendrift:2085: 2025-09-13 06:58:19.812677 - step 95 of 120 - 10000 active elements (0 deactivated)
15:18:38 INFO    opendrift:2085: 2025-09-13 07:08:19.812677 - step 96 of 120 - 10000 active elements (0 deactivated)
15:18:38 INFO    opendrift:2085: 2025-09-13 07:18:19.812677 - step 97 of 120 - 10000 active elements (0 deactivated)
15:18:38 INFO    opendrift:2085: 2025-09-13 07:28:19.812677 - step 98 of 120 - 10000 active elements (0 deactivated)
15:18:38 INFO    opendrift:2085: 2025-09-13 07:38:19.812677 - step 99 of 120 - 10000 active elements (0 deactivated)
15:18:39 INFO    opendrift:2085: 2025-09-13 07:48:19.812677 - step 100 of 120 - 10000 active elements (0 deactivated)
15:18:39 INFO    opendrift:2085: 2025-09-13 07:58:19.812677 - step 101 of 120 - 10000 active elements (0 deactivated)
15:18:39 INFO    opendrift:2085: 2025-09-13 08:08:19.812677 - step 102 of 120 - 10000 active elements (0 deactivated)
15:18:39 INFO    opendrift:2085: 2025-09-13 08:18:19.812677 - step 103 of 120 - 10000 active elements (0 deactivated)
15:18:39 INFO    opendrift:2085: 2025-09-13 08:28:19.812677 - step 104 of 120 - 10000 active elements (0 deactivated)
15:18:40 INFO    opendrift:2085: 2025-09-13 08:38:19.812677 - step 105 of 120 - 10000 active elements (0 deactivated)
15:18:40 INFO    opendrift:2085: 2025-09-13 08:48:19.812677 - step 106 of 120 - 10000 active elements (0 deactivated)
15:18:40 INFO    opendrift:2085: 2025-09-13 08:58:19.812677 - step 107 of 120 - 10000 active elements (0 deactivated)
15:18:40 INFO    opendrift:2085: 2025-09-13 09:08:19.812677 - step 108 of 120 - 10000 active elements (0 deactivated)
15:18:40 INFO    opendrift:2085: 2025-09-13 09:18:19.812677 - step 109 of 120 - 10000 active elements (0 deactivated)
15:18:41 INFO    opendrift:2085: 2025-09-13 09:28:19.812677 - step 110 of 120 - 10000 active elements (0 deactivated)
15:18:41 INFO    opendrift:2085: 2025-09-13 09:38:19.812677 - step 111 of 120 - 10000 active elements (0 deactivated)
15:18:41 INFO    opendrift:2085: 2025-09-13 09:48:19.812677 - step 112 of 120 - 10000 active elements (0 deactivated)
15:18:41 INFO    opendrift:2085: 2025-09-13 09:58:19.812677 - step 113 of 120 - 10000 active elements (0 deactivated)
15:18:41 INFO    opendrift:2085: 2025-09-13 10:08:19.812677 - step 114 of 120 - 10000 active elements (0 deactivated)
15:18:42 INFO    opendrift:2085: 2025-09-13 10:18:19.812677 - step 115 of 120 - 10000 active elements (0 deactivated)
15:18:42 INFO    opendrift:2085: 2025-09-13 10:28:19.812677 - step 116 of 120 - 10000 active elements (0 deactivated)
15:18:42 INFO    opendrift:2085: 2025-09-13 10:38:19.812677 - step 117 of 120 - 10000 active elements (0 deactivated)
15:18:42 INFO    opendrift:2085: 2025-09-13 10:48:19.812677 - step 118 of 120 - 10000 active elements (0 deactivated)
15:18:42 INFO    opendrift:2085: 2025-09-13 10:58:19.812677 - step 119 of 120 - 10000 active elements (0 deactivated)
15:18:43 INFO    opendrift:2085: 2025-09-13 11:08:19.812677 - step 120 of 120 - 10000 active elements (0 deactivated)
<xarray.Dataset> Size: 36MB
Dimensions:                                                                              (
                                                                                          trajectory: 10000,
                                                                                          time: 21)
Coordinates:
  * trajectory                                                                           (trajectory) int64 80kB ...
  * time                                                                                 (time) datetime64[ns] 168B ...
Data variables: (12/43)
    status                                                                               (trajectory, time) float32 840kB ...
    moving                                                                               (trajectory, time) float32 840kB ...
    age_seconds                                                                          (trajectory, time) float32 840kB ...
    origin_marker                                                                        (trajectory, time) float32 840kB ...
    lon                                                                                  (trajectory, time) float32 840kB ...
    lat                                                                                  (trajectory, time) float32 840kB ...
    ...                                                                                   ...
    sea_water_temperature                                                                (trajectory, time) float32 840kB ...
    sea_water_salinity                                                                   (trajectory, time) float32 840kB ...
    sea_floor_depth_below_sea_level                                                      (trajectory, time) float32 840kB ...
    ocean_vertical_diffusivity                                                           (trajectory, time) float32 840kB ...
    land_binary_mask                                                                     (trajectory, time) float32 840kB ...
    ocean_mixed_layer_thickness                                                          (trajectory, time) float32 840kB ...
Attributes: (12/165)
    Conventions:                                                             ...
    standard_name_vocabulary:                                                ...
    featureType:                                                             ...
    title:                                                                   ...
    summary:                                                                 ...
    keywords:                                                                ...
    ...                                                                                                               ...
    geospatial_lon_units:                                                    ...
    geospatial_lon_resolution:                                               ...
    runtime:                                                                 ...
    geospatial_vertical_min:                                                 ...
    geospatial_vertical_max:                                                 ...
    geospatial_vertical_positive:                                            ...


Comparison plots

o1.plot_oil_budget()
o2.plot_oil_budget()
  • GENERIC DIESEL (855.8 kg/m3) - 2025-09-12 15:17 to 2025-09-13 11:17
  • GENERIC DIESEL (855.8 kg/m3) - 2025-09-12 15:18 to 2025-09-13 11:18

Entrainment

b1 = o1.get_oil_budget()
b2 = o2.get_oil_budget()
plt.plot(b1['mass_surface'], '-r', linewidth=2,
            label='Surface, updated thickness')
plt.plot(b1['mass_submerged'], '--r', linewidth=2,
            label='Submerged, updated thickness')
plt.plot(b1['mass_evaporated'], '-.r', linewidth=2,
            label='Evaporated, updated thickness')
plt.plot(b2['mass_surface'], '-b', linewidth=2,
            label='Surface, constant thickness')
plt.plot(b2['mass_submerged'], '--b', linewidth=2,
            label='Submerged, constant thickness')
plt.plot(b2['mass_evaporated'], '-.b', linewidth=2,
            label='Evaporated, constant thickness')
plt.legend()
plt.xlabel('Time step')
plt.show()
example oil thickness

We see that with the updated film thickness, the droplets are getting gradually smaller

plt.plot(1e6*o1.result.diameter.median(dim='trajectory'))
plt.plot(1e6*o2.result.diameter.median(dim='trajectory'))
plt.legend(['With updated film thickness', 'With constant film thickness'])
plt.xlabel('Time step')
plt.ylabel('Median droplet diameter  [micrometer]')
plt.show()
example oil thickness

We see that oil film thickness has virtually no impact on horizontal drift

o1.animation(compare=o2, fast=True,
             legend=['Updated film thickness',
                     'Constant/default film thickness'])
15:18:45 WARNING opendrift:2444: Plotting fast. This will make your plots less accurate.
15:18:46 INFO    opendrift:4612: Saving animation to /root/project/docs/source/gallery/animations/example_oil_thickness_1.gif...
15:19:08 INFO    opendrift:3053: Time to make animation: 0:00:23.078417
../_images/example_oil_thickness_1.gif

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

Gallery generated by Sphinx-Gallery