Note

Go to the end to download the full example code

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')

14:02:46 INFO    opendrift.models.basemodel:539: OpenDriftSimulation initialised (version 1.10.7 / v1.10.6-119-g1da5bec)

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.utcnow())
o1.run(time_step=timestep, time_step_output=timestep_output,
       duration=duration)

14:02:46 INFO    opendrift.models.openoil.openoil:1462: setting oil_type to: GENERIC DIESEL
14:02:46 INFO    opendrift.models.openoil.adios.dirjs:71: Querying ADIOS database for oil: GENERIC DIESEL
14:02:46 INFO    opendrift.models.openoil.openoil:1652: Using density 841.2 and viscosity 6.3404392151148136e-06 of oiltype GENERIC DIESEL
14:02:46 INFO    opendrift.models.basemodel:2591: Fallback values will be used for the following variables which have no readers:
14:02:46 INFO    opendrift.models.basemodel:2594:       x_sea_water_velocity: 0.100000
14:02:46 INFO    opendrift.models.basemodel:2594:       y_sea_water_velocity: 0.000000
14:02:46 INFO    opendrift.models.basemodel:2594:       x_wind: 0.000000
14:02:46 INFO    opendrift.models.basemodel:2594:       y_wind: 7.000000
14:02:46 INFO    opendrift.models.basemodel:2594:       upward_sea_water_velocity: 0.000000
14:02:46 INFO    opendrift.models.basemodel:2594:       sea_surface_wave_significant_height: 0.000000
14:02:46 INFO    opendrift.models.basemodel:2594:       sea_surface_wave_stokes_drift_x_velocity: 0.000000
14:02:46 INFO    opendrift.models.basemodel:2594:       sea_surface_wave_stokes_drift_y_velocity: 0.300000
14:02:46 INFO    opendrift.models.basemodel:2594:       sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
14:02:46 INFO    opendrift.models.basemodel:2594:       sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
14:02:46 INFO    opendrift.models.basemodel:2594:       sea_ice_area_fraction: 0.000000
14:02:46 INFO    opendrift.models.basemodel:2594:       sea_ice_x_velocity: 0.000000
14:02:46 INFO    opendrift.models.basemodel:2594:       sea_ice_y_velocity: 0.000000
14:02:46 INFO    opendrift.models.basemodel:2594:       sea_water_temperature: 10.000000
14:02:46 INFO    opendrift.models.basemodel:2594:       sea_water_salinity: 34.000000
14:02:46 INFO    opendrift.models.basemodel:2594:       sea_floor_depth_below_sea_level: 10000.000000
14:02:46 INFO    opendrift.models.basemodel:2594:       ocean_vertical_diffusivity: 0.020000
14:02:46 INFO    opendrift.models.basemodel:2594:       ocean_mixed_layer_thickness: 50.000000
14:02:46 INFO    opendrift.models.basemodel:1689: No land reader added, making a temporary landmask reader
14:02:52 INFO    opendrift.models.basemodel:539: OpenDriftSimulation initialised (version 1.10.7 / v1.10.6-119-g1da5bec)
14:02:52 INFO    opendrift.models.basemodel:1713: All points are in ocean
14:02:52 INFO    opendrift.models.openoil.openoil:646: Oil-water surface tension is 0.029464 Nm
14:02:52 INFO    opendrift.models.basemodel:2882: 2023-09-01 14:02:46.313665 - step 1 of 120 - 10000 active elements (0 deactivated)
14:02:53 INFO    opendrift.models.basemodel:2882: 2023-09-01 14:12:46.313665 - step 2 of 120 - 10000 active elements (0 deactivated)
14:02:53 INFO    opendrift.models.basemodel:2882: 2023-09-01 14:22:46.313665 - step 3 of 120 - 10000 active elements (0 deactivated)
14:02:53 INFO    opendrift.models.basemodel:2882: 2023-09-01 14:32:46.313665 - step 4 of 120 - 10000 active elements (0 deactivated)
14:02:54 INFO    opendrift.models.basemodel:2882: 2023-09-01 14:42:46.313665 - step 5 of 120 - 10000 active elements (0 deactivated)
14:02:54 INFO    opendrift.models.basemodel:2882: 2023-09-01 14:52:46.313665 - step 6 of 120 - 10000 active elements (0 deactivated)
14:02:54 INFO    opendrift.models.basemodel:2882: 2023-09-01 15:02:46.313665 - step 7 of 120 - 10000 active elements (0 deactivated)
14:02:54 INFO    opendrift.models.basemodel:2882: 2023-09-01 15:12:46.313665 - step 8 of 120 - 10000 active elements (0 deactivated)
14:02:54 INFO    opendrift.models.basemodel:2882: 2023-09-01 15:22:46.313665 - step 9 of 120 - 10000 active elements (0 deactivated)
14:02:54 INFO    opendrift.models.basemodel:2882: 2023-09-01 15:32:46.313665 - step 10 of 120 - 10000 active elements (0 deactivated)
14:02:55 INFO    opendrift.models.basemodel:2882: 2023-09-01 15:42:46.313665 - step 11 of 120 - 10000 active elements (0 deactivated)
14:02:55 INFO    opendrift.models.basemodel:2882: 2023-09-01 15:52:46.313665 - step 12 of 120 - 10000 active elements (0 deactivated)
14:02:55 INFO    opendrift.models.basemodel:2882: 2023-09-01 16:02:46.313665 - step 13 of 120 - 10000 active elements (0 deactivated)
14:02:55 INFO    opendrift.models.basemodel:2882: 2023-09-01 16:12:46.313665 - step 14 of 120 - 10000 active elements (0 deactivated)
14:02:55 INFO    opendrift.models.basemodel:2882: 2023-09-01 16:22:46.313665 - step 15 of 120 - 10000 active elements (0 deactivated)
14:02:55 INFO    opendrift.models.basemodel:2882: 2023-09-01 16:32:46.313665 - step 16 of 120 - 10000 active elements (0 deactivated)
14:02:55 INFO    opendrift.models.basemodel:2882: 2023-09-01 16:42:46.313665 - step 17 of 120 - 10000 active elements (0 deactivated)
14:02:56 INFO    opendrift.models.basemodel:2882: 2023-09-01 16:52:46.313665 - step 18 of 120 - 10000 active elements (0 deactivated)
14:02:56 INFO    opendrift.models.basemodel:2882: 2023-09-01 17:02:46.313665 - step 19 of 120 - 10000 active elements (0 deactivated)
14:02:56 INFO    opendrift.models.basemodel:2882: 2023-09-01 17:12:46.313665 - step 20 of 120 - 10000 active elements (0 deactivated)
14:02:56 INFO    opendrift.models.basemodel:2882: 2023-09-01 17:22:46.313665 - step 21 of 120 - 10000 active elements (0 deactivated)
14:02:56 INFO    opendrift.models.basemodel:2882: 2023-09-01 17:32:46.313665 - step 22 of 120 - 10000 active elements (0 deactivated)
14:02:56 INFO    opendrift.models.basemodel:2882: 2023-09-01 17:42:46.313665 - step 23 of 120 - 10000 active elements (0 deactivated)
14:02:56 INFO    opendrift.models.basemodel:2882: 2023-09-01 17:52:46.313665 - step 24 of 120 - 10000 active elements (0 deactivated)
14:02:57 INFO    opendrift.models.basemodel:2882: 2023-09-01 18:02:46.313665 - step 25 of 120 - 10000 active elements (0 deactivated)
14:02:57 INFO    opendrift.models.basemodel:2882: 2023-09-01 18:12:46.313665 - step 26 of 120 - 10000 active elements (0 deactivated)
14:02:57 INFO    opendrift.models.basemodel:2882: 2023-09-01 18:22:46.313665 - step 27 of 120 - 10000 active elements (0 deactivated)
14:02:57 INFO    opendrift.models.basemodel:2882: 2023-09-01 18:32:46.313665 - step 28 of 120 - 10000 active elements (0 deactivated)
14:02:57 INFO    opendrift.models.basemodel:2882: 2023-09-01 18:42:46.313665 - step 29 of 120 - 10000 active elements (0 deactivated)
14:02:57 INFO    opendrift.models.basemodel:2882: 2023-09-01 18:52:46.313665 - step 30 of 120 - 10000 active elements (0 deactivated)
14:02:57 INFO    opendrift.models.basemodel:2882: 2023-09-01 19:02:46.313665 - step 31 of 120 - 10000 active elements (0 deactivated)
14:02:58 INFO    opendrift.models.basemodel:2882: 2023-09-01 19:12:46.313665 - step 32 of 120 - 10000 active elements (0 deactivated)
14:02:58 INFO    opendrift.models.basemodel:2882: 2023-09-01 19:22:46.313665 - step 33 of 120 - 10000 active elements (0 deactivated)
14:02:58 INFO    opendrift.models.basemodel:2882: 2023-09-01 19:32:46.313665 - step 34 of 120 - 10000 active elements (0 deactivated)
14:02:58 INFO    opendrift.models.basemodel:2882: 2023-09-01 19:42:46.313665 - step 35 of 120 - 10000 active elements (0 deactivated)
14:02:58 INFO    opendrift.models.basemodel:2882: 2023-09-01 19:52:46.313665 - step 36 of 120 - 10000 active elements (0 deactivated)
14:02:58 INFO    opendrift.models.basemodel:2882: 2023-09-01 20:02:46.313665 - step 37 of 120 - 10000 active elements (0 deactivated)
14:02:58 INFO    opendrift.models.basemodel:2882: 2023-09-01 20:12:46.313665 - step 38 of 120 - 10000 active elements (0 deactivated)
14:02:58 INFO    opendrift.models.basemodel:2882: 2023-09-01 20:22:46.313665 - step 39 of 120 - 10000 active elements (0 deactivated)
14:02:59 INFO    opendrift.models.basemodel:2882: 2023-09-01 20:32:46.313665 - step 40 of 120 - 10000 active elements (0 deactivated)
14:02:59 INFO    opendrift.models.basemodel:2882: 2023-09-01 20:42:46.313665 - step 41 of 120 - 10000 active elements (0 deactivated)
14:02:59 INFO    opendrift.models.basemodel:2882: 2023-09-01 20:52:46.313665 - step 42 of 120 - 10000 active elements (0 deactivated)
14:02:59 INFO    opendrift.models.basemodel:2882: 2023-09-01 21:02:46.313665 - step 43 of 120 - 10000 active elements (0 deactivated)
14:02:59 INFO    opendrift.models.basemodel:2882: 2023-09-01 21:12:46.313665 - step 44 of 120 - 10000 active elements (0 deactivated)
14:02:59 INFO    opendrift.models.basemodel:2882: 2023-09-01 21:22:46.313665 - step 45 of 120 - 10000 active elements (0 deactivated)
14:02:59 INFO    opendrift.models.basemodel:2882: 2023-09-01 21:32:46.313665 - step 46 of 120 - 10000 active elements (0 deactivated)
14:03:00 INFO    opendrift.models.basemodel:2882: 2023-09-01 21:42:46.313665 - step 47 of 120 - 10000 active elements (0 deactivated)
14:03:00 INFO    opendrift.models.basemodel:2882: 2023-09-01 21:52:46.313665 - step 48 of 120 - 10000 active elements (0 deactivated)
14:03:00 INFO    opendrift.models.basemodel:2882: 2023-09-01 22:02:46.313665 - step 49 of 120 - 10000 active elements (0 deactivated)
14:03:00 INFO    opendrift.models.basemodel:2882: 2023-09-01 22:12:46.313665 - step 50 of 120 - 10000 active elements (0 deactivated)
14:03:00 INFO    opendrift.models.basemodel:2882: 2023-09-01 22:22:46.313665 - step 51 of 120 - 10000 active elements (0 deactivated)
14:03:00 INFO    opendrift.models.basemodel:2882: 2023-09-01 22:32:46.313665 - step 52 of 120 - 10000 active elements (0 deactivated)
14:03:00 INFO    opendrift.models.basemodel:2882: 2023-09-01 22:42:46.313665 - step 53 of 120 - 10000 active elements (0 deactivated)
14:03:00 INFO    opendrift.models.basemodel:2882: 2023-09-01 22:52:46.313665 - step 54 of 120 - 10000 active elements (0 deactivated)
14:03:01 INFO    opendrift.models.basemodel:2882: 2023-09-01 23:02:46.313665 - step 55 of 120 - 10000 active elements (0 deactivated)
14:03:01 INFO    opendrift.models.basemodel:2882: 2023-09-01 23:12:46.313665 - step 56 of 120 - 10000 active elements (0 deactivated)
14:03:01 INFO    opendrift.models.basemodel:2882: 2023-09-01 23:22:46.313665 - step 57 of 120 - 10000 active elements (0 deactivated)
14:03:01 INFO    opendrift.models.basemodel:2882: 2023-09-01 23:32:46.313665 - step 58 of 120 - 10000 active elements (0 deactivated)
14:03:01 INFO    opendrift.models.basemodel:2882: 2023-09-01 23:42:46.313665 - step 59 of 120 - 10000 active elements (0 deactivated)
14:03:01 INFO    opendrift.models.basemodel:2882: 2023-09-01 23:52:46.313665 - step 60 of 120 - 10000 active elements (0 deactivated)
14:03:01 INFO    opendrift.models.basemodel:2882: 2023-09-02 00:02:46.313665 - step 61 of 120 - 10000 active elements (0 deactivated)
14:03:01 INFO    opendrift.models.basemodel:2882: 2023-09-02 00:12:46.313665 - step 62 of 120 - 10000 active elements (0 deactivated)
14:03:02 INFO    opendrift.models.basemodel:2882: 2023-09-02 00:22:46.313665 - step 63 of 120 - 10000 active elements (0 deactivated)
14:03:02 INFO    opendrift.models.basemodel:2882: 2023-09-02 00:32:46.313665 - step 64 of 120 - 10000 active elements (0 deactivated)
14:03:02 INFO    opendrift.models.basemodel:2882: 2023-09-02 00:42:46.313665 - step 65 of 120 - 10000 active elements (0 deactivated)
14:03:02 INFO    opendrift.models.basemodel:2882: 2023-09-02 00:52:46.313665 - step 66 of 120 - 10000 active elements (0 deactivated)
14:03:02 INFO    opendrift.models.basemodel:2882: 2023-09-02 01:02:46.313665 - step 67 of 120 - 10000 active elements (0 deactivated)
14:03:02 INFO    opendrift.models.basemodel:2882: 2023-09-02 01:12:46.313665 - step 68 of 120 - 10000 active elements (0 deactivated)
14:03:02 INFO    opendrift.models.basemodel:2882: 2023-09-02 01:22:46.313665 - step 69 of 120 - 10000 active elements (0 deactivated)
14:03:02 INFO    opendrift.models.basemodel:2882: 2023-09-02 01:32:46.313665 - step 70 of 120 - 10000 active elements (0 deactivated)
14:03:03 INFO    opendrift.models.basemodel:2882: 2023-09-02 01:42:46.313665 - step 71 of 120 - 10000 active elements (0 deactivated)
14:03:03 INFO    opendrift.models.basemodel:2882: 2023-09-02 01:52:46.313665 - step 72 of 120 - 10000 active elements (0 deactivated)
14:03:03 INFO    opendrift.models.basemodel:2882: 2023-09-02 02:02:46.313665 - step 73 of 120 - 10000 active elements (0 deactivated)
14:03:03 INFO    opendrift.models.basemodel:2882: 2023-09-02 02:12:46.313665 - step 74 of 120 - 10000 active elements (0 deactivated)
14:03:03 INFO    opendrift.models.basemodel:2882: 2023-09-02 02:22:46.313665 - step 75 of 120 - 10000 active elements (0 deactivated)
14:03:03 INFO    opendrift.models.basemodel:2882: 2023-09-02 02:32:46.313665 - step 76 of 120 - 10000 active elements (0 deactivated)
14:03:03 INFO    opendrift.models.basemodel:2882: 2023-09-02 02:42:46.313665 - step 77 of 120 - 10000 active elements (0 deactivated)
14:03:03 INFO    opendrift.models.basemodel:2882: 2023-09-02 02:52:46.313665 - step 78 of 120 - 10000 active elements (0 deactivated)
14:03:04 INFO    opendrift.models.basemodel:2882: 2023-09-02 03:02:46.313665 - step 79 of 120 - 10000 active elements (0 deactivated)
14:03:04 INFO    opendrift.models.basemodel:2882: 2023-09-02 03:12:46.313665 - step 80 of 120 - 10000 active elements (0 deactivated)
14:03:04 INFO    opendrift.models.basemodel:2882: 2023-09-02 03:22:46.313665 - step 81 of 120 - 10000 active elements (0 deactivated)
14:03:04 INFO    opendrift.models.basemodel:2882: 2023-09-02 03:32:46.313665 - step 82 of 120 - 10000 active elements (0 deactivated)
14:03:04 INFO    opendrift.models.basemodel:2882: 2023-09-02 03:42:46.313665 - step 83 of 120 - 10000 active elements (0 deactivated)
14:03:04 INFO    opendrift.models.basemodel:2882: 2023-09-02 03:52:46.313665 - step 84 of 120 - 10000 active elements (0 deactivated)
14:03:04 INFO    opendrift.models.basemodel:2882: 2023-09-02 04:02:46.313665 - step 85 of 120 - 10000 active elements (0 deactivated)
14:03:05 INFO    opendrift.models.basemodel:2882: 2023-09-02 04:12:46.313665 - step 86 of 120 - 10000 active elements (0 deactivated)
14:03:05 INFO    opendrift.models.basemodel:2882: 2023-09-02 04:22:46.313665 - step 87 of 120 - 10000 active elements (0 deactivated)
14:03:05 INFO    opendrift.models.basemodel:2882: 2023-09-02 04:32:46.313665 - step 88 of 120 - 10000 active elements (0 deactivated)
14:03:05 INFO    opendrift.models.basemodel:2882: 2023-09-02 04:42:46.313665 - step 89 of 120 - 10000 active elements (0 deactivated)
14:03:05 INFO    opendrift.models.basemodel:2882: 2023-09-02 04:52:46.313665 - step 90 of 120 - 10000 active elements (0 deactivated)
14:03:05 INFO    opendrift.models.basemodel:2882: 2023-09-02 05:02:46.313665 - step 91 of 120 - 10000 active elements (0 deactivated)
14:03:05 INFO    opendrift.models.basemodel:2882: 2023-09-02 05:12:46.313665 - step 92 of 120 - 10000 active elements (0 deactivated)
14:03:05 INFO    opendrift.models.basemodel:2882: 2023-09-02 05:22:46.313665 - step 93 of 120 - 10000 active elements (0 deactivated)
14:03:06 INFO    opendrift.models.basemodel:2882: 2023-09-02 05:32:46.313665 - step 94 of 120 - 10000 active elements (0 deactivated)
14:03:06 INFO    opendrift.models.basemodel:2882: 2023-09-02 05:42:46.313665 - step 95 of 120 - 10000 active elements (0 deactivated)
14:03:06 INFO    opendrift.models.basemodel:2882: 2023-09-02 05:52:46.313665 - step 96 of 120 - 10000 active elements (0 deactivated)
14:03:06 INFO    opendrift.models.basemodel:2882: 2023-09-02 06:02:46.313665 - step 97 of 120 - 10000 active elements (0 deactivated)
14:03:06 INFO    opendrift.models.basemodel:2882: 2023-09-02 06:12:46.313665 - step 98 of 120 - 10000 active elements (0 deactivated)
14:03:06 INFO    opendrift.models.basemodel:2882: 2023-09-02 06:22:46.313665 - step 99 of 120 - 10000 active elements (0 deactivated)
14:03:06 INFO    opendrift.models.basemodel:2882: 2023-09-02 06:32:46.313665 - step 100 of 120 - 10000 active elements (0 deactivated)
14:03:06 INFO    opendrift.models.basemodel:2882: 2023-09-02 06:42:46.313665 - step 101 of 120 - 10000 active elements (0 deactivated)
14:03:07 INFO    opendrift.models.basemodel:2882: 2023-09-02 06:52:46.313665 - step 102 of 120 - 10000 active elements (0 deactivated)
14:03:07 INFO    opendrift.models.basemodel:2882: 2023-09-02 07:02:46.313665 - step 103 of 120 - 10000 active elements (0 deactivated)
14:03:07 INFO    opendrift.models.basemodel:2882: 2023-09-02 07:12:46.313665 - step 104 of 120 - 10000 active elements (0 deactivated)
14:03:07 INFO    opendrift.models.basemodel:2882: 2023-09-02 07:22:46.313665 - step 105 of 120 - 10000 active elements (0 deactivated)
14:03:07 INFO    opendrift.models.basemodel:2882: 2023-09-02 07:32:46.313665 - step 106 of 120 - 10000 active elements (0 deactivated)
14:03:07 INFO    opendrift.models.basemodel:2882: 2023-09-02 07:42:46.313665 - step 107 of 120 - 10000 active elements (0 deactivated)
14:03:07 INFO    opendrift.models.basemodel:2882: 2023-09-02 07:52:46.313665 - step 108 of 120 - 10000 active elements (0 deactivated)
14:03:08 INFO    opendrift.models.basemodel:2882: 2023-09-02 08:02:46.313665 - step 109 of 120 - 10000 active elements (0 deactivated)
14:03:08 INFO    opendrift.models.basemodel:2882: 2023-09-02 08:12:46.313665 - step 110 of 120 - 10000 active elements (0 deactivated)
14:03:08 INFO    opendrift.models.basemodel:2882: 2023-09-02 08:22:46.313665 - step 111 of 120 - 10000 active elements (0 deactivated)
14:03:08 INFO    opendrift.models.basemodel:2882: 2023-09-02 08:32:46.313665 - step 112 of 120 - 10000 active elements (0 deactivated)
14:03:08 INFO    opendrift.models.basemodel:2882: 2023-09-02 08:42:46.313665 - step 113 of 120 - 10000 active elements (0 deactivated)
14:03:08 INFO    opendrift.models.basemodel:2882: 2023-09-02 08:52:46.313665 - step 114 of 120 - 10000 active elements (0 deactivated)
14:03:08 INFO    opendrift.models.basemodel:2882: 2023-09-02 09:02:46.313665 - step 115 of 120 - 10000 active elements (0 deactivated)
14:03:09 INFO    opendrift.models.basemodel:2882: 2023-09-02 09:12:46.313665 - step 116 of 120 - 10000 active elements (0 deactivated)
14:03:09 INFO    opendrift.models.basemodel:2882: 2023-09-02 09:22:46.313665 - step 117 of 120 - 10000 active elements (0 deactivated)
14:03:09 INFO    opendrift.models.basemodel:2882: 2023-09-02 09:32:46.313665 - step 118 of 120 - 10000 active elements (0 deactivated)
14:03:09 INFO    opendrift.models.basemodel:2882: 2023-09-02 09:42:46.313665 - step 119 of 120 - 10000 active elements (0 deactivated)
14:03:09 INFO    opendrift.models.basemodel:2882: 2023-09-02 09:52:46.313665 - step 120 of 120 - 10000 active elements (0 deactivated)

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)

