pvlib 0.9.4a1__py3-none-any.whl → 0.10.0__py3-none-any.whl

pvlib/ivtools/sde.py

@@ -10,16 +10,6 @@ import numpy as np
 from pvlib.ivtools.utils import _schumaker_qspline
 
 
-# set constant for numpy.linalg.lstsq parameter rcond
-# rcond=-1 for numpy<1.14, rcond=None for numpy>=1.14
-# TODO remove after minimum numpy version >= 1.14
-minor = int(np.__version__.split('.')[1])
-if minor < 14:
-    RCOND = -1
-else:
-    RCOND = None
-
-
 def fit_sandia_simple(voltage, current, v_oc=None, i_sc=None, v_mp_i_mp=None,
                       vlim=0.2, ilim=0.1):
     r"""
@@ -74,7 +64,8 @@ def fit_sandia_simple(voltage, current, v_oc=None, i_sc=None, v_mp_i_mp=None,
 
     Raises
     ------
-    RuntimeError if parameter extraction is not successful.
+    RuntimeError
+        if parameter extraction is not successful.
 
     Notes
     -----
@@ -116,7 +107,7 @@ def fit_sandia_simple(voltage, current, v_oc=None, i_sc=None, v_mp_i_mp=None,
     .. math::
 
         I &\approx \frac{I_{L}}{1 + G_{p} R_{s}}
-        - \frac{G_{p}}{1 + G_{p}R_{s}} V
+        - \frac{G_{p}}{1 + G_{p}R_{s}} V \\
         &= \beta_{0} + \beta_{1} V
 
     4. The exponential portion of the IV curve is defined by
@@ -134,8 +125,8 @@ def fit_sandia_simple(voltage, current, v_oc=None, i_sc=None, v_mp_i_mp=None,
     .. math::
 
         \log(\beta_{0} - \beta_{1} V - I)
-        &\approx \log(\frac{I_{0}}{1 + G_{p} R_{s}} + \frac{V}{nN_sV_{th}}
-        + \frac{I R_{s}}{nN_sV_{th}}) \\
+        &\approx \log(\frac{I_{0}}{1 + G_{p} R_{s}}) + \frac{V}{nN_sV_{th}}
+        + \frac{I R_{s}}{nN_sV_{th}} \\
         &= \beta_{2} + \beta_{3} V + \beta_{4} I
 
     6. Calculate values for ``IL, I0, Rs, Rsh,`` and ``nNsVth`` from the
@@ -149,7 +140,8 @@ def fit_sandia_simple(voltage, current, v_oc=None, i_sc=None, v_mp_i_mp=None,
        0 86758 909 4
     .. [2] C. B. Jones, C. W. Hansen, "Single Diode Parameter Extraction from
        In-Field Photovoltaic I-V Curves on a Single Board Computer", 46th IEEE
-       Photovoltaic Specialist Conference, Chicago, IL, 2019
+       Photovoltaic Specialist Conference, Chicago, IL, 2019.
+       :doi:`10.1109/PVSC40753.2019.8981330`
     """
 
     # If not provided, extract v_oc, i_sc, v_mp and i_mp from the IV curve data
@@ -228,7 +220,7 @@ def _sandia_beta3_beta4(voltage, current, beta0, beta1, ilim, i_sc):
     x = np.array([np.ones_like(voltage), voltage, current]).T
     # Select points where y > ilim * i_sc to regress log(y) onto x
     idx = (y > ilim * i_sc)
-    result = np.linalg.lstsq(x[idx], np.log(y[idx]), rcond=RCOND)
+    result = np.linalg.lstsq(x[idx], np.log(y[idx]), rcond=None)
     coef = result[0]
     beta3 = coef[1].item()
     beta4 = coef[2].item()
@@ -299,8 +291,10 @@ def _fit_sandia_cocontent(voltage, current, nsvth):
 
     Raises
     ------
-    ValueError if ``voltage`` and ``current`` are different lengths.
-    ValueError if ``len(voltage)`` < 6
+    ValueError
+        if ``voltage`` and ``current`` are different lengths.
+    ValueError
+        if ``len(voltage)`` < 6
 
     Notes
     -----
