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')
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()
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', fast=True) o.animation(filename='codegg_forecast.gif', fast=True) #%% # .. image:: /gallery/animations/codegg_forecast.gif o.plot_vertical_distribution()
o.add_reader([reader_norkyst, reader_arome]) # spawn NEA cod eggs at defined position and time time = reader_arome.start_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 # 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()
# 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()