Below is an idea, where the pvsystem.single_diode function is essentially the same as the current pvsystem.singlediode. Note the overloading of F and H vs. effective_irradiance and cell_temp, respectively, in the highest level method calc_model_dc of PVSystem, because if one is using a matched reference cell with a back-of-cell temperature sensor, then F and H are directly available, otherwise F and H would typically be first derived from SAPM's effective irradiance and cell_temp, like the other single diode models. This removes a lot of wrappers in PVSystem, but I haven't totally thought through the SAPM and PVWatts integration at this level or "all the things" into the higher level ModelChain.

class PVSystem(object):
...
    # This property would normally be set by PVSystem constructor
    self.model_dc = 'pvsyst_v6'  # or 'desoto' or 'campanelli' or 'sapm' or 'pvwatts'

    def calc_model_dc(self, overloaded_irradiance, overloaded_temperature):
        if self.model_dc ==  'pvsyst_v6':
            return sdm(get_params_pvsyst_v6(overloaded_irradiance, overloaded_temperature, **self.module_parameters))
        elif self.model_dc ==  'desoto':
            return sdm(get_params_desoto(overloaded_irradiance, overloaded_temperature, **self.module_parameters)
        elif self.model_dc ==  'campanelli':
            return sdm(get_params_campanelli(overloaded_irradiance, overloaded_temperature, **self.module_parameters)
        elif self.model_dc ==  'sapm':
            return sapm(overloaded_irradiance, overloaded_temperature, **self.module_parameters)
        elif self.model_dc ==  'pvwatts':
            return pvwatts_dc(overloaded_irradiance, overloaded_temperature, **self.module_parameters)
        else:
            # Assume custom case is given a function to call
            return self.model_dc(irr, temp, **self.module_parameters)

# These are module-level functions in pvsystem, without wrappers.
def sdm(...):
    # Essentially the same as current singlediode function (not sure about system scaling here)
    ...
    return ...

def get_params_pvsyst_v6(effective_irradiance, cell_temp, **kwargs):
    ...
    return IL, I0, Rs, Rsh, NsVth

def get_params_desoto(effective_irradiance, cell_temp, **kwargs):
    ...
    return IL, I0, Rs, Rsh, NsVth

def get_params_campanelli(F, H, **kwargs):
    ...
    return IL, I0, Rs, Rsh, NsVth

def sapm(self, effective_irradiance, temp_cell, **kwargs)
    ...
    return ...

def pvwatts_dc(g_poa_effective, temp_cell, **kwargs)
    ...
    return ...

Add PVsyst as option to ModelChain #487

Description

Metadata

Metadata

Assignees

Labels

Type

Fields

Projects

Milestone

Relationships

Development

Issue actions