@@ -444,7 +438,7 @@ def _cocontent_regress(v, i, voc, isc, cci):
     sx = np.vstack((s[:, 0], s[:, 1], s[:, 0] * s[:, 1], s[:, 0] * s[:, 0],
                     s[:, 1] * s[:, 1], col1)).T
 
-    gamma = np.linalg.lstsq(sx, scc, rcond=RCOND)[0]
+    gamma = np.linalg.lstsq(sx, scc, rcond=None)[0]
     # coefficients from regression in rotated coordinates
 
     # Principle components transformation steps
@@ -473,5 +467,5 @@ def _cocontent_regress(v, i, voc, isc, cci):
 
     # translate from coefficients in rotated space (gamma) to coefficients in
     # original coordinates (beta)
-    beta = np.linalg.lstsq(np.dot(mb, ma), gamma[0:5], rcond=RCOND)[0]
+    beta = np.linalg.lstsq(np.dot(mb, ma), gamma[0:5], rcond=None)[0]
     return beta

pvlib/ivtools/sdm.py

@@ -8,10 +8,9 @@ Function names should follow the pattern "fit_" + name of model + "_" +
 
 import numpy as np
 
-import scipy.constants
+from scipy import constants
 from scipy import optimize
 from scipy.special import lambertw
-from scipy.misc import derivative
 
 from pvlib.pvsystem import calcparams_pvsyst, singlediode, v_from_i
 from pvlib.singlediode import bishop88_mpp
@@ -19,6 +18,11 @@ from pvlib.singlediode import bishop88_mpp
 from pvlib.ivtools.utils import rectify_iv_curve, _numdiff
 from pvlib.ivtools.sde import _fit_sandia_cocontent
 
+from pvlib.tools import _first_order_centered_difference
+
+
+CONSTANTS = {'E0': 1000.0, 'T0': 25.0, 'k': constants.k, 'q': constants.e}
+
 
 def fit_cec_sam(celltype, v_mp, i_mp, v_oc, i_sc, alpha_sc, beta_voc,
                 gamma_pmp, cells_in_series, temp_ref=25):
