#lons=[-0.8, 0.4]; lats=[59.9, 59.95] # NorKyst border
#lons=[4.8, 5.1]; lats=[59.9, 59.95] # Western Norway coast
#lons=[13.11, 13.13]; lats=[67.81, 67.80] # Lofoten
#lons=[19.37, 19.33]; lats=[70.32, 70.34] # Troms
lons = [3.8, 3.82]
lats = [59.6, 59.61] # North Sea
o.seed_elements(
lon=lons,
lat=lats,
time=[datetime.now() - timedelta(hours=3),
datetime.now()],
number=1000,
radius=[0, 1000],
cone=True,
**args)
print(o)
o.run(duration=timedelta(hours=24),
time_step=timestep,
time_step_output=1800,
outfile='halo_test.nc')
print(o)
o.animation()
if case == 'oil':
o.plot()
o.plot_oil_budget()
o.plot_property('water_fraction')
o.animation_profile()
number=1000,
radius=100,
z=0,
oiltype='TIA JUANA LIGHT, OIL & GAS'
) #'EKOFISK BLEND, STATOIL' similar ent.
o3.run(duration=timedelta(hours=12), time_step=900, time_step_output=3600)
#%%
# Plotting and comparing
print('#######################')
print('Entrainment rate (heavy)', np.mean(o2.oil_wave_entrainment_rate()))
print('Entrainment rate (light)', np.mean(o3.oil_wave_entrainment_rate()))
print('Viscosity (heavy)', np.mean(o2.elements.viscosity))
print('Viscosity (light)', np.mean(o3.elements.viscosity))
print('Density (heavy)', np.mean(o2.elements.density))
print('Density (light)', np.mean(o3.elements.density))
print('#######################')
o2.plot_oil_budget()
o3.plot_oil_budget()
legend = ['TIA JUANA HEAVY', 'TIA JUANA LIGHT']
o2.animation_profile(compare=o3, legend=legend)
#%%
# .. image:: /gallery/animations/example_entrainment_rate_oil_types_0.gif
o2.animation(compare=o3, legend=legend, fast=True)
#%%
# .. image:: /gallery/animations/example_entrainment_rate_oil_types_1.gif
