Note

Go to the end to download the full example code.

Radionuclides

from opendrift import test_data_folder as tdf
from opendrift.readers import reader_netCDF_CF_generic
from opendrift.models.radionuclides import RadionuclideDrift
from datetime import timedelta, datetime
import numpy as np


o = RadionuclideDrift(loglevel=20)  # Set loglevel to 0 for debug information

# Norkyst
#reader_norkyst = reader_netCDF_CF_generic.Reader(tdf + '/14Jan2016_NorKyst_z_3d/NorKyst-800m_ZDEPTHS_his_00_3Dsubset.nc')
reader_norkyst = reader_netCDF_CF_generic.Reader('https://thredds.met.no/thredds/dodsC/fou-hi/norkystv3_800m_m00_be')

o.add_reader([reader_norkyst])

# Adjusting some configuration
o.set_config('drift:vertical_mixing', True)
#o.set_config('vertical_mixing:diffusivitymodel','constant')  # include settling without vertical turbulent mixing
o.set_config('vertical_mixing:diffusivitymodel','environment')  # apply vertical diffusivity from ocean model
# Vertical mixing requires fast time step
o.set_config('vertical_mixing:timestep', 600.) # seconds
o.set_config('environment:constant:horizontal_diffusivity', 10)

16:11:51 INFO    opendrift:579: OpenDriftSimulation initialised (version 1.14.11 / v1.14.12-75-gf9e7122)
16:11:51 INFO    opendrift.readers:67: Opening file with xr.open_dataset
16:11:55 INFO    opendrift.readers.reader_netCDF_CF_generic:340: Detected dimensions: {'x': 'X', 'y': 'Y', 'z': 'depth', 'time': 'time'}
16:11:55 INFO    opendrift.readers.basereader:175: Variable x_sea_water_velocity will be rotated from eastward_sea_water_velocity
16:11:55 INFO    opendrift.readers.basereader:175: Variable y_sea_water_velocity will be rotated from northward_sea_water_velocity
16:11:55 INFO    opendrift.readers.basereader:175: Variable x_wind will be rotated from eastward_wind
16:11:55 INFO    opendrift.readers.basereader:175: Variable y_wind will be rotated from northward_wind

o.set_config('radionuclide:particle_diameter',5.e-6)  # m
o.set_config('radionuclide:sediment:resuspension_depth',2.)
o.set_config('radionuclide:sediment:resuspension_depth_uncert',0.1)
o.set_config('radionuclide:sediment:resuspension_critvel',0.15)
o.set_config('radionuclide:isotope', '137Cs')
o.set_config('radionuclide:specie_setup','LMM + Rev')
# By default, radionuclides do not strand towards coastline
o.set_config('general:coastline_action', 'previous')
o.set_config('general:seafloor_action','lift_to_seafloor')
o.set_config('seed:LMM_fraction',.45)
o.set_config('seed:particle_fraction',.55)


# SEEDING
td=datetime.today()
time = datetime(td.year, td.month, td.day, 0)

lat = 59.825673
lon = 10.724714
ntraj=5000
iniz=np.random.rand(ntraj) * -10. # seeding the radionuclides in the upper 10m
o.seed_elements(lon, lat, z=iniz, radius=10, number=ntraj, time=time)

16:11:55 INFO    opendrift.models.radionuclides:338: All good, isotop: 137Cs, species: LMM + Rev
16:11:55 INFO    opendrift.models.radionuclides:273: Number of species: 3
16:11:55 INFO    opendrift.models.radionuclides:275:   0 LMM
16:11:55 INFO    opendrift.models.radionuclides:275:   1 Particle reversible
16:11:55 INFO    opendrift.models.radionuclides:275:   2 Sediment reversible
ISOTOPE 137Cs
16:11:55 INFO    opendrift.models.radionuclides:523: ISOTOPE 137Cs
16:11:55 INFO    opendrift.models.radionuclides:524: SPECIE SETUP LMM + Rev
16:11:55 INFO    opendrift.models.radionuclides:525: Kd: 4.0
16:11:55 INFO    opendrift.models.radionuclides:418: Initial speciation:
16:11:55 INFO    opendrift.models.radionuclides:420:      2238   0 LMM
16:11:55 INFO    opendrift.models.radionuclides:420:      2762   1 Particle reversible
16:11:55 INFO    opendrift.models.radionuclides:420:         0   2 Sediment reversible
16:11:55 INFO    opendrift.models.radionuclides:440: PARTICLE DIAMETER 5e-06
16:11:55 INFO    opendrift.models.radionuclides:452: Min: 4.57e-06, Max: 5.39e-06, Numer of particles seeded: 2762, at normal distribution
16:11:55 INFO    opendrift.models.basemodel.environment:203: Adding a global landmask from GSHHG
16:12:01 INFO    opendrift.models.basemodel.environment:227: Fallback values will be used for the following variables which have no readers:
16:12:01 INFO    opendrift.models.basemodel.environment:230:    ocean_mixed_layer_thickness: 50.000000
16:12:01 INFO    opendrift.models.basemodel.environment:230:    conc3: 0.001000

Running model

o.run(duration=timedelta(hours=48), time_step=1800, time_step_output=3600)

