docs: def f_solpos(times, latitude, longitude):
      """
      Calculate solar position for lat/long at times.
  
      :param times: Python :class:`datetime.datetime` object.
      :type times: list
      :param latitude: latitude [degrees]
      :type latitude: float
      :param longitude: longitude [degrees]
      :type longitude: float
      :returns: apparent zenith, azimuth
      """
      # pvlib converts Python datetime objects to pandas DatetimeIndex
      solpos = pvlib.solarposition.get_solarposition(times, latitude, longitude)
      # solpos is a pandas DataFrame, so unpack the desired values
      # return shape is (2, NOBS), so unc_wrapper sees 2 dependent variables
      return solpos['apparent_zenith'].values, solpos['azimuth'].values