Note
Go to the end to download the full example code.
Vietnam
from datetime import datetime, timedelta
from opendrift.models.leeway import Leeway
o = Leeway(loglevel=20) # Set loglevel to 0 for debug information
# Adding readers for global Thredds datasets:
# - Ocean forecast from global Hycom
# - Weather forecast from NOAA/NCEP
o.add_readers_from_list([
'https://tds.hycom.org/thredds/dodsC/FMRC_ESPC-D-V02_uv3z/FMRC_ESPC-D-V02_uv3z_best.ncd',
'https://pae-paha.pacioos.hawaii.edu/thredds/dodsC/ncep_global/NCEP_Global_Atmospheric_Model_best.ncd'])
# Seed some particles
object_type = 26 # 26 = Life-raft, no ballast
o.seed_elements(lon=107.8, lat=10.0, radius=1000, number=1000,
object_type=object_type, time=datetime.utcnow())
# Run model
o.run(duration=timedelta(days=3),
time_step=timedelta(hours=1),
time_step_output=timedelta(hours=3))
# Print and plot results
print(o)
o.plot(fast=True)
o.animation(fast=True)