number=500,
time=time,
diameter=0.0014,
neutral_buoyancy_salinity=31.25)
#%%
# Adjusting some configuration
o.set_config('processes:turbulentmixing', True)
#o.set_config('turbulentmixing:diffusivitymodel', 'windspeed_Sundby1983') # windspeed parameterization for eddy diffusivity
o.set_config('turbulentmixing:diffusivitymodel',
'environment') # use eddy diffusivity from ocean model
#%%
# Vertical mixing requires fast time step
o.set_config('turbulentmixing:timestep', 60.) # seconds
#%%
# Running model (until end of driver data)
o.run(steps=96, time_step=3600)
#%%
# Print and plot results
print(o)
o.plot(fast=True)
o.animation(fast=True)
#%%
# .. image:: /gallery/animations/example_codegg_0.gif
o.plot_vertical_distribution()
number=500,
time=time,
diameter=0.0014,
neutral_buoyancy_salinity=31.25)
o.seed_elements(13.,
67.8,
z=-40,
radius=2000,
number=500,
time=time,
diameter=0.0014,
neutral_buoyancy_salinity=31.25)
# Adjusting some configuration
o.set_config('processes:turbulentmixing', True)
#o.set_config('turbulentmixing:diffusivitymodel', 'windspeed_Sundby1983') # windspeed parameterization for eddy diffusivity
o.set_config('turbulentmixing:diffusivitymodel',
'environment') # use eddy diffusivity from ocean model
# Vertical mixing requires fast time step
o.set_config('turbulentmixing:timestep', 60.) # seconds
# Running model (until end of driver data)
o.run(steps=96, time_step=3600)
# Print and plot results
print(o)
o.plot(filename='codegg_forecast.png')
o.animation(filename='codegg_forecast.gif')
o.plot_vertical_distribution(filename='codegg_vertical.png')
# spawn NEA cod eggs at defined position and time
o.seed_elements(14. , 68.1, z=-40, radius=2000, number=500,
time=time, diameter=0.0014, neutral_buoyancy_salinity=31.25)
o.seed_elements(12.5, 68.5, z=-40, radius=2000, number=500,
time=time, diameter=0.0014, neutral_buoyancy_salinity=31.25)
# Adjusting some configuration
o.set_config('processes:turbulentmixing', True)
o.set_config('turbulentmixing:diffusivitymodel', 'windspeed_Sundby1983') # windspeed parameterization for eddy diffusivity
#o.set_config('turbulentmixing:diffusivitymodel', 'environment') # use eddy diffusivity from ocean model
# Vertical resolution and time step should be adjusted so to avoid getting
# output warnings like 'DEBUG: WARNING! some elements have p+q>1.'
o.set_config('turbulentmixing:timestep', 10.) # seconds
o.set_config('turbulentmixing:verticalresolution', 2.) # m
o.set_config('turbulentmixing:TSprofiles', True) # update T,S frequently in each inner loop for turbulent
o.set_config('turbulentmixing:TSprofiles', False) # False is faster
# Running model (until end of driver data)
o.run(steps=6*2, time_step=1800)
# Print and plot results
print o
o.plot()
o.animation()
o.plot_vertical_distribution()
