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)
15:35:44 INFO    opendrift:569: OpenDriftSimulation initialised (version 1.14.11 / v1.14.12-32-g22db1c0)
15:35:44 INFO    opendrift.readers:67: Opening file with xr.open_dataset
15:35:45 INFO    opendrift.readers.reader_netCDF_CF_generic:340: Detected dimensions: {'x': 'X', 'y': 'Y', 'z': 'depth', 'time': 'time'}
15:35:45 INFO    opendrift.readers.basereader:175: Variable x_sea_water_velocity will be rotated from eastward_sea_water_velocity
15:35:45 INFO    opendrift.readers.basereader:175: Variable y_sea_water_velocity will be rotated from northward_sea_water_velocity
15:35:45 INFO    opendrift.readers.basereader:175: Variable x_wind will be rotated from eastward_wind
15:35:45 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)
15:35:45 INFO    opendrift.models.radionuclides:338: All good, isotop: 137Cs, species: LMM + Rev
15:35:45 INFO    opendrift.models.radionuclides:273: Number of species: 3
15:35:45 INFO    opendrift.models.radionuclides:275:   0 LMM
15:35:45 INFO    opendrift.models.radionuclides:275:   1 Particle reversible
15:35:45 INFO    opendrift.models.radionuclides:275:   2 Sediment reversible
ISOTOPE 137Cs
15:35:45 INFO    opendrift.models.radionuclides:523: ISOTOPE 137Cs
15:35:45 INFO    opendrift.models.radionuclides:524: SPECIE SETUP LMM + Rev
15:35:45 INFO    opendrift.models.radionuclides:525: Kd: 4.0
15:35:45 INFO    opendrift.models.radionuclides:418: Initial speciation:
15:35:45 INFO    opendrift.models.radionuclides:420:      2238   0 LMM
15:35:45 INFO    opendrift.models.radionuclides:420:      2762   1 Particle reversible
15:35:45 INFO    opendrift.models.radionuclides:420:         0   2 Sediment reversible
15:35:45 INFO    opendrift.models.radionuclides:440: PARTICLE DIAMETER 5e-06
15:35:45 INFO    opendrift.models.radionuclides:452: Min: 4.63e-06, Max: 5.31e-06, Numer of particles seeded: 2762, at normal distribution
15:35:45 INFO    opendrift.models.basemodel.environment:203: Adding a global landmask from GSHHG
15:35:52 INFO    opendrift.models.basemodel.environment:227: Fallback values will be used for the following variables which have no readers:
15:35:52 INFO    opendrift.models.basemodel.environment:230:    ocean_mixed_layer_thickness: 50.000000
15:35:52 INFO    opendrift.models.basemodel.environment:230:    conc3: 0.001000

Running model

