Hello,
I am trying to plot VAD and VVP for the radar data. I couldn’t find the VVP and VAD retrieval Techniques using wradlib. Any help ?
A notebook will be really helpful
Hello,
I am trying to plot VAD and VVP for the radar data. I couldn’t find the VVP and VAD retrieval Techniques using wradlib. Any help ?
A notebook will be really helpful
@adithiyrnair Velocity Azimuth Display (VAD) and volume velocity processing (VVP) is currently not implemented in wradlib. VAD algorithm(s) are implemented in Py-ART, not sure if there is a full featured VVP, too. Maybe @mgrover1 or @rcjackson can give more guidance here.
Thanks for the information
Kai
We do have VADs in PyART, but not VVP. I don’t think I know of a package that uses VVP to get 2D winds that’s written in Python.
PyDDA doesn’t have VVD, but it can generate winds from a single doppler radar using the weak variational technique. I would just use u and v from the results you get from PyDDA.
I believe you’re referring to something like the below attached VVP. Also, IMD uses one hour radar data for plotting vertical wind profile .
I think it should be something like this, but I am doubtful about the barbs direction .
import matplotlib.pyplot as plt
import numpy as np
import pyart
radar = pyart.io.read_cfradial("goa_data/cfradial/polar_GOA210516024101-IMD-B.nc")
# Loop on all sweeps and compute VAD
zlevels = np.arange(100, 5000, 100) # height above radar
u_allsweeps = []
v_allsweeps = []
for idx in range(radar.nsweeps):
radar_1sweep = radar.extract_sweeps([idx])
vad = pyart.retrieve.vad_browning(
radar_1sweep, "VELH", z_want=zlevels
)
u_allsweeps.append(vad.u_wind)
v_allsweeps.append(vad.v_wind)
# Average U and V over all sweeps and compute magnitude and angle
u_avg = np.nanmean(np.array(u_allsweeps), axis=0)
v_avg = np.nanmean(np.array(v_allsweeps), axis=0)
orientation = np.rad2deg(np.arctan2(-u_avg, -v_avg)) % 360
speed = np.sqrt(u_avg**2 + v_avg**2)
# Display vertical profile of wind
fig, ax = plt.subplots(1, 3, sharey=True)
ax[0].plot(speed * 2, zlevels + radar.altitude["data"])
ax[1].plot(orientation, zlevels + radar.altitude["data"])
ax[0].set_xlabel("Wind speed [m/s]")
ax[1].set_xlabel("Wind direction [deg]")
ax[0].set_ylabel("Altitude [m]")
fig.suptitle("Wind profile obtained from VAD")
u = np.ma.masked_array(u_avg, mask=False)
v = np.ma.masked_array(v_avg, mask=False)
# this line should no longer cause an exception
nbarbs = (~u.mask).sum()
skip = max(1, int(nbarbs//5))
bloc = 0.5
ax[2].barbs((np.zeros(zlevels.shape)+bloc)[::skip]-0.5, zlevels[::skip]+radar.altitude['data'],
-u[::skip], -v[::skip],)
ax[2].axis("off")