Rice-Holmberg rain rate distribution estimation model

calculation of a rain rate cumulative distribution function following the RH approach

%function [P R]=RH_distrib(M,Mm,U)
%Function implementing the Rice Holmberg model for the estimation of a rain rate cumulative distribution function
%
%Inputs
%M: SCALAR. average annual rain accumulation, in mm
%Mm: SCALAR highest monthly accumulation in a period of 30 consecutive years, in
%mm
%U: SCALAR average number of thunderstorms in an average year
%override: boolean. 1 if Beta factor override (allows input of beta
%directly, if it is available)
%OUTPUTS
%P: vector of percentages of exceedance.
%R : vector of rain rate values exceeded during P percent of the time
%{P,R} are a cumulative distribution function useful for telecom system
%design at frequencies above 7 GHz in tropical regions, or in general, for
%systems affected by rain attenuation
%
%Example of usage: [P R]=RH_distrib(1669.2,400.1,125)
%will deliver the values for Medellin, Colombia in table 7 on the paper
%Emiliani, L., Agudelo, J., Gutierrez, E., Restrepo, J., Fradique-Mendez, C.
%“Development of rain-attenuation and rain-rate maps for satellite system design
%in the Ku and Ka bands in Colombia”.
%Antennas and Propagation Magazine, IEEE. Vol 46,
%Issue 6, Dec 2004 Page(s):54 - 68
%
%approach
%This mfile iterates over a sufficiently large vector of R (0.000001 to 400
%mm/hr) and finds P as per the model, then interpolates to standard ITU P
%values.