16:12:01 INFO    opendrift:1932: Storing previous values of element property lon because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
16:12:01 INFO    opendrift:1932: Storing previous values of element property lat because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
16:12:01 INFO    opendrift:968: Using existing reader for land_binary_mask to move elements to ocean
16:12:01 INFO    opendrift:1000: All points are in ocean
16:12:01 INFO    opendrift.models.radionuclides:213: Number of species: 3
16:12:01 INFO    opendrift.models.radionuclides:215:   0 LMM
16:12:01 INFO    opendrift.models.radionuclides:215:   1 Particle reversible
16:12:01 INFO    opendrift.models.radionuclides:215:   2 Sediment reversible
16:12:01 INFO    opendrift.models.radionuclides:218: Isotope: 137Cs
16:12:01 INFO    opendrift.models.radionuclides:219: specie setup: LMM + Rev
16:12:01 INFO    opendrift.models.radionuclides:221: nspecies: 3
16:12:01 INFO    opendrift.models.radionuclides:222: Transfer rates:
 [[0.00000e+00 4.64000e-08 4.34304e-03]
 [1.16000e-05 0.00000e+00 0.00000e+00]
 [1.16000e-06 0.00000e+00 0.00000e+00]]
16:12:01 INFO    opendrift.models.radionuclides:226: DEPTH_INTERVALS: [-25.0, -10.0, -5.0, -1.0]
16:12:01 INFO    opendrift:2248: 2026-11-12 00:00:00 - step 1 of 96 - 5000 active elements (0 deactivated)
16:12:12 INFO    opendrift.models.radionuclides:779: Number of transformations: 68
16:12:12 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:12:12 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2306), np.int64(2694), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:12:12 INFO    opendrift:2248: 2026-11-12 00:30:00 - step 2 of 96 - 5000 active elements (0 deactivated)
16:12:24 INFO    opendrift.models.radionuclides:779: Number of transformations: 51
16:12:24 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:12:24 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2357), np.int64(2643), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:12:24 INFO    opendrift:2248: 2026-11-12 01:00:00 - step 3 of 96 - 5000 active elements (0 deactivated)
16:12:24 INFO    opendrift.models.radionuclides:779: Number of transformations: 39
16:12:25 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:12:25 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2396), np.int64(2604), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:12:25 INFO    opendrift:2248: 2026-11-12 01:30:00 - step 4 of 96 - 5000 active elements (0 deactivated)
16:12:37 INFO    opendrift.models.radionuclides:779: Number of transformations: 64
16:12:37 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:12:37 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2456), np.int64(2544), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:12:37 INFO    opendrift:2248: 2026-11-12 02:00:00 - step 5 of 96 - 5000 active elements (0 deactivated)
16:12:37 INFO    opendrift.models.radionuclides:779: Number of transformations: 55
16:12:37 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:12:37 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2511), np.int64(2489), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:12:37 INFO    opendrift:2248: 2026-11-12 02:30:00 - step 6 of 96 - 5000 active elements (0 deactivated)
16:12:48 INFO    opendrift.models.radionuclides:779: Number of transformations: 50
16:12:48 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:12:48 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2561), np.int64(2439), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:12:48 INFO    opendrift:2248: 2026-11-12 03:00:00 - step 7 of 96 - 5000 active elements (0 deactivated)
16:12:48 INFO    opendrift.models.radionuclides:779: Number of transformations: 53
16:12:48 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:12:48 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2614), np.int64(2386), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:12:48 INFO    opendrift:2248: 2026-11-12 03:30:00 - step 8 of 96 - 5000 active elements (0 deactivated)
16:12:59 INFO    opendrift.models.radionuclides:779: Number of transformations: 59
16:12:59 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:12:59 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2671), np.int64(2329), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:12:59 INFO    opendrift:2248: 2026-11-12 04:00:00 - step 9 of 96 - 5000 active elements (0 deactivated)
16:12:59 INFO    opendrift.models.radionuclides:779: Number of transformations: 47
16:12:59 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:12:59 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2718), np.int64(2282), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:12:59 INFO    opendrift:2248: 2026-11-12 04:30:00 - step 10 of 96 - 5000 active elements (0 deactivated)
16:13:12 INFO    opendrift.models.radionuclides:779: Number of transformations: 54
16:13:12 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:13:12 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2772), np.int64(2228), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:13:12 INFO    opendrift:2248: 2026-11-12 05:00:00 - step 11 of 96 - 5000 active elements (0 deactivated)
16:13:12 INFO    opendrift.models.radionuclides:779: Number of transformations: 49
16:13:12 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:13:12 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2821), np.int64(2179), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:13:12 INFO    opendrift:2248: 2026-11-12 05:30:00 - step 12 of 96 - 5000 active elements (0 deactivated)
16:13:26 INFO    opendrift.models.radionuclides:779: Number of transformations: 56
16:13:26 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:13:26 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2875), np.int64(2125), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:13:26 INFO    opendrift:2248: 2026-11-12 06:00:00 - step 13 of 96 - 5000 active elements (0 deactivated)
16:13:26 INFO    opendrift.models.radionuclides:779: Number of transformations: 46
16:13:26 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:13:26 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2921), np.int64(2079), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:13:26 INFO    opendrift:2248: 2026-11-12 06:30:00 - step 14 of 96 - 5000 active elements (0 deactivated)
16:13:37 INFO    opendrift.models.radionuclides:779: Number of transformations: 30
16:13:37 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:13:37 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2951), np.int64(2049), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:13:37 INFO    opendrift:2248: 2026-11-12 07:00:00 - step 15 of 96 - 5000 active elements (0 deactivated)
16:13:37 INFO    opendrift.models.radionuclides:779: Number of transformations: 50
16:13:37 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:13:37 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2999), np.int64(2000), np.int64(1)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:13:37 INFO    opendrift:2248: 2026-11-12 07:30:00 - step 16 of 96 - 5000 active elements (0 deactivated)
16:13:48 INFO    opendrift.models.radionuclides:779: Number of transformations: 39
16:13:48 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:13:48 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3036), np.int64(1962), np.int64(2)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:13:48 INFO    opendrift:2248: 2026-11-12 08:00:00 - step 17 of 96 - 5000 active elements (0 deactivated)
16:13:48 INFO    opendrift.models.radionuclides:779: Number of transformations: 37
16:13:48 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:13:48 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3073), np.int64(1925), np.int64(2)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:13:48 INFO    opendrift:2248: 2026-11-12 08:30:00 - step 18 of 96 - 5000 active elements (0 deactivated)
16:14:01 INFO    opendrift.models.radionuclides:779: Number of transformations: 42
16:14:01 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:14:01 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3113), np.int64(1884), np.int64(3)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:14:01 INFO    opendrift:2248: 2026-11-12 09:00:00 - step 19 of 96 - 5000 active elements (0 deactivated)
16:14:01 INFO    opendrift.models.radionuclides:779: Number of transformations: 48
16:14:01 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:14:01 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3159), np.int64(1838), np.int64(3)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:14:01 INFO    opendrift:2248: 2026-11-12 09:30:00 - step 20 of 96 - 5000 active elements (0 deactivated)
16:14:12 INFO    opendrift.models.radionuclides:779: Number of transformations: 40
16:14:12 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:14:12 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3199), np.int64(1798), np.int64(3)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:14:12 INFO    opendrift:2248: 2026-11-12 10:00:00 - step 21 of 96 - 5000 active elements (0 deactivated)
16:14:12 INFO    opendrift.models.radionuclides:779: Number of transformations: 45
16:14:12 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:14:12 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3242), np.int64(1754), np.int64(4)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:14:12 INFO    opendrift:2248: 2026-11-12 10:30:00 - step 22 of 96 - 5000 active elements (0 deactivated)
16:14:24 INFO    opendrift.models.radionuclides:779: Number of transformations: 40
16:14:24 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:14:24 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3282), np.int64(1714), np.int64(4)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:14:24 INFO    opendrift:2248: 2026-11-12 11:00:00 - step 23 of 96 - 5000 active elements (0 deactivated)
16:14:24 INFO    opendrift.models.radionuclides:779: Number of transformations: 38
16:14:24 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:14:24 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3320), np.int64(1676), np.int64(4)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:14:24 INFO    opendrift:2248: 2026-11-12 11:30:00 - step 24 of 96 - 5000 active elements (0 deactivated)
16:14:35 INFO    opendrift.models.radionuclides:779: Number of transformations: 22
16:14:35 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:14:35 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3340), np.int64(1655), np.int64(5)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:14:35 INFO    opendrift:2248: 2026-11-12 12:00:00 - step 25 of 96 - 5000 active elements (0 deactivated)
16:14:35 INFO    opendrift.models.radionuclides:779: Number of transformations: 43
16:14:35 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:14:35 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3383), np.int64(1612), np.int64(5)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:14:35 INFO    opendrift:2248: 2026-11-12 12:30:00 - step 26 of 96 - 5000 active elements (0 deactivated)
16:14:47 INFO    opendrift.models.radionuclides:779: Number of transformations: 32
16:14:47 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:14:47 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3411), np.int64(1583), np.int64(6)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:14:47 INFO    opendrift:2248: 2026-11-12 13:00:00 - step 27 of 96 - 5000 active elements (0 deactivated)
16:14:47 INFO    opendrift.models.radionuclides:779: Number of transformations: 27
16:14:47 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:14:47 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3438), np.int64(1556), np.int64(6)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:14:47 INFO    opendrift:2248: 2026-11-12 13:30:00 - step 28 of 96 - 5000 active elements (0 deactivated)
16:15:00 INFO    opendrift.models.radionuclides:779: Number of transformations: 39
16:15:00 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:15:00 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3475), np.int64(1519), np.int64(6)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:15:00 INFO    opendrift:2248: 2026-11-12 14:00:00 - step 29 of 96 - 5000 active elements (0 deactivated)
16:15:00 INFO    opendrift.models.radionuclides:779: Number of transformations: 29
16:15:00 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:15:00 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3504), np.int64(1490), np.int64(6)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:15:00 INFO    opendrift:2248: 2026-11-12 14:30:00 - step 30 of 96 - 5000 active elements (0 deactivated)
16:15:11 INFO    opendrift.models.radionuclides:779: Number of transformations: 34
16:15:12 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:15:12 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3538), np.int64(1456), np.int64(6)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:15:12 INFO    opendrift:2248: 2026-11-12 15:00:00 - step 31 of 96 - 5000 active elements (0 deactivated)
16:15:12 INFO    opendrift.models.radionuclides:779: Number of transformations: 29
16:15:12 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:15:12 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3567), np.int64(1428), np.int64(5)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:15:12 INFO    opendrift:2248: 2026-11-12 15:30:00 - step 32 of 96 - 5000 active elements (0 deactivated)
16:15:23 INFO    opendrift.models.radionuclides:779: Number of transformations: 28
16:15:23 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:15:23 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3591), np.int64(1403), np.int64(6)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:15:23 INFO    opendrift:2248: 2026-11-12 16:00:00 - step 33 of 96 - 5000 active elements (0 deactivated)
16:15:24 INFO    opendrift.models.radionuclides:779: Number of transformations: 28
16:15:24 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:15:24 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3619), np.int64(1375), np.int64(6)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:15:24 INFO    opendrift:2248: 2026-11-12 16:30:00 - step 34 of 96 - 5000 active elements (0 deactivated)
16:16:51 INFO    opendrift.models.radionuclides:779: Number of transformations: 25
16:16:51 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:16:51 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3642), np.int64(1352), np.int64(6)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:16:51 INFO    opendrift:2248: 2026-11-12 17:00:00 - step 35 of 96 - 5000 active elements (0 deactivated)
16:16:52 INFO    opendrift.models.radionuclides:779: Number of transformations: 36
16:16:52 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:16:52 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3678), np.int64(1316), np.int64(6)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:16:52 INFO    opendrift:2248: 2026-11-12 17:30:00 - step 36 of 96 - 5000 active elements (0 deactivated)
16:17:03 INFO    opendrift.models.radionuclides:779: Number of transformations: 25
16:17:03 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:17:03 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3703), np.int64(1290), np.int64(7)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:17:03 INFO    opendrift:2248: 2026-11-12 18:00:00 - step 37 of 96 - 5000 active elements (0 deactivated)
16:17:03 INFO    opendrift.models.radionuclides:779: Number of transformations: 17
16:17:03 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:17:03 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3720), np.int64(1273), np.int64(7)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:17:03 INFO    opendrift:2248: 2026-11-12 18:30:00 - step 38 of 96 - 5000 active elements (0 deactivated)
16:17:15 INFO    opendrift.models.radionuclides:779: Number of transformations: 16
16:17:15 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:17:15 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3734), np.int64(1258), np.int64(8)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:17:15 INFO    opendrift:2248: 2026-11-12 19:00:00 - step 39 of 96 - 5000 active elements (0 deactivated)
16:17:15 INFO    opendrift.models.radionuclides:779: Number of transformations: 29
16:17:15 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:17:15 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3761), np.int64(1231), np.int64(8)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:17:15 INFO    opendrift:2248: 2026-11-12 19:30:00 - step 40 of 96 - 5000 active elements (0 deactivated)
16:17:27 INFO    opendrift.models.radionuclides:779: Number of transformations: 25
16:17:27 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:17:27 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3786), np.int64(1206), np.int64(8)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:17:27 INFO    opendrift:2248: 2026-11-12 20:00:00 - step 41 of 96 - 5000 active elements (0 deactivated)
16:17:27 INFO    opendrift.models.radionuclides:779: Number of transformations: 25
16:17:27 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:17:27 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3809), np.int64(1183), np.int64(8)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:17:28 INFO    opendrift:2248: 2026-11-12 20:30:00 - step 42 of 96 - 5000 active elements (0 deactivated)
16:17:40 INFO    opendrift.models.radionuclides:779: Number of transformations: 23
16:17:40 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:17:40 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3830), np.int64(1162), np.int64(8)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:17:40 INFO    opendrift:2248: 2026-11-12 21:00:00 - step 43 of 96 - 5000 active elements (0 deactivated)
16:17:40 INFO    opendrift.models.radionuclides:779: Number of transformations: 32
16:17:40 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:17:40 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3858), np.int64(1133), np.int64(9)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:17:40 INFO    opendrift:2248: 2026-11-12 21:30:00 - step 44 of 96 - 5000 active elements (0 deactivated)
16:17:53 INFO    opendrift.models.radionuclides:779: Number of transformations: 22
16:17:53 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:17:53 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3878), np.int64(1112), np.int64(10)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:17:53 INFO    opendrift:2248: 2026-11-12 22:00:00 - step 45 of 96 - 5000 active elements (0 deactivated)
16:17:53 INFO    opendrift.models.radionuclides:779: Number of transformations: 29
16:17:53 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:17:53 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3903), np.int64(1085), np.int64(12)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:17:53 INFO    opendrift:2248: 2026-11-12 22:30:00 - step 46 of 96 - 5000 active elements (0 deactivated)
16:18:05 INFO    opendrift.models.radionuclides:779: Number of transformations: 33
16:18:05 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:18:05 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3934), np.int64(1054), np.int64(12)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:18:05 INFO    opendrift:2248: 2026-11-12 23:00:00 - step 47 of 96 - 5000 active elements (0 deactivated)
16:18:05 INFO    opendrift.models.radionuclides:779: Number of transformations: 26
16:18:05 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:18:05 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3958), np.int64(1028), np.int64(14)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:18:05 INFO    opendrift:2248: 2026-11-12 23:30:00 - step 48 of 96 - 5000 active elements (0 deactivated)
16:18:18 INFO    opendrift.models.radionuclides:779: Number of transformations: 19
16:18:18 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:18:18 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3975), np.int64(1011), np.int64(14)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:18:18 INFO    opendrift:2248: 2026-11-13 00:00:00 - step 49 of 96 - 5000 active elements (0 deactivated)
16:18:18 INFO    opendrift.models.radionuclides:779: Number of transformations: 20
16:18:18 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:18:18 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3991), np.int64(993), np.int64(16)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:18:18 INFO    opendrift:2248: 2026-11-13 00:30:00 - step 50 of 96 - 5000 active elements (0 deactivated)
16:18:30 INFO    opendrift.models.radionuclides:779: Number of transformations: 20
16:18:30 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:18:30 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4009), np.int64(974), np.int64(17)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:18:30 INFO    opendrift:2248: 2026-11-13 01:00:00 - step 51 of 96 - 5000 active elements (0 deactivated)
16:18:31 INFO    opendrift.models.radionuclides:779: Number of transformations: 10
16:18:31 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:18:31 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4019), np.int64(964), np.int64(17)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:18:31 INFO    opendrift:2248: 2026-11-13 01:30:00 - step 52 of 96 - 5000 active elements (0 deactivated)
16:18:42 INFO    opendrift.models.radionuclides:779: Number of transformations: 26
16:18:42 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:18:42 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4045), np.int64(938), np.int64(17)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:18:42 INFO    opendrift:2248: 2026-11-13 02:00:00 - step 53 of 96 - 5000 active elements (0 deactivated)
16:18:42 INFO    opendrift.models.radionuclides:779: Number of transformations: 24
16:18:42 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:18:42 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4067), np.int64(916), np.int64(17)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:18:42 INFO    opendrift:2248: 2026-11-13 02:30:00 - step 54 of 96 - 5000 active elements (0 deactivated)
16:18:54 INFO    opendrift.models.radionuclides:779: Number of transformations: 12
16:18:54 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:18:54 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4079), np.int64(904), np.int64(17)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:18:54 INFO    opendrift:2248: 2026-11-13 03:00:00 - step 55 of 96 - 5000 active elements (0 deactivated)
16:18:54 INFO    opendrift.models.radionuclides:779: Number of transformations: 13
16:18:54 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:18:54 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4090), np.int64(892), np.int64(18)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:18:54 INFO    opendrift:2248: 2026-11-13 03:30:00 - step 56 of 96 - 5000 active elements (0 deactivated)
16:19:04 INFO    opendrift.models.radionuclides:779: Number of transformations: 24
16:19:04 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:19:04 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4112), np.int64(869), np.int64(19)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:19:04 INFO    opendrift:2248: 2026-11-13 04:00:00 - step 57 of 96 - 5000 active elements (0 deactivated)
16:19:04 INFO    opendrift.models.radionuclides:779: Number of transformations: 16
16:19:04 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:19:04 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4126), np.int64(854), np.int64(20)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:19:04 INFO    opendrift:2248: 2026-11-13 04:30:00 - step 58 of 96 - 5000 active elements (0 deactivated)
16:19:16 INFO    opendrift.models.radionuclides:779: Number of transformations: 26
16:19:16 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:19:16 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4144), np.int64(832), np.int64(24)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:19:16 INFO    opendrift:2248: 2026-11-13 05:00:00 - step 59 of 96 - 5000 active elements (0 deactivated)
16:19:16 INFO    opendrift.models.radionuclides:779: Number of transformations: 14
16:19:16 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:19:16 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4158), np.int64(818), np.int64(24)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:19:16 INFO    opendrift:2248: 2026-11-13 05:30:00 - step 60 of 96 - 5000 active elements (0 deactivated)
16:19:28 INFO    opendrift.models.radionuclides:779: Number of transformations: 19
16:19:28 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:19:28 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4175), np.int64(800), np.int64(25)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:19:28 INFO    opendrift:2248: 2026-11-13 06:00:00 - step 61 of 96 - 5000 active elements (0 deactivated)
16:19:28 INFO    opendrift.models.radionuclides:779: Number of transformations: 19
16:19:28 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:19:28 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4188), np.int64(784), np.int64(28)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:19:28 INFO    opendrift:2248: 2026-11-13 06:30:00 - step 62 of 96 - 5000 active elements (0 deactivated)
16:19:40 INFO    opendrift.models.radionuclides:779: Number of transformations: 11
16:19:40 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:19:40 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4197), np.int64(774), np.int64(29)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:19:40 INFO    opendrift:2248: 2026-11-13 07:00:00 - step 63 of 96 - 5000 active elements (0 deactivated)
16:19:40 INFO    opendrift.models.radionuclides:779: Number of transformations: 18
16:19:40 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:19:40 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4213), np.int64(758), np.int64(29)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:19:40 INFO    opendrift:2248: 2026-11-13 07:30:00 - step 64 of 96 - 5000 active elements (0 deactivated)
16:19:52 INFO    opendrift.models.radionuclides:779: Number of transformations: 11
16:19:52 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:19:52 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4224), np.int64(747), np.int64(29)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:19:52 INFO    opendrift:2248: 2026-11-13 08:00:00 - step 65 of 96 - 5000 active elements (0 deactivated)
16:19:52 INFO    opendrift.models.radionuclides:779: Number of transformations: 17
16:19:52 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:19:52 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4235), np.int64(733), np.int64(32)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:19:52 INFO    opendrift:2248: 2026-11-13 08:30:00 - step 66 of 96 - 5000 active elements (0 deactivated)
16:20:03 INFO    opendrift.models.radionuclides:779: Number of transformations: 22
16:20:03 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:20:03 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4251), np.int64(714), np.int64(35)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:20:03 INFO    opendrift:2248: 2026-11-13 09:00:00 - step 67 of 96 - 5000 active elements (0 deactivated)
16:20:04 INFO    opendrift.models.radionuclides:779: Number of transformations: 21
16:20:04 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:20:04 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4268), np.int64(697), np.int64(35)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:20:04 INFO    opendrift:2248: 2026-11-13 09:30:00 - step 68 of 96 - 5000 active elements (0 deactivated)
16:20:15 INFO    opendrift.models.radionuclides:779: Number of transformations: 13
16:20:15 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:20:15 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4277), np.int64(686), np.int64(37)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:20:15 INFO    opendrift:2248: 2026-11-13 10:00:00 - step 69 of 96 - 5000 active elements (0 deactivated)
16:20:15 INFO    opendrift.models.radionuclides:779: Number of transformations: 15
16:20:15 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:20:15 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4288), np.int64(673), np.int64(39)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:20:15 INFO    opendrift:2248: 2026-11-13 10:30:00 - step 70 of 96 - 5000 active elements (0 deactivated)
16:20:29 INFO    opendrift.models.radionuclides:779: Number of transformations: 12
16:20:29 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:20:29 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4298), np.int64(662), np.int64(40)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:20:29 INFO    opendrift:2248: 2026-11-13 11:00:00 - step 71 of 96 - 5000 active elements (0 deactivated)
16:20:29 INFO    opendrift.models.radionuclides:779: Number of transformations: 20
16:20:29 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:20:29 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4308), np.int64(649), np.int64(43)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:20:29 INFO    opendrift:2248: 2026-11-13 11:30:00 - step 72 of 96 - 5000 active elements (0 deactivated)
16:20:40 INFO    opendrift.models.radionuclides:779: Number of transformations: 17
16:20:40 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:20:40 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4319), np.int64(635), np.int64(46)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:20:40 INFO    opendrift:2248: 2026-11-13 12:00:00 - step 73 of 96 - 5000 active elements (0 deactivated)
16:20:40 INFO    opendrift.models.radionuclides:779: Number of transformations: 19
16:20:40 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:20:40 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4336), np.int64(618), np.int64(46)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:20:40 INFO    opendrift:2248: 2026-11-13 12:30:00 - step 74 of 96 - 5000 active elements (0 deactivated)
16:20:52 INFO    opendrift.models.radionuclides:779: Number of transformations: 22
16:20:52 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:20:52 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4350), np.int64(600), np.int64(50)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:20:52 INFO    opendrift:2248: 2026-11-13 13:00:00 - step 75 of 96 - 5000 active elements (0 deactivated)
16:20:52 INFO    opendrift.models.radionuclides:779: Number of transformations: 12
16:20:52 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:20:52 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4358), np.int64(590), np.int64(52)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:20:52 INFO    opendrift:2248: 2026-11-13 13:30:00 - step 76 of 96 - 5000 active elements (0 deactivated)
16:21:03 INFO    opendrift.models.radionuclides:779: Number of transformations: 18
16:21:03 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:21:03 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4372), np.int64(575), np.int64(53)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:21:03 INFO    opendrift:2248: 2026-11-13 14:00:00 - step 77 of 96 - 5000 active elements (0 deactivated)
16:21:03 INFO    opendrift.models.radionuclides:779: Number of transformations: 24
16:21:03 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:21:03 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4388), np.int64(556), np.int64(56)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:21:03 INFO    opendrift:2248: 2026-11-13 14:30:00 - step 78 of 96 - 5000 active elements (0 deactivated)
16:21:14 INFO    opendrift.models.radionuclides:779: Number of transformations: 6
16:21:14 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:21:14 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4390), np.int64(553), np.int64(57)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:21:14 INFO    opendrift:2248: 2026-11-13 15:00:00 - step 79 of 96 - 5000 active elements (0 deactivated)
16:21:14 INFO    opendrift.models.radionuclides:779: Number of transformations: 11
16:21:14 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:21:14 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4391), np.int64(547), np.int64(62)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:21:14 INFO    opendrift:2248: 2026-11-13 15:30:00 - step 80 of 96 - 5000 active elements (0 deactivated)
16:21:24 INFO    opendrift.models.radionuclides:779: Number of transformations: 12
16:21:24 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:21:24 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4393), np.int64(540), np.int64(67)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:21:24 INFO    opendrift:2248: 2026-11-13 16:00:00 - step 81 of 96 - 5000 active elements (0 deactivated)
16:21:25 INFO    opendrift.models.radionuclides:779: Number of transformations: 5
16:21:25 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:21:25 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4396), np.int64(537), np.int64(67)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:21:25 INFO    opendrift:2248: 2026-11-13 16:30:00 - step 82 of 96 - 5000 active elements (0 deactivated)
16:21:34 INFO    opendrift.models.radionuclides:779: Number of transformations: 16
16:21:34 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:21:34 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4404), np.int64(525), np.int64(71)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:21:34 INFO    opendrift:2248: 2026-11-13 17:00:00 - step 83 of 96 - 5000 active elements (0 deactivated)
16:21:35 INFO    opendrift.models.radionuclides:779: Number of transformations: 11
16:21:35 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:21:35 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4409), np.int64(517), np.int64(74)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:21:35 INFO    opendrift:2248: 2026-11-13 17:30:00 - step 84 of 96 - 5000 active elements (0 deactivated)
16:21:45 INFO    opendrift.models.radionuclides:779: Number of transformations: 11
16:21:45 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:21:45 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4414), np.int64(509), np.int64(77)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:21:45 INFO    opendrift:2248: 2026-11-13 18:00:00 - step 85 of 96 - 5000 active elements (0 deactivated)
16:21:45 INFO    opendrift.models.radionuclides:779: Number of transformations: 12
16:21:45 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:21:45 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4424), np.int64(497), np.int64(79)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:21:45 INFO    opendrift:2248: 2026-11-13 18:30:00 - step 86 of 96 - 5000 active elements (0 deactivated)
16:21:56 INFO    opendrift.models.radionuclides:779: Number of transformations: 8
16:21:56 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:21:56 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4428), np.int64(491), np.int64(81)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:21:56 INFO    opendrift:2248: 2026-11-13 19:00:00 - step 87 of 96 - 5000 active elements (0 deactivated)
16:21:56 INFO    opendrift.models.radionuclides:779: Number of transformations: 7
16:21:56 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:21:56 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4435), np.int64(484), np.int64(81)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:21:56 INFO    opendrift:2248: 2026-11-13 19:30:00 - step 88 of 96 - 5000 active elements (0 deactivated)
16:22:06 INFO    opendrift.models.radionuclides:779: Number of transformations: 12
16:22:06 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:22:06 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4441), np.int64(476), np.int64(83)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:22:06 INFO    opendrift:2248: 2026-11-13 20:00:00 - step 89 of 96 - 5000 active elements (0 deactivated)
16:22:06 INFO    opendrift.models.radionuclides:779: Number of transformations: 6
16:22:06 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:22:06 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4445), np.int64(472), np.int64(83)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:22:06 INFO    opendrift:2248: 2026-11-13 20:30:00 - step 90 of 96 - 5000 active elements (0 deactivated)
16:22:14 INFO    opendrift.models.radionuclides:779: Number of transformations: 12
16:22:14 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:22:14 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4455), np.int64(461), np.int64(84)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:22:14 INFO    opendrift:2248: 2026-11-13 21:00:00 - step 91 of 96 - 5000 active elements (0 deactivated)
16:22:15 INFO    opendrift.models.radionuclides:779: Number of transformations: 15
16:22:15 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:22:15 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4466), np.int64(448), np.int64(86)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:22:15 INFO    opendrift:2248: 2026-11-13 21:30:00 - step 92 of 96 - 5000 active elements (0 deactivated)
16:22:26 INFO    opendrift.models.radionuclides:779: Number of transformations: 16
16:22:26 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:22:26 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4468), np.int64(439), np.int64(93)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:22:26 INFO    opendrift:2248: 2026-11-13 22:00:00 - step 93 of 96 - 5000 active elements (0 deactivated)
16:22:26 INFO    opendrift.models.radionuclides:779: Number of transformations: 13
16:22:26 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:22:26 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4477), np.int64(429), np.int64(94)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:22:26 INFO    opendrift:2248: 2026-11-13 22:30:00 - step 94 of 96 - 5000 active elements (0 deactivated)
16:22:37 INFO    opendrift.models.radionuclides:779: Number of transformations: 13
16:22:37 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:22:37 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4486), np.int64(420), np.int64(94)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:22:37 INFO    opendrift:2248: 2026-11-13 23:00:00 - step 95 of 96 - 5000 active elements (0 deactivated)
16:22:37 INFO    opendrift.models.radionuclides:779: Number of transformations: 19
16:22:37 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:22:37 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4489), np.int64(411), np.int64(100)] ['LMM', 'Particle reversible', 'Sediment reversible']
16:22:37 INFO    opendrift:2248: 2026-11-13 23:30:00 - step 96 of 96 - 5000 active elements (0 deactivated)
16:22:47 INFO    opendrift.models.radionuclides:779: Number of transformations: 12
16:22:47 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
16:22:47 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4489), np.int64(406), np.int64(105)] ['LMM', 'Particle reversible', 'Sediment reversible']

