Note
Go to the end to download the full example code
Thredds resources for GUI
from datetime import datetime
from opendrift.models.oceandrift import OceanDrift
from opendrift.readers.reader_netCDF_CF_generic import Reader
o = OceanDrift(loglevel=0)
thredds_resources = open(o.test_data_folder()+'../../opendrift/scripts/data_sources.txt').readlines()
times = {}
14:03:50 DEBUG opendrift.models.basemodel:625: Adding 18 config items from basemodel
14:03:50 DEBUG opendrift.models.basemodel:625: Adding 5 config items from basemodel
14:03:50 DEBUG opendrift.models.basemodel:625: Adding 48 config items from basemodel
14:03:51 INFO opendrift.models.basemodel:539: OpenDriftSimulation initialised (version 1.10.7 / v1.10.6-119-g1da5bec)
14:03:51 DEBUG opendrift.models.basemodel:625: Adding 15 config items from oceandrift
14:03:51 DEBUG opendrift.models.basemodel:637: Overwriting config item seed:z
Open each thredds dataset to check contents and spatial coverage
for t in thredds_resources:
if t.startswith('http') and 'nrt.cmems' not in t:
start = datetime.now()
print('\n#%%\n%s\n' % t)
r = Reader(t)
print(r)
ts = str(datetime.now() - start)
times[t] = ts
print('Time to open reader: ', ts)
if r.global_coverage():
lscale = 'coarse'
else:
lscale = 'intermediate'
r.plot(lscale=lscale)
#%
https://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be
14:03:51 INFO opendrift.readers.reader_netCDF_CF_generic:92: Opening dataset: https://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be
14:03:52 DEBUG opendrift.readers.reader_netCDF_CF_generic:109: Finding coordinate variables.
14:03:52 DEBUG opendrift.readers.reader_netCDF_CF_generic:122: Parsing CF grid mapping dictionary: {'grid_mapping_name': 'polar_stereographic', 'straight_vertical_longitude_from_pole': 70.0, 'latitude_of_projection_origin': 90.0, 'standard_parallel': 60.0, 'false_easting': 3192800.0, 'false_northing': 1784000.0, 'semi_major_axis': 6378137.0, 'semi_minor_axis': 6356752.3142, 'proj4': '+proj=stere +lat_0=90 +lat_ts=60 +lon_0=70 +x_0=3192800 +y_0=1784000 +a=6378137 +b=6356752.3142 +units=m +no_defs +type=crs'}
/opt/conda/envs/opendrift/lib/python3.11/site-packages/pyproj/crs/crs.py:1286: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems
proj = self._crs.to_proj4(version=version)
14:03:52 DEBUG opendrift.readers.reader_netCDF_CF_generic:310: Skipped variables without standard_name: ['angle', 'tke', 'ubar', 'vbar']
14:03:52 DEBUG opendrift.readers.basereader.variables:593: Setting buffer size 25 for reader https://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be
, assuming a maximum average speed of 5 m/s and time span of 1:00:00
14:03:52 DEBUG opendrift.readers.basereader:176: Variable mapping: ['sea_floor_depth_below_sea_level'] -> ['land_binary_mask'] is not activated
14:03:52 DEBUG opendrift.readers.basereader.variables:548: Adding variable mapping: ['x_wind', 'y_wind'] -> wind_speed
14:03:52 DEBUG opendrift.readers.basereader.variables:548: Adding variable mapping: ['x_sea_water_velocity', 'y_sea_water_velocity'] -> sea_water_speed
===========================
Reader: https://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be
Projection:
+proj=stere +lat_0=90 +lat_ts=60 +lon_0=70 +x_0=3192800 +y_0=1784000 +a=6378137 +b=6356752.3142 +units=m +no_defs +type=crs
Coverage: [degrees]
xmin: 0.000000 xmax: 2080800.000000 step: 800 numx: 2602
ymin: 0.000000 ymax: 720800.000000 step: 800 numy: 902
Corners (lon, lat):
( -1.58, 58.50) ( 23.71, 75.32)
( 9.19, 55.91) ( 38.06, 70.03)
Vertical levels [m]:
[ -0. -3. -10. -15. -25. -50. -75. -100. -150. -200.
-250. -300. -500. -1000. -2000. -3000.]
Available time range:
start: 2017-02-20 00:00:00 end: 2023-09-06 00:00:00 step: 1:00:00
57337 times (3300 missing)
Variables:
sea_floor_depth_below_sea_level
latitude
longitude
x_wind
y_wind
sea_water_salinity
sea_water_temperature
x_sea_water_velocity
eastward_sea_water_velocity
y_sea_water_velocity
northward_sea_water_velocity
upward_sea_water_velocity
sea_surface_elevation
wind_speed - derived from ['x_wind', 'y_wind']
sea_water_speed - derived from ['x_sea_water_velocity', 'y_sea_water_velocity']
===========================
Time to open reader: 0:00:01.676639
#%
https://thredds.met.no/thredds/dodsC/mepslatest/meps_lagged_6_h_latest_2_5km_latest.nc
14:03:54 INFO opendrift.readers.reader_netCDF_CF_generic:92: Opening dataset: https://thredds.met.no/thredds/dodsC/mepslatest/meps_lagged_6_h_latest_2_5km_latest.nc
14:03:55 DEBUG opendrift.readers.reader_netCDF_CF_generic:109: Finding coordinate variables.
14:03:55 DEBUG opendrift.readers.reader_netCDF_CF_generic:211: 30 ensemble members available
14:03:55 DEBUG opendrift.readers.reader_netCDF_CF_generic:122: Parsing CF grid mapping dictionary: {'grid_mapping_name': 'lambert_conformal_conic', 'standard_parallel': array([63.3, 63.3]), 'longitude_of_central_meridian': 15.0, 'latitude_of_projection_origin': 63.3, 'earth_radius': 6371000.0}
/opt/conda/envs/opendrift/lib/python3.11/site-packages/pyproj/crs/crs.py:1286: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems
proj = self._crs.to_proj4(version=version)
14:03:55 DEBUG opendrift.readers.reader_netCDF_CF_generic:310: Skipped variables without standard_name: ['fog_area_fraction']
14:03:55 DEBUG opendrift.readers.basereader.variables:593: Setting buffer size 10 for reader https://thredds.met.no/thredds/dodsC/mepslatest/meps_lagged_6_h_latest_2_5km_latest.nc
, assuming a maximum average speed of 5 m/s and time span of 1:00:00
14:03:55 DEBUG opendrift.readers.basereader:176: Variable mapping: ['sea_floor_depth_below_sea_level'] -> ['land_binary_mask'] is not activated
14:03:55 DEBUG opendrift.readers.basereader.variables:548: Adding variable mapping: ['x_wind', 'y_wind'] -> wind_speed
===========================
Reader: https://thredds.met.no/thredds/dodsC/mepslatest/meps_lagged_6_h_latest_2_5km_latest.nc
Projection:
+proj=lcc +lat_0=63.3 +lon_0=15 +lat_1=63.3 +lat_2=63.3 +x_0=0 +y_0=0 +R=6371000 +units=m +no_defs +type=crs
Coverage: [degrees]
xmin: -1060084.000000 xmax: 1309916.000000 step: 2500 numx: 949
ymin: -1332517.875000 ymax: 1337482.125000 step: 2500 numy: 1069
Corners (lon, lat):
(-18.12, 72.76) ( 54.24, 71.58)
( 0.28, 50.32) ( 33.03, 49.77)
Vertical levels [m]:
Not specified
Available time range:
start: 2023-09-01 11:00:00 end: 2023-09-04 00:00:00 step: 1:00:00
62 times (0 missing)
Variables (30 ensemble members):
longitude
latitude
x_wind
y_wind
air_temperature
relative_humidity
cloud_area_fraction
air_pressure_at_sea_level
precipitation_amount
snowfall_amount
wind_speed_of_gust
wind_speed - derived from ['x_wind', 'y_wind']
===========================
Time to open reader: 0:00:00.941629
#%
https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
14:03:55 INFO opendrift.readers.reader_netCDF_CF_generic:92: Opening dataset: https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
14:03:57 DEBUG opendrift.readers.reader_netCDF_CF_generic:109: Finding coordinate variables.
14:03:57 DEBUG opendrift.readers.reader_netCDF_CF_generic:122: Parsing CF grid mapping dictionary: {'grid_mapping_name': 'lambert_conformal_conic', 'standard_parallel': array([77.5, 77.5], dtype=float32), 'longitude_of_central_meridian': -25.0, 'latitude_of_projection_origin': 77.5, 'earth_radius': 6371000.0, 'proj4': '+proj=lcc +lat_0=77.5 +lon_0=-25 +lat_1=77.5 +lat_2=77.5 +no_defs +R=6.371e+06'}
/opt/conda/envs/opendrift/lib/python3.11/site-packages/pyproj/crs/crs.py:1286: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems
proj = self._crs.to_proj4(version=version)
14:03:57 DEBUG opendrift.readers.reader_netCDF_CF_generic:310: Skipped variables without standard_name: ['angle']
14:03:57 DEBUG opendrift.readers.basereader.variables:593: Setting buffer size 10 for reader https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
, assuming a maximum average speed of 5 m/s and time span of 1:00:00
14:03:57 INFO opendrift.readers.basereader:166: Variable x_sea_water_velocity will be rotated from eastward_sea_water_velocity
14:03:57 INFO opendrift.readers.basereader:166: Variable y_sea_water_velocity will be rotated from northward_sea_water_velocity
14:03:57 DEBUG opendrift.readers.basereader:176: Variable mapping: ['sea_floor_depth_below_sea_level'] -> ['land_binary_mask'] is not activated
14:03:57 DEBUG opendrift.readers.basereader.variables:548: Adding variable mapping: ['x_wind', 'y_wind'] -> wind_speed
14:03:57 DEBUG opendrift.readers.basereader.variables:548: Adding variable mapping: ['sea_ice_x_velocity', 'sea_ice_y_velocity'] -> sea_ice_speed
14:03:57 DEBUG opendrift.readers.basereader.variables:548: Adding variable mapping: ['x_sea_water_velocity', 'y_sea_water_velocity'] -> sea_water_speed
===========================
Reader: https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
Projection:
+proj=lcc +lat_0=77.5 +lon_0=-25 +lat_1=77.5 +lat_2=77.5 +x_0=0 +y_0=0 +R=6371000 +units=m +no_defs +type=crs
Coverage: [degrees]
xmin: 278603.187500 xmax: 2123602.750000 step: 2500 numx: 738
ymin: -897931.625000 ymax: 1472067.625000 step: 2500 numy: 948
Corners (lon, lat):
( 79.56, 87.57) ( 68.83, 71.11)
(-17.96, 69.30) ( 18.62, 62.13)
Vertical levels [m]:
[ -0. -3. -10. -15. -25. -50. -75. -100. -150. -200.
-250. -300. -500. -1000. -2000. -3000.]
Available time range:
start: 2022-06-29 00:00:00 end: 2023-09-03 18:00:00 step: 1:00:00
10363 times (408 missing)
Variables:
sea_floor_depth_below_sea_level
latitude
longitude
sea_binary_mask
ocean_vertical_diffusivity
x_wind
y_wind
sea_ice_area_fraction
sea_ice_thickness
sea_ice_x_velocity
sea_ice_y_velocity
sea_water_salinity
sea_water_temperature
eastward_sea_water_velocity
northward_sea_water_velocity
upward_sea_water_velocity
sea_surface_elevation
x_sea_water_velocity
y_sea_water_velocity
wind_speed - derived from ['x_wind', 'y_wind']
sea_ice_speed - derived from ['sea_ice_x_velocity', 'sea_ice_y_velocity']
sea_water_speed - derived from ['x_sea_water_velocity', 'y_sea_water_velocity']
===========================
Time to open reader: 0:00:01.596576
#%
https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml
14:03:57 INFO opendrift.readers.reader_netCDF_CF_generic:92: Opening dataset: https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml
14:04:02 DEBUG opendrift.readers.reader_netCDF_CF_generic:109: Finding coordinate variables.
14:04:02 DEBUG opendrift.readers.reader_netCDF_CF_generic:122: Parsing CF grid mapping dictionary: {'grid_mapping_name': 'polar_stereographic', 'scale_factor_at_projection_origin': 1.0, 'straight_vertical_longitude_from_pole': -45.0, 'latitude_of_projection_origin': 90.0, 'earth_radius': 6371000.0, 'proj4': '+proj=stere +lon_0=-45 +lat_0=90 +k=1 +R=6371000 +no_defs'}
/opt/conda/envs/opendrift/lib/python3.11/site-packages/pyproj/crs/crs.py:1286: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems
proj = self._crs.to_proj4(version=version)
14:04:03 DEBUG opendrift.readers.basereader.variables:593: Setting buffer size 8 for reader https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml
, assuming a maximum average speed of 5 m/s and time span of 1:00:00
14:04:03 DEBUG opendrift.readers.basereader:176: Variable mapping: ['sea_floor_depth_below_sea_level'] -> ['land_binary_mask'] is not activated
14:04:03 DEBUG opendrift.readers.basereader.variables:548: Adding variable mapping: ['sea_water_speed', 'sea_water_to_direction'] -> x_sea_water_velocity
14:04:03 DEBUG opendrift.readers.basereader.variables:548: Adding variable mapping: ['sea_water_speed', 'sea_water_to_direction'] -> y_sea_water_velocity
14:04:03 DEBUG opendrift.readers.basereader.variables:548: Adding variable mapping: ['sea_surface_wave_stokes_drift_x_velocity', 'sea_surface_wave_stokes_drift_y_velocity'] -> sea_surface_wave_stokes_drift_speed
===========================
Reader: https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml
Projection:
+proj=stere +lat_0=90 +lon_0=-45 +k=1 +x_0=0 +y_0=0 +R=6371000 +units=m +no_defs +type=crs
Coverage: [degrees]
xmin: -3600000.000000 xmax: 3798000.000000 step: 3000 numx: 2467
ymin: -4300000.000000 ymax: 2798000.000000 step: 3000 numy: 2367
Corners (lon, lat):
(-172.86, 50.62) ( 81.38, 49.37)
(-84.94, 42.49) ( -3.55, 41.52)
Vertical levels [m]:
Not specified
Available time range:
start: 2017-11-01 00:00:00 end: 2023-09-10 12:00:00 step: 1:00:00
51349 times (0 missing)
Variables:
sea_floor_depth_below_sea_level
latitude
longitude
sea_water_speed
sea_water_to_direction
sea_ice_area_fraction
sea_ice_thickness
sea_surface_wave_significant_height
sea_surface_primary_swell_wave_significant_height
sea_surface_secondary_swell_wave_significant_height
sea_surface_wind_wave_significant_height
sea_surface_wave_from_direction
sea_surface_primary_swell_wave_from_direction
sea_surface_secondary_swell_wave_from_direction
sea_surface_wind_wave_from_direction
sea_surface_wave_from_direction_at_variance_spectral_density_maximum
sea_surface_wave_stokes_drift_x_velocity
sea_surface_wave_stokes_drift_y_velocity
sea_surface_primary_swell_wave_mean_period
sea_surface_secondary_swell_wave_mean_period
sea_surface_wind_wave_mean_period
sea_surface_wave_mean_period_from_variance_spectral_density_second_frequency_moment
sea_surface_wave_mean_period_from_variance_spectral_density_inverse_frequency_moment
sea_surface_wave_period_at_variance_spectral_density_maximum
x_sea_water_velocity - derived from ['sea_water_speed', 'sea_water_to_direction']
y_sea_water_velocity - derived from ['sea_water_speed', 'sea_water_to_direction']
sea_surface_wave_stokes_drift_speed - derived from ['sea_surface_wave_stokes_drift_x_velocity', 'sea_surface_wave_stokes_drift_y_velocity']
===========================
Time to open reader: 0:00:05.479299
#%
https://thredds.met.no/thredds/dodsC/cmems/topaz6/dataset-topaz6-arc-15min-3km-be.ncml
14:04:03 INFO opendrift.readers.reader_netCDF_CF_generic:92: Opening dataset: https://thredds.met.no/thredds/dodsC/cmems/topaz6/dataset-topaz6-arc-15min-3km-be.ncml
14:04:06 DEBUG opendrift.readers.reader_netCDF_CF_generic:109: Finding coordinate variables.
14:04:06 DEBUG opendrift.readers.reader_netCDF_CF_generic:122: Parsing CF grid mapping dictionary: {'grid_mapping_name': 'polar_stereographic', 'latitude_of_projection_origin': 90.0, 'longitude_of_projection_origin': -45.0, 'scale_factor_at_projection_origin': 1.0, 'straight_vertical_longitude_from_pole': -45.0, 'earth_redius': 6378273.0, 'proj4': '+proj=stere +lon_0=-45 +lat_0=90 +k=1 +R=6378273 +no_defs', 'false_easting': 0.0, 'false_northing': 0.0}
/opt/conda/envs/opendrift/lib/python3.11/site-packages/pyproj/crs/crs.py:1286: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems
proj = self._crs.to_proj4(version=version)
14:04:07 DEBUG opendrift.readers.basereader.variables:593: Setting buffer size 4 for reader https://thredds.met.no/thredds/dodsC/cmems/topaz6/dataset-topaz6-arc-15min-3km-be.ncml
, assuming a maximum average speed of 5 m/s and time span of 0:15:00
14:04:07 DEBUG opendrift.readers.basereader:176: Variable mapping: ['sea_floor_depth_below_sea_level'] -> ['land_binary_mask'] is not activated
14:04:07 DEBUG opendrift.readers.basereader.variables:548: Adding variable mapping: ['x_sea_water_velocity', 'y_sea_water_velocity'] -> sea_water_speed
===========================
Reader: https://thredds.met.no/thredds/dodsC/cmems/topaz6/dataset-topaz6-arc-15min-3km-be.ncml
Projection:
+proj=stere +lat_0=90 +lon_0=-45 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs +type=crs
Coverage: [degrees]
xmin: -3600000.000000 xmax: 3797999.954224 step: 2999.88 numx: 2467
ymin: -4300000.000000 ymax: 2797999.954224 step: 2999.88 numy: 2367
Corners (lon, lat):
(-172.86, 50.73) ( 81.38, 49.48)
(-84.94, 42.59) ( -3.55, 41.62)
Vertical levels [m]:
Not specified
Available time range:
start: 2018-01-01 00:00:00 end: 2023-09-09 23:45:00 step: 0:15:00
199488 times (0 missing)
Variables:
longitude
latitude
sea_floor_depth_below_sea_level
x_sea_water_velocity
y_sea_water_velocity
sea_surface_height_above_geoid
sea_water_speed - derived from ['x_sea_water_velocity', 'y_sea_water_velocity']
===========================
Time to open reader: 0:00:04.472325
#%
https://thredds.met.no/thredds/dodsC/aromearcticlatest/latest/arome_arctic_lagged_12_h_latest_2_5km_latest.nc
14:04:07 INFO opendrift.readers.reader_netCDF_CF_generic:92: Opening dataset: https://thredds.met.no/thredds/dodsC/aromearcticlatest/latest/arome_arctic_lagged_12_h_latest_2_5km_latest.nc
14:04:08 DEBUG opendrift.readers.reader_netCDF_CF_generic:109: Finding coordinate variables.
14:04:08 DEBUG opendrift.readers.reader_netCDF_CF_generic:211: 4 ensemble members available
14:04:08 DEBUG opendrift.readers.reader_netCDF_CF_generic:122: Parsing CF grid mapping dictionary: {'grid_mapping_name': 'lambert_conformal_conic', 'standard_parallel': array([77.5, 77.5]), 'longitude_of_central_meridian': -25.0, 'latitude_of_projection_origin': 77.5, 'earth_radius': 6371000.0}
/opt/conda/envs/opendrift/lib/python3.11/site-packages/pyproj/crs/crs.py:1286: UserWarning: You will likely lose important projection information when converting to a PROJ string from another format. See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems
proj = self._crs.to_proj4(version=version)
14:04:08 DEBUG opendrift.readers.reader_netCDF_CF_generic:310: Skipped variables without standard_name: ['fog_area_fraction']
14:04:08 DEBUG opendrift.readers.basereader.variables:593: Setting buffer size 10 for reader https://thredds.met.no/thredds/dodsC/aromearcticlatest/latest/arome_arctic_lagged_12_h_latest_2_5km_latest.nc
, assuming a maximum average speed of 5 m/s and time span of 1:00:00
14:04:08 DEBUG opendrift.readers.basereader:176: Variable mapping: ['sea_floor_depth_below_sea_level'] -> ['land_binary_mask'] is not activated
14:04:08 DEBUG opendrift.readers.basereader.variables:548: Adding variable mapping: ['x_wind', 'y_wind'] -> wind_speed
===========================
Reader: https://thredds.met.no/thredds/dodsC/aromearcticlatest/latest/arome_arctic_lagged_12_h_latest_2_5km_latest.nc
Projection:
+proj=lcc +lat_0=77.5 +lon_0=-25 +lat_1=77.5 +lat_2=77.5 +x_0=0 +y_0=0 +R=6371000 +units=m +no_defs +type=crs
Coverage: [degrees]
xmin: 278620.937500 xmax: 2123621.000000 step: 2500 numx: 739
ymin: -897985.687500 ymax: 1472014.250000 step: 2500 numy: 949
Corners (lon, lat):
( 79.55, 87.57) ( 68.83, 71.11)
(-17.96, 69.30) ( 18.62, 62.13)
Vertical levels [m]:
Not specified
Available time range:
start: 2023-09-01 09:00:00 end: 2023-09-03 18:00:00 step: 1:00:00
58 times (0 missing)
Variables (4 ensemble members):
longitude
latitude
x_wind
y_wind
air_temperature
relative_humidity
cloud_area_fraction
air_pressure_at_sea_level
precipitation_amount
snowfall_amount
wind_speed_of_gust
wind_speed - derived from ['x_wind', 'y_wind']
===========================
Time to open reader: 0:00:00.984705
#%
https://pae-paha.pacioos.hawaii.edu/thredds/dodsC/ncep_global/NCEP_Global_Atmospheric_Model_best.ncd
14:04:08 INFO opendrift.readers.reader_netCDF_CF_generic:92: Opening dataset: https://pae-paha.pacioos.hawaii.edu/thredds/dodsC/ncep_global/NCEP_Global_Atmospheric_Model_best.ncd
14:04:09 DEBUG opendrift.readers.reader_netCDF_CF_generic:109: Finding coordinate variables.
14:04:09 DEBUG opendrift.readers.reader_netCDF_CF_generic:234: Lon and lat are 1D arrays - using as projection coordinates
14:04:09 DEBUG opendrift.readers.basereader.variables:593: Setting buffer size 3 for reader https://pae-paha.pacioos.hawaii.edu/thredds/dodsC/ncep_global/NCEP_Global_Atmospheric_Model_best.ncd
, assuming a maximum average speed of 5 m/s and time span of 3:00:00
14:04:09 INFO opendrift.readers.basereader:166: Variable x_wind will be rotated from eastward_wind
14:04:09 INFO opendrift.readers.basereader:166: Variable y_wind will be rotated from northward_wind
14:04:09 DEBUG opendrift.readers.basereader:176: Variable mapping: ['sea_floor_depth_below_sea_level'] -> ['land_binary_mask'] is not activated
14:04:09 DEBUG opendrift.readers.basereader.variables:548: Adding variable mapping: ['x_wind', 'y_wind'] -> wind_speed
===========================
Reader: https://pae-paha.pacioos.hawaii.edu/thredds/dodsC/ncep_global/NCEP_Global_Atmospheric_Model_best.ncd
Projection:
+proj=latlong
Coverage: [degrees]
xmin: 0.000000 xmax: 359.500000 step: 0.5 numx: 720
ymin: -90.000000 ymax: 90.000000 step: 0.5 numy: 361
Corners (lon, lat):
( 0.00, 90.00) (359.50, 90.00)
( 0.00, -90.00) (359.50, -90.00)
Vertical levels [m]:
Not specified
Available time range:
start: 2022-12-01 12:00:00 end: 2023-09-08 12:00:00 step: 3:00:00
2249 times (0 missing)
Variables:
air_pressure_at_sea_level
surface_downwelling_longwave_flux
surface_downwelling_shortwave_flux
precipitation_flux
relative_humidity
surface_temperature
air_temperature
eastward_wind
surface_upwelling_longwave_flux
surface_upwelling_shortwave_flux
northward_wind
x_wind
y_wind
wind_speed - derived from ['x_wind', 'y_wind']
===========================
Time to open reader: 0:00:01.212201
#%
https://tds.hycom.org/thredds/dodsC/GLBy0.08/latest
14:04:09 INFO opendrift.readers.reader_netCDF_CF_generic:92: Opening dataset: https://tds.hycom.org/thredds/dodsC/GLBy0.08/latest
14:04:10 DEBUG opendrift.readers.reader_netCDF_CF_generic:109: Finding coordinate variables.
14:04:10 INFO opendrift.readers.reader_netCDF_CF_generic:253: Grid coordinates are detected, but proj4 string not given: assuming latlong
14:04:10 DEBUG opendrift.readers.basereader.variables:593: Setting buffer size 9 for reader https://tds.hycom.org/thredds/dodsC/GLBy0.08/latest
, assuming a maximum average speed of 5 m/s and time span of 3:00:00
14:04:10 INFO opendrift.readers.basereader:166: Variable x_sea_water_velocity will be rotated from eastward_sea_water_velocity
14:04:10 INFO opendrift.readers.basereader:166: Variable y_sea_water_velocity will be rotated from northward_sea_water_velocity
14:04:10 DEBUG opendrift.readers.basereader:176: Variable mapping: ['sea_floor_depth_below_sea_level'] -> ['land_binary_mask'] is not activated
14:04:10 DEBUG opendrift.readers.basereader.variables:548: Adding variable mapping: ['x_sea_water_velocity', 'y_sea_water_velocity'] -> sea_water_speed
===========================
Reader: https://tds.hycom.org/thredds/dodsC/GLBy0.08/latest
Projection:
+proj=latlong
Coverage: [degrees]
xmin: 0.000000 xmax: 359.920044 step: 0.0799561 numx: 4502
ymin: -80.000000 ymax: 90.000000 step: 0.0400009 numy: 4250
Corners (lon, lat):
( 0.00, 90.00) (359.92, 90.00)
( 0.00, -80.00) (359.92, -80.00)
Vertical levels [m]:
[ -0. -2. -4. -6. -8. -10. -12. -15. -20. -25.
-30. -35. -40. -45. -50. -60. -70. -80. -90. -100.
-125. -150. -200. -250. -300. -350. -400. -500. -600. -700.
-800. -900. -1000. -1250. -1500. -2000. -2500. -3000. -4000. -5000.]
Available time range:
start: 2023-08-26 12:00:00 end: 2023-09-08 00:00:00 step: 3:00:00
101 times (0 missing)
Variables:
sea_surface_elevation
eastward_sea_water_velocity
eastward_sea_water_velocity_at_bottom
northward_sea_water_velocity
northward_sea_water_velocity_at_bottom
sea_water_temperature
sea_water_temperature_at_bottom
sea_water_salinity
sea_water_salinity_at_bottom
x_sea_water_velocity
y_sea_water_velocity
sea_water_speed - derived from ['x_sea_water_velocity', 'y_sea_water_velocity']
===========================
Time to open reader: 0:00:00.429052
#%
https://tds.hycom.org/thredds/dodsC/GLBy0.08/expt_93.0/uv3z
14:04:10 INFO opendrift.readers.reader_netCDF_CF_generic:92: Opening dataset: https://tds.hycom.org/thredds/dodsC/GLBy0.08/expt_93.0/uv3z
14:10:15 DEBUG opendrift.readers.reader_netCDF_CF_generic:109: Finding coordinate variables.
14:10:15 INFO opendrift.readers.reader_netCDF_CF_generic:253: Grid coordinates are detected, but proj4 string not given: assuming latlong
14:10:15 DEBUG opendrift.readers.basereader.variables:593: Setting buffer size 9 for reader https://tds.hycom.org/thredds/dodsC/GLBy0.08/expt_93.0/uv3z
, assuming a maximum average speed of 5 m/s and time span of 3:00:00
14:10:15 INFO opendrift.readers.basereader:166: Variable x_sea_water_velocity will be rotated from eastward_sea_water_velocity
14:10:15 INFO opendrift.readers.basereader:166: Variable y_sea_water_velocity will be rotated from northward_sea_water_velocity
14:10:15 DEBUG opendrift.readers.basereader:176: Variable mapping: ['sea_floor_depth_below_sea_level'] -> ['land_binary_mask'] is not activated
14:10:15 DEBUG opendrift.readers.basereader.variables:548: Adding variable mapping: ['x_sea_water_velocity', 'y_sea_water_velocity'] -> sea_water_speed
===========================
Reader: https://tds.hycom.org/thredds/dodsC/GLBy0.08/expt_93.0/uv3z
Projection:
+proj=latlong
Coverage: [degrees]
xmin: 0.000000 xmax: 359.920044 step: 0.0799561 numx: 4502
ymin: -80.000000 ymax: 90.000000 step: 0.0400009 numy: 4250
Corners (lon, lat):
( 0.00, 90.00) (359.92, 90.00)
( 0.00, -80.00) (359.92, -80.00)
Vertical levels [m]:
[ -0. -2. -4. -6. -8. -10. -12. -15. -20. -25.
-30. -35. -40. -45. -50. -60. -70. -80. -90. -100.
-125. -150. -200. -250. -300. -350. -400. -500. -600. -700.
-800. -900. -1000. -1250. -1500. -2000. -2500. -3000. -4000. -5000.]
Available time range:
start: 2018-12-04 12:00:00 end: 2023-09-01 09:00:00 step: 3:00:00
13856 times (7 missing)
Variables:
eastward_sea_water_velocity
eastward_sea_water_velocity_at_bottom
northward_sea_water_velocity
northward_sea_water_velocity_at_bottom
x_sea_water_velocity
y_sea_water_velocity
sea_water_speed - derived from ['x_sea_water_velocity', 'y_sea_water_velocity']
===========================
Time to open reader: 0:06:05.038971
Summary of times to open each dataset:
for t, time in times.items():
print(time, t)
0:00:01.676639 https://thredds.met.no/thredds/dodsC/sea/norkyst800m/1h/aggregate_be
0:00:00.941629 https://thredds.met.no/thredds/dodsC/mepslatest/meps_lagged_6_h_latest_2_5km_latest.nc
0:00:01.596576 https://thredds.met.no/thredds/dodsC/fou-hi/barents_eps_zdepth_be
0:00:05.479299 https://thredds.met.no/thredds/dodsC/cmems/mywavewam3km/dataset-wam-arctic-1hr3km-be.ncml
0:00:04.472325 https://thredds.met.no/thredds/dodsC/cmems/topaz6/dataset-topaz6-arc-15min-3km-be.ncml
0:00:00.984705 https://thredds.met.no/thredds/dodsC/aromearcticlatest/latest/arome_arctic_lagged_12_h_latest_2_5km_latest.nc
0:00:01.212201 https://pae-paha.pacioos.hawaii.edu/thredds/dodsC/ncep_global/NCEP_Global_Atmospheric_Model_best.ncd
0:00:00.429052 https://tds.hycom.org/thredds/dodsC/GLBy0.08/latest
0:06:05.038971 https://tds.hycom.org/thredds/dodsC/GLBy0.08/expt_93.0/uv3z
Total running time of the script: (7 minutes 59.210 seconds)