pyart.retrieve.est_wind_vel#

pyart.retrieve.est_wind_vel(radar, vert_proj=False, vel_field=None, wind_field=None)[source]#

Estimates wind velocity. Projects the radial wind component to the horizontal or vertical of the azimuth plane. It assumes that the orthogonal component is negligible.

The horizontal wind component is given by:

v = v_r*cos(el)-v_el*sin(el)+v_az

where:

v_r is the radial wind component (measured by the radar) v_el is the perpendicular wind component in the azimuth plane. v_az is the horizontal component perpendicular to the radial direction and the azimuth plane el is the elevation

The horizontal wind component in the azimuth plane is given by:

v_h = v_r*cos(el) - v_el*sin(el)

which since we do not know v_el we assume:

v_h ~ v_r*cos(el)

This assumption holds for small elevation angles

The vertical wind component in the azimuth plane is given by:

v_h = v_r*sin(el) - v_el*cos(el)

which since we do not know v_el we assume:

v_h ~ v_r*sin(el)

This assumption holds for angles close to 90 deg

Parameters:

  • radar (Radar) – Radar object

  • vert_proj (Boolean) – If true estimates the vertical projection, otherwise the horizontal

  • vel_field (str) – name of the velocity field

  • wind_field (str) – name of the velocity field

Returns:

wind (dict) – Field dictionary containing the estimated wind velocity