add series and parallel multipliers, aoi modifier

Author: wholmgren

pvlib/modelchain.py

@@ -386,6 +386,12 @@ def run_model(self, times, irradiance=None, weather=None):
                                    self.airmass['airmass_absolute'],
                                    self.aoi)
 
+        voltages = ['v_mp', 'v_oc']
+        self.dc[voltages] *= self.system.series_modules
+        currents = ['i_mp', 'i_sc', 'i_x', 'i_xx']
+        self.dc[currents] *= self.system.parallel_modules
+        self.dc['p_mp'] = self.dc['v_mp'] * self.dc['i_mp']
+
         self.ac = self.system.snlinverter(self.dc['v_mp'], self.dc['p_mp'])
 
         return self
@@ -425,13 +431,19 @@ def run_model(self, times, irradiance=None, weather=None):
 
         self.prepare_inputs(times, irradiance, weather)
 
-        self.temps = self.system.sapm_celltemp(self.total_irrad['poa_global'],
-                                               self.weather['wind_speed'],
-                                               self.weather['temp_air'])
+        self.aoi_mod = self.system.ashraeiam(self.aoi).fillna(0)
+        self.total_irrad['poa_global_aoi'] = (
+            self.total_irrad['poa_direct'] * self.aoi_mod +
+            self.total_irrad['poa_diffuse'])
+
+        self.temps = self.system.sapm_celltemp(
+            self.total_irrad['poa_global_aoi'],
+            self.weather['wind_speed'],
+            self.weather['temp_air'])
 
         (photocurrent, saturation_current, resistance_series,
          resistance_shunt, nNsVth) = (
-            self.system.calcparams_desoto(self.total_irrad['poa_global'],
+            self.system.calcparams_desoto(self.total_irrad['poa_global_aoi'],
                                           self.temps['temp_cell']))
 
         self.desoto = (photocurrent, saturation_current, resistance_series,
@@ -441,6 +453,14 @@ def run_model(self, times, irradiance=None, weather=None):
             photocurrent, saturation_current, resistance_series,
             resistance_shunt, nNsVth)
 
+        self.dc = self.dc.fillna(0)
+
+        voltages = ['v_mp', 'v_oc']
+        self.dc[voltages] *= self.system.series_modules
+        currents = ['i_mp', 'i_sc', 'i_x', 'i_xx']
+        self.dc[currents] *= self.system.parallel_modules
+        self.dc['p_mp'] = self.dc['v_mp'] * self.dc['i_mp']
+
         self.ac = self.system.snlinverter(self.dc['v_mp'], self.dc['p_mp'])
 
         return self