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.

1 Like

Thanks for the information
Kai

1 Like

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.

1 Like

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

# Loop on all sweeps and compute VAD
zlevels = np.arange(100, 5000, 100)  # height above radar
u_allsweeps = []
v_allsweeps = []

)

# 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.plot(speed * 2, zlevels + radar.altitude["data"])
ax.set_xlabel("Wind speed [m/s]")
ax.set_xlabel("Wind direction [deg]")
ax.set_ylabel("Altitude [m]")
# this line should no longer cause an exception This is it thank you  