Note
Go to the end to download the full example code.
Macondo
from datetime import datetime, timedelta
from opendrift.readers import reader_netCDF_CF_generic
from opendrift.models.openoil import OpenOil
This is a highly simplified 2D simulation, for illustration purpose only
o = OpenOil(loglevel=20) # Set loglevel to 0 for debug information
12:50:39 INFO opendrift.models.basemodel:515: OpenDriftSimulation initialised (version 1.11.13 / v1.11.13-99-gd2132d3)
For this datasource which does not contain standard_name, we impose a variable mapping
reader_globcurrent = None
try:
reader_globcurrent = reader_netCDF_CF_generic.Reader(
'https://tds0.ifremer.fr/thredds/dodsC/GLOBCURRENT-L4-CUREUL_HS-ALT_SUM-V03.0',
standard_name_mapping={'eastward_eulerian_current_velocity': 'x_sea_water_velocity',
'northward_eulerian_current_velocity': 'y_sea_water_velocity'})
except:
print('Thredds server not available, cannot run example')
if reader_globcurrent is not None:
reader_oceanwind = reader_netCDF_CF_generic.Reader('https://tds0.ifremer.fr/thredds/dodsC/CERSAT-GLO-CLIM_WIND_L4-OBS_FULL_TIME_SERIE')
# Add readers
o.add_reader([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*30)]
o.seed_elements(lon, lat, radius=0, number=15000, time=time, oil_type='LIGHT LOUISIANNA SWEET, BP')
# Run model
print(o)
o.run(duration=timedelta(days=30),
time_step=timedelta(hours=3),
time_step_output=timedelta(hours=6))
# Print and plot results
print(o)
o.plot(fast=True)
o.plot_oil_budget()
o.animation(fast=True)
o.animation(fast=True, frames = range(0, 120, 4), fps = 1)
12:50:39 INFO opendrift.readers.reader_netCDF_CF_generic:102: Opening dataset: https://tds0.ifremer.fr/thredds/dodsC/GLOBCURRENT-L4-CUREUL_HS-ALT_SUM-V03.0
Thredds server not available, cannot run example
Total running time of the script: (0 minutes 8.125 seconds)