- wradlib.qual.beam_block_frac(th, bh, a)[source]#
Partial beam blockage fraction.
Code was migrated from nguy/PyRadarMet.
From Bech et al. (2003), Eqn 2 and Appendix
numpy.ndarray) – Terrain height [m]
numpy.ndarray) – Beam height [m]
numpy.ndarray) – Half power beam radius [m]
numpy.ndarray) – Partial beam blockage fraction [unitless]
>>> pbb = beam_block_frac(th,bh,a)
See Beam Blockage Calculation using a DEM.
This procedure uses a simplified interception function where no vertical gradient of refractivity is considered. Other algorithms treat this more thoroughly. However, this is accurate in most cases other than the super-refractive case.
See the the half_power_radius function to calculate variable a.
The heights must be the same units!