<xarray.Dataset> Size: 27MB
Dimensions:                          (trajectory: 5000, time: 49)
Coordinates:
  * trajectory                       (trajectory) int64 40kB 0 1 2 ... 4998 4999
  * time                             (time) datetime64[us] 392B 2026-11-12 .....
Data variables: (12/28)
    status                           (trajectory, time) float32 980kB 0.0 ......
    moving                           (trajectory, time) float32 980kB 1.0 ......
    age_seconds                      (trajectory, time) float32 980kB 0.0 ......
    origin_marker                    (trajectory, time) float32 980kB 0.0 ......
    lon                              (trajectory, time) float32 980kB 10.72 ....
    lat                              (trajectory, time) float32 980kB 59.83 ....
    ...                               ...
    ocean_mixed_layer_thickness      (trajectory, time) float32 980kB 50.0 .....
    sea_water_temperature            (trajectory, time) float32 980kB 8.91 .....
    sea_water_salinity               (trajectory, time) float32 980kB 29.49 ....
    horizontal_diffusivity           (trajectory, time) float32 980kB 10.0 .....
    upward_sea_water_velocity        (trajectory, time) float32 980kB -3.271e...
    conc3                            (trajectory, time) float32 980kB 0.001 ....
Attributes: (12/155)
    Conventions:                                                  CF-1.11, AC...
    standard_name_vocabulary:                                     CF Standard...
    featureType:                                                  trajectory
    title:                                                        OpenDrift t...
    summary:                                                      Output from...
    keywords:                                                     trajectory,...
    ...                                                           ...
    geospatial_lon_units:                                         degrees_east
    geospatial_lon_resolution:                                    point
    runtime:                                                      0:10:56.615468
    geospatial_vertical_min:                                      -46.81189
    geospatial_vertical_max:                                      0.0
    geospatial_vertical_positive:                                 up


