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