14:03:09 WARNING opendrift.models.basemodel:3226: Plotting fast. This will make your plots less accurate.
14:03:11 INFO    opendrift.models.basemodel:5350: Saving animation to /root/project/docs/source/gallery/animations/example_oil_thickness_0.gif...
14:03:16 INFO    opendrift.models.basemodel:3848: Time to make animation: 0:00:06.751197
../_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.utcnow())
o2.run(time_step=timestep, time_step_output=timestep_output,
       duration=duration)

14:03:16 INFO    opendrift.models.basemodel:539: OpenDriftSimulation initialised (version 1.10.7 / v1.10.6-119-g1da5bec)
14:03:16 INFO    opendrift.models.openoil.openoil:1462: setting oil_type to: GENERIC DIESEL
14:03:16 INFO    opendrift.models.openoil.adios.dirjs:71: Querying ADIOS database for oil: GENERIC DIESEL
14:03:16 INFO    opendrift.models.openoil.openoil:1652: Using density 841.2 and viscosity 6.3404392151148136e-06 of oiltype GENERIC DIESEL
14:03:16 INFO    opendrift.models.basemodel:2591: Fallback values will be used for the following variables which have no readers:
14:03:16 INFO    opendrift.models.basemodel:2594:       x_sea_water_velocity: 0.100000
14:03:16 INFO    opendrift.models.basemodel:2594:       y_sea_water_velocity: 0.000000
14:03:16 INFO    opendrift.models.basemodel:2594:       x_wind: 0.000000
14:03:16 INFO    opendrift.models.basemodel:2594:       y_wind: 7.000000
14:03:16 INFO    opendrift.models.basemodel:2594:       upward_sea_water_velocity: 0.000000
14:03:16 INFO    opendrift.models.basemodel:2594:       sea_surface_wave_significant_height: 0.000000
14:03:16 INFO    opendrift.models.basemodel:2594:       sea_surface_wave_stokes_drift_x_velocity: 0.000000
14:03:16 INFO    opendrift.models.basemodel:2594:       sea_surface_wave_stokes_drift_y_velocity: 0.300000
14:03:16 INFO    opendrift.models.basemodel:2594:       sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
14:03:16 INFO    opendrift.models.basemodel:2594:       sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
14:03:16 INFO    opendrift.models.basemodel:2594:       sea_ice_area_fraction: 0.000000
14:03:16 INFO    opendrift.models.basemodel:2594:       sea_ice_x_velocity: 0.000000
14:03:16 INFO    opendrift.models.basemodel:2594:       sea_ice_y_velocity: 0.000000
14:03:16 INFO    opendrift.models.basemodel:2594:       sea_water_temperature: 10.000000
14:03:16 INFO    opendrift.models.basemodel:2594:       sea_water_salinity: 34.000000
14:03:16 INFO    opendrift.models.basemodel:2594:       sea_floor_depth_below_sea_level: 10000.000000
14:03:16 INFO    opendrift.models.basemodel:2594:       ocean_vertical_diffusivity: 0.020000
14:03:16 INFO    opendrift.models.basemodel:2594:       ocean_mixed_layer_thickness: 50.000000
14:03:17 INFO    opendrift.models.basemodel:1689: No land reader added, making a temporary landmask reader
14:03:17 INFO    opendrift.models.basemodel:539: OpenDriftSimulation initialised (version 1.10.7 / v1.10.6-119-g1da5bec)
14:03:17 INFO    opendrift.models.basemodel:1713: All points are in ocean
14:03:17 INFO    opendrift.models.openoil.openoil:646: Oil-water surface tension is 0.029464 Nm
14:03:17 INFO    opendrift.models.basemodel:2882: 2023-09-01 14:03:16.850642 - step 1 of 120 - 10000 active elements (0 deactivated)
14:03:17 INFO    opendrift.models.basemodel:2882: 2023-09-01 14:13:16.850642 - step 2 of 120 - 10000 active elements (0 deactivated)
14:03:17 INFO    opendrift.models.basemodel:2882: 2023-09-01 14:23:16.850642 - step 3 of 120 - 10000 active elements (0 deactivated)
14:03:18 INFO    opendrift.models.basemodel:2882: 2023-09-01 14:33:16.850642 - step 4 of 120 - 10000 active elements (0 deactivated)
14:03:18 INFO    opendrift.models.basemodel:2882: 2023-09-01 14:43:16.850642 - step 5 of 120 - 10000 active elements (0 deactivated)
14:03:18 INFO    opendrift.models.basemodel:2882: 2023-09-01 14:53:16.850642 - step 6 of 120 - 10000 active elements (0 deactivated)
14:03:18 INFO    opendrift.models.basemodel:2882: 2023-09-01 15:03:16.850642 - step 7 of 120 - 10000 active elements (0 deactivated)
14:03:19 INFO    opendrift.models.basemodel:2882: 2023-09-01 15:13:16.850642 - step 8 of 120 - 10000 active elements (0 deactivated)
14:03:19 INFO    opendrift.models.basemodel:2882: 2023-09-01 15:23:16.850642 - step 9 of 120 - 10000 active elements (0 deactivated)
14:03:19 INFO    opendrift.models.basemodel:2882: 2023-09-01 15:33:16.850642 - step 10 of 120 - 10000 active elements (0 deactivated)
14:03:19 INFO    opendrift.models.basemodel:2882: 2023-09-01 15:43:16.850642 - step 11 of 120 - 10000 active elements (0 deactivated)
14:03:19 INFO    opendrift.models.basemodel:2882: 2023-09-01 15:53:16.850642 - step 12 of 120 - 10000 active elements (0 deactivated)
14:03:20 INFO    opendrift.models.basemodel:2882: 2023-09-01 16:03:16.850642 - step 13 of 120 - 10000 active elements (0 deactivated)
14:03:20 INFO    opendrift.models.basemodel:2882: 2023-09-01 16:13:16.850642 - step 14 of 120 - 10000 active elements (0 deactivated)
14:03:20 INFO    opendrift.models.basemodel:2882: 2023-09-01 16:23:16.850642 - step 15 of 120 - 10000 active elements (0 deactivated)
14:03:20 INFO    opendrift.models.basemodel:2882: 2023-09-01 16:33:16.850642 - step 16 of 120 - 10000 active elements (0 deactivated)
14:03:20 INFO    opendrift.models.basemodel:2882: 2023-09-01 16:43:16.850642 - step 17 of 120 - 10000 active elements (0 deactivated)
14:03:21 INFO    opendrift.models.basemodel:2882: 2023-09-01 16:53:16.850642 - step 18 of 120 - 10000 active elements (0 deactivated)
14:03:21 INFO    opendrift.models.basemodel:2882: 2023-09-01 17:03:16.850642 - step 19 of 120 - 10000 active elements (0 deactivated)
14:03:21 INFO    opendrift.models.basemodel:2882: 2023-09-01 17:13:16.850642 - step 20 of 120 - 10000 active elements (0 deactivated)
14:03:21 INFO    opendrift.models.basemodel:2882: 2023-09-01 17:23:16.850642 - step 21 of 120 - 10000 active elements (0 deactivated)
14:03:21 INFO    opendrift.models.basemodel:2882: 2023-09-01 17:33:16.850642 - step 22 of 120 - 10000 active elements (0 deactivated)
14:03:21 INFO    opendrift.models.basemodel:2882: 2023-09-01 17:43:16.850642 - step 23 of 120 - 10000 active elements (0 deactivated)
14:03:22 INFO    opendrift.models.basemodel:2882: 2023-09-01 17:53:16.850642 - step 24 of 120 - 10000 active elements (0 deactivated)
14:03:22 INFO    opendrift.models.basemodel:2882: 2023-09-01 18:03:16.850642 - step 25 of 120 - 10000 active elements (0 deactivated)
14:03:22 INFO    opendrift.models.basemodel:2882: 2023-09-01 18:13:16.850642 - step 26 of 120 - 10000 active elements (0 deactivated)
14:03:22 INFO    opendrift.models.basemodel:2882: 2023-09-01 18:23:16.850642 - step 27 of 120 - 10000 active elements (0 deactivated)
14:03:22 INFO    opendrift.models.basemodel:2882: 2023-09-01 18:33:16.850642 - step 28 of 120 - 10000 active elements (0 deactivated)
14:03:23 INFO    opendrift.models.basemodel:2882: 2023-09-01 18:43:16.850642 - step 29 of 120 - 10000 active elements (0 deactivated)
14:03:23 INFO    opendrift.models.basemodel:2882: 2023-09-01 18:53:16.850642 - step 30 of 120 - 10000 active elements (0 deactivated)
14:03:23 INFO    opendrift.models.basemodel:2882: 2023-09-01 19:03:16.850642 - step 31 of 120 - 10000 active elements (0 deactivated)
14:03:23 INFO    opendrift.models.basemodel:2882: 2023-09-01 19:13:16.850642 - step 32 of 120 - 10000 active elements (0 deactivated)
14:03:23 INFO    opendrift.models.basemodel:2882: 2023-09-01 19:23:16.850642 - step 33 of 120 - 10000 active elements (0 deactivated)
14:03:23 INFO    opendrift.models.basemodel:2882: 2023-09-01 19:33:16.850642 - step 34 of 120 - 10000 active elements (0 deactivated)
14:03:24 INFO    opendrift.models.basemodel:2882: 2023-09-01 19:43:16.850642 - step 35 of 120 - 10000 active elements (0 deactivated)
14:03:24 INFO    opendrift.models.basemodel:2882: 2023-09-01 19:53:16.850642 - step 36 of 120 - 10000 active elements (0 deactivated)
14:03:24 INFO    opendrift.models.basemodel:2882: 2023-09-01 20:03:16.850642 - step 37 of 120 - 10000 active elements (0 deactivated)
14:03:24 INFO    opendrift.models.basemodel:2882: 2023-09-01 20:13:16.850642 - step 38 of 120 - 10000 active elements (0 deactivated)
14:03:24 INFO    opendrift.models.basemodel:2882: 2023-09-01 20:23:16.850642 - step 39 of 120 - 10000 active elements (0 deactivated)
14:03:25 INFO    opendrift.models.basemodel:2882: 2023-09-01 20:33:16.850642 - step 40 of 120 - 10000 active elements (0 deactivated)
14:03:25 INFO    opendrift.models.basemodel:2882: 2023-09-01 20:43:16.850642 - step 41 of 120 - 10000 active elements (0 deactivated)
14:03:25 INFO    opendrift.models.basemodel:2882: 2023-09-01 20:53:16.850642 - step 42 of 120 - 10000 active elements (0 deactivated)
14:03:25 INFO    opendrift.models.basemodel:2882: 2023-09-01 21:03:16.850642 - step 43 of 120 - 10000 active elements (0 deactivated)
14:03:25 INFO    opendrift.models.basemodel:2882: 2023-09-01 21:13:16.850642 - step 44 of 120 - 10000 active elements (0 deactivated)
14:03:26 INFO    opendrift.models.basemodel:2882: 2023-09-01 21:23:16.850642 - step 45 of 120 - 10000 active elements (0 deactivated)
14:03:26 INFO    opendrift.models.basemodel:2882: 2023-09-01 21:33:16.850642 - step 46 of 120 - 10000 active elements (0 deactivated)
14:03:26 INFO    opendrift.models.basemodel:2882: 2023-09-01 21:43:16.850642 - step 47 of 120 - 10000 active elements (0 deactivated)
14:03:26 INFO    opendrift.models.basemodel:2882: 2023-09-01 21:53:16.850642 - step 48 of 120 - 10000 active elements (0 deactivated)
14:03:26 INFO    opendrift.models.basemodel:2882: 2023-09-01 22:03:16.850642 - step 49 of 120 - 10000 active elements (0 deactivated)
14:03:27 INFO    opendrift.models.basemodel:2882: 2023-09-01 22:13:16.850642 - step 50 of 120 - 10000 active elements (0 deactivated)
14:03:27 INFO    opendrift.models.basemodel:2882: 2023-09-01 22:23:16.850642 - step 51 of 120 - 10000 active elements (0 deactivated)
14:03:27 INFO    opendrift.models.basemodel:2882: 2023-09-01 22:33:16.850642 - step 52 of 120 - 10000 active elements (0 deactivated)
14:03:27 INFO    opendrift.models.basemodel:2882: 2023-09-01 22:43:16.850642 - step 53 of 120 - 10000 active elements (0 deactivated)
14:03:27 INFO    opendrift.models.basemodel:2882: 2023-09-01 22:53:16.850642 - step 54 of 120 - 10000 active elements (0 deactivated)
14:03:27 INFO    opendrift.models.basemodel:2882: 2023-09-01 23:03:16.850642 - step 55 of 120 - 10000 active elements (0 deactivated)
14:03:28 INFO    opendrift.models.basemodel:2882: 2023-09-01 23:13:16.850642 - step 56 of 120 - 10000 active elements (0 deactivated)
14:03:28 INFO    opendrift.models.basemodel:2882: 2023-09-01 23:23:16.850642 - step 57 of 120 - 10000 active elements (0 deactivated)
14:03:28 INFO    opendrift.models.basemodel:2882: 2023-09-01 23:33:16.850642 - step 58 of 120 - 10000 active elements (0 deactivated)
14:03:28 INFO    opendrift.models.basemodel:2882: 2023-09-01 23:43:16.850642 - step 59 of 120 - 10000 active elements (0 deactivated)
14:03:28 INFO    opendrift.models.basemodel:2882: 2023-09-01 23:53:16.850642 - step 60 of 120 - 10000 active elements (0 deactivated)
14:03:29 INFO    opendrift.models.basemodel:2882: 2023-09-02 00:03:16.850642 - step 61 of 120 - 10000 active elements (0 deactivated)
14:03:29 INFO    opendrift.models.basemodel:2882: 2023-09-02 00:13:16.850642 - step 62 of 120 - 10000 active elements (0 deactivated)
14:03:29 INFO    opendrift.models.basemodel:2882: 2023-09-02 00:23:16.850642 - step 63 of 120 - 10000 active elements (0 deactivated)
14:03:29 INFO    opendrift.models.basemodel:2882: 2023-09-02 00:33:16.850642 - step 64 of 120 - 10000 active elements (0 deactivated)
14:03:29 INFO    opendrift.models.basemodel:2882: 2023-09-02 00:43:16.850642 - step 65 of 120 - 10000 active elements (0 deactivated)
14:03:29 INFO    opendrift.models.basemodel:2882: 2023-09-02 00:53:16.850642 - step 66 of 120 - 10000 active elements (0 deactivated)
14:03:30 INFO    opendrift.models.basemodel:2882: 2023-09-02 01:03:16.850642 - step 67 of 120 - 10000 active elements (0 deactivated)
14:03:30 INFO    opendrift.models.basemodel:2882: 2023-09-02 01:13:16.850642 - step 68 of 120 - 10000 active elements (0 deactivated)
14:03:30 INFO    opendrift.models.basemodel:2882: 2023-09-02 01:23:16.850642 - step 69 of 120 - 10000 active elements (0 deactivated)
14:03:30 INFO    opendrift.models.basemodel:2882: 2023-09-02 01:33:16.850642 - step 70 of 120 - 10000 active elements (0 deactivated)
14:03:30 INFO    opendrift.models.basemodel:2882: 2023-09-02 01:43:16.850642 - step 71 of 120 - 10000 active elements (0 deactivated)
14:03:31 INFO    opendrift.models.basemodel:2882: 2023-09-02 01:53:16.850642 - step 72 of 120 - 10000 active elements (0 deactivated)
14:03:31 INFO    opendrift.models.basemodel:2882: 2023-09-02 02:03:16.850642 - step 73 of 120 - 10000 active elements (0 deactivated)
14:03:31 INFO    opendrift.models.basemodel:2882: 2023-09-02 02:13:16.850642 - step 74 of 120 - 10000 active elements (0 deactivated)
14:03:31 INFO    opendrift.models.basemodel:2882: 2023-09-02 02:23:16.850642 - step 75 of 120 - 10000 active elements (0 deactivated)
14:03:31 INFO    opendrift.models.basemodel:2882: 2023-09-02 02:33:16.850642 - step 76 of 120 - 10000 active elements (0 deactivated)
14:03:31 INFO    opendrift.models.basemodel:2882: 2023-09-02 02:43:16.850642 - step 77 of 120 - 10000 active elements (0 deactivated)
14:03:32 INFO    opendrift.models.basemodel:2882: 2023-09-02 02:53:16.850642 - step 78 of 120 - 10000 active elements (0 deactivated)
14:03:32 INFO    opendrift.models.basemodel:2882: 2023-09-02 03:03:16.850642 - step 79 of 120 - 10000 active elements (0 deactivated)
14:03:32 INFO    opendrift.models.basemodel:2882: 2023-09-02 03:13:16.850642 - step 80 of 120 - 10000 active elements (0 deactivated)
14:03:32 INFO    opendrift.models.basemodel:2882: 2023-09-02 03:23:16.850642 - step 81 of 120 - 10000 active elements (0 deactivated)
14:03:32 INFO    opendrift.models.basemodel:2882: 2023-09-02 03:33:16.850642 - step 82 of 120 - 10000 active elements (0 deactivated)
14:03:32 INFO    opendrift.models.basemodel:2882: 2023-09-02 03:43:16.850642 - step 83 of 120 - 10000 active elements (0 deactivated)
14:03:33 INFO    opendrift.models.basemodel:2882: 2023-09-02 03:53:16.850642 - step 84 of 120 - 10000 active elements (0 deactivated)
14:03:33 INFO    opendrift.models.basemodel:2882: 2023-09-02 04:03:16.850642 - step 85 of 120 - 10000 active elements (0 deactivated)
14:03:33 INFO    opendrift.models.basemodel:2882: 2023-09-02 04:13:16.850642 - step 86 of 120 - 10000 active elements (0 deactivated)
14:03:33 INFO    opendrift.models.basemodel:2882: 2023-09-02 04:23:16.850642 - step 87 of 120 - 10000 active elements (0 deactivated)
14:03:33 INFO    opendrift.models.basemodel:2882: 2023-09-02 04:33:16.850642 - step 88 of 120 - 10000 active elements (0 deactivated)
14:03:33 INFO    opendrift.models.basemodel:2882: 2023-09-02 04:43:16.850642 - step 89 of 120 - 10000 active elements (0 deactivated)
14:03:34 INFO    opendrift.models.basemodel:2882: 2023-09-02 04:53:16.850642 - step 90 of 120 - 10000 active elements (0 deactivated)
14:03:34 INFO    opendrift.models.basemodel:2882: 2023-09-02 05:03:16.850642 - step 91 of 120 - 10000 active elements (0 deactivated)
14:03:34 INFO    opendrift.models.basemodel:2882: 2023-09-02 05:13:16.850642 - step 92 of 120 - 10000 active elements (0 deactivated)
14:03:34 INFO    opendrift.models.basemodel:2882: 2023-09-02 05:23:16.850642 - step 93 of 120 - 10000 active elements (0 deactivated)
14:03:34 INFO    opendrift.models.basemodel:2882: 2023-09-02 05:33:16.850642 - step 94 of 120 - 10000 active elements (0 deactivated)
14:03:34 INFO    opendrift.models.basemodel:2882: 2023-09-02 05:43:16.850642 - step 95 of 120 - 10000 active elements (0 deactivated)
14:03:35 INFO    opendrift.models.basemodel:2882: 2023-09-02 05:53:16.850642 - step 96 of 120 - 10000 active elements (0 deactivated)
14:03:35 INFO    opendrift.models.basemodel:2882: 2023-09-02 06:03:16.850642 - step 97 of 120 - 10000 active elements (0 deactivated)
14:03:35 INFO    opendrift.models.basemodel:2882: 2023-09-02 06:13:16.850642 - step 98 of 120 - 10000 active elements (0 deactivated)
14:03:35 INFO    opendrift.models.basemodel:2882: 2023-09-02 06:23:16.850642 - step 99 of 120 - 10000 active elements (0 deactivated)
14:03:35 INFO    opendrift.models.basemodel:2882: 2023-09-02 06:33:16.850642 - step 100 of 120 - 10000 active elements (0 deactivated)
14:03:35 INFO    opendrift.models.basemodel:2882: 2023-09-02 06:43:16.850642 - step 101 of 120 - 10000 active elements (0 deactivated)
14:03:36 INFO    opendrift.models.basemodel:2882: 2023-09-02 06:53:16.850642 - step 102 of 120 - 10000 active elements (0 deactivated)
14:03:36 INFO    opendrift.models.basemodel:2882: 2023-09-02 07:03:16.850642 - step 103 of 120 - 10000 active elements (0 deactivated)
14:03:36 INFO    opendrift.models.basemodel:2882: 2023-09-02 07:13:16.850642 - step 104 of 120 - 10000 active elements (0 deactivated)
14:03:36 INFO    opendrift.models.basemodel:2882: 2023-09-02 07:23:16.850642 - step 105 of 120 - 10000 active elements (0 deactivated)
14:03:36 INFO    opendrift.models.basemodel:2882: 2023-09-02 07:33:16.850642 - step 106 of 120 - 10000 active elements (0 deactivated)
14:03:36 INFO    opendrift.models.basemodel:2882: 2023-09-02 07:43:16.850642 - step 107 of 120 - 10000 active elements (0 deactivated)
14:03:37 INFO    opendrift.models.basemodel:2882: 2023-09-02 07:53:16.850642 - step 108 of 120 - 10000 active elements (0 deactivated)
14:03:37 INFO    opendrift.models.basemodel:2882: 2023-09-02 08:03:16.850642 - step 109 of 120 - 10000 active elements (0 deactivated)
14:03:37 INFO    opendrift.models.basemodel:2882: 2023-09-02 08:13:16.850642 - step 110 of 120 - 10000 active elements (0 deactivated)
14:03:37 INFO    opendrift.models.basemodel:2882: 2023-09-02 08:23:16.850642 - step 111 of 120 - 10000 active elements (0 deactivated)
14:03:37 INFO    opendrift.models.basemodel:2882: 2023-09-02 08:33:16.850642 - step 112 of 120 - 10000 active elements (0 deactivated)
14:03:37 INFO    opendrift.models.basemodel:2882: 2023-09-02 08:43:16.850642 - step 113 of 120 - 10000 active elements (0 deactivated)
14:03:37 INFO    opendrift.models.basemodel:2882: 2023-09-02 08:53:16.850642 - step 114 of 120 - 10000 active elements (0 deactivated)
14:03:38 INFO    opendrift.models.basemodel:2882: 2023-09-02 09:03:16.850642 - step 115 of 120 - 10000 active elements (0 deactivated)
14:03:38 INFO    opendrift.models.basemodel:2882: 2023-09-02 09:13:16.850642 - step 116 of 120 - 10000 active elements (0 deactivated)
14:03:38 INFO    opendrift.models.basemodel:2882: 2023-09-02 09:23:16.850642 - step 117 of 120 - 10000 active elements (0 deactivated)
14:03:38 INFO    opendrift.models.basemodel:2882: 2023-09-02 09:33:16.850642 - step 118 of 120 - 10000 active elements (0 deactivated)
14:03:38 INFO    opendrift.models.basemodel:2882: 2023-09-02 09:43:16.850642 - step 119 of 120 - 10000 active elements (0 deactivated)
14:03:38 INFO    opendrift.models.basemodel:2882: 2023-09-02 09:53:16.850642 - step 120 of 120 - 10000 active elements (0 deactivated)

