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)
t1 = datetime.now()
t2 = t1 + timedelta(hours=6)
number = 10000
outfile = 'simulation.nc'  # Raw simulation output
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)

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.set_config('drift:horizontal_diffusivity', 10)
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)
13:59:44 DEBUG   opendrift.models.basemodel:625: Adding 18 config items from basemodel
13:59:44 DEBUG   opendrift.models.basemodel:625: Adding 5 config items from basemodel
13:59:44 DEBUG   opendrift.models.basemodel:625: Adding 48 config items from basemodel
13:59:44 INFO    opendrift.models.basemodel:539: OpenDriftSimulation initialised (version 1.10.7 / v1.10.6-119-g1da5bec)
13:59:44 DEBUG   opendrift.models.basemodel:625: Adding 15 config items from oceandrift
13:59:44 DEBUG   opendrift.models.basemodel:637:   Overwriting config item seed:z
13:59:44 DEBUG   opendrift.export.io_netcdf:216: Importing with Xarray from simulation.nc
13:59:44 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)
13:59:44 INFO    opendrift.models.basemodel:4667:       calculating for origin_marker 0...
13:59:45 INFO    opendrift.models.basemodel:4667:       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 2023-09-01 13:59 to 2023-09-02 13:59 UTC (49 steps)
13:59:45 DEBUG   opendrift.models.basemodel:3180: Setting up map: corners=None, fast=True, lscale=None
13:59:45 WARNING opendrift.models.basemodel:3226: Plotting fast. This will make your plots less accurate.
1.0009216427327419 ALPHA

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

First seeding

../_images/example_huge_output_0.gif

Second seeding

../_images/example_huge_output_1.gif

Total running time of the script: (5 minutes 0.762 seconds)

Gallery generated by Sphinx-Gallery