diff --git a/docs/sphinx/source/timetimezones.rst b/docs/sphinx/source/timetimezones.rst
index dc5f91f565..235bd4776c 100644
--- a/docs/sphinx/source/timetimezones.rst
+++ b/docs/sphinx/source/timetimezones.rst
@@ -278,7 +278,7 @@ Let's first examine how pvlib handles time when it imports a TMY3 file.
 
 The metadata has a ``'TZ'`` key with a value of ``-9.0``. This is the
 UTC offset in hours in which the data has been recorded. The
-:py:func:`~pvlib.tmy.readtmy3` function read the data in the file,
+:py:func:`~pvlib.iotools.read_tmy3` function read the data in the file,
 created a :py:class:`~pandas.DataFrame` with that data, and then
 localized the DataFrame's index to have this fixed offset. Here, we
 print just a few of the rows and columns of the large dataframe.
@@ -289,7 +289,7 @@ print just a few of the rows and columns of the large dataframe.
 
     tmy3_data.loc[tmy3_data.index[0:3], ['GHI', 'DNI', 'AOD']]
 
-The :py:func:`~pvlib.tmy.readtmy2` function also returns a DataFrame
+The :py:func:`~pvlib.iotools.read_tmy2` function also returns a DataFrame
 with a localized DatetimeIndex.
 
 Solar position
diff --git a/pvlib/ivtools/sde.py b/pvlib/ivtools/sde.py
index 1042a19736..dad1add22f 100644
--- a/pvlib/ivtools/sde.py
+++ b/pvlib/ivtools/sde.py
@@ -90,7 +90,7 @@ def fit_sandia_simple(voltage, current, v_oc=None, i_sc=None, v_mp_i_mp=None,
         I = I_{L} - I_{0} (\exp \frac{V + I R_{s}}{nNsVth} - 1)
         - \frac{V + I R_{s}}{R_{sh}}
 
-    See :py:func:`pvsystem.singlediode` for definition of the parameters.
+    See :py:func:`pvlib.pvsystem.singlediode` for definition of the parameters.
 
     The extraction method [2]_ proceeds in six steps.
 
diff --git a/pvlib/location.py b/pvlib/location.py
index 93d54fcd9e..97767a5801 100644
--- a/pvlib/location.py
+++ b/pvlib/location.py
@@ -164,7 +164,7 @@ def from_epw(cls, metadata, data=None, **kwargs):
     def get_solarposition(self, times, pressure=None, temperature=12,
                           **kwargs):
         """
-        Uses the :py:func:`solarposition.get_solarposition` function
+        Uses the :py:func:`pvlib.solarposition.get_solarposition` function
         to calculate the solar zenith, azimuth, etc. at this location.
 
         Parameters
@@ -173,11 +173,11 @@ def get_solarposition(self, times, pressure=None, temperature=12,
             Must be localized or UTC will be assumed.
         pressure : None, float, or array-like, default None
             If None, pressure will be calculated using
-            :py:func:`atmosphere.alt2pres` and ``self.altitude``.
+            :py:func:`pvlib.atmosphere.alt2pres` and ``self.altitude``.
         temperature : None, float, or array-like, default 12
 
         kwargs
-            passed to :py:func:`solarposition.get_solarposition`
+            passed to :py:func:`pvlib.solarposition.get_solarposition`
 
         Returns
         -------
diff --git a/pvlib/pvsystem.py b/pvlib/pvsystem.py
index eb288978d4..6bb89f34a3 100644
--- a/pvlib/pvsystem.py
+++ b/pvlib/pvsystem.py
@@ -627,7 +627,7 @@ def sapm(self, effective_irradiance, temp_cell):
     @deprecated('0.9', alternative='PVSystem.get_cell_temperature',
                 removal='0.10.0')
     def sapm_celltemp(self, poa_global, temp_air, wind_speed):
-        """Uses :py:func:`temperature.sapm_cell` to calculate cell
+        """Uses :py:func:`pvlib.temperature.sapm_cell` to calculate cell
         temperatures.
 
         Parameters
@@ -720,7 +720,7 @@ def sapm_effective_irradiance(self, poa_direct, poa_diffuse,
     @deprecated('0.9', alternative='PVSystem.get_cell_temperature',
                 removal='0.10.0')
     def pvsyst_celltemp(self, poa_global, temp_air, wind_speed=1.0):
-        """Uses :py:func:`temperature.pvsyst_cell` to calculate cell
+        """Uses :py:func:`pvlib.temperature.pvsyst_cell` to calculate cell
         temperature.
 
         Parameters
@@ -756,7 +756,7 @@ def pvsyst_celltemp(self, poa_global, temp_air, wind_speed=1.0):
                 removal='0.10.0')
     def faiman_celltemp(self, poa_global, temp_air, wind_speed=1.0):
         """
-        Use :py:func:`temperature.faiman` to calculate cell temperature.
+        Use :py:func:`pvlib.temperature.faiman` to calculate cell temperature.
 
         Parameters
         ----------
@@ -791,7 +791,7 @@ def faiman_celltemp(self, poa_global, temp_air, wind_speed=1.0):
                 removal='0.10.0')
     def fuentes_celltemp(self, poa_global, temp_air, wind_speed):
         """
-        Use :py:func:`temperature.fuentes` to calculate cell temperature.
+        Use :py:func:`pvlib.temperature.fuentes` to calculate cell temperature.
 
         Parameters
         ----------
@@ -834,7 +834,8 @@ def fuentes_celltemp(self, poa_global, temp_air, wind_speed):
     def noct_sam_celltemp(self, poa_global, temp_air, wind_speed,
                           effective_irradiance=None):
         """
-        Use :py:func:`temperature.noct_sam` to calculate cell temperature.
+        Use :py:func:`pvlib.temperature.noct_sam` to calculate cell
+        temperature.
 
         Parameters
         ----------
@@ -871,7 +872,7 @@ def noct_sam_celltemp(self, poa_global, temp_air, wind_speed,
     @_unwrap_single_value
     def first_solar_spectral_loss(self, pw, airmass_absolute):
         """
-        Use the :py:func:`first_solar_spectral_correction` function to
+        Use :py:func:`pvlib.atmosphere.first_solar_spectral_correction` to
         calculate the spectral loss modifier. The model coefficients are
         specific to the module's cell type, and are determined by searching
         for one of the following keys in self.module_parameters (in order):