def calc_rain_z_zdr(dz, zdr, a=6.7e-3, b=0.927, c=-0.343):
"""
dz = Reflectivity (dBZ)
zdr = Differential Reflectivity (dB)
a, b, c = Adjustable function parameters
"""
return a * linearize(dz)**b * 10.0**(c * zdr)
def calc_rain_zr(dz, a=300.0, b=1.4):
"""
dz = Reflectivity (dBZ), returns rainfall rate (mm h^-1)
Form Z = a * R**b
"""
return (linearize(dz) / a)**(1.0 / b)