From 4143e2ae93ee6f4970c566c7823c141cc42a4cf0 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Thu, 18 Dec 2025 23:11:55 +0000
Subject: [PATCH 01/18] Adds faiman_rad and ross models to
 get_cell_temperature().

---
 pvlib/pvsystem.py | 25 +++++++++++++++++++++----
 1 file changed, 21 insertions(+), 4 deletions(-)

diff --git a/pvlib/pvsystem.py b/pvlib/pvsystem.py
index fe85359b99..76f53b3e74 100644
--- a/pvlib/pvsystem.py
+++ b/pvlib/pvsystem.py
@@ -1217,11 +1217,12 @@ def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
             Ambient dry bulb temperature [C]
 
         wind_speed : numeric
-            Wind speed [m/s]
+            Wind speed [m/s], although can be ``None`` for ``'ross'`` model
 
         model : str
             Supported models include ``'sapm'``, ``'pvsyst'``,
-            ``'faiman'``, ``'fuentes'``, and ``'noct_sam'``
+            ``'faiman'``, ``'faiman_rad'``, ``'fuentes'``, ``'noct_sam'``,
+            and ``'ross'``
 
         effective_irradiance : numeric, optional
             The irradiance that is converted to photocurrent in W/m^2.
@@ -1267,6 +1268,12 @@ def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
             required = tuple()
             optional = _build_kwargs(['u0', 'u1'],
                                      self.temperature_model_parameters)
+        elif model == 'faiman_rad':
+            func = temperature.faiman_rad
+            required = ()
+            optional = _build_kwargs(['ir_down','u0','u1',
+                                      'sky_view','emissivity'],
+                                     self.temperature_model_parameters)
         elif model == 'fuentes':
             func = temperature.fuentes
             required = _build_tcell_args(['noct_installed'])
@@ -1283,11 +1290,21 @@ def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
             optional = _build_kwargs(['transmittance_absorptance',
                                       'array_height', 'mount_standoff'],
                                      self.temperature_model_parameters)
+        elif model == 'ross':
+            func = temperature.ross
+            required = ()
+            # either noct or k must be defined
+            optional = _build_kwargs(['noct','k'],
+                                      self.temperature_model_parameters)
         else:
             raise ValueError(f'{model} is not a valid cell temperature model')
 
-        temperature_cell = func(poa_global, temp_air, wind_speed,
-                                *required, **optional)
+        if model == 'ross':
+            temperature_cell = func(poa_global, temp_air,
+                                    *required, **optional)
+        else:
+            temperature_cell = func(poa_global, temp_air,  wind_speed,
+                                    *required, **optional)
         return temperature_cell
 
     def dc_ohms_from_percent(self):

From 70bb9483bb890433b3afb2142d84258afaaf5f2e Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Thu, 18 Dec 2025 23:24:47 +0000
Subject: [PATCH 02/18] Addresses flake8 formatting issues.

---
 pvlib/pvsystem.py | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/pvlib/pvsystem.py b/pvlib/pvsystem.py
index 76f53b3e74..98e30a168a 100644
--- a/pvlib/pvsystem.py
+++ b/pvlib/pvsystem.py
@@ -1271,8 +1271,8 @@ def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
         elif model == 'faiman_rad':
             func = temperature.faiman_rad
             required = ()
