schemes, using an idealised (analytical) eddy current field.
Double gyre current field from
http://shaddenlab.berkeley.edu/uploads/LCS-tutorial/examples.html
"""
from datetime import datetime, timedelta
from opendrift.readers import reader_double_gyre
from opendrift.models.oceandrift import OceanDrift
o = OceanDrift(loglevel=20) # Set loglevel to 0 for debug information
o.fallback_values['land_binary_mask'] = 0
o.set_config('drift:scheme', 'runge-kutta4')
double_gyre = reader_double_gyre.Reader(epsilon=.25, omega=0.628, A=0.1)
print(double_gyre)
o.add_reader(double_gyre)
x = [.9]
y = [.5]
lon, lat = double_gyre.xy2lonlat(x, y)
o.seed_elements(lon,
lat,
radius=.1,
number=1000,
time=double_gyre.initial_time)
o.run(duration=timedelta(seconds=10), time_step=0.1)
from datetime import timedelta
from opendrift.readers import reader_double_gyre
from opendrift.models.oceandrift import OceanDrift
from src.double_gyre.dispersion_utils import exportParticlesToDataset
# PARAMETERS
A = 0.1
epsilon = 0.25
NB_PARTS = 25000
altSeed = 3
o = OceanDrift(loglevel=20) # Set loglevel to 0 for debug information
o.fallback_values["land_binary_mask"] = 0
o.set_config("drift:scheme", "runge-kutta4")
double_gyre = reader_double_gyre.Reader(epsilon=epsilon, omega=6.28 / 10, A=A)
print(double_gyre)
o.add_reader(double_gyre)
x = [0.95]
y = [0.5]
lon, lat = double_gyre.xy2lonlat(x, y)
o.seed_elements(
lon,
lat,
radius=0.1,
number=NB_PARTS + altSeed,
time=[double_gyre.initial_time, double_gyre.initial_time + timedelta(seconds=50)],
)
