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)
07:40:45 DEBUG   opendrift.config:168: Adding 18 config items from __init__
07:40:45 DEBUG   opendrift.config:178:   Overwriting config item readers:max_number_of_fails
07:40:45 DEBUG   opendrift.config:168: Adding 5 config items from __init__
07:40:45 INFO    opendrift.models.basemodel:515: OpenDriftSimulation initialised (version 1.11.13 / v1.11.13)
07:40:45 DEBUG   opendrift.config:168: Adding 15 config items from oceandrift
07:40:45 DEBUG   opendrift.config:178:   Overwriting config item seed:z
07:40:45 DEBUG   opendrift.export.io_netcdf:226: Importing with Xarray from simulation.nc
07:40:45 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)
07:40:46 INFO    opendrift.models.basemodel:3892:       calculating for origin_marker 0...
07:40:46 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-07-24 07:40 to 2024-07-25 07:40 UTC (49 steps)
07:40:46 DEBUG   opendrift.models.basemodel:2365: Setting up map: corners=None, fast=True, lscale=None
07:40:46 WARNING opendrift.models.basemodel:2411: Plotting fast. This will make your plots less accurate.

(<GeoAxes: title={'center': 'OpenDrift - OceanDrift\n2024-07-24 07:40 to 2024-07-25 07:40 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)
07:40:48 DEBUG   opendrift.models.basemodel:2365: Setting up map: corners=[4.0, 6, 59.5, 61], fast=False, lscale=None
07:40:49 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)
07:41:08 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:09 DEBUG   opendrift.models.basemodel:3044: Saving animation..
07:41:09 INFO    opendrift.models.basemodel:4608: Saving animation to /root/project/docs/source/gallery/animations/example_huge_output_0.gif...
07:41:09 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:10 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:12 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:13 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:15 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:17 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:18 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:20 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:22 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:23 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:25 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:26 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:28 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:29 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:31 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:32 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:34 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:35 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:37 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:38 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:39 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:41 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:42 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:44 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:45 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:46 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:48 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:49 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:51 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:52 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:54 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:55 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:56 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:58 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:41:59 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:00 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:02 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:03 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:05 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:06 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:08 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:09 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:11 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:12 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:14 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:15 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:17 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:18 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:20 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:21 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:23 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:24 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:26 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:27 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:29 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:31 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:32 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:34 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:35 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:37 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:39 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:40 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:42 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:43 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:45 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:47 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:48 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:50 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:52 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:53 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:55 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:56 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:58 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:42:59 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:01 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:03 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:05 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:07 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:09 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:10 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:12 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:13 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:15 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:16 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:18 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:19 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:21 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:22 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:24 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:25 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:27 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:28 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:30 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:31 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:33 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:34 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:36 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:37 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:40 DEBUG   opendrift.models.basemodel:4646: MPLBACKEND = agg
07:43:40 DEBUG   opendrift.models.basemodel:4647: DISPLAY = None
07:43:40 DEBUG   opendrift.models.basemodel:4648: Time to save animation: 0:02:31.202722
07:43:40 INFO    opendrift.models.basemodel:3037: Time to make animation: 0:02:51.763378
07:43:40 DEBUG   opendrift.models.basemodel:2365: Setting up map: corners=[4.0, 6, 59.5, 61], fast=False, lscale=None
07:43:40 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)
07:43:42 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:43 DEBUG   opendrift.models.basemodel:3044: Saving animation..
07:43:44 INFO    opendrift.models.basemodel:4608: Saving animation to /root/project/docs/source/gallery/animations/example_huge_output_1.gif...
07:43:44 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:45 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:47 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:48 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:50 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:51 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:53 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:54 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:55 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:57 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:43:58 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:00 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:01 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:02 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:04 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:05 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:07 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:08 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:10 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:12 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:13 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:15 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:17 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:18 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:20 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:21 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:22 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:24 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:25 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:26 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:28 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:29 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:31 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:32 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:34 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:35 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:36 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:38 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:39 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:41 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:42 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:44 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:45 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:47 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:48 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:50 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:51 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:53 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:54 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:56 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:57 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:44:59 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:01 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:02 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:04 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:06 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:08 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:09 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:11 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:12 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:14 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:16 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:18 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:19 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:21 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:23 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:24 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:26 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:27 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:29 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:30 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:32 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:34 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:35 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:37 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:38 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:40 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:41 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:43 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:44 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:45 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:47 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:49 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:50 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:52 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:53 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:54 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:56 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:57 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:45:59 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:46:00 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:46:01 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:46:03 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:46:04 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:46:06 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:46:07 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:46:09 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:46:10 DEBUG   opendrift.readers.reader_global_landmask:83: Adding GSHHG shapes from cartopy, scale: h, extent: (3.999999999999999, 6.0, 59.5, 60.99999999999999)..
07:46:12 DEBUG   opendrift.models.basemodel:4646: MPLBACKEND = agg
07:46:12 DEBUG   opendrift.models.basemodel:4647: DISPLAY = None
07:46:12 DEBUG   opendrift.models.basemodel:4648: Time to save animation: 0:02:28.882306
07:46:12 INFO    opendrift.models.basemodel:3037: Time to make animation: 0:02:32.289053

First seeding

../_images/example_huge_output_0.gif

Second seeding

../_images/example_huge_output_1.gif

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

Gallery generated by Sphinx-Gallery