Note
Click here to download the full example code
Example on integrating radar and HRRR dataΒΆ
This is an example of how to retrieve winds in Hurricane Florence. In this example, we use data from 2 NEXRAD radars as well as from the HRRR to retrieve the winds
Grided netCDF files are downloadable from: https://drive.google.com/drive/u/1/folders/1pcQxWRJV78xuJePTZnlXPPpMe1qut0ie

Author: Robert C. Jackson
import urllib
import pyart
import pydda
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import numpy as np
hrrr_url = ('https://pando-rgw01.chpc.utah.edu/hrrr/prs/20180914/' +
'hrrr.t06z.wrfprsf00.grib2')
urllib.request.urlretrieve(hrrr_url, 'test.grib2')
grid_mhx = pyart.io.read_grid(u'grid_mhx.nc')
grid_ltx = pyart.io.read_grid(u'grid_ltx.nc')
grid_mhx = pydda.constraints.add_hrrr_constraint_to_grid(grid_mhx,
u'test.grib2')
u_init, v_init, w_init = pydda.initialization.make_constant_wind_field(
grid_mhx, (0.0, 0.0, 0.0))
out_grids = pydda.retrieval.get_dd_wind_field(
[grid_mhx, grid_ltx], u_init, v_init, w_init, Co=0.1, Cm=1000.0, Cmod=1e-3,
mask_outside_opt=True, vel_name='corrected_velocity',
model_fields=["hrrr"])
fig = plt.figure(figsize=(25, 15))
ax = plt.axes(projection=ccrs.PlateCarree())
ax = pydda.vis.plot_horiz_xsection_barbs_map(
out_grids, ax=ax, bg_grid_no=-1, level=1, barb_spacing_x_km=20.0,
barb_spacing_y_km=20.0, cmap='pyart_HomeyerRainbow')
ax.set_xticks(np.arange(-80, -75, 0.5))
ax.set_yticks(np.arange(33., 35.5, 0.5))
plt.title(out_grids[0].time['units'][13:] + ' winds at 0.5 km')
plt.show()
Total running time of the script: ( 0 minutes 0.000 seconds)