diff --git a/pyspectral/rayleigh.py b/pyspectral/rayleigh.py
index <HASH>..<HASH> 100644
--- a/pyspectral/rayleigh.py
+++ b/pyspectral/rayleigh.py
@@ -174,6 +174,8 @@ class Rayleigh(object):
         clip_angle = np.rad2deg(np.arccos(1. / 25))
         sun_zenith = np.clip(np.asarray(sun_zenith), 0, clip_angle)
         sunzsec = 1. / np.cos(np.deg2rad(sun_zenith))
+        clip_angle = np.rad2deg(np.arccos(1. / 3.))
+        sat_zenith = np.clip(np.asarray(sat_zenith), 0, clip_angle)
         satzsec = 1. / np.cos(np.deg2rad(np.asarray(sat_zenith)))
 
         shape = sun_zenith.shape
Clip satellite-zenith angles outside range