Print and plot results

print(o)
print('Final speciation:')
for isp,sp in enumerate(o.name_species):
    print ('{:32}: {:>6}'.format(sp,sum(o.elements.specie==isp)))

print('Number of transformations:')
for isp in range(o.nspecies):
    print('{}'.format(['{:>9}'.format(np.int32(item)) for item in o.ntransformations[isp,:]]) )

o.animation(color='specie',
            vmin=0,vmax=o.nspecies-1,
            colorbar=False,
            legend=[o.specie_num2name(i) for i in range(o.nspecies)],
            )

===========================
--------------------
Reader performance:
--------------------
https://thredds.met.no/thredds/dodsC/fou-hi/norkystv3_800m_m00_be
 0:10:36.5  total
 0:00:00.0  preparing
 0:10:34.9  reading
 0:00:05.2  interpolation
 0:00:00.0  interpolation_time
 0:00:01.5  rotating vectors
 0:00:00.0  masking
--------------------
constant_reader
 0:00:00.0  total
 0:00:00.0  preparing
 0:00:00.0  reading
 0:00:00.0  masking
--------------------
global_landmask
 0:00:00.9  total
 0:00:00.0  preparing
 0:00:00.8  reading
 0:00:00.0  masking
--------------------
Performance:
10:56.6 total time
   10.2 configuration
    0.1 preparing main loop
      0.0 moving elements to ocean
