code: def sun_zenith_angle(utc_time, lon, lat):
      
      return np.rad2deg(np.arccos(cos_zen(utc_time, lon, lat)))
docs: Sun-zenith angle for *lon*, *lat* at *utc_time*.
      lon,lat in degrees.
      The angle returned is given in degrees