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')
12:59:13 INFO    opendrift.models.basemodel:515: OpenDriftSimulation initialised (version 1.11.13 / v1.11.13-99-gd2132d3)

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)
12:59:13 INFO    opendrift.models.openoil.adios.dirjs:90: Querying ADIOS database for oil: GENERIC DIESEL
12:59:13 INFO    opendrift.models.openoil.openoil:1719: Using density 841.2 and viscosity 6.3404392151148136e-06 of oiltype GENERIC DIESEL
12:59:13 INFO    opendrift.models.basemodel.environment:218: Adding a dynamical landmask with max. priority based on assumed maximum speed of 1.3 m/s. Adding a customised landmask may be faster...
12:59:17 INFO    opendrift.models.basemodel.environment:245: Fallback values will be used for the following variables which have no readers:
12:59:17 INFO    opendrift.models.basemodel.environment:248:    x_sea_water_velocity: 0.100000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    y_sea_water_velocity: 0.000000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    x_wind: 0.000000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    y_wind: 7.000000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    sea_surface_height: 0.000000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    upward_sea_water_velocity: 0.000000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    sea_surface_wave_significant_height: 0.000000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    sea_surface_wave_stokes_drift_y_velocity: 0.300000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    sea_ice_area_fraction: 0.000000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    sea_ice_x_velocity: 0.000000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    sea_ice_y_velocity: 0.000000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    sea_water_temperature: 10.000000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    sea_water_salinity: 34.000000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    sea_floor_depth_below_sea_level: 10000.000000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    ocean_vertical_diffusivity: 0.020000
12:59:17 INFO    opendrift.models.basemodel.environment:248:    ocean_mixed_layer_thickness: 50.000000
12:59:18 INFO    opendrift.models.basemodel:936: Using existing reader for land_binary_mask
12:59:18 INFO    opendrift.models.basemodel:947: All points are in ocean
12:59:18 INFO    opendrift.models.openoil.openoil:687: Oil-water surface tension is 0.029464 Nm
12:59:18 INFO    opendrift.models.openoil.openoil:700: Max water fraction not available for GENERIC DIESEL, using default
12:59:18 INFO    opendrift.models.basemodel:2038: 2024-11-11 12:59:13.351819 - step 1 of 120 - 10000 active elements (0 deactivated)
12:59:19 INFO    opendrift.models.basemodel:2038: 2024-11-11 13:09:13.351819 - step 2 of 120 - 10000 active elements (0 deactivated)
12:59:19 INFO    opendrift.models.basemodel:2038: 2024-11-11 13:19:13.351819 - step 3 of 120 - 10000 active elements (0 deactivated)
12:59:19 INFO    opendrift.models.basemodel:2038: 2024-11-11 13:29:13.351819 - step 4 of 120 - 10000 active elements (0 deactivated)
12:59:19 INFO    opendrift.models.basemodel:2038: 2024-11-11 13:39:13.351819 - step 5 of 120 - 10000 active elements (0 deactivated)
12:59:19 INFO    opendrift.models.basemodel:2038: 2024-11-11 13:49:13.351819 - step 6 of 120 - 10000 active elements (0 deactivated)
12:59:20 INFO    opendrift.models.basemodel:2038: 2024-11-11 13:59:13.351819 - step 7 of 120 - 10000 active elements (0 deactivated)
12:59:20 INFO    opendrift.models.basemodel:2038: 2024-11-11 14:09:13.351819 - step 8 of 120 - 10000 active elements (0 deactivated)
12:59:20 INFO    opendrift.models.basemodel:2038: 2024-11-11 14:19:13.351819 - step 9 of 120 - 10000 active elements (0 deactivated)
12:59:20 INFO    opendrift.models.basemodel:2038: 2024-11-11 14:29:13.351819 - step 10 of 120 - 10000 active elements (0 deactivated)
12:59:20 INFO    opendrift.models.basemodel:2038: 2024-11-11 14:39:13.351819 - step 11 of 120 - 10000 active elements (0 deactivated)
12:59:20 INFO    opendrift.models.basemodel:2038: 2024-11-11 14:49:13.351819 - step 12 of 120 - 10000 active elements (0 deactivated)
12:59:21 INFO    opendrift.models.basemodel:2038: 2024-11-11 14:59:13.351819 - step 13 of 120 - 10000 active elements (0 deactivated)
12:59:21 INFO    opendrift.models.basemodel:2038: 2024-11-11 15:09:13.351819 - step 14 of 120 - 10000 active elements (0 deactivated)
12:59:21 INFO    opendrift.models.basemodel:2038: 2024-11-11 15:19:13.351819 - step 15 of 120 - 10000 active elements (0 deactivated)
12:59:21 INFO    opendrift.models.basemodel:2038: 2024-11-11 15:29:13.351819 - step 16 of 120 - 10000 active elements (0 deactivated)
12:59:21 INFO    opendrift.models.basemodel:2038: 2024-11-11 15:39:13.351819 - step 17 of 120 - 10000 active elements (0 deactivated)
12:59:21 INFO    opendrift.models.basemodel:2038: 2024-11-11 15:49:13.351819 - step 18 of 120 - 10000 active elements (0 deactivated)
12:59:22 INFO    opendrift.models.basemodel:2038: 2024-11-11 15:59:13.351819 - step 19 of 120 - 10000 active elements (0 deactivated)
12:59:22 INFO    opendrift.models.basemodel:2038: 2024-11-11 16:09:13.351819 - step 20 of 120 - 10000 active elements (0 deactivated)
12:59:22 INFO    opendrift.models.basemodel:2038: 2024-11-11 16:19:13.351819 - step 21 of 120 - 10000 active elements (0 deactivated)
12:59:22 INFO    opendrift.models.basemodel:2038: 2024-11-11 16:29:13.351819 - step 22 of 120 - 10000 active elements (0 deactivated)
12:59:22 INFO    opendrift.models.basemodel:2038: 2024-11-11 16:39:13.351819 - step 23 of 120 - 10000 active elements (0 deactivated)
12:59:22 INFO    opendrift.models.basemodel:2038: 2024-11-11 16:49:13.351819 - step 24 of 120 - 10000 active elements (0 deactivated)
12:59:23 INFO    opendrift.models.basemodel:2038: 2024-11-11 16:59:13.351819 - step 25 of 120 - 10000 active elements (0 deactivated)
12:59:23 INFO    opendrift.models.basemodel:2038: 2024-11-11 17:09:13.351819 - step 26 of 120 - 10000 active elements (0 deactivated)
12:59:23 INFO    opendrift.models.basemodel:2038: 2024-11-11 17:19:13.351819 - step 27 of 120 - 10000 active elements (0 deactivated)
12:59:23 INFO    opendrift.models.basemodel:2038: 2024-11-11 17:29:13.351819 - step 28 of 120 - 10000 active elements (0 deactivated)
12:59:23 INFO    opendrift.models.basemodel:2038: 2024-11-11 17:39:13.351819 - step 29 of 120 - 10000 active elements (0 deactivated)
12:59:23 INFO    opendrift.models.basemodel:2038: 2024-11-11 17:49:13.351819 - step 30 of 120 - 10000 active elements (0 deactivated)
12:59:24 INFO    opendrift.models.basemodel:2038: 2024-11-11 17:59:13.351819 - step 31 of 120 - 10000 active elements (0 deactivated)
12:59:24 INFO    opendrift.models.basemodel:2038: 2024-11-11 18:09:13.351819 - step 32 of 120 - 10000 active elements (0 deactivated)
12:59:24 INFO    opendrift.models.basemodel:2038: 2024-11-11 18:19:13.351819 - step 33 of 120 - 10000 active elements (0 deactivated)
12:59:24 INFO    opendrift.models.basemodel:2038: 2024-11-11 18:29:13.351819 - step 34 of 120 - 10000 active elements (0 deactivated)
12:59:24 INFO    opendrift.models.basemodel:2038: 2024-11-11 18:39:13.351819 - step 35 of 120 - 10000 active elements (0 deactivated)
12:59:24 INFO    opendrift.models.basemodel:2038: 2024-11-11 18:49:13.351819 - step 36 of 120 - 10000 active elements (0 deactivated)
12:59:24 INFO    opendrift.models.basemodel:2038: 2024-11-11 18:59:13.351819 - step 37 of 120 - 10000 active elements (0 deactivated)
12:59:25 INFO    opendrift.models.basemodel:2038: 2024-11-11 19:09:13.351819 - step 38 of 120 - 10000 active elements (0 deactivated)
12:59:25 INFO    opendrift.models.basemodel:2038: 2024-11-11 19:19:13.351819 - step 39 of 120 - 10000 active elements (0 deactivated)
12:59:25 INFO    opendrift.models.basemodel:2038: 2024-11-11 19:29:13.351819 - step 40 of 120 - 10000 active elements (0 deactivated)
12:59:25 INFO    opendrift.models.basemodel:2038: 2024-11-11 19:39:13.351819 - step 41 of 120 - 10000 active elements (0 deactivated)
12:59:25 INFO    opendrift.models.basemodel:2038: 2024-11-11 19:49:13.351819 - step 42 of 120 - 10000 active elements (0 deactivated)
12:59:25 INFO    opendrift.models.basemodel:2038: 2024-11-11 19:59:13.351819 - step 43 of 120 - 10000 active elements (0 deactivated)
12:59:26 INFO    opendrift.models.basemodel:2038: 2024-11-11 20:09:13.351819 - step 44 of 120 - 10000 active elements (0 deactivated)
12:59:26 INFO    opendrift.models.basemodel:2038: 2024-11-11 20:19:13.351819 - step 45 of 120 - 10000 active elements (0 deactivated)
12:59:26 INFO    opendrift.models.basemodel:2038: 2024-11-11 20:29:13.351819 - step 46 of 120 - 10000 active elements (0 deactivated)
12:59:26 INFO    opendrift.models.basemodel:2038: 2024-11-11 20:39:13.351819 - step 47 of 120 - 10000 active elements (0 deactivated)
12:59:26 INFO    opendrift.models.basemodel:2038: 2024-11-11 20:49:13.351819 - step 48 of 120 - 10000 active elements (0 deactivated)
12:59:26 INFO    opendrift.models.basemodel:2038: 2024-11-11 20:59:13.351819 - step 49 of 120 - 10000 active elements (0 deactivated)
12:59:27 INFO    opendrift.models.basemodel:2038: 2024-11-11 21:09:13.351819 - step 50 of 120 - 10000 active elements (0 deactivated)
12:59:27 INFO    opendrift.models.basemodel:2038: 2024-11-11 21:19:13.351819 - step 51 of 120 - 10000 active elements (0 deactivated)
12:59:27 INFO    opendrift.models.basemodel:2038: 2024-11-11 21:29:13.351819 - step 52 of 120 - 10000 active elements (0 deactivated)
12:59:27 INFO    opendrift.models.basemodel:2038: 2024-11-11 21:39:13.351819 - step 53 of 120 - 10000 active elements (0 deactivated)
12:59:27 INFO    opendrift.models.basemodel:2038: 2024-11-11 21:49:13.351819 - step 54 of 120 - 10000 active elements (0 deactivated)
12:59:28 INFO    opendrift.models.basemodel:2038: 2024-11-11 21:59:13.351819 - step 55 of 120 - 10000 active elements (0 deactivated)
12:59:28 INFO    opendrift.models.basemodel:2038: 2024-11-11 22:09:13.351819 - step 56 of 120 - 10000 active elements (0 deactivated)
12:59:28 INFO    opendrift.models.basemodel:2038: 2024-11-11 22:19:13.351819 - step 57 of 120 - 10000 active elements (0 deactivated)
12:59:28 INFO    opendrift.models.basemodel:2038: 2024-11-11 22:29:13.351819 - step 58 of 120 - 10000 active elements (0 deactivated)
12:59:28 INFO    opendrift.models.basemodel:2038: 2024-11-11 22:39:13.351819 - step 59 of 120 - 10000 active elements (0 deactivated)
12:59:29 INFO    opendrift.models.basemodel:2038: 2024-11-11 22:49:13.351819 - step 60 of 120 - 10000 active elements (0 deactivated)
12:59:29 INFO    opendrift.models.basemodel:2038: 2024-11-11 22:59:13.351819 - step 61 of 120 - 10000 active elements (0 deactivated)
12:59:29 INFO    opendrift.models.basemodel:2038: 2024-11-11 23:09:13.351819 - step 62 of 120 - 10000 active elements (0 deactivated)
12:59:29 INFO    opendrift.models.basemodel:2038: 2024-11-11 23:19:13.351819 - step 63 of 120 - 10000 active elements (0 deactivated)
12:59:29 INFO    opendrift.models.basemodel:2038: 2024-11-11 23:29:13.351819 - step 64 of 120 - 10000 active elements (0 deactivated)
12:59:29 INFO    opendrift.models.basemodel:2038: 2024-11-11 23:39:13.351819 - step 65 of 120 - 10000 active elements (0 deactivated)
12:59:30 INFO    opendrift.models.basemodel:2038: 2024-11-11 23:49:13.351819 - step 66 of 120 - 10000 active elements (0 deactivated)
12:59:30 INFO    opendrift.models.basemodel:2038: 2024-11-11 23:59:13.351819 - step 67 of 120 - 10000 active elements (0 deactivated)
12:59:30 INFO    opendrift.models.basemodel:2038: 2024-11-12 00:09:13.351819 - step 68 of 120 - 10000 active elements (0 deactivated)
12:59:30 INFO    opendrift.models.basemodel:2038: 2024-11-12 00:19:13.351819 - step 69 of 120 - 10000 active elements (0 deactivated)
12:59:30 INFO    opendrift.models.basemodel:2038: 2024-11-12 00:29:13.351819 - step 70 of 120 - 10000 active elements (0 deactivated)
12:59:30 INFO    opendrift.models.basemodel:2038: 2024-11-12 00:39:13.351819 - step 71 of 120 - 10000 active elements (0 deactivated)
12:59:30 INFO    opendrift.models.basemodel:2038: 2024-11-12 00:49:13.351819 - step 72 of 120 - 10000 active elements (0 deactivated)
12:59:31 INFO    opendrift.models.basemodel:2038: 2024-11-12 00:59:13.351819 - step 73 of 120 - 10000 active elements (0 deactivated)
12:59:31 INFO    opendrift.models.basemodel:2038: 2024-11-12 01:09:13.351819 - step 74 of 120 - 10000 active elements (0 deactivated)
12:59:31 INFO    opendrift.models.basemodel:2038: 2024-11-12 01:19:13.351819 - step 75 of 120 - 10000 active elements (0 deactivated)
12:59:31 INFO    opendrift.models.basemodel:2038: 2024-11-12 01:29:13.351819 - step 76 of 120 - 10000 active elements (0 deactivated)
12:59:31 INFO    opendrift.models.basemodel:2038: 2024-11-12 01:39:13.351819 - step 77 of 120 - 10000 active elements (0 deactivated)
12:59:31 INFO    opendrift.models.basemodel:2038: 2024-11-12 01:49:13.351819 - step 78 of 120 - 10000 active elements (0 deactivated)
12:59:32 INFO    opendrift.models.basemodel:2038: 2024-11-12 01:59:13.351819 - step 79 of 120 - 10000 active elements (0 deactivated)
12:59:32 INFO    opendrift.models.basemodel:2038: 2024-11-12 02:09:13.351819 - step 80 of 120 - 10000 active elements (0 deactivated)
12:59:32 INFO    opendrift.models.basemodel:2038: 2024-11-12 02:19:13.351819 - step 81 of 120 - 10000 active elements (0 deactivated)
12:59:32 INFO    opendrift.models.basemodel:2038: 2024-11-12 02:29:13.351819 - step 82 of 120 - 10000 active elements (0 deactivated)
12:59:32 INFO    opendrift.models.basemodel:2038: 2024-11-12 02:39:13.351819 - step 83 of 120 - 10000 active elements (0 deactivated)
12:59:32 INFO    opendrift.models.basemodel:2038: 2024-11-12 02:49:13.351819 - step 84 of 120 - 10000 active elements (0 deactivated)
12:59:33 INFO    opendrift.models.basemodel:2038: 2024-11-12 02:59:13.351819 - step 85 of 120 - 10000 active elements (0 deactivated)
12:59:33 INFO    opendrift.models.basemodel:2038: 2024-11-12 03:09:13.351819 - step 86 of 120 - 10000 active elements (0 deactivated)
12:59:33 INFO    opendrift.models.basemodel:2038: 2024-11-12 03:19:13.351819 - step 87 of 120 - 10000 active elements (0 deactivated)
12:59:33 INFO    opendrift.models.basemodel:2038: 2024-11-12 03:29:13.351819 - step 88 of 120 - 10000 active elements (0 deactivated)
12:59:33 INFO    opendrift.models.basemodel:2038: 2024-11-12 03:39:13.351819 - step 89 of 120 - 10000 active elements (0 deactivated)
12:59:34 INFO    opendrift.models.basemodel:2038: 2024-11-12 03:49:13.351819 - step 90 of 120 - 10000 active elements (0 deactivated)
12:59:34 INFO    opendrift.models.basemodel:2038: 2024-11-12 03:59:13.351819 - step 91 of 120 - 10000 active elements (0 deactivated)
12:59:34 INFO    opendrift.models.basemodel:2038: 2024-11-12 04:09:13.351819 - step 92 of 120 - 10000 active elements (0 deactivated)
12:59:34 INFO    opendrift.models.basemodel:2038: 2024-11-12 04:19:13.351819 - step 93 of 120 - 10000 active elements (0 deactivated)
12:59:34 INFO    opendrift.models.basemodel:2038: 2024-11-12 04:29:13.351819 - step 94 of 120 - 10000 active elements (0 deactivated)
12:59:35 INFO    opendrift.models.basemodel:2038: 2024-11-12 04:39:13.351819 - step 95 of 120 - 10000 active elements (0 deactivated)
12:59:35 INFO    opendrift.models.basemodel:2038: 2024-11-12 04:49:13.351819 - step 96 of 120 - 10000 active elements (0 deactivated)
12:59:35 INFO    opendrift.models.basemodel:2038: 2024-11-12 04:59:13.351819 - step 97 of 120 - 10000 active elements (0 deactivated)
12:59:35 INFO    opendrift.models.basemodel:2038: 2024-11-12 05:09:13.351819 - step 98 of 120 - 10000 active elements (0 deactivated)
12:59:35 INFO    opendrift.models.basemodel:2038: 2024-11-12 05:19:13.351819 - step 99 of 120 - 10000 active elements (0 deactivated)
12:59:35 INFO    opendrift.models.basemodel:2038: 2024-11-12 05:29:13.351819 - step 100 of 120 - 10000 active elements (0 deactivated)
12:59:36 INFO    opendrift.models.basemodel:2038: 2024-11-12 05:39:13.351819 - step 101 of 120 - 10000 active elements (0 deactivated)
12:59:36 INFO    opendrift.models.basemodel:2038: 2024-11-12 05:49:13.351819 - step 102 of 120 - 10000 active elements (0 deactivated)
12:59:36 INFO    opendrift.models.basemodel:2038: 2024-11-12 05:59:13.351819 - step 103 of 120 - 10000 active elements (0 deactivated)
12:59:36 INFO    opendrift.models.basemodel:2038: 2024-11-12 06:09:13.351819 - step 104 of 120 - 10000 active elements (0 deactivated)
12:59:36 INFO    opendrift.models.basemodel:2038: 2024-11-12 06:19:13.351819 - step 105 of 120 - 10000 active elements (0 deactivated)
12:59:37 INFO    opendrift.models.basemodel:2038: 2024-11-12 06:29:13.351819 - step 106 of 120 - 10000 active elements (0 deactivated)
12:59:37 INFO    opendrift.models.basemodel:2038: 2024-11-12 06:39:13.351819 - step 107 of 120 - 10000 active elements (0 deactivated)
12:59:37 INFO    opendrift.models.basemodel:2038: 2024-11-12 06:49:13.351819 - step 108 of 120 - 10000 active elements (0 deactivated)
12:59:37 INFO    opendrift.models.basemodel:2038: 2024-11-12 06:59:13.351819 - step 109 of 120 - 10000 active elements (0 deactivated)
12:59:37 INFO    opendrift.models.basemodel:2038: 2024-11-12 07:09:13.351819 - step 110 of 120 - 10000 active elements (0 deactivated)
12:59:37 INFO    opendrift.models.basemodel:2038: 2024-11-12 07:19:13.351819 - step 111 of 120 - 10000 active elements (0 deactivated)
12:59:38 INFO    opendrift.models.basemodel:2038: 2024-11-12 07:29:13.351819 - step 112 of 120 - 10000 active elements (0 deactivated)
12:59:38 INFO    opendrift.models.basemodel:2038: 2024-11-12 07:39:13.351819 - step 113 of 120 - 10000 active elements (0 deactivated)
12:59:38 INFO    opendrift.models.basemodel:2038: 2024-11-12 07:49:13.351819 - step 114 of 120 - 10000 active elements (0 deactivated)
12:59:38 INFO    opendrift.models.basemodel:2038: 2024-11-12 07:59:13.351819 - step 115 of 120 - 10000 active elements (0 deactivated)
12:59:38 INFO    opendrift.models.basemodel:2038: 2024-11-12 08:09:13.351819 - step 116 of 120 - 10000 active elements (0 deactivated)
12:59:39 INFO    opendrift.models.basemodel:2038: 2024-11-12 08:19:13.351819 - step 117 of 120 - 10000 active elements (0 deactivated)
12:59:39 INFO    opendrift.models.basemodel:2038: 2024-11-12 08:29:13.351819 - step 118 of 120 - 10000 active elements (0 deactivated)
12:59:39 INFO    opendrift.models.basemodel:2038: 2024-11-12 08:39:13.351819 - step 119 of 120 - 10000 active elements (0 deactivated)
12:59:39 INFO    opendrift.models.basemodel:2038: 2024-11-12 08:49:13.351819 - 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)
12:59:40 WARNING opendrift.models.basemodel:2411: Plotting fast. This will make your plots less accurate.
12:59:43 INFO    opendrift.models.basemodel:4613: Saving animation to /root/project/docs/source/gallery/animations/example_oil_thickness_0.gif...
13:00:00 INFO    opendrift.models.basemodel:3037: Time to make animation: 0:00:20.941283
../_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)
13:00:01 INFO    opendrift.models.basemodel:515: OpenDriftSimulation initialised (version 1.11.13 / v1.11.13-99-gd2132d3)
13:00:01 INFO    opendrift.models.openoil.adios.dirjs:90: Querying ADIOS database for oil: GENERIC DIESEL
13:00:01 INFO    opendrift.models.openoil.openoil:1719: Using density 841.2 and viscosity 6.3404392151148136e-06 of oiltype GENERIC DIESEL
13:00:01 INFO    opendrift.models.basemodel.environment:218: Adding a dynamical landmask with max. priority based on assumed maximum speed of 1.3 m/s. Adding a customised landmask may be faster...
13:00:01 INFO    opendrift.models.basemodel.environment:245: Fallback values will be used for the following variables which have no readers:
13:00:01 INFO    opendrift.models.basemodel.environment:248:    x_sea_water_velocity: 0.100000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    y_sea_water_velocity: 0.000000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    x_wind: 0.000000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    y_wind: 7.000000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    sea_surface_height: 0.000000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    upward_sea_water_velocity: 0.000000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    sea_surface_wave_significant_height: 0.000000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    sea_surface_wave_stokes_drift_x_velocity: 0.000000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    sea_surface_wave_stokes_drift_y_velocity: 0.300000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    sea_surface_wave_period_at_variance_spectral_density_maximum: 0.000000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment: 0.000000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    sea_ice_area_fraction: 0.000000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    sea_ice_x_velocity: 0.000000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    sea_ice_y_velocity: 0.000000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    sea_water_temperature: 10.000000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    sea_water_salinity: 34.000000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    sea_floor_depth_below_sea_level: 10000.000000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    ocean_vertical_diffusivity: 0.020000
13:00:01 INFO    opendrift.models.basemodel.environment:248:    ocean_mixed_layer_thickness: 50.000000
13:00:01 INFO    opendrift.models.basemodel:936: Using existing reader for land_binary_mask
13:00:01 INFO    opendrift.models.basemodel:947: All points are in ocean
13:00:01 INFO    opendrift.models.openoil.openoil:687: Oil-water surface tension is 0.029464 Nm
13:00:01 INFO    opendrift.models.openoil.openoil:700: Max water fraction not available for GENERIC DIESEL, using default
13:00:01 INFO    opendrift.models.basemodel:2038: 2024-11-11 13:00:01.048759 - step 1 of 120 - 10000 active elements (0 deactivated)
13:00:01 INFO    opendrift.models.basemodel:2038: 2024-11-11 13:10:01.048759 - step 2 of 120 - 10000 active elements (0 deactivated)
13:00:01 INFO    opendrift.models.basemodel:2038: 2024-11-11 13:20:01.048759 - step 3 of 120 - 10000 active elements (0 deactivated)
13:00:02 INFO    opendrift.models.basemodel:2038: 2024-11-11 13:30:01.048759 - step 4 of 120 - 10000 active elements (0 deactivated)
13:00:02 INFO    opendrift.models.basemodel:2038: 2024-11-11 13:40:01.048759 - step 5 of 120 - 10000 active elements (0 deactivated)
13:00:02 INFO    opendrift.models.basemodel:2038: 2024-11-11 13:50:01.048759 - step 6 of 120 - 10000 active elements (0 deactivated)
13:00:02 INFO    opendrift.models.basemodel:2038: 2024-11-11 14:00:01.048759 - step 7 of 120 - 10000 active elements (0 deactivated)
13:00:03 INFO    opendrift.models.basemodel:2038: 2024-11-11 14:10:01.048759 - step 8 of 120 - 10000 active elements (0 deactivated)
13:00:03 INFO    opendrift.models.basemodel:2038: 2024-11-11 14:20:01.048759 - step 9 of 120 - 10000 active elements (0 deactivated)
13:00:03 INFO    opendrift.models.basemodel:2038: 2024-11-11 14:30:01.048759 - step 10 of 120 - 10000 active elements (0 deactivated)
13:00:03 INFO    opendrift.models.basemodel:2038: 2024-11-11 14:40:01.048759 - step 11 of 120 - 10000 active elements (0 deactivated)
13:00:04 INFO    opendrift.models.basemodel:2038: 2024-11-11 14:50:01.048759 - step 12 of 120 - 10000 active elements (0 deactivated)
13:00:04 INFO    opendrift.models.basemodel:2038: 2024-11-11 15:00:01.048759 - step 13 of 120 - 10000 active elements (0 deactivated)
13:00:04 INFO    opendrift.models.basemodel:2038: 2024-11-11 15:10:01.048759 - step 14 of 120 - 10000 active elements (0 deactivated)
13:00:04 INFO    opendrift.models.basemodel:2038: 2024-11-11 15:20:01.048759 - step 15 of 120 - 10000 active elements (0 deactivated)
13:00:05 INFO    opendrift.models.basemodel:2038: 2024-11-11 15:30:01.048759 - step 16 of 120 - 10000 active elements (0 deactivated)
13:00:05 INFO    opendrift.models.basemodel:2038: 2024-11-11 15:40:01.048759 - step 17 of 120 - 10000 active elements (0 deactivated)
13:00:05 INFO    opendrift.models.basemodel:2038: 2024-11-11 15:50:01.048759 - step 18 of 120 - 10000 active elements (0 deactivated)
13:00:05 INFO    opendrift.models.basemodel:2038: 2024-11-11 16:00:01.048759 - step 19 of 120 - 10000 active elements (0 deactivated)
13:00:06 INFO    opendrift.models.basemodel:2038: 2024-11-11 16:10:01.048759 - step 20 of 120 - 10000 active elements (0 deactivated)
13:00:06 INFO    opendrift.models.basemodel:2038: 2024-11-11 16:20:01.048759 - step 21 of 120 - 10000 active elements (0 deactivated)
13:00:06 INFO    opendrift.models.basemodel:2038: 2024-11-11 16:30:01.048759 - step 22 of 120 - 10000 active elements (0 deactivated)
13:00:06 INFO    opendrift.models.basemodel:2038: 2024-11-11 16:40:01.048759 - step 23 of 120 - 10000 active elements (0 deactivated)
13:00:06 INFO    opendrift.models.basemodel:2038: 2024-11-11 16:50:01.048759 - step 24 of 120 - 10000 active elements (0 deactivated)
13:00:07 INFO    opendrift.models.basemodel:2038: 2024-11-11 17:00:01.048759 - step 25 of 120 - 10000 active elements (0 deactivated)
13:00:07 INFO    opendrift.models.basemodel:2038: 2024-11-11 17:10:01.048759 - step 26 of 120 - 10000 active elements (0 deactivated)
13:00:07 INFO    opendrift.models.basemodel:2038: 2024-11-11 17:20:01.048759 - step 27 of 120 - 10000 active elements (0 deactivated)
13:00:07 INFO    opendrift.models.basemodel:2038: 2024-11-11 17:30:01.048759 - step 28 of 120 - 10000 active elements (0 deactivated)
13:00:08 INFO    opendrift.models.basemodel:2038: 2024-11-11 17:40:01.048759 - step 29 of 120 - 10000 active elements (0 deactivated)
13:00:08 INFO    opendrift.models.basemodel:2038: 2024-11-11 17:50:01.048759 - step 30 of 120 - 10000 active elements (0 deactivated)
13:00:08 INFO    opendrift.models.basemodel:2038: 2024-11-11 18:00:01.048759 - step 31 of 120 - 10000 active elements (0 deactivated)
13:00:08 INFO    opendrift.models.basemodel:2038: 2024-11-11 18:10:01.048759 - step 32 of 120 - 10000 active elements (0 deactivated)
13:00:08 INFO    opendrift.models.basemodel:2038: 2024-11-11 18:20:01.048759 - step 33 of 120 - 10000 active elements (0 deactivated)
13:00:09 INFO    opendrift.models.basemodel:2038: 2024-11-11 18:30:01.048759 - step 34 of 120 - 10000 active elements (0 deactivated)
13:00:09 INFO    opendrift.models.basemodel:2038: 2024-11-11 18:40:01.048759 - step 35 of 120 - 10000 active elements (0 deactivated)
13:00:09 INFO    opendrift.models.basemodel:2038: 2024-11-11 18:50:01.048759 - step 36 of 120 - 10000 active elements (0 deactivated)
13:00:09 INFO    opendrift.models.basemodel:2038: 2024-11-11 19:00:01.048759 - step 37 of 120 - 10000 active elements (0 deactivated)
13:00:09 INFO    opendrift.models.basemodel:2038: 2024-11-11 19:10:01.048759 - step 38 of 120 - 10000 active elements (0 deactivated)
13:00:10 INFO    opendrift.models.basemodel:2038: 2024-11-11 19:20:01.048759 - step 39 of 120 - 10000 active elements (0 deactivated)
13:00:10 INFO    opendrift.models.basemodel:2038: 2024-11-11 19:30:01.048759 - step 40 of 120 - 10000 active elements (0 deactivated)
13:00:10 INFO    opendrift.models.basemodel:2038: 2024-11-11 19:40:01.048759 - step 41 of 120 - 10000 active elements (0 deactivated)
13:00:10 INFO    opendrift.models.basemodel:2038: 2024-11-11 19:50:01.048759 - step 42 of 120 - 10000 active elements (0 deactivated)
13:00:10 INFO    opendrift.models.basemodel:2038: 2024-11-11 20:00:01.048759 - step 43 of 120 - 10000 active elements (0 deactivated)
13:00:11 INFO    opendrift.models.basemodel:2038: 2024-11-11 20:10:01.048759 - step 44 of 120 - 10000 active elements (0 deactivated)
13:00:11 INFO    opendrift.models.basemodel:2038: 2024-11-11 20:20:01.048759 - step 45 of 120 - 10000 active elements (0 deactivated)
13:00:11 INFO    opendrift.models.basemodel:2038: 2024-11-11 20:30:01.048759 - step 46 of 120 - 10000 active elements (0 deactivated)
13:00:11 INFO    opendrift.models.basemodel:2038: 2024-11-11 20:40:01.048759 - step 47 of 120 - 10000 active elements (0 deactivated)
13:00:11 INFO    opendrift.models.basemodel:2038: 2024-11-11 20:50:01.048759 - step 48 of 120 - 10000 active elements (0 deactivated)
13:00:11 INFO    opendrift.models.basemodel:2038: 2024-11-11 21:00:01.048759 - step 49 of 120 - 10000 active elements (0 deactivated)
13:00:12 INFO    opendrift.models.basemodel:2038: 2024-11-11 21:10:01.048759 - step 50 of 120 - 10000 active elements (0 deactivated)
13:00:12 INFO    opendrift.models.basemodel:2038: 2024-11-11 21:20:01.048759 - step 51 of 120 - 10000 active elements (0 deactivated)
13:00:12 INFO    opendrift.models.basemodel:2038: 2024-11-11 21:30:01.048759 - step 52 of 120 - 10000 active elements (0 deactivated)
13:00:12 INFO    opendrift.models.basemodel:2038: 2024-11-11 21:40:01.048759 - step 53 of 120 - 10000 active elements (0 deactivated)
13:00:13 INFO    opendrift.models.basemodel:2038: 2024-11-11 21:50:01.048759 - step 54 of 120 - 10000 active elements (0 deactivated)
13:00:13 INFO    opendrift.models.basemodel:2038: 2024-11-11 22:00:01.048759 - step 55 of 120 - 10000 active elements (0 deactivated)
13:00:13 INFO    opendrift.models.basemodel:2038: 2024-11-11 22:10:01.048759 - step 56 of 120 - 10000 active elements (0 deactivated)
13:00:13 INFO    opendrift.models.basemodel:2038: 2024-11-11 22:20:01.048759 - step 57 of 120 - 10000 active elements (0 deactivated)
13:00:13 INFO    opendrift.models.basemodel:2038: 2024-11-11 22:30:01.048759 - step 58 of 120 - 10000 active elements (0 deactivated)
13:00:14 INFO    opendrift.models.basemodel:2038: 2024-11-11 22:40:01.048759 - step 59 of 120 - 10000 active elements (0 deactivated)
13:00:14 INFO    opendrift.models.basemodel:2038: 2024-11-11 22:50:01.048759 - step 60 of 120 - 10000 active elements (0 deactivated)
13:00:14 INFO    opendrift.models.basemodel:2038: 2024-11-11 23:00:01.048759 - step 61 of 120 - 10000 active elements (0 deactivated)
13:00:14 INFO    opendrift.models.basemodel:2038: 2024-11-11 23:10:01.048759 - step 62 of 120 - 10000 active elements (0 deactivated)
13:00:14 INFO    opendrift.models.basemodel:2038: 2024-11-11 23:20:01.048759 - step 63 of 120 - 10000 active elements (0 deactivated)
13:00:14 INFO    opendrift.models.basemodel:2038: 2024-11-11 23:30:01.048759 - step 64 of 120 - 10000 active elements (0 deactivated)
13:00:15 INFO    opendrift.models.basemodel:2038: 2024-11-11 23:40:01.048759 - step 65 of 120 - 10000 active elements (0 deactivated)
13:00:15 INFO    opendrift.models.basemodel:2038: 2024-11-11 23:50:01.048759 - step 66 of 120 - 10000 active elements (0 deactivated)
13:00:15 INFO    opendrift.models.basemodel:2038: 2024-11-12 00:00:01.048759 - step 67 of 120 - 10000 active elements (0 deactivated)
13:00:15 INFO    opendrift.models.basemodel:2038: 2024-11-12 00:10:01.048759 - step 68 of 120 - 10000 active elements (0 deactivated)
13:00:15 INFO    opendrift.models.basemodel:2038: 2024-11-12 00:20:01.048759 - step 69 of 120 - 10000 active elements (0 deactivated)
13:00:16 INFO    opendrift.models.basemodel:2038: 2024-11-12 00:30:01.048759 - step 70 of 120 - 10000 active elements (0 deactivated)
13:00:16 INFO    opendrift.models.basemodel:2038: 2024-11-12 00:40:01.048759 - step 71 of 120 - 10000 active elements (0 deactivated)
13:00:16 INFO    opendrift.models.basemodel:2038: 2024-11-12 00:50:01.048759 - step 72 of 120 - 10000 active elements (0 deactivated)
13:00:16 INFO    opendrift.models.basemodel:2038: 2024-11-12 01:00:01.048759 - step 73 of 120 - 10000 active elements (0 deactivated)
13:00:16 INFO    opendrift.models.basemodel:2038: 2024-11-12 01:10:01.048759 - step 74 of 120 - 10000 active elements (0 deactivated)
13:00:16 INFO    opendrift.models.basemodel:2038: 2024-11-12 01:20:01.048759 - step 75 of 120 - 10000 active elements (0 deactivated)
13:00:17 INFO    opendrift.models.basemodel:2038: 2024-11-12 01:30:01.048759 - step 76 of 120 - 10000 active elements (0 deactivated)
13:00:17 INFO    opendrift.models.basemodel:2038: 2024-11-12 01:40:01.048759 - step 77 of 120 - 10000 active elements (0 deactivated)
13:00:17 INFO    opendrift.models.basemodel:2038: 2024-11-12 01:50:01.048759 - step 78 of 120 - 10000 active elements (0 deactivated)
13:00:17 INFO    opendrift.models.basemodel:2038: 2024-11-12 02:00:01.048759 - step 79 of 120 - 10000 active elements (0 deactivated)
13:00:17 INFO    opendrift.models.basemodel:2038: 2024-11-12 02:10:01.048759 - step 80 of 120 - 10000 active elements (0 deactivated)
13:00:17 INFO    opendrift.models.basemodel:2038: 2024-11-12 02:20:01.048759 - step 81 of 120 - 10000 active elements (0 deactivated)
13:00:18 INFO    opendrift.models.basemodel:2038: 2024-11-12 02:30:01.048759 - step 82 of 120 - 10000 active elements (0 deactivated)
13:00:18 INFO    opendrift.models.basemodel:2038: 2024-11-12 02:40:01.048759 - step 83 of 120 - 10000 active elements (0 deactivated)
13:00:18 INFO    opendrift.models.basemodel:2038: 2024-11-12 02:50:01.048759 - step 84 of 120 - 10000 active elements (0 deactivated)
13:00:18 INFO    opendrift.models.basemodel:2038: 2024-11-12 03:00:01.048759 - step 85 of 120 - 10000 active elements (0 deactivated)
13:00:18 INFO    opendrift.models.basemodel:2038: 2024-11-12 03:10:01.048759 - step 86 of 120 - 10000 active elements (0 deactivated)
13:00:19 INFO    opendrift.models.basemodel:2038: 2024-11-12 03:20:01.048759 - step 87 of 120 - 10000 active elements (0 deactivated)
13:00:19 INFO    opendrift.models.basemodel:2038: 2024-11-12 03:30:01.048759 - step 88 of 120 - 10000 active elements (0 deactivated)
13:00:19 INFO    opendrift.models.basemodel:2038: 2024-11-12 03:40:01.048759 - step 89 of 120 - 10000 active elements (0 deactivated)
13:00:19 INFO    opendrift.models.basemodel:2038: 2024-11-12 03:50:01.048759 - step 90 of 120 - 10000 active elements (0 deactivated)
13:00:19 INFO    opendrift.models.basemodel:2038: 2024-11-12 04:00:01.048759 - step 91 of 120 - 10000 active elements (0 deactivated)
13:00:19 INFO    opendrift.models.basemodel:2038: 2024-11-12 04:10:01.048759 - step 92 of 120 - 10000 active elements (0 deactivated)
13:00:20 INFO    opendrift.models.basemodel:2038: 2024-11-12 04:20:01.048759 - step 93 of 120 - 10000 active elements (0 deactivated)
13:00:20 INFO    opendrift.models.basemodel:2038: 2024-11-12 04:30:01.048759 - step 94 of 120 - 10000 active elements (0 deactivated)
13:00:20 INFO    opendrift.models.basemodel:2038: 2024-11-12 04:40:01.048759 - step 95 of 120 - 10000 active elements (0 deactivated)
13:00:20 INFO    opendrift.models.basemodel:2038: 2024-11-12 04:50:01.048759 - step 96 of 120 - 10000 active elements (0 deactivated)
13:00:20 INFO    opendrift.models.basemodel:2038: 2024-11-12 05:00:01.048759 - step 97 of 120 - 10000 active elements (0 deactivated)
13:00:20 INFO    opendrift.models.basemodel:2038: 2024-11-12 05:10:01.048759 - step 98 of 120 - 10000 active elements (0 deactivated)
13:00:21 INFO    opendrift.models.basemodel:2038: 2024-11-12 05:20:01.048759 - step 99 of 120 - 10000 active elements (0 deactivated)
13:00:21 INFO    opendrift.models.basemodel:2038: 2024-11-12 05:30:01.048759 - step 100 of 120 - 10000 active elements (0 deactivated)
13:00:21 INFO    opendrift.models.basemodel:2038: 2024-11-12 05:40:01.048759 - step 101 of 120 - 10000 active elements (0 deactivated)
13:00:21 INFO    opendrift.models.basemodel:2038: 2024-11-12 05:50:01.048759 - step 102 of 120 - 10000 active elements (0 deactivated)
13:00:21 INFO    opendrift.models.basemodel:2038: 2024-11-12 06:00:01.048759 - step 103 of 120 - 10000 active elements (0 deactivated)
13:00:21 INFO    opendrift.models.basemodel:2038: 2024-11-12 06:10:01.048759 - step 104 of 120 - 10000 active elements (0 deactivated)
13:00:22 INFO    opendrift.models.basemodel:2038: 2024-11-12 06:20:01.048759 - step 105 of 120 - 10000 active elements (0 deactivated)
13:00:22 INFO    opendrift.models.basemodel:2038: 2024-11-12 06:30:01.048759 - step 106 of 120 - 10000 active elements (0 deactivated)
13:00:22 INFO    opendrift.models.basemodel:2038: 2024-11-12 06:40:01.048759 - step 107 of 120 - 10000 active elements (0 deactivated)
13:00:22 INFO    opendrift.models.basemodel:2038: 2024-11-12 06:50:01.048759 - step 108 of 120 - 10000 active elements (0 deactivated)
13:00:22 INFO    opendrift.models.basemodel:2038: 2024-11-12 07:00:01.048759 - step 109 of 120 - 10000 active elements (0 deactivated)
13:00:23 INFO    opendrift.models.basemodel:2038: 2024-11-12 07:10:01.048759 - step 110 of 120 - 10000 active elements (0 deactivated)
13:00:23 INFO    opendrift.models.basemodel:2038: 2024-11-12 07:20:01.048759 - step 111 of 120 - 10000 active elements (0 deactivated)
13:00:23 INFO    opendrift.models.basemodel:2038: 2024-11-12 07:30:01.048759 - step 112 of 120 - 10000 active elements (0 deactivated)
13:00:23 INFO    opendrift.models.basemodel:2038: 2024-11-12 07:40:01.048759 - step 113 of 120 - 10000 active elements (0 deactivated)
13:00:23 INFO    opendrift.models.basemodel:2038: 2024-11-12 07:50:01.048759 - step 114 of 120 - 10000 active elements (0 deactivated)
13:00:23 INFO    opendrift.models.basemodel:2038: 2024-11-12 08:00:01.048759 - step 115 of 120 - 10000 active elements (0 deactivated)
13:00:24 INFO    opendrift.models.basemodel:2038: 2024-11-12 08:10:01.048759 - step 116 of 120 - 10000 active elements (0 deactivated)
13:00:24 INFO    opendrift.models.basemodel:2038: 2024-11-12 08:20:01.048759 - step 117 of 120 - 10000 active elements (0 deactivated)
13:00:24 INFO    opendrift.models.basemodel:2038: 2024-11-12 08:30:01.048759 - step 118 of 120 - 10000 active elements (0 deactivated)
13:00:24 INFO    opendrift.models.basemodel:2038: 2024-11-12 08:40:01.048759 - step 119 of 120 - 10000 active elements (0 deactivated)
13:00:24 INFO    opendrift.models.basemodel:2038: 2024-11-12 08:50:01.048759 - 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) - 2024-11-11 12:59 to 2024-11-12 08:59
  • GENERIC DIESEL (841.2 kg/m3) - 2024-11-11 13:00 to 2024-11-12 09:00

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'])
13:00:35 WARNING opendrift.models.basemodel:2411: Plotting fast. This will make your plots less accurate.
/opt/conda/envs/opendrift/lib/python3.11/site-packages/cartopy/mpl/geoaxes.py:1692: UserWarning: No data for colormapping provided via 'c'. Parameters 'cmap' will be ignored
  result = super().scatter(*args, **kwargs)
13:00:37 INFO    opendrift.models.basemodel:4613: Saving animation to /root/project/docs/source/gallery/animations/example_oil_thickness_1.gif...
13:00:58 INFO    opendrift.models.basemodel:3037: Time to make animation: 0:00:23.729973
../_images/example_oil_thickness_1.gif

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

Gallery generated by Sphinx-Gallery