Comparison plots

o1.plot_oil_budget()
o2.plot_oil_budget()
  • GENERIC DIESEL (841.2 kg/m3) - 2023-09-01 14:02 to 2023-09-02 10:02
  • GENERIC DIESEL (841.2 kg/m3) - 2023-09-01 14:03 to 2023-09-02 10:03

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

r1 = o1.get_property('diameter')[0]
r2 = o2.get_property('diameter')[0]
plt.plot(np.median(r1*1e6, 1))
plt.plot(np.median(r2*1e6, 1))
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
/opt/conda/envs/opendrift/lib/python3.11/site-packages/numpy/core/fromnumeric.py:771: UserWarning: Warning: 'partition' will ignore the 'mask' of the MaskedArray.
  a.partition(kth, axis=axis, kind=kind, order=order)

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'])

14:03:47 WARNING opendrift.models.basemodel:3226: Plotting fast. This will make your plots less accurate.
/opt/conda/envs/opendrift/lib/python3.11/site-packages/cartopy/mpl/geoaxes.py:1696: UserWarning: No data for colormapping provided via 'c'. Parameters 'cmap' will be ignored
  result = super().scatter(*args, **kwargs)
14:03:48 INFO    opendrift.models.basemodel:5350: Saving animation to /root/project/docs/source/gallery/animations/example_oil_thickness_1.gif...
14:03:54 INFO    opendrift.models.basemodel:3848: Time to make animation: 0:00:07.658758
../_images/example_oil_thickness_1.gif

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

Gallery generated by Sphinx-Gallery