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)
08:34:18 DEBUG   opendrift.config:161: Adding 16 config items from __init__
08:34:18 DEBUG   opendrift.config:171:   Overwriting config item readers:max_number_of_fails
08:34:18 DEBUG   opendrift.config:161: Adding 5 config items from __init__
08:34:18 INFO    opendrift.models.basemodel:533: OpenDriftSimulation initialised (version 1.11.2 / v1.11.2-45-g6369924)
08:34:18 DEBUG   opendrift.config:161: Adding 15 config items from oceandrift
08:34:18 DEBUG   opendrift.config:171:   Overwriting config item seed:z
08:34:18 DEBUG   opendrift.export.io_netcdf:226: Importing with Xarray from simulation.nc
08:34:18 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)
08:34:19 INFO    opendrift.models.basemodel:3854:       calculating for origin_marker 0...
08:34:19 INFO    opendrift.models.basemodel:3854:       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-05-23 08:33 to 2024-05-24 08:33 UTC (49 steps)
08:34:19 DEBUG   opendrift.models.basemodel:2333: Setting up map: corners=None, fast=True, lscale=None
08:34:19 WARNING opendrift.models.basemodel:2379: Plotting fast. This will make your plots less accurate.

(<GeoAxes: title={'center': 'OpenDrift - OceanDrift\n2024-05-23 08:33 to 2024-05-24 08:33 UTC (49 steps)'}>, <Figure size 856.033x1100 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)
08:34:22 DEBUG   opendrift.models.basemodel:2333: Setting up map: corners=[4.0, 6, 59.5, 61], fast=False, lscale=None
08:34:22 DEBUG   opendrift.readers.reader_global_landmask:83: 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:1683: UserWarning: No data for colormapping provided via 'c'. Parameters 'cmap', 'vmin', 'vmax' will be ignored
  result = super().scatter(*args, **kwargs)
08:34:40 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:34:42 DEBUG   opendrift.models.basemodel:3010: Saving animation..
08:34:42 INFO    opendrift.models.basemodel:4570: Saving animation to /root/project/docs/source/gallery/animations/example_huge_output_0.gif...
08:34:42 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:34:44 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:34:46 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:34:47 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:34:49 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:34:50 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:34:51 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:34:53 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:34:54 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:34:56 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:34:57 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:34:58 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:00 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:01 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:02 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:04 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:05 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:07 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:08 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:09 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:11 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:12 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:14 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:15 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:17 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:18 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:20 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:21 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:23 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:24 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:26 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:27 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:29 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:30 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:31 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:33 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:34 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:36 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:37 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:38 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:40 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:41 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:43 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:45 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:47 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:49 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:51 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:52 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:54 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:56 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:35:58 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:00 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:02 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:04 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:07 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:08 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:10 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:11 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:13 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:15 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:16 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:18 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:21 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:22 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:25 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:27 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:29 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:31 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:33 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:35 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:38 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:40 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:42 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:44 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:47 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:49 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:52 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:54 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:56 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:36:58 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:00 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:03 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:05 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:07 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:08 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:10 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:12 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:14 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:15 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:17 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:19 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:21 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:24 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:26 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:28 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:30 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:32 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:34 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:38 DEBUG   opendrift.models.basemodel:4608: MPLBACKEND = agg
08:37:38 DEBUG   opendrift.models.basemodel:4609: DISPLAY = None
08:37:38 DEBUG   opendrift.models.basemodel:4610: Time to save animation: 0:02:55.558562
08:37:38 INFO    opendrift.models.basemodel:3003: Time to make animation: 0:03:15.539569
08:37:38 DEBUG   opendrift.models.basemodel:2333: Setting up map: corners=[4.0, 6, 59.5, 61], fast=False, lscale=None
08:37:38 DEBUG   opendrift.readers.reader_global_landmask:83: 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:1683: UserWarning: No data for colormapping provided via 'c'. Parameters 'cmap', 'vmin', 'vmax' will be ignored
  result = super().scatter(*args, **kwargs)
08:37:40 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:42 DEBUG   opendrift.models.basemodel:3010: Saving animation..
08:37:43 INFO    opendrift.models.basemodel:4570: Saving animation to /root/project/docs/source/gallery/animations/example_huge_output_1.gif...
08:37:43 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:45 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:47 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:49 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:51 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:53 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:56 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:58 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:37:59 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:01 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:03 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:04 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:05 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:07 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:08 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:09 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:11 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:12 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:14 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:15 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:17 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:18 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:20 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:21 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:22 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:24 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:25 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:26 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:28 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:29 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:31 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:32 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:33 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:35 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:36 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:38 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:39 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:41 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:42 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:44 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:45 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:47 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:48 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:49 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:51 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:52 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:54 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:55 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:56 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:58 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:38:59 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:00 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:02 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:04 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:05 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:06 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:08 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:09 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:11 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:12 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:14 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:15 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:16 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:18 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:19 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:20 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:22 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:23 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:25 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:26 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:28 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:29 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:31 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:32 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:33 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:35 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:36 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:37 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:39 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:40 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:42 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:43 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:45 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:46 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:47 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:49 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:50 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:51 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:53 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:54 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:56 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:57 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:39:58 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:40:00 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:40:01 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:40:03 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:40:04 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:40:05 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
08:40:08 DEBUG   opendrift.models.basemodel:4608: MPLBACKEND = agg
08:40:08 DEBUG   opendrift.models.basemodel:4609: DISPLAY = None
08:40:08 DEBUG   opendrift.models.basemodel:4610: Time to save animation: 0:02:25.428993
08:40:08 INFO    opendrift.models.basemodel:3003: Time to make animation: 0:02:30.268990

First seeding

../_images/example_huge_output_0.gif

Second seeding

../_images/example_huge_output_1.gif

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

Gallery generated by Sphinx-Gallery