-            optional = _build_kwargs(['ir_down','u0','u1',
-                                      'sky_view','emissivity'],
+            optional = _build_kwargs(['ir_down', 'u0', 'u1',
+                                      'sky_view', 'emissivity'],
                                      self.temperature_model_parameters)
         elif model == 'fuentes':
             func = temperature.fuentes
@@ -1294,8 +1294,8 @@ def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
             func = temperature.ross
             required = ()
             # either noct or k must be defined
-            optional = _build_kwargs(['noct','k'],
-                                      self.temperature_model_parameters)
+            optional = _build_kwargs(['noct', 'k'],
+                                     self.temperature_model_parameters)
         else:
             raise ValueError(f'{model} is not a valid cell temperature model')
 

From acc43b45c60bade542ee36c53418c846468e5434 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Sun, 21 Dec 2025 10:44:12 +0000
Subject: [PATCH 03/18] Updates contribution md file.

---
 docs/sphinx/source/whatsnew/v0.13.2.rst | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/docs/sphinx/source/whatsnew/v0.13.2.rst b/docs/sphinx/source/whatsnew/v0.13.2.rst
index b6170a6053..85aab28275 100644
--- a/docs/sphinx/source/whatsnew/v0.13.2.rst
+++ b/docs/sphinx/source/whatsnew/v0.13.2.rst
@@ -46,6 +46,8 @@ Enhancements
   MERRA-2 reanalysis data. (:pull:`2572`)
 * Add :py:func:`~pvlib.spectrum.spectral_factor_polo`, a function for estimating
   spectral mismatch factors for vertical PV façades. (:issue:`2406`, :pull:`2491`)
+* Includes `ross` and `faiman_rad` in the allowed models within 
+  `pvlib.pvsystem.PVSystem.get_cell_temperature` (:issue:`2625`, :pull:`2631`)
 
 Documentation
 ~~~~~~~~~~~~~
@@ -87,3 +89,4 @@ Contributors
 * Anton Driesse (:ghuser:`adriesse`)
 * Rajiv Daxini (:ghuser:`RDaxini`)
 * Kevin Anderson (:ghuser:`kandersolar`)
+* Rodrigo Amaro e Silva (:ghuser:`ramaroesilva`)
\ No newline at end of file

From 69c7f434a207d17972795e930aeeed8c72d06f59 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Sun, 21 Dec 2025 13:52:52 +0000
Subject: [PATCH 04/18] Adds test for ross option, placeholder for faiman_rad.

---
 tests/test_pvsystem.py | 34 ++++++++++++++++++++++++++++++++++
 1 file changed, 34 insertions(+)

diff --git a/tests/test_pvsystem.py b/tests/test_pvsystem.py
index 4fbd782e65..852e84d602 100644
--- a/tests/test_pvsystem.py
+++ b/tests/test_pvsystem.py
@@ -494,6 +494,40 @@ def test_PVSystem_faiman_celltemp(mocker):
     assert_allclose(out, 56.4, atol=1e-1)
 
 
+def test_PVSystem_faiman_rad_celltemp(mocker):
+    # default values (need to account for a certain ir_down and adjust u0, u1)
+    ir_down = None
+    u0, u1 = 25.0, 6.84
+    sky_view = 1.0
+    emissivity = 0.88
+    
+    temp_model_params = {'ir_down': ir_down, 'u0': u0, 'u1': u1,
+                         'sky_view': sky_view, 'emissivity': emissivity}
+    system = pvsystem.PVSystem(temperature_model_parameters=temp_model_params)
+    mocker.spy(temperature, 'faiman_rad')
+    temps = 25
+    irrads = 1000
+    winds = 1
+    out = system.get_cell_temperature(irrads, temps, winds, model='faiman_rad')
+    temperature.faiman_rad.assert_called_once_with(irrads, temps, winds, u0, u1,
+                                                   sky_view, emissivity)
+    assert_allclose(out, 56.4, atol=1e-1)
+
+
+def test_PVSystem_ross_celltemp(mocker):
+    # example value
+    k = 0.0208 # free-standing system
+    
+    temp_model_params = {'k': k}
+    system = pvsystem.PVSystem(temperature_model_parameters=temp_model_params)
+    mocker.spy(temperature, 'ross')
+    temps = 25
+    irrads = 1000
+    out = system.get_cell_temperature(irrads, temps, model='ross')
+    temperature.faiman_rad.assert_called_once_with(irrads, temps, k)
+    assert_allclose(out, 45.8, atol=1e-1)
+
+
 def test_PVSystem_noct_celltemp(mocker):
     poa_global, temp_air, wind_speed, noct, module_efficiency = (
         1000., 25., 1., 45., 0.2)

From 5998725732f7af84c7778f41d483604f6155c837 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Sun, 21 Dec 2025 13:55:57 +0000
Subject: [PATCH 05/18] Corrects flake8 issues.

---
 tests/test_pvsystem.py | 9 +++++----
 1 file changed, 5 insertions(+), 4 deletions(-)

diff --git a/tests/test_pvsystem.py b/tests/test_pvsystem.py
index 852e84d602..c103f9cff5 100644
--- a/tests/test_pvsystem.py
+++ b/tests/test_pvsystem.py
@@ -500,7 +500,7 @@ def test_PVSystem_faiman_rad_celltemp(mocker):
     u0, u1 = 25.0, 6.84
     sky_view = 1.0
     emissivity = 0.88
-    
+
     temp_model_params = {'ir_down': ir_down, 'u0': u0, 'u1': u1,
                          'sky_view': sky_view, 'emissivity': emissivity}
     system = pvsystem.PVSystem(temperature_model_parameters=temp_model_params)
@@ -509,15 +509,16 @@ def test_PVSystem_faiman_rad_celltemp(mocker):
     irrads = 1000
     winds = 1
     out = system.get_cell_temperature(irrads, temps, winds, model='faiman_rad')
-    temperature.faiman_rad.assert_called_once_with(irrads, temps, winds, u0, u1,
+    temperature.faiman_rad.assert_called_once_with(irrads, temps, winds,
+                                                   u0, u1,
                                                    sky_view, emissivity)
     assert_allclose(out, 56.4, atol=1e-1)
 
 
 def test_PVSystem_ross_celltemp(mocker):
     # example value
-    k = 0.0208 # free-standing system
-    
+    k = 0.0208  # free-standing system
+
     temp_model_params = {'k': k}
     system = pvsystem.PVSystem(temperature_model_parameters=temp_model_params)
     mocker.spy(temperature, 'ross')

From 61e8163be1a86c9bd455eedff6b040be242e6b83 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Sun, 21 Dec 2025 14:06:06 +0000
Subject: [PATCH 06/18] Adds missing wind_speed=None when testing
 get_cell_temperature for ross.

---
 tests/test_pvsystem.py | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/tests/test_pvsystem.py b/tests/test_pvsystem.py
index c103f9cff5..a135932631 100644
--- a/tests/test_pvsystem.py
+++ b/tests/test_pvsystem.py
@@ -524,7 +524,8 @@ def test_PVSystem_ross_celltemp(mocker):
     mocker.spy(temperature, 'ross')
     temps = 25
     irrads = 1000
-    out = system.get_cell_temperature(irrads, temps, model='ross')
+    winds = None
+    out = system.get_cell_temperature(irrads, temps, winds, model='ross')
     temperature.faiman_rad.assert_called_once_with(irrads, temps, k)
     assert_allclose(out, 45.8, atol=1e-1)
 

From 22f94a786d8b0056d1e6839bbaafd2359f1d09b6 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Sun, 21 Dec 2025 15:34:56 +0000
Subject: [PATCH 07/18] Corrects wrong function naming.

---
 tests/test_pvsystem.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/tests/test_pvsystem.py b/tests/test_pvsystem.py
index a135932631..a4e7dcc285 100644
--- a/tests/test_pvsystem.py
+++ b/tests/test_pvsystem.py
@@ -526,7 +526,7 @@ def test_PVSystem_ross_celltemp(mocker):
     irrads = 1000
     winds = None
     out = system.get_cell_temperature(irrads, temps, winds, model='ross')
-    temperature.faiman_rad.assert_called_once_with(irrads, temps, k)
+    temperature.ross.assert_called_once_with(irrads, temps, k)
     assert_allclose(out, 45.8, atol=1e-1)
 
 

From aaa52e32255b7e004c4dffca22aebd51b6226ce9 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Sun, 21 Dec 2025 15:45:58 +0000
Subject: [PATCH 08/18] Adds missing input to functions within tests.

---
 tests/test_pvsystem.py | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/tests/test_pvsystem.py b/tests/test_pvsystem.py
index a4e7dcc285..f9619b79db 100644
--- a/tests/test_pvsystem.py
+++ b/tests/test_pvsystem.py
@@ -510,7 +510,7 @@ def test_PVSystem_faiman_rad_celltemp(mocker):
     winds = 1
     out = system.get_cell_temperature(irrads, temps, winds, model='faiman_rad')
     temperature.faiman_rad.assert_called_once_with(irrads, temps, winds,
-                                                   u0, u1,
+                                                   ir_down, u0, u1,
                                                    sky_view, emissivity)
     assert_allclose(out, 56.4, atol=1e-1)
 
@@ -518,6 +518,7 @@ def test_PVSystem_faiman_rad_celltemp(mocker):
 def test_PVSystem_ross_celltemp(mocker):
     # example value
     k = 0.0208  # free-standing system
+    noct = None # function expects only one of two, the other as None
 
     temp_model_params = {'k': k}
     system = pvsystem.PVSystem(temperature_model_parameters=temp_model_params)
@@ -526,7 +527,7 @@ def test_PVSystem_ross_celltemp(mocker):
     irrads = 1000
     winds = None
     out = system.get_cell_temperature(irrads, temps, winds, model='ross')
-    temperature.ross.assert_called_once_with(irrads, temps, k)
+    temperature.ross.assert_called_once_with(irrads, temps, noct, k)
     assert_allclose(out, 45.8, atol=1e-1)
 
 

From c6d8ab9601cee7cb1cb579f652452d81d73b6368 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Sun, 21 Dec 2025 15:53:25 +0000
Subject: [PATCH 09/18] Corrects test + flake8 issue.

---
 tests/test_pvsystem.py | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/tests/test_pvsystem.py b/tests/test_pvsystem.py
index f9619b79db..54a7b3759b 100644
--- a/tests/test_pvsystem.py
+++ b/tests/test_pvsystem.py
@@ -518,7 +518,7 @@ def test_PVSystem_faiman_rad_celltemp(mocker):
 def test_PVSystem_ross_celltemp(mocker):
     # example value
     k = 0.0208  # free-standing system
-    noct = None # function expects only one of two, the other as None
+    noct = None  # function expects only one of two, the other as None
 
     temp_model_params = {'k': k}
     system = pvsystem.PVSystem(temperature_model_parameters=temp_model_params)
@@ -527,7 +527,7 @@ def test_PVSystem_ross_celltemp(mocker):
     irrads = 1000
     winds = None
     out = system.get_cell_temperature(irrads, temps, winds, model='ross')
-    temperature.ross.assert_called_once_with(irrads, temps, noct, k)
+    temperature.ross.assert_called_once_with(irrads, temps, k=k)
     assert_allclose(out, 45.8, atol=1e-1)
 
 

From 78eb4c3e64a0f2ae8beb012212c462074ab828a9 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Sun, 21 Dec 2025 15:56:00 +0000
Subject: [PATCH 10/18] Minor update in documentation + flake8 issue.

---
 tests/test_pvsystem.py | 5 ++---
 1 file changed, 2 insertions(+), 3 deletions(-)

diff --git a/tests/test_pvsystem.py b/tests/test_pvsystem.py
index 54a7b3759b..c373216c43 100644
--- a/tests/test_pvsystem.py
+++ b/tests/test_pvsystem.py
@@ -516,10 +516,9 @@ def test_PVSystem_faiman_rad_celltemp(mocker):
 
 
 def test_PVSystem_ross_celltemp(mocker):
-    # example value
+    # example value (could use equivalent noct as alternative input)
     k = 0.0208  # free-standing system
-    noct = None  # function expects only one of two, the other as None
-
+  
     temp_model_params = {'k': k}
     system = pvsystem.PVSystem(temperature_model_parameters=temp_model_params)
     mocker.spy(temperature, 'ross')

From 402b53caab321d11e0d68d2764c54b0e2b063568 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Sun, 21 Dec 2025 15:57:12 +0000
Subject: [PATCH 11/18] Another flake8 issue handled.

---
 tests/test_pvsystem.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/tests/test_pvsystem.py b/tests/test_pvsystem.py
index c373216c43..3445253370 100644
--- a/tests/test_pvsystem.py
+++ b/tests/test_pvsystem.py
@@ -518,7 +518,7 @@ def test_PVSystem_faiman_rad_celltemp(mocker):
 def test_PVSystem_ross_celltemp(mocker):
     # example value (could use equivalent noct as alternative input)
     k = 0.0208  # free-standing system
-  
+
     temp_model_params = {'k': k}
     system = pvsystem.PVSystem(temperature_model_parameters=temp_model_params)
     mocker.spy(temperature, 'ross')

From 9461e013af32d00cab0aa3d23df609e111111344 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Mon, 22 Dec 2025 08:09:02 +0000
Subject: [PATCH 12/18] Improve faiman_rad test and completes
 get_cell_temperature doc.

---
 pvlib/pvsystem.py      | 8 +++++---
 tests/test_pvsystem.py | 9 +++++----
 2 files changed, 10 insertions(+), 7 deletions(-)

diff --git a/pvlib/pvsystem.py b/pvlib/pvsystem.py
index 98e30a168a..22da2f5f2a 100644
--- a/pvlib/pvsystem.py
+++ b/pvlib/pvsystem.py
@@ -431,7 +431,8 @@ def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
 
         model : str
             Supported models include ``'sapm'``, ``'pvsyst'``,
-            ``'faiman'``, ``'fuentes'``, and ``'noct_sam'``
+            ``'faiman'``, ``'faiman_rad'``, ``'fuentes'``, ``'noct_sam'``,
+            and ``'ross'``
 
         effective_irradiance : numeric or tuple of numeric, optional
             The irradiance that is converted to photocurrent in W/m^2.
@@ -1236,8 +1237,9 @@ def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
         See Also
         --------
         pvlib.temperature.sapm_cell, pvlib.temperature.pvsyst_cell,
-        pvlib.temperature.faiman, pvlib.temperature.fuentes,
-        pvlib.temperature.noct_sam
+        pvlib.temperature.faiman, pvlib.temperature.faiman_rad,
+        pvlib.temperature.fuentes, pvlib.temperature.noct_sam,
+        pvlib.temperature.ross
 
         Notes
         -----
diff --git a/tests/test_pvsystem.py b/tests/test_pvsystem.py
index 3445253370..1cf4f87ed3 100644
--- a/tests/test_pvsystem.py
+++ b/tests/test_pvsystem.py
@@ -495,9 +495,10 @@ def test_PVSystem_faiman_celltemp(mocker):
 
 
 def test_PVSystem_faiman_rad_celltemp(mocker):
-    # default values (need to account for a certain ir_down and adjust u0, u1)
-    ir_down = None
-    u0, u1 = 25.0, 6.84
+    ir_down = 50  # arbitrary value
+    # default values, u0 and u1 being adjusted in same proportion as in
+    # https://www.osti.gov/servlets/purl/1884890/ (not suggested, just example)
+    u0,u1 = 25.0*0.86, 6.84*0.88 
     sky_view = 1.0
     emissivity = 0.88
 
@@ -512,7 +513,7 @@ def test_PVSystem_faiman_rad_celltemp(mocker):
     temperature.faiman_rad.assert_called_once_with(irrads, temps, winds,
                                                    ir_down, u0, u1,
                                                    sky_view, emissivity)
-    assert_allclose(out, 56.4, atol=1e-1)
+    assert_allclose(out, 48.6, atol=1e-1)
 
 
 def test_PVSystem_ross_celltemp(mocker):

From 90350669fb02d0b7fe64394169c8cb25bb4a9427 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Mon, 22 Dec 2025 08:12:05 +0000
Subject: [PATCH 13/18] Addresses flake8.

---
 tests/test_pvsystem.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/tests/test_pvsystem.py b/tests/test_pvsystem.py
index 1cf4f87ed3..109db1859c 100644
--- a/tests/test_pvsystem.py
+++ b/tests/test_pvsystem.py
@@ -498,7 +498,7 @@ def test_PVSystem_faiman_rad_celltemp(mocker):
     ir_down = 50  # arbitrary value
     # default values, u0 and u1 being adjusted in same proportion as in
     # https://www.osti.gov/servlets/purl/1884890/ (not suggested, just example)
-    u0,u1 = 25.0*0.86, 6.84*0.88 
+    u0, u1 = 25.0*0.86, 6.84*0.88
     sky_view = 1.0
     emissivity = 0.88
 

From 6f5bc0d1efdb54bea50e146369ea2c970078cb45 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Mon, 22 Dec 2025 16:04:44 +0000
Subject: [PATCH 14/18] Update whatsnew doc.

---
 docs/sphinx/source/whatsnew/v0.13.2.rst | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/docs/sphinx/source/whatsnew/v0.13.2.rst b/docs/sphinx/source/whatsnew/v0.13.2.rst
index 85aab28275..72dce0efb0 100644
--- a/docs/sphinx/source/whatsnew/v0.13.2.rst
+++ b/docs/sphinx/source/whatsnew/v0.13.2.rst
@@ -46,8 +46,8 @@ Enhancements
   MERRA-2 reanalysis data. (:pull:`2572`)
 * Add :py:func:`~pvlib.spectrum.spectral_factor_polo`, a function for estimating
   spectral mismatch factors for vertical PV façades. (:issue:`2406`, :pull:`2491`)
-* Includes `ross` and `faiman_rad` in the allowed models within 
-  `pvlib.pvsystem.PVSystem.get_cell_temperature` (:issue:`2625`, :pull:`2631`)
+* Include `ross` and `faiman_rad` in the allowed models within 
+  :py:meth:`pvlib.pvsystem.PVSystem.get_cell_temperature` (:issue:`2625`, :pull:`2631`)
 
 Documentation
 ~~~~~~~~~~~~~

From a536cb27745f38047bbd4d613a6c192c42166bc7 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Mon, 22 Dec 2025 16:23:47 +0000
Subject: [PATCH 15/18] Defines ir_down as weather variable Within faiman_rad.

---
 pvlib/pvsystem.py      | 23 +++++++++++++++++------
 tests/test_pvsystem.py |  7 ++++---
 2 files changed, 21 insertions(+), 9 deletions(-)

diff --git a/pvlib/pvsystem.py b/pvlib/pvsystem.py
index 22da2f5f2a..cfc4f0d6b4 100644
--- a/pvlib/pvsystem.py
+++ b/pvlib/pvsystem.py
@@ -414,7 +414,7 @@ def get_iam(self, aoi, iam_model='physical'):
 
     @_unwrap_single_value
     def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
-                             effective_irradiance=None):
+                             effective_irradiance=None, ir_down=None):
         """
         Determine cell temperature using the method specified by ``model``.
 
@@ -437,6 +437,10 @@ def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
         effective_irradiance : numeric or tuple of numeric, optional
             The irradiance that is converted to photocurrent in W/m^2.
             Only used for some models.
+            
+        ir_down: numeric, optional
+            Downwelling infrared radiation from the sky, measured on a 
+            horizontal surface in W/m^2. Only used in ``'faiman_rad'`` model.
 
         Returns
         -------
@@ -460,14 +464,16 @@ def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
         # Not used for all models, but Array.get_cell_temperature handles it
         effective_irradiance = self._validate_per_array(effective_irradiance,
                                                         system_wide=True)
+        ir_down = self._validate_per_array(ir_down, system_wide=True)
 
         return tuple(
             array.get_cell_temperature(poa_global, temp_air, wind_speed,
-                                       model, effective_irradiance)
-            for array, poa_global, temp_air, wind_speed, effective_irradiance
+                                       model, effective_irradiance, ir_down)
+            for array, poa_global, temp_air, wind_speed, effective_irradiance,
+            ir_down
             in zip(
                 self.arrays, poa_global, temp_air, wind_speed,
-                effective_irradiance
+                effective_irradiance, ir_down
             )
         )
 
@@ -1205,7 +1211,7 @@ def get_iam(self, aoi, iam_model='physical'):
             raise ValueError(model + ' is not a valid IAM model')
 
     def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
-                             effective_irradiance=None):
+                             effective_irradiance=None, ir_down=None):
         """
         Determine cell temperature using the method specified by ``model``.
 
