Analysing huge output files

import os
from datetime import datetime, timedelta
import opendrift
from opendrift.models.oceandrift import OceanDrift
from opendrift.readers import reader_oscillating

First make a simulation with two seedings, marked by origin_marker

o = OceanDrift(loglevel=50)
o.set_config('drift:horizontal_diffusivity', 10)
t1 = datetime.now()
t2 = t1 + timedelta(hours=6)
number = 10000
outfile = 'simulation.nc'  # Raw simulation output
reader_x = reader_oscillating.Reader('x_sea_water_velocity',
                amplitude=1, zero_time=t1)
reader_y = reader_oscillating.Reader('y_sea_water_velocity',
                amplitude=1, zero_time=t2)
o.add_reader([reader_x, reader_y])
o.seed_elements(time=t1, lon=4, lat=60, number=number,
                origin_marker=0)
o.seed_elements(time=[t1, t2], lon=4.2, lat=60.4, number=number,
                origin_marker=1)

o.run(duration=timedelta(hours=24),
      time_step=900, time_step_output=1800, outfile=outfile)
/opt/conda/envs/opendrift/lib/python3.11/site-packages/numpy/ma/core.py:467: RuntimeWarning: invalid value encountered in cast
  fill_value = np.array(fill_value, copy=False, dtype=ndtype)

Opening the output file lazily with Xarray. This will work even if the file is too large to fit in memory, as it will read and process data chuck-by-chunk directly from file using Dask. (See also example_river_runoff.py)

oa = opendrift.open_xarray(outfile)
20:11:41 DEBUG   opendrift.config:168: Adding 18 config items from __init__
20:11:41 DEBUG   opendrift.config:178:   Overwriting config item readers:max_number_of_fails
20:11:41 DEBUG   opendrift.config:168: Adding 5 config items from __init__
20:11:41 INFO    opendrift.models.basemodel:515: OpenDriftSimulation initialised (version 1.12.0 / v1.12.0-26-g390e945)
20:11:41 DEBUG   opendrift.config:168: Adding 15 config items from oceandrift
20:11:41 DEBUG   opendrift.config:178:   Overwriting config item seed:z
20:11:41 DEBUG   opendrift.export.io_netcdf:239: Importing with Xarray from simulation.nc
20:11:41 INFO    opendrift:118: Returning <class 'opendrift.models.oceandrift.OceanDrift'> object

Calculating histogram The histogram may be stored/cached to a netCDF file for later re-use, as the calculation may be time consuming for huge output files.

h = oa.get_histogram(pixelsize_m=500)
20:11:42 INFO    opendrift.models.basemodel:3892:       calculating for origin_marker 0...
20:11:42 INFO    opendrift.models.basemodel:3892:       calculating for origin_marker 1...

Plot the cumulative coverage of first seeding (origin_marker=0)

b=h.isel(origin_marker=0).sum(dim='time')
oa.plot(background=b.where(b>0), fast=True, show_elements=False, vmin=0, vmax=1000, clabel='First seeding')
OpenDrift - OceanDrift 2024-12-12 20:11 to 2024-12-13 20:11 UTC (49 steps)
20:11:42 DEBUG   opendrift.models.basemodel:2364: Setting up map: corners=None, fast=True, lscale=None
20:11:42 WARNING opendrift.models.basemodel:2410: Plotting fast. This will make your plots less accurate.

(<GeoAxes: title={'center': 'OpenDrift - OceanDrift\n2024-12-12 20:11 to 2024-12-13 20:11 UTC (49 steps)'}>, <Figure size 858.237x1100 with 2 Axes>)

Making two animations, for each of the two seedings / origin_markers. The calculated density fields may be stored/cached to a netCDF file for later re-use, as their calculation may be time consuming for huge output files. Note that other analysis/plotting methods are not yet adapted to datasets opened lazily with open_xarray

for om in [0, 1]:
    background=h.isel(origin_marker=om)
    oa.animation(background=background.where(background>0), bgalpha=1,
                corners=[4.0, 6, 59.5, 61], fast=False, show_elements=False, vmin=0, vmax=200)

# Cleaning up
os.remove(outfile)
20:11:44 DEBUG   opendrift.models.basemodel:2364: Setting up map: corners=[4.0, 6, 59.5, 61], fast=False, lscale=None
20:11:45 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
/opt/conda/envs/opendrift/lib/python3.11/site-packages/cartopy/mpl/geoaxes.py:1692: UserWarning: No data for colormapping provided via 'c'. Parameters 'cmap', 'vmin', 'vmax' will be ignored
  result = super().scatter(*args, **kwargs)
