- wradlib.qual.beam_block_frac(th, bh, a)¶
Partial beam blockage fraction.
Code was migrated from https://github.com/nguy/PyRadarMet.
From Bech et al. (2003), Eqn 2 and Appendix
th (float |
numpy:numpy.ndarrayof floats) – Terrain height [m]
bh (float |
numpy:numpy.ndarrayof floats) – Beam height [m]
a (float |
numpy:numpy.ndarrayof floats) – Half power beam radius [m]
pbb (float) – Partial beam blockage fraction [unitless]
>>> pbb = beam_block_frac(th,bh,a)
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!