o.seed_elements(3.3,
60.,
radius=3000,
time=reader_current.start_time,
water_line_length=100,
keel_depth=90,
number=100)
print('Starting free run .../n')
print('Start time: ' + str(o.start_time))
#######################
# Running model
#######################
o.run(time_step=3600, steps=steps)
#########################
# Print and plot results
#########################
o.plot()
o.animation(filename='example_det.gif')
#%%
# .. image:: /gallery/animations/example_det.gif
print('############## Latitudes:', o.history['lat'])
print('############## Longitudes:', o.history['lon'])
# Seeding elements
#
# Icebergs are moved with the ocean current as per Barker et al (2004),
# in addition to a fraction of the wind speed (wind_drift_factor).
# This factor depends on the properties of the elements.
# Default empirical values are:
# - Wind drift fraction: 0.018 (1.8 %) (Garret 1985)
# - Iceberg size: Keel dept = 60m
# Waterline length = 90.5m
# NB! Iceberg size is irrelevant for current_reader with 1D z-profile
o.seed_elements(3.3, 61.3, radius=3000, number=500,
time=reader_current.start_time)
#%%
# Run model
print('Starting free run .../n')
print('Start time: ' + str(o.start_time))
o.run(time_step=3600, steps=steps)
#%%
# Print and plot results
o.plot(fast=True)
o.animation(fast=True)
#%%
# .. image:: /gallery/animations/example_openberg_stat_0.gif
# - Wind drift fraction: 0.018 (1.8 %) (Garret 1985)
# - Iceberg size: Keel dept = 60m
# Waterline length = 90.5m
o.seed_elements(3.3,
60.,
radius=3000,
time=reader_current.start_time,
water_line_length=100,
keel_depth=90,
number=100)
print('Starting free run .../n')
print('Start time: ' + str(o.start_time))
#######################
# Running model
#######################
o.run(time_step=3600, steps=steps)
#########################
# Print and plot results
#########################
o.plot(filename='example_det.pdf')
o.animation(filename='example_det.mp4')
print('############## Latitudes:', o.history['lat'])
print('############## Longitudes:', o.history['lon'])
# Icebergs are moved with the ocean current as per Barker et al (2004),
# in addition to a fraction of the wind speed (wind_drift_factor).
# This factor depends on the properties of the elements.
# Default empirical values are:
# - Wind drift fraction: 0.018 (1.8 %) (Garret 1985)
# - Iceberg size: Keel dept = 60m
# Waterline length = 90.5m
# NB! Iceberg size is irrelevant for current_reader with 1D z-profile
o.seed_elements(3.3,
61.3,
radius=3000,
number=10,
time=reader_current.start_time)
#%%
# Run model
print('Starting free run .../n')
print('Start time: ' + str(o.start_time))
o.run(time_step=3600, steps=steps)
#%%
# Print and plot results
o.plot()
o.animation()
#%%
# .. image:: /gallery/animations/example_openberg_stat_0.gif
# This factor depends on the properties of the elements.
# Default empirical values are:
# - Wind drift fraction: 0.018 (1.8 %) (Garret 1985)
# - Iceberg size: Keel dept = 60m
# Waterline length = 90.5m
# NB! Iceberg size is irrelevant for current_reader with 1D z-profile
o.seed_elements(3.3,
61.3,
radius=3000,
number=10,
time=reader_current.start_time)
#######################
# Run model
#######################
print('Starting free run .../n')
print('Start time: ' + str(o.start_time))
o.run(time_step=3600, steps=steps)
#########################
# Print and plot results
#########################
o.plot(filename='example_stat.pdf')
o.animation(filename='example_stat.mp4')
#
# print('############## Latitudes:', o.history['lat'])
# print('############## Longitudes:',o.history['lon'])