def edr_short_range(sw, rng, theta, gs):
"""
For gate spacing > range * beamwidth
sw (spectrum width) in m/s,
rng (range) in km,
theta (beamwidth) in deg,
gs (gate spacing) in km
"""
beta = gs * RNG_MULT / RRV_SCALING_FACTOR
alpha = rng * RNG_MULT * np.deg2rad(theta) / RRV_SCALING_FACTOR
z = 1.0 - (alpha**2 / beta**2) # alpha always <= beta
series = hypergeometric_gaussian(-0.3333, 2.0, 2.5, z)
edr = sw**3 * (1.0/beta) * (CONSTANT * series)**(-1.5)
return edr**0.33333
def edr_short_range(sw, rng, theta, gs):
"""
For gate spacing > range * beamwidth
sw (spectrum width) in m/s,
rng (range) in km,
theta (beamwidth) in deg,
gs (gate spacing) in km
"""
beta = gs * RNG_MULT / RRV_SCALING_FACTOR
alpha = rng * RNG_MULT * np.deg2rad(theta) / RRV_SCALING_FACTOR
z = 1.0 - (alpha**2 / beta**2) # alpha always <= beta
series = hypergeometric_gaussian(-0.3333, 2.0, 2.5, z)
edr = sw**3 * (1.0/beta) * (CONSTANT * series)**(-1.5)
return edr**0.33333