added main placeholder test and temporarily took out module area in short wave heat transfer calculation

Author: stephenjkaplan

pvlib/temperature.py

@@ -796,7 +796,7 @@ def hayes(poa_effective, temp_air, wind_speed, module_efficiency, module_area,
     dt = pd.to_timedelta(to_offset(freq)).seconds
 
     # radiation (from sun)
-    q_short_wave_radiation = module_area * poa_effective
+    q_short_wave_radiation = poa_effective # module_area * poa_effective
 
     # converted electrical energy
     p_out = module_efficiency * module_area * poa_effective

pvlib/tests/test_temperature.py

@@ -5,7 +5,7 @@
 from conftest import DATA_DIR, assert_series_equal
 from numpy.testing import assert_allclose
 
-from pvlib import temperature
+from pvlib import temperature, location, irradiance, iam
 
 
 @pytest.fixture
@@ -206,5 +206,56 @@ def test__calculate_radiative_heat():
 
 def test_hayes():
     # TODO placeholder until final model is validated
-    pass
 
+    # simulate psm3 data
+    data_psm3 = [
+        {'Latitude': 39.66, 'Longitude': -105.207},
+        pd.DataFrame(
+            data={
+                'DNI': [0, 163, 133, 189],
+                'DHI': [0, 4, 12, 16],
+                'GHI': [0, 7, 16, 25],
+                'Temperature': [-13.2, -13.1, -13.1, -13],
+                'Wind Speed': [1.6, 1.7, 1.7, 1.7]
+            },
+            index=pd.date_range('2019-01-01 07:25:00',
+                                '2019-01-01 07:40:00',
+                                freq='5min')
+        )
+    ]
+
+    # data preparation
+    module_tilt = 30
+    module_azimuth = 180
+    site = location.Location(
+        latitude=data_psm3[0]['Latitude'],
+        longitude=data_psm3[0]['Longitude'],
+        tz='MST'
+    )
+    solar_position = site.get_solarposition(times=data_psm3[1].index)
+    poa_global = irradiance.get_total_irradiance(
+        surface_tilt=module_tilt,
+        surface_azimuth=module_azimuth,
+        dni=data_psm3[1]['DNI'],
+        ghi=data_psm3[1]['GHI'],
+        dhi=data_psm3[1]['DHI'],
+        solar_zenith=solar_position['apparent_zenith'],
+        solar_azimuth=solar_position['azimuth']
+    )['poa_global']
+    temp_air = data_psm3[1]['Temperature']
+    wind_speed = data_psm3[1]['Wind Speed']
+
+    # 1. Calculate module temp with new model
+    aoi = irradiance.aoi(module_tilt, module_azimuth,
+                         solar_position['zenith'],
+                         solar_position['azimuth'])
+    poa_effective = poa_global.multiply(iam.ashrae(aoi))
+    module_efficiency = 0.176
+    module_area = 2.47  # m^2
+    module_weight = 34.5
+    tmod_hayes = temperature.hayes(poa_effective, temp_air, wind_speed,
+                                   module_efficiency, module_area,
+                                   module_weight, module_tilt)
+
+    assert [round(t, 2) for t in tmod_hayes.values] == \
+           [-13.20, -7.81, -8.98, -9.85]