#!/usr/bin/env python
"""
Radionuclides
=============
"""
from opendrift.readers import reader_netCDF_CF_generic, reader_ROMS_native
from opendrift.models.radionuclides import RadionuclideDrift
from datetime import timedelta, datetime
import numpy as np
o = RadionuclideDrift(loglevel=20,
seed=0) # Set loglevel to 0 for debug information
# Norkyst
#reader_norkyst = reader_netCDF_CF_generic.Reader(o.test_data_folder() + '/14Jan2016_NorKyst_z_3d/NorKyst-800m_ZDEPTHS_his_00_3Dsubset.nc')
reader_norkyst = reader_netCDF_CF_generic.Reader(
'https://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')
o.add_reader([reader_norkyst])
# Adjusting some configuration
o.set_config('drift:vertical_mixing', True)
#o.set_config('environment:fallback:ocean_vertical_diffusivity', 0)
#o.set_config('vertical_mixing:diffusivitymodel','constant') # include settling without vertical turbulent mixing
#o.set_config('vertical_mixing:diffusivitymodel','environment') # include settling without vertical turbulent mixing
# Vertical mixing requires fast time step
o.set_config('vertical_mixing:timestep', 600.) # seconds
#%%
# Activate the desired species
#!/usr/bin/env python
"""
Radionuclides
=============
"""
from opendrift.readers import reader_netCDF_CF_generic, reader_ROMS_native
from opendrift.models.radionuclides import RadionuclideDrift
from datetime import timedelta, datetime
import numpy as np
from numpy import dtype
o = RadionuclideDrift(loglevel=20,
seed=0) # Set loglevel to 0 for debug information
# Norkyst
#reader_norkyst = reader_netCDF_CF_generic.Reader(o.test_data_folder() + '/14Jan2016_NorKyst_z_3d/NorKyst-800m_ZDEPTHS_his_00_3Dsubset.nc')
reader_norkyst = reader_netCDF_CF_generic.Reader(
'http://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')
o.add_reader([reader_norkyst])
td = datetime.today()
time = datetime(td.year, td.month, td.day, 0)
#latseed= 61.2; lonseed= 4.3 # Sognesjen
latseed = 59.0
lonseed = 10.75 # Hvaler/Koster
ntraj = 2000
iniz = np.random.rand(
#!/usr/bin/env python
from opendrift.readers import reader_netCDF_CF_generic, reader_ROMS_native
from opendrift.models.radionuclides import RadionuclideDrift
from datetime import timedelta, datetime
import numpy as np
from numpy import dtype
o = RadionuclideDrift(loglevel=20,
seed=0) # Set loglevel to 0 for debug information
# Arome
#reader_arome = reader_netCDF_CF_generic.Reader(o.test_data_folder() + '16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc')
#reader_arome = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/meps25files/meps_det_extracted_2_5km_latest.nc')
# Norkyst
#reader_norkyst = reader_netCDF_CF_generic.Reader(o.test_data_folder() + '/14Jan2016_NorKyst_z_3d/NorKyst-800m_ZDEPTHS_his_00_3Dsubset.nc')
reader_norkyst = reader_netCDF_CF_generic.Reader(
'http://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')
#reader_bokna = reader_ROMS_native.Reader('/home/magnes/projects/MARTINI/data/norfjords_160m_his.nc4_2015110?01-2015110?00')
#o.add_reader([reader_norkyst, reader_arome])
o.add_reader([reader_norkyst])
#o.add_reader([reader_bokna])
# release radionuclides at defined position and time
#time = reader_norkyst.start_time
td = datetime.today()
time = datetime(td.year, td.month, td.day, 0)
#time = datetime(2015,11,2,1) # Boknafj
def test_radio_nuclide(test_data):
o = RadionuclideDrift(loglevel=20,
seed=0) # Set loglevel to 0 for debug information
# Norkyst
reader_norkyst = reader_netCDF_CF_generic.Reader(
test_data +
'/14Jan2016_NorKyst_z_3d/NorKyst-800m_ZDEPTHS_his_00_3Dsubset.nc')
o.add_reader([reader_norkyst])
# Adjusting some configuration
o.set_config('drift:vertical_mixing', True)
o.set_config(
'vertical_mixing:diffusivitymodel',
'environment') # include settling without vertical turbulent mixing
# Vertical mixing requires fast time step
o.set_config('vertical_mixing:timestep', 600.) # seconds
o.set_config('drift:horizontal_diffusivity', 10)
o.set_config('radionuclide:particle_diameter', 5.e-6) # m
o.set_config('radionuclide:transformations:Kd', 2.e0) # (m3/kg)
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:transfer_setup', 'Bokna_137Cs')
o.set_config('general:coastline_action', 'previous')
o.set_config('seed:LMM_fraction', .45)
o.set_config('seed:particle_fraction', .55)
time = reader_norkyst.start_time
latseed = 62.0
lonseed = 4.3 # Bergen (?)
ntraj = 5000
iniz = np.random.rand(
ntraj) * -10. # seeding the radionuclides in the upper 10m
o.seed_elements(
lonseed,
latseed,
z=iniz,
radius=1000,
number=ntraj,
time=time,
)
o.run(steps=4, time_step=1800, time_step_output=3600)
#!/usr/bin/env python
"""
Radionuclides
=============
"""
from opendrift.readers import reader_netCDF_CF_generic, reader_ROMS_native
from opendrift.models.radionuclides import RadionuclideDrift
from datetime import timedelta, datetime
import numpy as np
o = RadionuclideDrift(loglevel=0,
seed=0) # Set loglevel to 0 for debug information
# Norkyst
#reader_norkyst = reader_netCDF_CF_generic.Reader(o.test_data_folder() + '/14Jan2016_NorKyst_z_3d/NorKyst-800m_ZDEPTHS_his_00_3Dsubset.nc')
reader_norkyst = reader_netCDF_CF_generic.Reader(
'https://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')
o.add_reader([reader_norkyst])
# Adjusting some configuration
o.set_config('drift:vertical_mixing', True)
#o.set_config('environment:constant:ocean_vertical_diffusivity', 0)
#o.set_config('vertical_mixing:diffusivitymodel','constant') # include settling without vertical turbulent mixing
o.set_config(
'vertical_mixing:diffusivitymodel',
'environment') # include settling without vertical turbulent mixing
# Vertical mixing requires fast time step
o.set_config('vertical_mixing:timestep', 600.) # seconds
o.set_config('drift:horizontal_diffusivity', 10)