class TestArray(unittest.TestCase):
"""Tests for (non-scalar) LagrangianArray"""
def setUp(self):
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])
def test_seed(self):
"""Test seeding"""
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)
class TestArray(unittest.TestCase):
"""Tests for (non-scalar) LagrangianArray"""
def setUp(self):
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])
def test_seed(self):
"""Test seeding"""
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)
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
# Arome
#reader_arome = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/arome25/arome_metcoop_default2_5km_latest.nc')
reader_arome = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc')
# Norkyst
#reader_norkyst = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')
reader_norkyst = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc')
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(llcrnrlon=2, llcrnrlat=59,
urcrnrlon=6, urcrnrlat=62,
resolution='i', projection='merc')
o.add_reader([reader_basemap, reader_norkyst, reader_arome])
# Seeding some particles
lon = [3.6, 5.1]; lat = [61., 59.6]; # Outside Bergen
# Seed elements along cone, e.g. ship track with
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
# HYCOM
#reader_hycom = reader_netCDF_CF_generic.Reader('http://tds.hycom.org/thredds/dodsC/GLBu0.08/expt_19.1/2010/3hrly')
#print reader_hycom
reader_globcurrent = reader_netCDF_CF_generic.Reader('http://tds0.ifremer.fr/thredds/dodsC/CLS-L4-CUREUL_HS-ALT_SUM-V01.0_FULL_TIME_SERIE')
# OceanWind
reader_oceanwind = reader_netCDF_CF_generic.Reader('http://www.ncdc.noaa.gov/thredds/dodsC/oceanwinds6hr')
#print reader_oceanwind
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(llcrnrlon=100, llcrnrlat=5,
urcrnrlon=120, urcrnrlat=15, resolution='h')
# Add readers
print 'adding...'
o.add_reader([reader_globcurrent, reader_oceanwind, reader_basemap])
# Seed some particles
lat=10.228248; lon=106.973337
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
# Arome
#reader_arome = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/arome25/arome_metcoop_default2_5km_latest.nc')
reader_arome = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc')
# Norkyst
#reader_norkyst = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')
reader_norkyst = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc')
# WAM10
#reader_wam10 = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/sea/wam/wam10_be') # 6 hourly aggregates two months back
#reader_wam10 = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/fou-hi/wam10h/g10kmwave.nc')
# GlobCurrent
#reader_globcurrent = reader_netCDF_CF_generic.Reader('http://tds0.ifremer.fr/thredds/dodsC/CLS-L4-CUREUL_HS-ALT_SUM-V01.0_FULL_TIME_SERIE') # Total
#reader_globcurrent = reader_netCDF_CF_generic.Reader('http://tds0.ifremer.fr/thredds/dodsC/GC_MOD_TIDE_GLO_010_FES2012_FULL_TIME_SERIE') # FES Tidal
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(llcrnrlon=4,
llcrnrlat=59.8,
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
# Arome
reader_arome = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc')
#reader_arome = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/arome25/arome_metcoop_default2_5km_latest.nc')
# Norkyst
reader_norkyst = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc')
#reader_norkyst = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(
llcrnrlon=4, llcrnrlat=59.7,
urcrnrlon=7, urcrnrlat=61.5,
resolution='h', projection='merc')
o.add_reader([reader_basemap, reader_norkyst, reader_arome])
# Seeding some particles
lon = 4.9; lat = 60.0; # Outside Bergen
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
#reader_hycom = reader_netCDF_CF_generic.Reader('http://tds.hycom.org/thredds/dodsC/GLBu0.08/expt_19.1/2010/3hrly')
#print reader_hycom0
reader_globcurrent = reader_netCDF_CF_generic.Reader('http://tds0.ifremer.fr/thredds/dodsC/CLS-L4-CUREUL_HS-ALT_SUM-V01.0_FULL_TIME_SERIE') # Total
reader_oceanwind = reader_netCDF_CF_generic.Reader('http://www.ncdc.noaa.gov/thredds/dodsC/oceanwinds6hr')
#print reader_oceanwind
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(
llcrnrlon=-94, llcrnrlat=20,
urcrnrlon=-80, urcrnrlat=32, resolution='i')
# Add readers
#o.add_reader([reader_basemap, reader_hycom, reader_oceanwind])
o.add_reader([reader_basemap, reader_globcurrent, reader_oceanwind])
#o.add_reader([reader_basemap, reader_globcurrent])
# Seed some particles
lon = -88.387161; lat = 28.736669 # Macondo location
starttime = datetime(2010, 4, 21, 6, 0, 0) # 4 hours after explosion
time = [starttime, starttime + timedelta(hours=24*40)]
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=20) # Set loglevel to 0 for debug information
# Arome
reader_arome = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc')
# Norkyst
reader_norkyst = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc')
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(llcrnrlon=3, llcrnrlat=59.5,
urcrnrlon=7, urcrnrlat=62,
resolution='h', projection='merc')
o.add_reader([reader_basemap, reader_norkyst, reader_arome])
# Seeding some particles
lon = 4.3; lat = 60.0; # Outside Bergen
time = [reader_arome.start_time,
reader_arome.start_time + timedelta(hours=30)]
# Seed oil elements at defined position and time
#!/usr/bin/env python from datetime import datetime, timedelta import numpy as np #from models.oceandrift import OceanDrift from models.openoil import OpenOil #o = OceanDrift(loglevel=50) o = OpenOil(loglevel=50) o.config['processes']['diffusion'] = False # We do not care about landmask or current for this seeding demonstration, # so we simple specify fallback_values instead of adding any readers o.fallback_values['land_binary_mask'] = 0 o.fallback_values['x_sea_water_velocity'] = 0 o.fallback_values['y_sea_water_velocity'] = 0 time=datetime(2016, 1, 20, 12, 30, 0) print '\n' + '='*70 print 'Seeding a single element at a point:' print 'o.seed_elements(lon=4, lat=60, time=time)' print '='*70 o.seed_elements(lon=4, lat=60, time=time) o.run(steps=1) o.plot(buffer=1) print '\n' + '='*70
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
# Arome
reader_arome = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc')
#reader_arome = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/arome25/arome_metcoop_default2_5km_latest.nc')
# Norkyst
reader_norkyst = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc')
#reader_norkyst = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(llcrnrlon=4,
llcrnrlat=59.7,
urcrnrlon=7,
urcrnrlat=61.5,
resolution='h',
projection='merc')
o.add_reader([reader_basemap, reader_norkyst, reader_arome])
# Seeding some particles
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
# Arome
reader_arome = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc')
# Norkyst
reader_norkyst = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc')
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(llcrnrlon=3.5,
llcrnrlat=60.4,
urcrnrlon=7,
urcrnrlat=61.8,
resolution='h',
projection='merc')
o.add_reader([reader_basemap, reader_norkyst, reader_arome])
############################################################
# Seed oil particles within contour detected from satellite
############################################################
#!/usr/bin/env python from datetime import datetime, timedelta import numpy as np #from models.oceandrift import OceanDrift from models.openoil import OpenOil #o = OceanDrift(loglevel=50) o = OpenOil(loglevel=50) o.config['processes']['diffusion'] = False # We do not care about landmask or current for this seeding demonstration, # so we simple specify fallback_values instead of adding any readers o.fallback_values['land_binary_mask'] = 0 o.fallback_values['x_sea_water_velocity'] = 0 o.fallback_values['y_sea_water_velocity'] = 0 time = datetime(2016, 1, 20, 12, 30, 0) print '\n' + '=' * 70 print 'Seeding a single element at a point:' print 'o.seed_elements(lon=4, lat=60, time=time)' print '=' * 70 o.seed_elements(lon=4, lat=60, time=time) o.run(steps=1) o.plot(buffer=1) print '\n' + '=' * 70 print 'Seeding 100 elements within a radius of 1000 m:' print 'o.seed_elements(lon=4, lat=60, number=100, radius=1000, time=time)'
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
o.io_import_file('openoil.nc')
print o
o.plot(buffer=.1)
o.plot_property('mass_oil')
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
# Arome
#reader_arome = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/arome25/arome_metcoop_default2_5km_latest.nc')
# Nordc4
reader_nordic4 = reader_netCDF_CF_generic.Reader(
'http://thredds.met.no/thredds/dodsC/fou-hi/nordic4km-1h/Nordic-4km_SURF_1h_avg_00.nc'
)
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(llcrnrlon=2,
llcrnrlat=78,
urcrnrlon=22,
urcrnrlat=82.3,
resolution='i',
projection='merc')
o.add_reader([reader_basemap, reader_nordic4])
o.fallback_values['y_wind'] = 35 # Northerly wind, to get stranding in ice
# Seeding some particles
lon = 8
#!/usr/bin/env python
import os
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
if not os.path.exists('openoil.nc'):
raise ValueError('Please run example.py first to generate a '
'netCDF file to be imported.')
o.io_import_file('openoil.nc')
print o
o.plot(buffer=.1)
o.plot_property('mass_oil')
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
# Arome
reader_arome = reader_netCDF_CF_generic.Reader(
'http://thredds.met.no/thredds/dodsC/arome25/arome_metcoop_default2_5km_latest.nc'
)
#reader_arome = reader_netCDF_CF_generic.Reader('/disk1/data/opendrift_testdata/arome.nc')
# Norkyst
reader_norkyst = reader_netCDF_CF_generic.Reader(
'http://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')
#reader_norkyst = reader_netCDF_CF_generic.Reader('/disk1/data/opendrift_testdata/norkyst800.nc')
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(llcrnrlon=15.4,
llcrnrlat=68.6,
urcrnrlon=20.0,
urcrnrlat=70.5,
resolution='h',
projection='merc')
o.add_reader([reader_basemap, reader_norkyst, reader_arome])
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
# Arome
#reader_arome = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/arome25/arome_metcoop_default2_5km_latest.nc')
# Nordc4
reader_nordic4 = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/fou-hi/nordic4km-1h/Nordic-4km_SURF_1h_avg_00.nc')
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(
llcrnrlon=2, llcrnrlat=78,
urcrnrlon=22, urcrnrlat=82.3,
resolution='i', projection='merc')
o.add_reader([reader_basemap, reader_nordic4])
o.fallback_values['y_wind'] = 35 # Northerly wind, to get stranding in ice
# Seeding some particles
lon = 8; lat = 79.5; # West of Svalbard
#time = datetime(2015, 9, 22, 6, 0, 0)
#time = [reader_nordic4.start_time,
# reader_nordic4.start_time + timedelta(hours=30)]
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
# Arome
reader_arome = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/arome25/arome_metcoop_default2_5km_latest.nc')
# Norkyst
reader_norkyst = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(llcrnrlon=4, llcrnrlat=59.8,
urcrnrlon=6, urcrnrlat=61,
resolution='h', projection='merc')
o.add_reader([reader_basemap, reader_norkyst, reader_arome])
# Seeding some particles
lon = 4.8; lat = 60.0; # Outside Bergen
time = None
time = reader_arome.start_time
# Seed elements at defined position and time
import numpy as np
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=20) # Set loglevel to 0 for debug information
# Arome
reader_arome = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc')
#reader_arome = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/arome25/arome_metcoop_default2_5km_latest.nc')
# Norkyst
reader_norkyst = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc')
#reader_norkyst = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(llcrnrlon=3, llcrnrlat=59.5,
urcrnrlon=7, urcrnrlat=62,
resolution='h', projection='merc')
o.add_reader([reader_basemap, reader_norkyst, reader_arome])
# Seeding some particles
lon = 4.2; lat = 60.0; # Outside Bergen
time = [reader_arome.start_time,
reader_arome.start_time + timedelta(hours=30)]
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
# HYCOM
#reader_hycom = reader_netCDF_CF_generic.Reader('http://tds.hycom.org/thredds/dodsC/GLBu0.08/expt_19.1/2010/3hrly')
#print reader_hycom
reader_globcurrent = reader_netCDF_CF_generic.Reader('http://tds0.ifremer.fr/thredds/dodsC/CLS-L4-CUREUL_HS-ALT_SUM-V02.0_FULL_TIME_SERIE')
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(llcrnrlon=100, llcrnrlat=5,
urcrnrlon=120, urcrnrlat=15, resolution='h')
# OceanWind
try:
reader_oceanwind = reader_netCDF_CF_generic.Reader(
'http://www.ncdc.noaa.gov/thredds/dodsC/oceanwinds6hr')
print reader_oceanwind
o.add_reader([reader_globcurrent, reader_oceanwind, reader_basemap])
except:
o.add_reader([reader_globcurrent, reader_basemap])
# Seed some particles
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
# Arome
#reader_arome = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/arome25/arome_metcoop_default2_5km_latest.nc')
reader_arome = reader_netCDF_CF_generic.Reader('test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc')
# Norkyst
#reader_norkyst = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')
reader_norkyst = reader_netCDF_CF_generic.Reader('test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc')
# WAM10
#reader_wam10 = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/sea/wam/wam10_be') # 6 hourly aggregates two months back
#reader_wam10 = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/fou-hi/wam10h/g10kmwave.nc')
# GlobCurrent
#reader_globcurrent = reader_netCDF_CF_generic.Reader('http://tds0.ifremer.fr/thredds/dodsC/CLS-L4-CUREUL_HS-ALT_SUM-V01.0_FULL_TIME_SERIE') # Total
#reader_globcurrent = reader_netCDF_CF_generic.Reader('http://tds0.ifremer.fr/thredds/dodsC/GC_MOD_TIDE_GLO_010_FES2012_FULL_TIME_SERIE') # FES Tidal
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(llcrnrlon=4, llcrnrlat=59.8,
urcrnrlon=6, urcrnrlat=61,
resolution='h', projection='merc')
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
reader_hycom = reader_netCDF_CF_generic.Reader('http://tds.hycom.org/thredds/dodsC/GLBu0.08/expt_19.1/2010/3hrly')
#print reader_hycom0
#reader_globcurrent = reader_netCDF_CF_generic.Reader('http://tds0.ifremer.fr/thredds/dodsC/CLS-L4-CUREUL_HS-ALT_SUM-V01.0_FULL_TIME_SERIE') # Total
reader_oceanwind = reader_netCDF_CF_generic.Reader('http://www.ncdc.noaa.gov/thredds/dodsC/oceanwinds6hr')
#print reader_oceanwind
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(
llcrnrlon=-94, llcrnrlat=20,
urcrnrlon=-80, urcrnrlat=32, resolution='i')
# Add readers
o.add_reader([reader_basemap, reader_hycom, reader_oceanwind])
#o.add_reader([reader_basemap, reader_globcurrent, reader_oceanwind])
# Seed some particles
lon = -88.387161; lat = 28.736669 # Macondo location
starttime = datetime(2010, 4, 21, 6, 0, 0) # 4 hours after explosion
time = [starttime, starttime + timedelta(hours=24*40)]
o.seed_elements(lon, lat, radius=0, number=5000, time=time)
def setUp(self):
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])
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=20) # Set loglevel to 0 for debug information
# Arome
reader_arome = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/arome_subset_16Nov2015.nc')
#reader_arome = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/arome25/arome_metcoop_default2_5km_latest.nc')
# Norkyst
reader_norkyst = reader_netCDF_CF_generic.Reader(
'test_data/16Nov2015_NorKyst_z_surface/norkyst800_subset_16Nov2015.nc')
#reader_norkyst = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(llcrnrlon=3,
llcrnrlat=59.5,
urcrnrlon=7,
urcrnrlat=62,
resolution='h',
projection='merc')
o.add_reader([reader_basemap, reader_norkyst, reader_arome])
# Seeding some particles
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
reader_hycom = reader_netCDF_CF_generic.Reader(
'http://tds.hycom.org/thredds/dodsC/GLBu0.08/expt_19.1/2010/3hrly')
#print reader_hycom0
#reader_globcurrent = reader_netCDF_CF_generic.Reader('http://tds0.ifremer.fr/thredds/dodsC/CLS-L4-CUREUL_HS-ALT_SUM-V01.0_FULL_TIME_SERIE') # Total
reader_oceanwind = reader_netCDF_CF_generic.Reader(
'http://www.ncdc.noaa.gov/thredds/dodsC/oceanwinds6hr')
#print reader_oceanwind
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(llcrnrlon=-94,
llcrnrlat=20,
urcrnrlon=-80,
urcrnrlat=32,
resolution='i')
# Add readers
o.add_reader([reader_basemap, reader_hycom, reader_oceanwind])
#o.add_reader([reader_basemap, reader_globcurrent, reader_oceanwind])
# Seed some particles
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
# Arome
reader_arome = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/arome25/arome_metcoop_default2_5km_latest.nc')
#reader_arome = reader_netCDF_CF_generic.Reader('/disk1/data/opendrift_testdata/arome.nc')
# Norkyst
reader_norkyst = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')
#reader_norkyst = reader_netCDF_CF_generic.Reader('/disk1/data/opendrift_testdata/norkyst800.nc')
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(
llcrnrlon=15.4, llcrnrlat=68.6,
urcrnrlon=20.0, urcrnrlat=70.5,
resolution='h', projection='merc')
o.add_reader([reader_basemap, reader_norkyst, reader_arome])
# Seeding some particles
latstart = 68.988911
lonstart = 16.040701
#latstart = 69.477754
def setUp(self):
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])
#!/usr/bin/env python
from datetime import datetime, timedelta
from readers import reader_basemap_landmask
from readers import reader_netCDF_CF_generic
from models.openoil import OpenOil
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
# Arome
reader_arome = reader_netCDF_CF_generic.Reader(
'http://thredds.met.no/thredds/dodsC/arome25/arome_metcoop_default2_5km_latest.nc'
)
# Norkyst
reader_norkyst = reader_netCDF_CF_generic.Reader(
'http://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be')
# Landmask (Basemap)
reader_basemap = reader_basemap_landmask.Reader(llcrnrlon=4,
llcrnrlat=59.8,
urcrnrlon=6,
urcrnrlat=61,
resolution='h',
projection='merc')
o.add_reader([reader_basemap, reader_norkyst, reader_arome])
# Seeding some particles
lon = 4.8
# You should have received a copy of the GNU General Public License
# along with OpenDrift. If not, see <http://www.gnu.org/licenses/>.
#
# Copyright 2015, Knut-Frode Dagestad, MET Norway
import sys
import argparse
import numpy as np
sys.path.append("..")
from models.openoil import OpenOil
if __name__ == '__main__':
parser = argparse.ArgumentParser()
parser.add_argument('filename',
help='<OpenDriftaoutput filename (netCDF)>')
parser.add_argument('-b',
dest='buffer',
default=1.0,
help='Buffer around plot in degrees lon/lat.')
args = parser.parse_args()
o = OpenOil(loglevel=0) # Set loglevel to 0 for debug information
print o
print type(np.array(args.buffer))
o.io_import_file(args.filename)
o.plot(buffer=np.float(args.buffer))