def make_OceanDrift_object(self):
self.o = OceanDrift(loglevel=30)
self.fake_eddy = reader_ArtificialOceanEddy.Reader(2, 62)
self.o.use_block = False
self.reader_landmask = reader_global_landmask.Reader(
extent=[-1.5, 59, 7, 64])
self.o.add_reader([self.fake_eddy, self.reader_landmask])
def make_OceanDrift_object(self):
self.o = OceanDrift(loglevel=30)
self.fake_eddy = reader_ArtificialOceanEddy.Reader(2, 62)
self.o.use_block = False
self.reader_basemap = reader_basemap_landmask.Reader(
llcrnrlon=-1.5, llcrnrlat=59,
urcrnrlon=7, urcrnrlat=64, resolution='c')
self.o.add_reader([self.fake_eddy, self.reader_basemap])
def test_seed(self):
"""Test seeding"""
self.o = OpenOil()
self.fake_eddy = reader_ArtificialOceanEddy.Reader(2, 62)
self.fake_eddy.start_time = datetime(2015, 1, 1)
self.o.add_reader([self.fake_eddy])
self.o.seed_elements(lon=4, lat=60, number=100,
time=self.fake_eddy.start_time)
self.assertEqual(len(self.o.elements), 0)
self.assertEqual(len(self.o.elements_deactivated), 0)
self.assertEqual(len(self.o.elements_scheduled), 100)
Double gyre current field from
http://shaddenlab.berkeley.edu/uploads/LCS-tutorial/examples.html
"""
import numpy as np
from datetime import datetime, timedelta
from opendrift.readers import reader_global_landmask
from opendrift.readers import reader_double_gyre
from opendrift.readers import reader_ArtificialOceanEddy
from opendrift.models.oceandrift import OceanDrift
o = OceanDrift(loglevel=0) # Set loglevel to 0 for debug information
fake_eddy = reader_ArtificialOceanEddy.Reader(2, 62)
reader_landmask = reader_global_landmask.Reader(llcrnrlon=-1.5,
llcrnrlat=59,
urcrnrlon=7,
urcrnrlat=64)
lon = 2.0
lat = 63.0
# Close to Station M
runs = []
leg = []
i = 0
for scheme in ['euler', 'runge-kutta', 'runge-kutta4']:
for time_step in [1800, 3600 * 3]:
