Python WindBlow Examples

Example #1

 def test_windblow(self):
     o = WindBlow(loglevel=20)
     reader_arome = reader_netCDF_CF_generic.Reader(o.test_data_folder() + '2Feb2016_Nordic_sigma_3d/AROME_MetCoOp_00_DEF.nc_20160202_subset')
     o.add_reader([reader_arome])
     lat = 67.711251; lon = 13.556971  # Lofoten
     o.seed_elements(lon, lat, radius=5000, number=1000,
                     time=reader_arome.start_time)
     o.run(steps=48*4, time_step=900)
     self.assertAlmostEqual(o.elements.lon.max(), 17.54, 2)

Example #2

File: test_models.py Project: knutfrode/opendrift

 def test_windblow(self):
     o = WindBlow(loglevel=0)
     reader_arome = reader_netCDF_CF_generic.Reader(o.test_data_folder() + '2Feb2016_Nordic_sigma_3d/AROME_MetCoOp_00_DEF.nc_20160202_subset')
     o.add_reader([reader_arome])
     lat = 67.711251; lon = 13.556971  # Lofoten
     o.seed_elements(lon, lat, radius=5000, number=1000,
                     time=reader_arome.start_time)
     o.run(steps=48*4, time_step=900)
     self.assertAlmostEqual(o.elements.lon.max(), 17.54, 2)

Example #3

File: example_windblow.py Project: ValV/drifting

#!/usr/bin/env python
"""
Wind blow model
==================================
"""

from opendrift.readers import reader_netCDF_CF_generic
from opendrift.models.windblow import WindBlow

o = WindBlow(loglevel=20)  # Set loglevel to 0 for debug information

#%%
# Example of elements blowing with the wind, also over land

#reader_arome = reader_netCDF_CF_generic.Reader('https://thredds.met.no/thredds/dodsC/mepslatest/meps_lagged_6_h_latest_2_5km_latest.nc')
reader_arome = reader_netCDF_CF_generic.Reader(
    o.test_data_folder() +
    '2Feb2016_Nordic_sigma_3d/AROME_MetCoOp_00_DEF.nc_20160202_subset')

o.add_reader([reader_arome])

#%%
# Seeding some particles
lat = 68.5
lon = 16.0  # Lofoten
o.seed_elements(lon,
                lat,
                radius=5000,
                number=1000,
                time=reader_arome.start_time)

Example #4

#!/usr/bin/env python
"""
Wind blow model
==================================
"""

from opendrift.readers import reader_netCDF_CF_generic
from opendrift.models.windblow import WindBlow

o = WindBlow(loglevel=20)  # Set loglevel to 0 for debug information

#%%
# Example of elements blowing with the wind, also over land

#reader_arome = reader_netCDF_CF_generic.Reader('https://thredds.met.no/thredds/dodsC/mepslatest/meps_lagged_6_h_latest_2_5km_latest.nc')
reader_arome = reader_netCDF_CF_generic.Reader(o.test_data_folder() +
    '2Feb2016_Nordic_sigma_3d/AROME_MetCoOp_00_DEF_20160202_subset.nc')

o.add_reader([reader_arome])

#%%
# Seeding some particles
lat = 68.5; lon = 16.0  # Lofoten
o.seed_elements(lon, lat, radius=5000, number=1000,
                time=reader_arome.start_time)

#%%
# Running model for 48 hours
o.run(steps=48*4, time_step=900, time_step_output=3600)

#%%

Example #5

File: example_windblow.py Project: p1heidary/opendrift

#!/usr/bin/env python

from opendrift.readers import reader_netCDF_CF_generic
from opendrift.models.windblow import WindBlow

o = WindBlow(loglevel=0)  # Set loglevel to 0 for debug information

##############################################################
# Example of elements blowing with the wind, also over land
##############################################################

#reader_arome = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/meps25files/meps_det_extracted_2_5km_latest.nc')
reader_arome = reader_netCDF_CF_generic.Reader(o.test_data_folder() + 
    '2Feb2016_Nordic_sigma_3d/AROME_MetCoOp_00_DEF.nc_20160202_subset')

o.add_reader([reader_arome])

# Seeding some particles
lat = 68.5; lon = 16.0  # Lofoten
o.seed_elements(lon, lat, radius=5000, number=1000,
                time=reader_arome.start_time)

# Running model
o.run(steps=48*4, time_step=900)

# Print and plot results
print o
o.plot(buffer=.5)
o.animation()

Example #6

#!/usr/bin/env python
"""
Wind blow model
==================================
"""

from opendrift.readers import reader_netCDF_CF_generic
from opendrift.models.windblow import WindBlow

o = WindBlow(loglevel=0)  # Set loglevel to 0 for debug information

##############################################################
# Example of elements blowing with the wind, also over land
##############################################################

#reader_arome = reader_netCDF_CF_generic.Reader('http://thredds.met.no/thredds/dodsC/meps25files/meps_det_extracted_2_5km_latest.nc')
reader_arome = reader_netCDF_CF_generic.Reader(o.test_data_folder() +
    '2Feb2016_Nordic_sigma_3d/AROME_MetCoOp_00_DEF.nc_20160202_subset')

o.add_reader([reader_arome])

# Seeding some particles
lat = 68.5; lon = 16.0  # Lofoten
o.seed_elements(lon, lat, radius=5000, number=1000,
                time=reader_arome.start_time)

# Running model
o.run(steps=48*4, time_step=900)

# Print and plot results
print(o)