@@ -1229,6 +1235,10 @@ def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
             The irradiance that is converted to photocurrent in W/m^2.
             Only used for some models.
 
+        ir_down: numeric, optional
+            Downwelling infrared radiation from the sky, measured on a 
+            horizontal surface in W/m^2. Only used in ``'faiman_rad'`` model.
+
         Returns
         -------
         numeric
@@ -1271,7 +1281,8 @@ def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
             optional = _build_kwargs(['u0', 'u1'],
                                      self.temperature_model_parameters)
         elif model == 'faiman_rad':
-            func = temperature.faiman_rad
+            func = functools.partial(temperature.faiman_rad,
+                                     ir_down=ir_down)
             required = ()
             optional = _build_kwargs(['ir_down', 'u0', 'u1',
                                       'sky_view', 'emissivity'],
diff --git a/tests/test_pvsystem.py b/tests/test_pvsystem.py
index 109db1859c..50e338fa00 100644
--- a/tests/test_pvsystem.py
+++ b/tests/test_pvsystem.py
@@ -502,14 +502,15 @@ def test_PVSystem_faiman_rad_celltemp(mocker):
     sky_view = 1.0
     emissivity = 0.88
 
-    temp_model_params = {'ir_down': ir_down, 'u0': u0, 'u1': u1,
-                         'sky_view': sky_view, 'emissivity': emissivity}
+    temp_model_params = {'u0': u0, 'u1': u1, 'sky_view': sky_view,
+                         'emissivity': emissivity}
     system = pvsystem.PVSystem(temperature_model_parameters=temp_model_params)
     mocker.spy(temperature, 'faiman_rad')
     temps = 25
     irrads = 1000
     winds = 1
-    out = system.get_cell_temperature(irrads, temps, winds, model='faiman_rad')
+    out = system.get_cell_temperature(irrads, temps, winds, ir_down,
+                                      model='faiman_rad')
     temperature.faiman_rad.assert_called_once_with(irrads, temps, winds,
                                                    ir_down, u0, u1,
                                                    sky_view, emissivity)

From bc625789412e211ab1a1ec7d28e75fd22d736777 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Mon, 22 Dec 2025 16:29:06 +0000
Subject: [PATCH 16/18] Handles flake8 issues.

---
 pvlib/pvsystem.py | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/pvlib/pvsystem.py b/pvlib/pvsystem.py
index cfc4f0d6b4..d307c57613 100644
--- a/pvlib/pvsystem.py
+++ b/pvlib/pvsystem.py
@@ -437,9 +437,9 @@ def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
         effective_irradiance : numeric or tuple of numeric, optional
             The irradiance that is converted to photocurrent in W/m^2.
             Only used for some models.
-            
+   
         ir_down: numeric, optional
-            Downwelling infrared radiation from the sky, measured on a 
+            Downwelling infrared radiation from the sky, measured on a
             horizontal surface in W/m^2. Only used in ``'faiman_rad'`` model.
 
         Returns
@@ -1236,7 +1236,7 @@ def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
             Only used for some models.
 
         ir_down: numeric, optional
-            Downwelling infrared radiation from the sky, measured on a 
+            Downwelling infrared radiation from the sky, measured on a
             horizontal surface in W/m^2. Only used in ``'faiman_rad'`` model.
 
         Returns

From 7612e56f9a5b28ab82fb5b909241de26f77923e2 Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Mon, 22 Dec 2025 22:55:29 +0000
Subject: [PATCH 17/18] Corrects bug in test.

---
 tests/test_pvsystem.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/tests/test_pvsystem.py b/tests/test_pvsystem.py
index 50e338fa00..ab74566c92 100644
--- a/tests/test_pvsystem.py
+++ b/tests/test_pvsystem.py
@@ -509,7 +509,7 @@ def test_PVSystem_faiman_rad_celltemp(mocker):
     temps = 25
     irrads = 1000
     winds = 1
-    out = system.get_cell_temperature(irrads, temps, winds, ir_down,
+    out = system.get_cell_temperature(irrads, temps, winds, ir_down=ir_down,
                                       model='faiman_rad')
     temperature.faiman_rad.assert_called_once_with(irrads, temps, winds,
                                                    ir_down, u0, u1,

From 534e2f6f43cd825ffce5d03993d9af2123667acd Mon Sep 17 00:00:00 2001
From: ramaroesilva <rasilva@fc.ul.pt>
Date: Mon, 22 Dec 2025 23:00:26 +0000
Subject: [PATCH 18/18] Handles flake8.

---
 pvlib/pvsystem.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/pvlib/pvsystem.py b/pvlib/pvsystem.py
index d307c57613..bf4dcb685a 100644
--- a/pvlib/pvsystem.py
+++ b/pvlib/pvsystem.py
@@ -437,7 +437,7 @@ def get_cell_temperature(self, poa_global, temp_air, wind_speed, model,
         effective_irradiance : numeric or tuple of numeric, optional
             The irradiance that is converted to photocurrent in W/m^2.
             Only used for some models.
-   
+
         ir_down: numeric, optional
             Downwelling infrared radiation from the sky, measured on a
             horizontal surface in W/m^2. Only used in ``'faiman_rad'`` model.