@@ -72,9 +76,10 @@ def fit_cec_sam(celltype, v_mp, i_mp, v_oc, i_sc, alpha_sc, beta_voc,
 
     Raises
     ------
-        ImportError if NREL-PySAM is not installed.
-
-        RuntimeError if parameter extraction is not successful.
+    ImportError
+        if NREL-PySAM is not installed.
+    RuntimeError
+        if parameter extraction is not successful.
 
     Notes
     -----
@@ -90,7 +95,7 @@ def fit_cec_sam(celltype, v_mp, i_mp, v_oc, i_sc, alpha_sc, beta_voc,
     ----------
     .. [1] A. Dobos, "An Improved Coefficient Calculator for the California
        Energy Commission 6 Parameter Photovoltaic Module Model", Journal of
-       Solar Energy Engineering, vol 134, 2012.
+       Solar Energy Engineering, vol 134, 2012. :doi:`10.1115/1.4005759`
     """
 
     try:
@@ -200,11 +205,11 @@ def fit_desoto(v_mp, i_mp, v_oc, i_sc, alpha_sc, beta_voc, cells_in_series,
     ----------
     .. [1] W. De Soto et al., "Improvement and validation of a model for
        photovoltaic array performance", Solar Energy, vol 80, pp. 78-88,
-       2006.
+       2006. :doi:`10.1016/j.solener.2005.06.010`
     """
 
     # Constants
-    k = scipy.constants.value('Boltzmann constant in eV/K')
+    k = constants.value('Boltzmann constant in eV/K')  # in eV/K
     Tref = temp_ref + 273.15  # [K]
 
     # initial guesses of variables for computing convergence:
@@ -340,9 +345,9 @@ def fit_pvsyst_sandia(ivcurves, specs, const=None, maxiter=5, eps1=1.e-3):
         T0 : float
             cell temperature at STC, default 25 [C]
         k : float
-            1.38066E-23 J/K (Boltzmann's constant)
+            Boltzmann's constant [J/K]
         q : float
-            1.60218E-19 Coulomb (elementary charge)
+            elementary charge [Coulomb]
 
     maxiter : int, default 5
         input that sets the maximum number of iterations for the parameter
@@ -401,6 +406,7 @@ def fit_pvsyst_sandia(ivcurves, specs, const=None, maxiter=5, eps1=1.e-3):
     .. [1] K. Sauer, T. Roessler, C. W. Hansen, Modeling the Irradiance and
        Temperature Dependence of Photovoltaic Modules in PVsyst, IEEE Journal
        of Photovoltaics v5(1), January 2015.
+       :doi:`10.1109/JPHOTOV.2014.2364133`
     .. [2] A. Mermoud, PV Modules modeling, Presentation at the 2nd PV
        Performance Modeling Workshop, Santa Clara, CA, May 2013
     .. [3] A. Mermoud, T. Lejeuene, Performance Assessment of a Simulation
@@ -408,16 +414,18 @@ def fit_pvsyst_sandia(ivcurves, specs, const=None, maxiter=5, eps1=1.e-3):
        Photovoltaic Solar Energy Conference, Valencia, Spain, Sept. 2010
     .. [4] C. Hansen, Estimating Parameters for the PVsyst Version 6
        Photovoltaic Module Performance Model, Sandia National Laboratories
-       Report SAND2015-8598
+       Report SAND2015-8598. :doi:`10.2172/1223058`
     .. [5] C. Hansen, Parameter Estimation for Single Diode Models of
-       Photovoltaic Modules, Sandia National Laboratories Report SAND2015-2065
+       Photovoltaic Modules, Sandia National Laboratories Report SAND2015-2065.
+       :doi:`10.2172/1177157`
     .. [6] C. Hansen, Estimation of Parameters for Single Diode Models using
         Measured IV Curves, Proc. of the 39th IEEE PVSC, June 2013.
+        :doi:`10.1109/PVSC.2013.6744135`
     .. [7] PVLib MATLAB https://github.com/sandialabs/MATLAB_PV_LIB
     """
 
     if const is None:
-        const = {'E0': 1000.0, 'T0': 25.0, 'k': 1.38066e-23, 'q': 1.60218e-19}
+        const = CONSTANTS
 
     ee = ivcurves['ee']
     tc = ivcurves['tc']
@@ -520,9 +528,9 @@ def fit_desoto_sandia(ivcurves, specs, const=None, maxiter=5, eps1=1.e-3):
         T0 : float
             cell temperature at STC, default 25 [C]
         k : float
-            1.38066E-23 J/K (Boltzmann's constant)
+            Boltzmann's constant [J/K]
         q : float
-            1.60218E-19 Coulomb (elementary charge)
+            elementary charge [Coulomb]
 
     maxiter : int, default 5
         input that sets the maximum number of iterations for the parameter
@@ -570,16 +578,18 @@ def fit_desoto_sandia(ivcurves, specs, const=None, maxiter=5, eps1=1.e-3):
     ----------
     .. [1] W. De Soto et al., "Improvement and validation of a model for
        photovoltaic array performance", Solar Energy, vol 80, pp. 78-88,
-       2006.
+       2006. :doi:`10.1016/j.solener.2005.06.010`
     .. [2] C. Hansen, Parameter Estimation for Single Diode Models of
-       Photovoltaic Modules, Sandia National Laboratories Report SAND2015-2065
+       Photovoltaic Modules, Sandia National Laboratories Report SAND2015-2065.
+       :doi:`10.2172/1177157`
     .. [3] C. Hansen, Estimation of Parameters for Single Diode Models using
         Measured IV Curves, Proc. of the 39th IEEE PVSC, June 2013.
+        :doi:`10.1109/PVSC.2013.6744135`
     .. [4] PVLib MATLAB https://github.com/sandialabs/MATLAB_PV_LIB
     """
 
     if const is None:
-        const = {'E0': 1000.0, 'T0': 25.0, 'k': 1.38066e-23, 'q': 1.60218e-19}
+        const = CONSTANTS
 
     ee = ivcurves['ee']
     tc = ivcurves['tc']
@@ -938,7 +948,7 @@ def _update_io(voc, iph, io, rs, rsh, nnsvth):
 
     while maxerr > eps and k < niter:
         # Predict Voc
-        pvoc = v_from_i(rsh, rs, nnsvth, 0., tio, iph)
+        pvoc = v_from_i(0., iph, tio, rs, rsh, nnsvth)
 
         # Difference in Voc
         dvoc = pvoc - voc
@@ -1341,5 +1351,6 @@ def pvsyst_temperature_coeff(alpha_sc, gamma_ref, mu_gamma, I_L_ref, I_o_ref,
             I_o_ref, R_sh_ref, R_sh_0, R_s, cells_in_series, R_sh_exp, EgRef,
             temp_ref)
     pmp = maxp(temp_ref, *args)
-    gamma_pdc = derivative(maxp, temp_ref, args=args)
+    gamma_pdc = _first_order_centered_difference(maxp, x0=temp_ref, args=args)
+
     return gamma_pdc / pmp

pvlib/ivtools/utils.py

@@ -6,6 +6,7 @@ working with IV curves, or fitting equations to IV curve data.
 
 import numpy as np
 import pandas as pd
+from numpy.polynomial.polynomial import Polynomial as Poly
 
 
 # A small number used to decide when a slope is equivalent to zero
@@ -149,12 +150,12 @@ def rectify_iv_curve(voltage, current, decimals=None):
     ``rectify_iv_curve`` ensures that the IV curve lies in the first quadrant
     of the (voltage, current) plane. The returned IV curve:
 
-      * increases in voltage
-      * contains no negative current or voltage values
-      * contains no NaNs
-      * contains no points with duplicate voltage values. Where voltage
-        values are repeated, a single data point is substituted with current
-        equal to the average of current at duplicated voltages.
+    * increases in voltage
+    * contains no negative current or voltage values
+    * contains no NaNs
+    * contains no points with duplicate voltage values. Where voltage
+      values are repeated, a single data point is substituted with current
+      equal to the average of current at duplicated voltages.
     """
 
     df = pd.DataFrame(data=np.vstack((voltage, current)).T, columns=['v', 'i'])
@@ -423,3 +424,123 @@ def _schumaker_qspline(x, y):
     yhat = tmp2[:, 4]
     kflag = tmp2[:, 5]
     return t, c, yhat, kflag
+
+
+def astm_e1036(v, i, imax_limits=(0.75, 1.15), vmax_limits=(0.75, 1.15),
+               voc_points=3, isc_points=3, mp_fit_order=4):
+    '''
+    Extract photovoltaic IV parameters according to ASTM E1036. Assumes that
+    the power producing portion of the curve is in the first quadrant.
+
+    Parameters
+    ----------
+    v : array-like
+        Voltage points
+    i : array-like
+        Current points
+    imax_limits : tuple, default (0.75, 1.15)
+        Two-element tuple (low, high) specifying the fraction of estimated
+        Imp within which to fit a polynomial for max power calculation
+    vmax_limits : tuple, default (0.75, 1.15)
+        Two-element tuple (low, high) specifying the fraction of estimated
+        Vmp within which to fit a polynomial for max power calculation
+    voc_points : int, default 3
+        The number of points near open circuit to use for linear fit
+        and Voc calculation
+    isc_points : int, default 3
+        The number of points near short circuit to use for linear fit and
+        Isc calculation
+    mp_fit_order : int, default 4
+        The order of the polynomial fit of power vs. voltage near maximum
+        power
+
+
+    Returns
+    -------
+    dict
+        Results. The IV parameters are given by the keys 'voc', 'isc',
+        'vmp', 'imp', 'pmp', and 'ff'. The key 'mp_fit' gives the numpy
+        Polynomial object for the fit of power vs voltage near maximum
+        power.
+
+    References
+    ----------
+    .. [1] Standard Test Methods for Electrical Performance of Nonconcentrator
+       Terrestrial Photovoltaic Modules and Arrays Using Reference Cells,
+       ASTM E1036-15(2019), :doi:`10.1520/E1036-15R19`
+    '''
+
+    # Adapted from https://github.com/NREL/iv_params
+    # Copyright (c) 2022, Alliance for Sustainable Energy, LLC
+    # All rights reserved.
+
+    df = pd.DataFrame()
+    df['v'] = v
+    df['i'] = i
+    df['p'] = df['v'] * df['i']
+
+    # determine if we can use voc and isc estimates
+    i_min_ind = df['i'].abs().idxmin()
+    v_min_ind = df['v'].abs().idxmin()
+    voc_est = df['v'][i_min_ind]
+    isc_est = df['i'][v_min_ind]
+
+    # accept the estimates if they are close enough
+    # if not, perform a linear fit
+    if abs(df['i'][i_min_ind]) <= isc_est * 0.001:
+        voc = voc_est
+    else:
+        df['i_abs'] = df['i'].abs()
+        voc_df = df.nsmallest(voc_points, 'i_abs')
+        voc_fit = Poly.fit(voc_df['i'], voc_df['v'], 1)
+        voc = voc_fit(0)
+
+    if abs(df['v'][v_min_ind]) <= voc_est * 0.005:
+        isc = isc_est
+    else:
+        df['v_abs'] = df['v'].abs()
+        isc_df = df.nsmallest(isc_points, 'v_abs')
+        isc_fit = Poly.fit(isc_df['v'], isc_df['i'], 1)
+        isc = isc_fit(0)
+
+    # estimate max power point
+    max_index = df['p'].idxmax()
+    mp_est = df.loc[max_index]
+
+    # filter around max power
+    mask = (
+        (df['i'] >= imax_limits[0] * mp_est['i']) &
+        (df['i'] <= imax_limits[1] * mp_est['i']) &
+        (df['v'] >= vmax_limits[0] * mp_est['v']) &
+        (df['v'] <= vmax_limits[1] * mp_est['v'])
+    )
+    filtered = df[mask]
+
+    # fit polynomial and find max
+    mp_fit = Poly.fit(filtered['v'], filtered['p'], mp_fit_order)
+    # Note that this root finding procedure differs from
+    # the suggestion in the standard
+    roots = mp_fit.deriv().roots()
+    # only consider real roots
+    roots = roots.real[abs(roots.imag) < 1e-5]
+    # only consider roots in the relevant part of the domain
+    roots = roots[(roots < filtered['v'].max()) &
+                  (roots > filtered['v'].min())]
+    vmp = roots[np.argmax(mp_fit(roots))]
+    pmp = mp_fit(vmp)
+    # Imp isn't mentioned for update in the
+    # standard, but this seems to be in the intended spirit
+    imp = pmp / vmp
+
+    ff = pmp / (voc * isc)
+
+    result = {}
+    result['voc'] = voc
+    result['isc'] = isc
+    result['vmp'] = vmp
+    result['imp'] = imp
+    result['pmp'] = pmp
+    result['ff'] = ff
+    result['mp_fit'] = mp_fit
+
+    return result

pvlib/location.py

@@ -331,8 +331,9 @@ class Location:
         method : str, default 'pyephem'
             'pyephem', 'spa', or 'geometric'
 
-        kwargs are passed to the relevant functions. See
-        solarposition.sun_rise_set_transit_<method> for details.
+        kwargs :
+            Passed to the relevant functions. See
+            solarposition.sun_rise_set_transit_<method> for details.
 
         Returns
         -------

pvlib/modelchain.py

@@ -14,8 +14,7 @@ from dataclasses import dataclass, field
 from typing import Union, Tuple, Optional, TypeVar
 
 from pvlib import (atmosphere, clearsky, inverter, pvsystem, solarposition,
-                   temperature, tools)
-from pvlib.tracking import SingleAxisTracker
+                   temperature)
 import pvlib.irradiance  # avoid name conflict with full import
 from pvlib.pvsystem import _DC_MODEL_PARAMS
 from pvlib._deprecation import pvlibDeprecationWarning
@@ -254,6 +253,33 @@ def get_orientation(strategy, **kwargs):
     return surface_tilt, surface_azimuth
 
 
+def _getmcattr(self, attr):
+    """
+    Helper for __repr__ methods, needed to avoid recursion in property
+    lookups
+    """
+    out = getattr(self, attr)
+    try:
+        out = out.__name__
+    except AttributeError:
+        pass
+    return out
+
+
+def _mcr_repr(obj):
+    '''
+    Helper for ModelChainResult.__repr__
+    '''
+    if isinstance(obj, tuple):
+        return "Tuple (" + ", ".join([_mcr_repr(o) for o in obj]) + ")"
+    if isinstance(obj, pd.DataFrame):
+        return "DataFrame ({} rows x {} columns)".format(*obj.shape)
+    if isinstance(obj, pd.Series):
+        return "Series (length {})".format(len(obj))
+    # scalar, None, other?
+    return repr(obj)
+
+    
 # Type for fields that vary between arrays
 T = TypeVar('T')
 
@@ -385,6 +411,33 @@ class ModelChainResult:
             value = self._result_type(value)
         super().__setattr__(key, value)
 
+    def __repr__(self):
+        mc_attrs = dir(self)
+
+        def _head(obj):
+            try:
+                return obj[:3]
+            except:
+                return obj
+
+        if type(self.dc) is tuple:
+            num_arrays = len(self.dc)
+        else:
+            num_arrays = 1
+
+        desc1 = ('=== ModelChainResult === \n')
+        desc2 = (f'Number of Arrays: {num_arrays} \n')
+        attr = 'times'
+        desc3 = ('times (first 3)\n' +
+                 f'{_head(_getmcattr(self, attr))}' +
+                 '\n')
+        lines = []
+        for attr in mc_attrs:
+            if not (attr.startswith('_') or attr=='times'):
+                lines.append(f' {attr}: ' + _mcr_repr(getattr(self, attr)))
+        desc4 = '\n'.join(lines)
+        return (desc1 + desc2 + desc3 + desc4)
+
 
 class ModelChain:
     """
@@ -465,13 +518,6 @@ class ModelChain:
         Name of ModelChain instance.
     """
 
-    # list of deprecated attributes
-    _deprecated_attrs = ['solar_position', 'airmass', 'total_irrad',
-                         'aoi', 'aoi_modifier', 'spectral_modifier',
-                         'cell_temperature', 'effective_irradiance',
-                         'dc', 'ac', 'diode_params', 'tracking',
-                         'weather', 'times', 'losses']
-
     def __init__(self, system, location,
                  clearsky_model='ineichen',
                  transposition_model='haydavies',
@@ -503,26 +549,6 @@ class ModelChain:
 
         self.results = ModelChainResult()
 
-    def __getattr__(self, key):
-        if key in ModelChain._deprecated_attrs:
-            msg = f'ModelChain.{key} is deprecated and will' \
-                  f' be removed in v0.10. Use' \
-                  f' ModelChain.results.{key} instead'
-            warnings.warn(msg, pvlibDeprecationWarning)
-            return getattr(self.results, key)
-        # __getattr__ is only called if __getattribute__ fails.
-        # In that case we should check if key is a deprecated attribute,
-        # and fail with an AttributeError if it is not.
-        raise AttributeError
-
-    def __setattr__(self, key, value):
-        if key in ModelChain._deprecated_attrs:
-            msg = f'ModelChain.{key} is deprecated from v0.9. Use' \
-                  f' ModelChain.results.{key} instead'
-            warnings.warn(msg, pvlibDeprecationWarning)
-            setattr(self.results, key, value)
-        else:
-            super().__setattr__(key, value)
 
     @classmethod
     def with_pvwatts(cls, system, location,
@@ -560,7 +586,7 @@ class ModelChain:
         Examples
         --------
         >>> module_parameters = dict(gamma_pdc=-0.003, pdc0=4500)
-        >>> inverter_parameters = dict(pac0=4000)
+        >>> inverter_parameters = dict(pdc0=4000)
         >>> tparams = TEMPERATURE_MODEL_PARAMETERS['sapm']['open_rack_glass_glass']
         >>> system = PVSystem(surface_tilt=30, surface_azimuth=180,
         ...     module_parameters=module_parameters,
@@ -678,18 +704,8 @@ class ModelChain:
             'airmass_model', 'dc_model', 'ac_model', 'aoi_model',
             'spectral_model', 'temperature_model', 'losses_model'
         ]
-
-        def getmcattr(self, attr):
-            """needed to avoid recursion in property lookups"""
-            out = getattr(self, attr)
-            try:
-                out = out.__name__
-            except AttributeError:
-                pass
-            return out
-
         return ('ModelChain: \n  ' + '\n  '.join(
-            f'{attr}: {getmcattr(self, attr)}' for attr in attrs))
+            f'{attr}: {_getmcattr(self, attr)}' for attr in attrs))
 
     @property
     def dc_model(self):
@@ -1535,27 +1551,11 @@ class ModelChain:
         self._prep_inputs_solar_pos(weather)
         self._prep_inputs_airmass()
         self._prep_inputs_albedo(weather)
+        self._prep_inputs_fixed()
 
-        # PVSystem.get_irradiance and SingleAxisTracker.get_irradiance
-        # and PVSystem.get_aoi and SingleAxisTracker.get_aoi
-        # have different method signatures. Use partial to handle
-        # the differences.
-        if isinstance(self.system, SingleAxisTracker):
-            self._prep_inputs_tracking()
-            get_irradiance = partial(
-                self.system.get_irradiance,
-                self.results.tracking['surface_tilt'],
-                self.results.tracking['surface_azimuth'],
-                self.results.solar_position['apparent_zenith'],
-                self.results.solar_position['azimuth'])
-        else:
-            self._prep_inputs_fixed()
-            get_irradiance = partial(
-                self.system.get_irradiance,
-                self.results.solar_position['apparent_zenith'],
-                self.results.solar_position['azimuth'])
-
-        self.results.total_irrad = get_irradiance(
+        self.results.total_irrad = self.system.get_irradiance(
+            self.results.solar_position['apparent_zenith'],
+            self.results.solar_position['azimuth'],
             _tuple_from_dfs(self.results.weather, 'dni'),
             _tuple_from_dfs(self.results.weather, 'ghi'),
             _tuple_from_dfs(self.results.weather, 'dhi'),
@@ -1636,10 +1636,7 @@ class ModelChain:
         self._prep_inputs_solar_pos(data)
         self._prep_inputs_airmass()
 
-        if isinstance(self.system, SingleAxisTracker):
-            self._prep_inputs_tracking()
-        else:
-            self._prep_inputs_fixed()
+        self._prep_inputs_fixed()
 
         return self
 
@@ -1686,7 +1683,7 @@ class ModelChain:
             self.temperature_model()
         return self
 
-    def _prepare_temperature(self, data=None):
+    def _prepare_temperature(self, data):
         """
         Sets cell_temperature using inputs in data and the specified
         temperature model.
@@ -1699,7 +1696,7 @@ class ModelChain:
 
         Parameters
         ----------
-        data : DataFrame, default None
+        data : DataFrame
             May contain columns ``'cell_temperature'`` or
             ``'module_temperaure'``.
 
@@ -1878,13 +1875,13 @@ class ModelChain:
 
         return self
 
-    def _run_from_effective_irrad(self, data=None):
+    def _run_from_effective_irrad(self, data):
         """
         Executes the temperature, DC, losses and AC models.
 
         Parameters
         ----------
-        data : DataFrame, or tuple of DataFrame, default None
+        data : DataFrame, or tuple of DataFrame
             If optional column ``'cell_temperature'`` is provided, these values
             are used instead of `temperature_model`. If optional column
             `module_temperature` is provided, `temperature_model` must be
@@ -1907,7 +1904,7 @@ class ModelChain:
 
         return self
 
-    def run_model_from_effective_irradiance(self, data=None):
+    def run_model_from_effective_irradiance(self, data):
         """
         Run the model starting with effective irradiance in the plane of array.