Exercice

The goal of this exercice is to do a QPE for every sweep and perform a weighted average of precipitation aloft to get an estimation of the precipitation intensity at the ground. The following steps will need to be taken:

Load all radar files in /data/question_pyart_meteoswiss and merge them into one single radar object
Perform attenuation correction of ZH, you can use a constant freezing level height of 2700 m.
Estimate the QPE with a a polynomial Z-R relation.
Compute a CAPPI of the resulting radar estimate rain rate from 500 to 8000 m above the radar using a vertical resolution of 100 m and a horizontal resolution of 500 m.
Using numpy, perform a weighted average of all CAPPI levels using the weights

$w(z) = exp(-0.5\cdot z)$, where $z$ is the height above the radar in km.

Finally display the resulting QPE on its latitude/longitude grid using pcolormesh (matplotlib).

The solution can be found here