replace equation with new function for test purposes

Author: Uwe Krien

pvlib/modelchain.py

@@ -645,9 +645,9 @@ def complete_irradiance(self, times=None, weather=None):
         warnings.warn(wrn_txt, UserWarning)
         if {'ghi', 'dhi'} <= icolumns and 'dni' not in icolumns:
             logging.debug('Estimate dni from ghi and dhi')
-            self.weather.loc[:, 'dni'] = (
-                (self.weather.loc[:, 'ghi'] - self.weather.loc[:, 'dhi']) /
-                tools.cosd(self.solar_position.loc[:, 'zenith']))
+            self.weather.loc[:, 'dni'] = pvlib.irradiance.dni(
+                self.weather.loc[:, 'ghi'], self.weather.loc[:, 'dhi'],
+                self.solar_position.loc[:, 'zenith'])
         elif {'dni', 'dhi'} <= icolumns and 'ghi' not in icolumns:
             logging.debug('Estimate ghi from dni and dhi')
             self.weather.loc[:, 'ghi'] = (