o.run(duration=timedelta(hours=48), time_step=1800, time_step_output=3600)
15:35:52 INFO    opendrift:1921: Storing previous values of element property lon because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
15:35:52 INFO    opendrift:1921: Storing previous values of element property lat because of condition (('general:coastline_action', 'in', ['stranding', 'previous']), 'or', ('general:seafloor_action', 'in', ['previous']))
15:35:52 INFO    opendrift:958: Using existing reader for land_binary_mask to move elements to ocean
15:35:52 INFO    opendrift:989: All points are in ocean
15:35:52 INFO    opendrift.models.radionuclides:213: Number of species: 3
15:35:52 INFO    opendrift.models.radionuclides:215:   0 LMM
15:35:52 INFO    opendrift.models.radionuclides:215:   1 Particle reversible
15:35:52 INFO    opendrift.models.radionuclides:215:   2 Sediment reversible
15:35:52 INFO    opendrift.models.radionuclides:218: Isotope: 137Cs
15:35:52 INFO    opendrift.models.radionuclides:219: specie setup: LMM + Rev
15:35:52 INFO    opendrift.models.radionuclides:221: nspecies: 3
15:35:52 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]]
15:35:52 INFO    opendrift.models.radionuclides:226: DEPTH_INTERVALS: [-25.0, -10.0, -5.0, -1.0]
15:35:52 INFO    opendrift:2237: 2026-10-05 00:00:00 - step 1 of 96 - 5000 active elements (0 deactivated)
15:35:54 INFO    opendrift.models.radionuclides:779: Number of transformations: 68
15:35:54 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:35:54 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2306), np.int64(2694), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:35:54 INFO    opendrift:2237: 2026-10-05 00:30:00 - step 2 of 96 - 5000 active elements (0 deactivated)
15:35:59 INFO    opendrift.models.radionuclides:779: Number of transformations: 51
15:35:59 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:35:59 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2357), np.int64(2643), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:35:59 INFO    opendrift:2237: 2026-10-05 01:00:00 - step 3 of 96 - 5000 active elements (0 deactivated)
15:35:59 INFO    opendrift.models.radionuclides:779: Number of transformations: 39
15:35:59 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:35:59 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2396), np.int64(2604), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:35:59 INFO    opendrift:2237: 2026-10-05 01:30:00 - step 4 of 96 - 5000 active elements (0 deactivated)
15:36:02 INFO    opendrift.models.radionuclides:779: Number of transformations: 64
15:36:02 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:02 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2456), np.int64(2544), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:02 INFO    opendrift:2237: 2026-10-05 02:00:00 - step 5 of 96 - 5000 active elements (0 deactivated)
15:36:02 INFO    opendrift.models.radionuclides:779: Number of transformations: 55
15:36:02 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:02 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2511), np.int64(2489), np.int64(0)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:02 INFO    opendrift:2237: 2026-10-05 02:30:00 - step 6 of 96 - 5000 active elements (0 deactivated)
15:36:05 INFO    opendrift.models.radionuclides:779: Number of transformations: 51
15:36:05 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:05 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2560), np.int64(2439), np.int64(1)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:05 INFO    opendrift:2237: 2026-10-05 03:00:00 - step 7 of 96 - 5000 active elements (0 deactivated)
15:36:06 INFO    opendrift.models.radionuclides:779: Number of transformations: 45
15:36:06 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:06 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2605), np.int64(2394), np.int64(1)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:06 INFO    opendrift:2237: 2026-10-05 03:30:00 - step 8 of 96 - 5000 active elements (0 deactivated)
15:36:08 INFO    opendrift.models.radionuclides:779: Number of transformations: 38
15:36:08 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:08 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2643), np.int64(2356), np.int64(1)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:08 INFO    opendrift:2237: 2026-10-05 04:00:00 - step 9 of 96 - 5000 active elements (0 deactivated)
15:36:08 INFO    opendrift.models.radionuclides:779: Number of transformations: 52
15:36:08 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:08 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2695), np.int64(2304), np.int64(1)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:08 INFO    opendrift:2237: 2026-10-05 04:30:00 - step 10 of 96 - 5000 active elements (0 deactivated)
15:36:11 INFO    opendrift.models.radionuclides:779: Number of transformations: 49
15:36:11 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:11 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2744), np.int64(2255), np.int64(1)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:11 INFO    opendrift:2237: 2026-10-05 05:00:00 - step 11 of 96 - 5000 active elements (0 deactivated)
15:36:11 INFO    opendrift.models.radionuclides:779: Number of transformations: 52
15:36:11 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:11 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2794), np.int64(2205), np.int64(1)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:11 INFO    opendrift:2237: 2026-10-05 05:30:00 - step 12 of 96 - 5000 active elements (0 deactivated)
15:36:14 INFO    opendrift.models.radionuclides:779: Number of transformations: 50
15:36:14 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:14 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2844), np.int64(2155), np.int64(1)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:14 INFO    opendrift:2237: 2026-10-05 06:00:00 - step 13 of 96 - 5000 active elements (0 deactivated)
15:36:14 INFO    opendrift.models.radionuclides:779: Number of transformations: 48
15:36:14 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:14 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2892), np.int64(2107), np.int64(1)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:14 INFO    opendrift:2237: 2026-10-05 06:30:00 - step 14 of 96 - 5000 active elements (0 deactivated)
15:36:18 INFO    opendrift.models.radionuclides:779: Number of transformations: 40
15:36:18 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:18 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2932), np.int64(2067), np.int64(1)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:18 INFO    opendrift:2237: 2026-10-05 07:00:00 - step 15 of 96 - 5000 active elements (0 deactivated)
15:36:18 INFO    opendrift.models.radionuclides:779: Number of transformations: 57
15:36:18 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:18 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(2985), np.int64(2012), np.int64(3)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:18 INFO    opendrift:2237: 2026-10-05 07:30:00 - step 16 of 96 - 5000 active elements (0 deactivated)
15:36:22 INFO    opendrift.models.radionuclides:779: Number of transformations: 40
15:36:22 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:22 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3025), np.int64(1972), np.int64(3)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:22 INFO    opendrift:2237: 2026-10-05 08:00:00 - step 17 of 96 - 5000 active elements (0 deactivated)
15:36:23 INFO    opendrift.models.radionuclides:779: Number of transformations: 39
15:36:23 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:23 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3064), np.int64(1933), np.int64(3)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:23 INFO    opendrift:2237: 2026-10-05 08:30:00 - step 18 of 96 - 5000 active elements (0 deactivated)
15:36:26 INFO    opendrift.models.radionuclides:779: Number of transformations: 30
15:36:26 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:26 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3094), np.int64(1903), np.int64(3)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:26 INFO    opendrift:2237: 2026-10-05 09:00:00 - step 19 of 96 - 5000 active elements (0 deactivated)
15:36:26 INFO    opendrift.models.radionuclides:779: Number of transformations: 43
15:36:26 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:26 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3137), np.int64(1860), np.int64(3)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:26 INFO    opendrift:2237: 2026-10-05 09:30:00 - step 20 of 96 - 5000 active elements (0 deactivated)
15:36:29 INFO    opendrift.models.radionuclides:779: Number of transformations: 39
15:36:29 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:29 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3172), np.int64(1824), np.int64(4)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:29 INFO    opendrift:2237: 2026-10-05 10:00:00 - step 21 of 96 - 5000 active elements (0 deactivated)
15:36:29 INFO    opendrift.models.radionuclides:779: Number of transformations: 46
15:36:29 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:29 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3218), np.int64(1778), np.int64(4)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:29 INFO    opendrift:2237: 2026-10-05 10:30:00 - step 22 of 96 - 5000 active elements (0 deactivated)
15:36:32 INFO    opendrift.models.radionuclides:779: Number of transformations: 45
15:36:32 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:32 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3259), np.int64(1735), np.int64(6)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:32 INFO    opendrift:2237: 2026-10-05 11:00:00 - step 23 of 96 - 5000 active elements (0 deactivated)
15:36:32 INFO    opendrift.models.radionuclides:779: Number of transformations: 33
15:36:32 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:32 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3292), np.int64(1702), np.int64(6)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:32 INFO    opendrift:2237: 2026-10-05 11:30:00 - step 24 of 96 - 5000 active elements (0 deactivated)
15:36:34 INFO    opendrift.models.radionuclides:779: Number of transformations: 32
15:36:34 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:34 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3320), np.int64(1672), np.int64(8)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:34 INFO    opendrift:2237: 2026-10-05 12:00:00 - step 25 of 96 - 5000 active elements (0 deactivated)
15:36:34 INFO    opendrift.models.radionuclides:779: Number of transformations: 39
15:36:34 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:34 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3359), np.int64(1633), np.int64(8)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:34 INFO    opendrift:2237: 2026-10-05 12:30:00 - step 26 of 96 - 5000 active elements (0 deactivated)
15:36:37 INFO    opendrift.models.radionuclides:779: Number of transformations: 38
15:36:37 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:37 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3395), np.int64(1596), np.int64(9)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:37 INFO    opendrift:2237: 2026-10-05 13:00:00 - step 27 of 96 - 5000 active elements (0 deactivated)
15:36:37 INFO    opendrift.models.radionuclides:779: Number of transformations: 48
15:36:37 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:37 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3437), np.int64(1553), np.int64(10)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:37 INFO    opendrift:2237: 2026-10-05 13:30:00 - step 28 of 96 - 5000 active elements (0 deactivated)
15:36:40 INFO    opendrift.models.radionuclides:779: Number of transformations: 31
15:36:40 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:40 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3464), np.int64(1525), np.int64(11)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:40 INFO    opendrift:2237: 2026-10-05 14:00:00 - step 29 of 96 - 5000 active elements (0 deactivated)
15:36:40 INFO    opendrift.models.radionuclides:779: Number of transformations: 33
15:36:40 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:40 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3493), np.int64(1495), np.int64(12)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:40 INFO    opendrift:2237: 2026-10-05 14:30:00 - step 30 of 96 - 5000 active elements (0 deactivated)
15:36:44 INFO    opendrift.models.radionuclides:779: Number of transformations: 33
15:36:44 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:44 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3520), np.int64(1465), np.int64(15)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:44 INFO    opendrift:2237: 2026-10-05 15:00:00 - step 31 of 96 - 5000 active elements (0 deactivated)
15:36:44 INFO    opendrift.models.radionuclides:779: Number of transformations: 41
15:36:44 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:44 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3555), np.int64(1427), np.int64(18)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:44 INFO    opendrift:2237: 2026-10-05 15:30:00 - step 32 of 96 - 5000 active elements (0 deactivated)
15:36:48 INFO    opendrift.models.radionuclides:779: Number of transformations: 25
15:36:48 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:48 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3578), np.int64(1403), np.int64(19)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:48 INFO    opendrift:2237: 2026-10-05 16:00:00 - step 33 of 96 - 5000 active elements (0 deactivated)
15:36:49 INFO    opendrift.models.radionuclides:779: Number of transformations: 37
15:36:49 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:49 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3615), np.int64(1366), np.int64(19)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:49 INFO    opendrift:2237: 2026-10-05 16:30:00 - step 34 of 96 - 5000 active elements (0 deactivated)
15:36:52 INFO    opendrift.models.radionuclides:779: Number of transformations: 32
15:36:52 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:52 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3647), np.int64(1335), np.int64(18)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:52 INFO    opendrift:2237: 2026-10-05 17:00:00 - step 35 of 96 - 5000 active elements (0 deactivated)
15:36:52 INFO    opendrift.models.radionuclides:779: Number of transformations: 32
15:36:52 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:52 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3677), np.int64(1304), np.int64(19)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:52 INFO    opendrift:2237: 2026-10-05 17:30:00 - step 36 of 96 - 5000 active elements (0 deactivated)
15:36:55 INFO    opendrift.models.radionuclides:779: Number of transformations: 32
15:36:55 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:55 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3703), np.int64(1275), np.int64(22)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:55 INFO    opendrift:2237: 2026-10-05 18:00:00 - step 37 of 96 - 5000 active elements (0 deactivated)
15:36:55 INFO    opendrift.models.radionuclides:779: Number of transformations: 37
15:36:55 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:55 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3728), np.int64(1244), np.int64(28)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:55 INFO    opendrift:2237: 2026-10-05 18:30:00 - step 38 of 96 - 5000 active elements (0 deactivated)
15:36:58 INFO    opendrift.models.radionuclides:779: Number of transformations: 24
15:36:58 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:58 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3748), np.int64(1222), np.int64(30)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:58 INFO    opendrift:2237: 2026-10-05 19:00:00 - step 39 of 96 - 5000 active elements (0 deactivated)
15:36:58 INFO    opendrift.models.radionuclides:779: Number of transformations: 30
15:36:58 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:36:58 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3768), np.int64(1198), np.int64(34)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:36:58 INFO    opendrift:2237: 2026-10-05 19:30:00 - step 40 of 96 - 5000 active elements (0 deactivated)
15:37:01 INFO    opendrift.models.radionuclides:779: Number of transformations: 27
15:37:01 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:01 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3791), np.int64(1173), np.int64(36)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:01 INFO    opendrift:2237: 2026-10-05 20:00:00 - step 41 of 96 - 5000 active elements (0 deactivated)
15:37:01 INFO    opendrift.models.radionuclides:779: Number of transformations: 22
15:37:01 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:01 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3807), np.int64(1154), np.int64(39)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:01 INFO    opendrift:2237: 2026-10-05 20:30:00 - step 42 of 96 - 5000 active elements (0 deactivated)
15:37:03 INFO    opendrift.models.radionuclides:779: Number of transformations: 22
15:37:03 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:03 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3819), np.int64(1138), np.int64(43)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:03 INFO    opendrift:2237: 2026-10-05 21:00:00 - step 43 of 96 - 5000 active elements (0 deactivated)
15:37:03 INFO    opendrift.models.radionuclides:779: Number of transformations: 26
15:37:03 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:03 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3843), np.int64(1113), np.int64(44)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:03 INFO    opendrift:2237: 2026-10-05 21:30:00 - step 44 of 96 - 5000 active elements (0 deactivated)
15:37:05 INFO    opendrift.models.radionuclides:779: Number of transformations: 35
15:37:05 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:05 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3866), np.int64(1085), np.int64(49)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:05 INFO    opendrift:2237: 2026-10-05 22:00:00 - step 45 of 96 - 5000 active elements (0 deactivated)
15:37:06 INFO    opendrift.models.radionuclides:779: Number of transformations: 26
15:37:06 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:06 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3880), np.int64(1065), np.int64(55)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:06 INFO    opendrift:2237: 2026-10-05 22:30:00 - step 46 of 96 - 5000 active elements (0 deactivated)
15:37:08 INFO    opendrift.models.radionuclides:779: Number of transformations: 25
15:37:08 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:08 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3903), np.int64(1041), np.int64(56)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:08 INFO    opendrift:2237: 2026-10-05 23:00:00 - step 47 of 96 - 5000 active elements (0 deactivated)
15:37:08 INFO    opendrift.models.radionuclides:779: Number of transformations: 31
15:37:08 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:08 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3922), np.int64(1016), np.int64(62)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:08 INFO    opendrift:2237: 2026-10-05 23:30:00 - step 48 of 96 - 5000 active elements (0 deactivated)
15:37:12 INFO    opendrift.models.radionuclides:779: Number of transformations: 24
15:37:12 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:12 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3936), np.int64(998), np.int64(66)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:12 INFO    opendrift:2237: 2026-10-06 00:00:00 - step 49 of 96 - 5000 active elements (0 deactivated)
15:37:12 INFO    opendrift.models.radionuclides:779: Number of transformations: 23
15:37:12 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:12 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3939), np.int64(985), np.int64(76)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:12 INFO    opendrift:2237: 2026-10-06 00:30:00 - step 50 of 96 - 5000 active elements (0 deactivated)
15:37:14 INFO    opendrift.models.radionuclides:779: Number of transformations: 20
15:37:14 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:14 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3951), np.int64(969), np.int64(80)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:14 INFO    opendrift:2237: 2026-10-06 01:00:00 - step 51 of 96 - 5000 active elements (0 deactivated)
15:37:14 INFO    opendrift.models.radionuclides:779: Number of transformations: 27
15:37:14 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:14 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3966), np.int64(948), np.int64(86)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:14 INFO    opendrift:2237: 2026-10-06 01:30:00 - step 52 of 96 - 5000 active elements (0 deactivated)
15:37:16 INFO    opendrift.models.radionuclides:779: Number of transformations: 18
15:37:17 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:17 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3974), np.int64(935), np.int64(91)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:17 INFO    opendrift:2237: 2026-10-06 02:00:00 - step 53 of 96 - 5000 active elements (0 deactivated)
15:37:17 INFO    opendrift.models.radionuclides:779: Number of transformations: 30
15:37:17 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:17 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(3992), np.int64(912), np.int64(96)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:17 INFO    opendrift:2237: 2026-10-06 02:30:00 - step 54 of 96 - 5000 active elements (0 deactivated)
15:37:19 INFO    opendrift.models.radionuclides:779: Number of transformations: 31
15:37:19 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:19 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4001), np.int64(893), np.int64(106)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:19 INFO    opendrift:2237: 2026-10-06 03:00:00 - step 55 of 96 - 5000 active elements (0 deactivated)
15:37:19 INFO    opendrift.models.radionuclides:779: Number of transformations: 35
15:37:19 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:19 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4020), np.int64(866), np.int64(114)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:19 INFO    opendrift:2237: 2026-10-06 03:30:00 - step 56 of 96 - 5000 active elements (0 deactivated)
15:37:21 INFO    opendrift.models.radionuclides:779: Number of transformations: 32
15:37:21 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:21 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4038), np.int64(841), np.int64(121)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:21 INFO    opendrift:2237: 2026-10-06 04:00:00 - step 57 of 96 - 5000 active elements (0 deactivated)
15:37:21 INFO    opendrift.models.radionuclides:779: Number of transformations: 21
15:37:21 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:21 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4047), np.int64(826), np.int64(127)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:21 INFO    opendrift:2237: 2026-10-06 04:30:00 - step 58 of 96 - 5000 active elements (0 deactivated)
15:37:24 INFO    opendrift.models.radionuclides:779: Number of transformations: 20
15:37:24 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:24 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4049), np.int64(815), np.int64(136)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:24 INFO    opendrift:2237: 2026-10-06 05:00:00 - step 59 of 96 - 5000 active elements (0 deactivated)
15:37:24 INFO    opendrift.models.radionuclides:779: Number of transformations: 24
15:37:24 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:24 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4039), np.int64(808), np.int64(153)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:24 INFO    opendrift:2237: 2026-10-06 05:30:00 - step 60 of 96 - 5000 active elements (0 deactivated)
15:37:26 INFO    opendrift.models.radionuclides:779: Number of transformations: 30
15:37:26 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:26 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4049), np.int64(788), np.int64(163)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:26 INFO    opendrift:2237: 2026-10-06 06:00:00 - step 61 of 96 - 5000 active elements (0 deactivated)
15:37:26 INFO    opendrift.models.radionuclides:779: Number of transformations: 27
15:37:26 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:26 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4058), np.int64(770), np.int64(172)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:26 INFO    opendrift:2237: 2026-10-06 06:30:00 - step 62 of 96 - 5000 active elements (0 deactivated)
15:37:28 INFO    opendrift.models.radionuclides:779: Number of transformations: 26
15:37:28 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:28 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4056), np.int64(758), np.int64(186)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:28 INFO    opendrift:2237: 2026-10-06 07:00:00 - step 63 of 96 - 5000 active elements (0 deactivated)
15:37:28 INFO    opendrift.models.radionuclides:779: Number of transformations: 17
15:37:28 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:28 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4065), np.int64(746), np.int64(189)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:28 INFO    opendrift:2237: 2026-10-06 07:30:00 - step 64 of 96 - 5000 active elements (0 deactivated)
15:37:30 INFO    opendrift.models.radionuclides:779: Number of transformations: 25
15:37:30 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:30 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4078), np.int64(727), np.int64(195)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:30 INFO    opendrift:2237: 2026-10-06 08:00:00 - step 65 of 96 - 5000 active elements (0 deactivated)
15:37:30 INFO    opendrift.models.radionuclides:779: Number of transformations: 18
15:37:30 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:30 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4086), np.int64(714), np.int64(200)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:30 INFO    opendrift:2237: 2026-10-06 08:30:00 - step 66 of 96 - 5000 active elements (0 deactivated)
15:37:32 INFO    opendrift.models.radionuclides:779: Number of transformations: 24
15:37:32 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:32 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4098), np.int64(696), np.int64(206)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:32 INFO    opendrift:2237: 2026-10-06 09:00:00 - step 67 of 96 - 5000 active elements (0 deactivated)
15:37:32 INFO    opendrift.models.radionuclides:779: Number of transformations: 25
15:37:32 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:32 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4111), np.int64(679), np.int64(210)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:32 INFO    opendrift:2237: 2026-10-06 09:30:00 - step 68 of 96 - 5000 active elements (0 deactivated)
15:37:35 INFO    opendrift.models.radionuclides:779: Number of transformations: 16
15:37:35 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:35 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4109), np.int64(672), np.int64(219)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:35 INFO    opendrift:2237: 2026-10-06 10:00:00 - step 69 of 96 - 5000 active elements (0 deactivated)
15:37:35 INFO    opendrift.models.radionuclides:779: Number of transformations: 16
15:37:35 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:35 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4111), np.int64(663), np.int64(226)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:35 INFO    opendrift:2237: 2026-10-06 10:30:00 - step 70 of 96 - 5000 active elements (0 deactivated)
15:37:37 INFO    opendrift.models.radionuclides:779: Number of transformations: 26
15:37:37 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:37 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4119), np.int64(646), np.int64(235)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:37 INFO    opendrift:2237: 2026-10-06 11:00:00 - step 71 of 96 - 5000 active elements (0 deactivated)
15:37:37 INFO    opendrift.models.radionuclides:779: Number of transformations: 25
15:37:37 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:37 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4134), np.int64(628), np.int64(238)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:37 INFO    opendrift:2237: 2026-10-06 11:30:00 - step 72 of 96 - 5000 active elements (0 deactivated)
15:37:39 INFO    opendrift.models.radionuclides:779: Number of transformations: 16
15:37:39 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:39 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4142), np.int64(616), np.int64(242)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:39 INFO    opendrift:2237: 2026-10-06 12:00:00 - step 73 of 96 - 5000 active elements (0 deactivated)
15:37:39 INFO    opendrift.models.radionuclides:779: Number of transformations: 17
15:37:39 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:39 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4153), np.int64(604), np.int64(243)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:39 INFO    opendrift:2237: 2026-10-06 12:30:00 - step 74 of 96 - 5000 active elements (0 deactivated)
15:37:41 INFO    opendrift.models.radionuclides:779: Number of transformations: 30
15:37:41 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:41 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4165), np.int64(583), np.int64(252)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:41 INFO    opendrift:2237: 2026-10-06 13:00:00 - step 75 of 96 - 5000 active elements (0 deactivated)
15:37:42 INFO    opendrift.models.radionuclides:779: Number of transformations: 16
15:37:42 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:42 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4167), np.int64(574), np.int64(259)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:42 INFO    opendrift:2237: 2026-10-06 13:30:00 - step 76 of 96 - 5000 active elements (0 deactivated)
15:37:44 INFO    opendrift.models.radionuclides:779: Number of transformations: 18
15:37:44 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:44 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4167), np.int64(567), np.int64(266)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:44 INFO    opendrift:2237: 2026-10-06 14:00:00 - step 77 of 96 - 5000 active elements (0 deactivated)
15:37:44 INFO    opendrift.models.radionuclides:779: Number of transformations: 23
15:37:44 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:44 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4176), np.int64(553), np.int64(271)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:44 INFO    opendrift:2237: 2026-10-06 14:30:00 - step 78 of 96 - 5000 active elements (0 deactivated)
15:37:47 INFO    opendrift.models.radionuclides:779: Number of transformations: 21
15:37:47 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:48 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4179), np.int64(542), np.int64(279)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:48 INFO    opendrift:2237: 2026-10-06 15:00:00 - step 79 of 96 - 5000 active elements (0 deactivated)
15:37:48 INFO    opendrift.models.radionuclides:779: Number of transformations: 16
15:37:48 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:48 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4179), np.int64(535), np.int64(286)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:48 INFO    opendrift:2237: 2026-10-06 15:30:00 - step 80 of 96 - 5000 active elements (0 deactivated)
15:37:50 INFO    opendrift.models.radionuclides:779: Number of transformations: 28
15:37:50 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:50 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4183), np.int64(519), np.int64(298)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:50 INFO    opendrift:2237: 2026-10-06 16:00:00 - step 81 of 96 - 5000 active elements (0 deactivated)
15:37:50 INFO    opendrift.models.radionuclides:779: Number of transformations: 24
15:37:50 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:50 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4191), np.int64(506), np.int64(303)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:50 INFO    opendrift:2237: 2026-10-06 16:30:00 - step 82 of 96 - 5000 active elements (0 deactivated)
15:37:52 INFO    opendrift.models.radionuclides:779: Number of transformations: 18
15:37:52 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:52 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4193), np.int64(500), np.int64(307)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:52 INFO    opendrift:2237: 2026-10-06 17:00:00 - step 83 of 96 - 5000 active elements (0 deactivated)
15:37:52 INFO    opendrift.models.radionuclides:779: Number of transformations: 25
15:37:52 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:52 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4192), np.int64(488), np.int64(320)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:52 INFO    opendrift:2237: 2026-10-06 17:30:00 - step 84 of 96 - 5000 active elements (0 deactivated)
15:37:54 INFO    opendrift.models.radionuclides:779: Number of transformations: 19
15:37:54 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:54 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4195), np.int64(478), np.int64(327)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:54 INFO    opendrift:2237: 2026-10-06 18:00:00 - step 85 of 96 - 5000 active elements (0 deactivated)
15:37:54 INFO    opendrift.models.radionuclides:779: Number of transformations: 19
15:37:54 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:54 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4206), np.int64(464), np.int64(330)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:54 INFO    opendrift:2237: 2026-10-06 18:30:00 - step 86 of 96 - 5000 active elements (0 deactivated)
15:37:56 INFO    opendrift.models.radionuclides:779: Number of transformations: 25
15:37:56 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:56 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4207), np.int64(451), np.int64(342)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:56 INFO    opendrift:2237: 2026-10-06 19:00:00 - step 87 of 96 - 5000 active elements (0 deactivated)
15:37:57 INFO    opendrift.models.radionuclides:779: Number of transformations: 19
15:37:57 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:57 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4200), np.int64(445), np.int64(355)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:57 INFO    opendrift:2237: 2026-10-06 19:30:00 - step 88 of 96 - 5000 active elements (0 deactivated)
15:37:58 INFO    opendrift.models.radionuclides:779: Number of transformations: 13
15:37:58 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:58 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4201), np.int64(439), np.int64(360)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:58 INFO    opendrift:2237: 2026-10-06 20:00:00 - step 89 of 96 - 5000 active elements (0 deactivated)
15:37:59 INFO    opendrift.models.radionuclides:779: Number of transformations: 16
15:37:59 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:37:59 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4199), np.int64(434), np.int64(367)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:37:59 INFO    opendrift:2237: 2026-10-06 20:30:00 - step 90 of 96 - 5000 active elements (0 deactivated)
15:38:01 INFO    opendrift.models.radionuclides:779: Number of transformations: 19
15:38:01 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:38:01 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4190), np.int64(429), np.int64(381)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:38:01 INFO    opendrift:2237: 2026-10-06 21:00:00 - step 91 of 96 - 5000 active elements (0 deactivated)
15:38:01 INFO    opendrift.models.radionuclides:779: Number of transformations: 29
15:38:01 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:38:01 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4189), np.int64(418), np.int64(393)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:38:01 INFO    opendrift:2237: 2026-10-06 21:30:00 - step 92 of 96 - 5000 active elements (0 deactivated)
15:38:03 INFO    opendrift.models.radionuclides:779: Number of transformations: 22
15:38:03 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:38:03 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4179), np.int64(412), np.int64(409)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:38:03 INFO    opendrift:2237: 2026-10-06 22:00:00 - step 93 of 96 - 5000 active elements (0 deactivated)
15:38:03 INFO    opendrift.models.radionuclides:779: Number of transformations: 25
15:38:03 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:38:03 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4174), np.int64(403), np.int64(423)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:38:03 INFO    opendrift:2237: 2026-10-06 22:30:00 - step 94 of 96 - 5000 active elements (0 deactivated)
15:38:05 INFO    opendrift.models.radionuclides:779: Number of transformations: 13
15:38:05 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:38:05 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4175), np.int64(399), np.int64(426)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:38:05 INFO    opendrift:2237: 2026-10-06 23:00:00 - step 95 of 96 - 5000 active elements (0 deactivated)
15:38:05 INFO    opendrift.models.radionuclides:779: Number of transformations: 22
15:38:05 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:38:05 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4171), np.int64(392), np.int64(437)] ['LMM', 'Particle reversible', 'Sediment reversible']
15:38:05 INFO    opendrift:2237: 2026-10-06 23:30:00 - step 96 of 96 - 5000 active elements (0 deactivated)
15:38:07 INFO    opendrift.models.radionuclides:779: Number of transformations: 19
15:38:07 INFO    opendrift.models.radionuclides:972: Number of resuspended particles: 0
15:38:07 INFO    opendrift.models.radionuclides:1026: Speciation: [np.int64(4176), np.int64(381), np.int64(443)] ['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[ns] 392B 2026-10-05 .....
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 14.89 ....
    sea_water_salinity               (trajectory, time) float32 980kB 30.48 ....
    horizontal_diffusivity           (trajectory, time) float32 980kB 10.0 .....
    upward_sea_water_velocity        (trajectory, time) float32 980kB -9.992e...
    conc3                            (trajectory, time) float32 980kB 0.001 ....
Attributes: (12/154)
    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:02:23.605696
    geospatial_vertical_min:                                      -51.761597
    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:02:04.9  total
 0:00:00.0  preparing
 0:02:03.1  reading
 0:00:06.0  interpolation
 0:00:00.1  interpolation_time
 0:00:01.7  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:01.0  total
 0:00:00.0  preparing
 0:00:01.0  reading
 0:00:00.0  masking
--------------------
Performance:
 2:23.6 total time
    8.0 configuration
    0.0 preparing main loop
      0.0 moving elements to ocean
 2:15.4 main loop
      1.7 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-10-05 00:00:00 UTC
        Present: 2026-10-07 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                             :   4176
Particle reversible             :    381
Sediment reversible             :    443
Number of transformations:
['        0', '       24', '      469']
['     2403', '        0', '        2']
['       28', '        0', '        0']
15:38:17 INFO    opendrift:4823: Saving animation to /root/project/docs/source/gallery/animations/example_radionuclides_0.gif...
15:38:33 INFO    opendrift:3241: Time to make animation: 0:00:26.053081
../_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-10-07T00:00:00   Mean windspeed: 1.8 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()
15:38:34 INFO    opendrift:4823: Saving animation to /root/project/docs/source/gallery/animations/example_radionuclides_1.gif...
15:38:42 INFO    opendrift:3459: Time to make animation: 0:00:08.821050
../_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: (3 minutes 7.366 seconds)

Gallery generated by Sphinx-Gallery