10:46.2 main loop
      1.5 updating elements
        0.5 vertical mixing
    0.0 cleaning up
--------------------
===========================
Model:  RadionuclideDrift     (OpenDrift version 1.14.11)
        5000 active Radionuclide particles  (0 deactivated, 0 scheduled)
-------------------
Environment variables:
  -----
  ocean_vertical_diffusivity
  sea_floor_depth_below_sea_level
  sea_surface_height
  sea_water_salinity
  sea_water_temperature
  upward_sea_water_velocity
  x_sea_water_velocity
  x_wind
  y_sea_water_velocity
  y_wind
     1) https://thredds.met.no/thredds/dodsC/fou-hi/norkystv3_800m_m00_be
  -----
  horizontal_diffusivity
     1) constant_reader
  -----
  land_binary_mask
     1) global_landmask
  -----
Readers not added for the following variables:
  conc3
  ocean_mixed_layer_thickness

Discarded readers:

Time:
        Start: 2026-11-12 00:00:00 UTC
        Present: 2026-11-14 00:00:00 UTC
        Calculation steps: 96 * 0:30:00 - total time: 2 days, 0:00:00
        Output steps: 49 * 1:00:00
===========================

Final speciation:
LMM                             :   4489
Particle reversible             :    406
Sediment reversible             :    105
Number of transformations:
['        0', '       28', '      109']
['     2381', '        0', '        3']
['        7', '        0', '        0']
16:22:56 INFO    opendrift:4899: Saving animation to /root/project/docs/source/gallery/animations/example_radionuclides_0.gif...
16:23:11 INFO    opendrift:3307: Time to make animation: 0:00:23.367949
../_images/example_radionuclides_0.gif 
o.plot_vertical_distribution()
#o.plot_property('specie')
o.animation_profile(color='specie',
            vmin=0,vmax=o.nspecies-1,
            legend=[o.specie_num2name(i) for i in range(o.nspecies)],
            legend_loc =3,
#            markersize=10
            )