20:12:04 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:06 DEBUG   opendrift.models.basemodel:3044: Saving animation..
20:12:06 INFO    opendrift.models.basemodel:4613: Saving animation to /root/project/docs/source/gallery/animations/example_huge_output_0.gif...
20:12:06 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:07 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:09 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:11 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:12 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:14 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:16 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:17 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:19 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:21 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:23 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:24 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:26 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:27 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:29 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:31 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:33 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:34 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:36 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:38 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:40 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:42 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:44 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:45 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:48 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:50 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:52 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:53 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:55 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:56 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:12:58 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:00 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:01 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:03 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:05 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:06 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:08 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:09 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:11 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:12 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:14 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:15 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:17 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:18 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:20 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:21 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:23 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:24 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:26 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:27 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:29 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:30 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:32 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:33 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:34 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:36 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:37 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:39 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:40 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:42 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:43 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:45 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:47 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:48 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:49 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:51 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:52 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:54 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:55 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:57 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:13:58 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:00 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:01 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:03 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:04 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:05 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:07 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:09 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:11 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:12 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:14 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:15 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:17 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:18 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:20 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:22 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:24 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:25 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:27 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:28 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:30 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:31 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:33 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:34 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:36 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:38 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:40 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:41 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:44 DEBUG   opendrift.models.basemodel:4651: MPLBACKEND = agg
20:14:44 DEBUG   opendrift.models.basemodel:4652: DISPLAY = None
20:14:44 DEBUG   opendrift.models.basemodel:4653: Time to save animation: 0:02:37.930307
20:14:44 INFO    opendrift.models.basemodel:3037: Time to make animation: 0:02:59.182404
20:14:44 DEBUG   opendrift.models.basemodel:2364: Setting up map: corners=[4.0, 6, 59.5, 61], fast=False, lscale=None
20:14:44 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
/opt/conda/envs/opendrift/lib/python3.11/site-packages/cartopy/mpl/geoaxes.py:1692: UserWarning: No data for colormapping provided via 'c'. Parameters 'cmap', 'vmin', 'vmax' will be ignored
  result = super().scatter(*args, **kwargs)
20:14:46 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:47 DEBUG   opendrift.models.basemodel:3044: Saving animation..
20:14:47 INFO    opendrift.models.basemodel:4613: Saving animation to /root/project/docs/source/gallery/animations/example_huge_output_1.gif...
20:14:48 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:49 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:51 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:52 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:54 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:55 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:57 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:14:59 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:01 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:02 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:04 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:05 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:07 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:08 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:10 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:11 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:13 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:14 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:16 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:17 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:19 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:20 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:22 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:23 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:25 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:27 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:28 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:30 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:31 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:33 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:35 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:36 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:38 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:39 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:41 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:42 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:44 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:45 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:47 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:48 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:50 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:51 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:53 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:54 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:55 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:57 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:15:58 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:00 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:01 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:03 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:04 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:06 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:07 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:09 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:10 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:12 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:13 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:15 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:17 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:18 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:19 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:21 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:23 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:24 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:26 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:27 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:29 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:30 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:32 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:34 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:36 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:38 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:39 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:41 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:43 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:45 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:47 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:48 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:50 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:51 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:53 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:54 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:56 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:57 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:16:59 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:17:01 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:17:03 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:17:04 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:17:06 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:17:07 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:17:09 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:17:11 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:17:12 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:17:14 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:17:16 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:17:17 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:17:19 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:17:20 DEBUG   opendrift.readers.reader_global_landmask:78: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
20:17:23 DEBUG   opendrift.models.basemodel:4651: MPLBACKEND = agg
20:17:23 DEBUG   opendrift.models.basemodel:4652: DISPLAY = None
20:17:23 DEBUG   opendrift.models.basemodel:4653: Time to save animation: 0:02:35.468985
20:17:23 INFO    opendrift.models.basemodel:3037: Time to make animation: 0:02:39.345384

First seeding

../_images/example_huge_output_0.gif

Second seeding

../_images/example_huge_output_1.gif

Total running time of the script: (6 minutes 19.204 seconds)

Gallery generated by Sphinx-Gallery