2026-11-14T00:00:00 Mean windspeed: 1.0 m/s
/root/project/opendrift/models/oceandrift.py:664: UserWarning: This figure includes Axes that are not compatible with tight_layout, so results might be incorrect.
  tslider = Slider(sliderax, 'Timestep', 0, steps_output-1,
/root/project/opendrift/models/oceandrift.py:706: UserWarning: This figure includes Axes that are not compatible with tight_layout, so results might be incorrect.
  fig.canvas.draw_idle()
16:23:11 INFO    opendrift:4899: Saving animation to /root/project/docs/source/gallery/animations/example_radionuclides_1.gif...
16:23:19 INFO    opendrift:3525: Time to make animation: 0:00:08.232413
../_images/example_radionuclides_1.gif 
# # Postprocessing: write to concentration netcdf file
o.write_netcdf_radionuclide_density_map(‘radio_conc.nc’, pixelsize_m=500.,

zlevels=[-2.], activity_unit=’Bq’, horizontal_smoothing=True, smoothing_cells=1, time_avg_conc=True, deltat=2., # hours llcrnrlon=4.4, llcrnrlat=59.9, urcrnrlon=4.8, urcrnrlat=60.2,

)

Total running time of the script: (11 minutes 36.231 seconds)

Gallery generated by Sphinx-Gallery