Skip to content

move ModelChain.weather and times to ModelChainResult #1197

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 7 commits into from
Mar 18, 2021
Merged

move ModelChain.weather and times to ModelChainResult #1197

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
5cf1a7e
move ModelChain.weather and times to ModelChainResult
wholmgren Mar 13, 2021
b328c87
stickler
wholmgren Mar 13, 2021
79a6b8a
Update pvlib/tests/test_modelchain.py
wholmgren Mar 14, 2021
2712e32
update docs
wholmgren Mar 14, 2021
e33fa10
Merge branch 'wxtimestomcresults' of github.com:wholmgren/pvlib-pytho…
wholmgren Mar 14, 2021
793a1f2
Merge remote-tracking branch 'pvlib/master' into wxtimestomcresults
wholmgren Mar 16, 2021
2681e1a
Merge remote-tracking branch 'pvlib/master' into wxtimestomcresults
wholmgren Mar 16, 2021
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
7 changes: 5 additions & 2 deletions docs/sphinx/source/modelchain.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -420,7 +420,10 @@ function if you wanted to.

def pvusa_mc_wrapper(mc):
# calculate the dc power and assign it to mc.dc
mc.dc = pvusa(mc.total_irrad['poa_global'], mc.weather['wind_speed'], mc.weather['temp_air'],
# in the future, need to explicitly iterate over system.arrays
# https://github.com/pvlib/pvlib-python/issues/1115
mc.dc = pvusa(mc.results.total_irrad['poa_global'],
mc.results.weather['wind_speed'], mc.results.weather['temp_air'],
mc.system.module_parameters['a'], mc.system.module_parameters['b'],
mc.system.module_parameters['c'], mc.system.module_parameters['d'])

Expand All @@ -436,7 +439,7 @@ function if you wanted to.

def no_loss_temperature(mc):
# keep it simple
mc.cell_temperature = mc.weather['temp_air']
mc.cell_temperature = mc.results.weather['temp_air']
return mc


Expand Down
2 changes: 2 additions & 0 deletions docs/sphinx/source/whatsnew/v0.9.0.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -64,6 +64,8 @@ Deprecations
* ``ModelChain.spectral_modifier``
* ``ModelChain.total_irrad``
* ``ModelChain.tracking``
* ``ModelChain.weather``
* ``ModelChain.times``

Enhancements
~~~~~~~~~~~~
Expand Down
121 changes: 68 additions & 53 deletions pvlib/modelchain.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -256,7 +256,7 @@ class ModelChainResult:
_per_array_fields = {'total_irrad', 'aoi', 'aoi_modifier',
'spectral_modifier', 'cell_temperature',
'effective_irradiance', 'dc', 'diode_params',
'dc_ohmic_losses'}
'dc_ohmic_losses', 'weather'}

# system-level information
solar_position: Optional[pd.DataFrame] = field(default=None)
Expand All @@ -278,6 +278,9 @@ class ModelChainResult:
diode_params: Optional[PerArray[pd.DataFrame]] = field(default=None)
dc_ohmic_losses: Optional[PerArray[pd.Series]] = field(default=None)

weather: Optional[PerArray[pd.DataFrame]] = None
times: Optional[pd.DatetimeIndex] = None

def _result_type(self, value):
"""Coerce `value` to the correct type according to
``self._singleton_tuples``."""
Expand Down Expand Up @@ -375,7 +378,8 @@ class ModelChain:
_deprecated_attrs = ['solar_position', 'airmass', 'total_irrad',
'aoi', 'aoi_modifier', 'spectral_modifier',
'cell_temperature', 'effective_irradiance',
'dc', 'ac', 'diode_params', 'tracking']
'dc', 'ac', 'diode_params', 'tracking',
'weather', 'times']

def __init__(self, system, location,
clearsky_model='ineichen',
Expand Down Expand Up @@ -406,9 +410,6 @@ def __init__(self, system, location,
self.dc_ohmic_model = dc_ohmic_model
self.losses_model = losses_model

self.weather = None
self.times = None

self.results = ModelChainResult()

def __getattr__(self, key):
Expand Down Expand Up @@ -908,7 +909,7 @@ def infer_spectral_model(self):

def first_solar_spectral_loss(self):
self.results.spectral_modifier = self.system.first_solar_spectral_loss(
_tuple_from_dfs(self.weather, 'precipitable_water'),
_tuple_from_dfs(self.results.weather, 'precipitable_water'),
self.results.airmass['airmass_absolute']
)
return self
Expand Down Expand Up @@ -1005,8 +1006,8 @@ def _set_celltemp(self, model):

poa = _irrad_for_celltemp(self.results.total_irrad,
self.results.effective_irradiance)
temp_air = _tuple_from_dfs(self.weather, 'temp_air')
wind_speed = _tuple_from_dfs(self.weather, 'wind_speed')
temp_air = _tuple_from_dfs(self.results.weather, 'temp_air')
wind_speed = _tuple_from_dfs(self.results.weather, 'wind_speed')
kwargs = {}
if model == self.system.noct_sam_celltemp:
kwargs['effective_irradiance'] = self.results.effective_irradiance
Expand Down Expand Up @@ -1161,7 +1162,7 @@ def complete_irradiance(self, weather):

Notes
-----
Assigns attributes: ``weather``
Assigns attributes to ``results``: ``times``, ``weather``

Examples
--------
Expand All @@ -1174,7 +1175,7 @@ def complete_irradiance(self, weather):
>>> # my_weather containing 'dhi' and 'ghi'.
>>> mc = ModelChain(my_system, my_location) # doctest: +SKIP
>>> mc.complete_irradiance(my_weather) # doctest: +SKIP
>>> mc.run_model(mc.weather) # doctest: +SKIP
>>> mc.run_model(mc.results.weather) # doctest: +SKIP

>>> # my_weather containing 'dhi', 'ghi' and 'dni'.
>>> mc = ModelChain(my_system, my_location) # doctest: +SKIP
Expand All @@ -1184,16 +1185,16 @@ def complete_irradiance(self, weather):
self._check_multiple_input(weather)
# Don't use ModelChain._assign_weather() here because it adds
# temperature and wind-speed columns which we do not need here.
self.weather = _copy(weather)
self.results.weather = _copy(weather)
self._assign_times()
self.results.solar_position = self.location.get_solarposition(
self.times, method=self.solar_position_method)
self.results.times, method=self.solar_position_method)

if isinstance(weather, tuple):
for w in self.weather:
for w in self.results.weather:
self._complete_irradiance(w)
else:
self._complete_irradiance(self.weather)
self._complete_irradiance(self.results.weather)

return self

Expand Down Expand Up @@ -1238,7 +1239,7 @@ def _prep_inputs_solar_pos(self, weather):
pass

self.results.solar_position = self.location.get_solarposition(
self.times, method=self.solar_position_method,
self.results.times, method=self.solar_position_method,
**kwargs)
return self

Expand Down Expand Up @@ -1301,15 +1302,23 @@ def _verify(data, index=None):
for (i, array_data) in enumerate(data):
_verify(array_data, i)

def _configure_results(self):
"""Configure the type used for per-array fields in ModelChainResult.
def _configure_results(self, per_array_data):
"""Configure the type used for per-array fields in
ModelChainResult.

Must be called after ``self.weather`` has been assigned. If
``self.weather`` is a tuple and the number of arrays in the system
is 1, then per-array results are stored as length-1 tuples.
If ``per_array_data`` is True and the number of arrays in the
system is 1, then per-array results are stored as length-1
tuples. This overrides the PVSystem defaults of unpacking a 1
length tuple into a singleton.

Parameters
----------
per_array_data : bool
If input data is provided for each array, pass True. If a
single input data is provided for all arrays, pass False.
"""
self.results._singleton_tuples = (
self.system.num_arrays == 1 and isinstance(self.weather, tuple)
self.system.num_arrays == 1 and per_array_data
)

def _assign_weather(self, data):
Expand All @@ -1321,13 +1330,13 @@ def _build_weather(data):
if weather.get('temp_air') is None:
weather['temp_air'] = 20
return weather
if not isinstance(data, tuple):
self.weather = _build_weather(data)
if isinstance(data, tuple):
weather = tuple(_build_weather(wx) for wx in data)
self._configure_results(per_array_data=True)
else:
self.weather = tuple(
_build_weather(weather) for weather in data
)
self._configure_results()
weather = _build_weather(data)
self._configure_results(per_array_data=False)
self.results.weather = weather
self._assign_times()
return self

Expand All @@ -1344,18 +1353,20 @@ def _build_irrad(data):
return self

def _assign_times(self):
"""Assign self.times according the the index of self.weather.

If there are multiple DataFrames in self.weather then the index
of the first one is assigned. It is assumed that the indices of
each data frame in `weather` are the same. This can be verified
by calling :py:func:`_all_same_index` or
:py:meth:`self._check_multiple_weather` before calling this method.
"""Assign self.results.times according the the index of
self.results.weather.

If there are multiple DataFrames in self.results.weather then
the index of the first one is assigned. It is assumed that the
indices of each DataFrame in self.results.weather are the same.
This can be verified by calling :py:func:`_all_same_index` or
:py:meth:`self._check_multiple_weather` before calling this
method.
"""
if isinstance(self.weather, tuple):
self.times = self.weather[0].index
if isinstance(self.results.weather, tuple):
self.results.times = self.results.weather[0].index
else:
self.times = self.weather.index
self.results.times = self.results.weather.index

def prepare_inputs(self, weather):
"""
Expand Down Expand Up @@ -1386,8 +1397,8 @@ def prepare_inputs(self, weather):

Notes
-----
Assigns attributes: ``weather``, ``solar_position``, ``airmass``,
``total_irrad``, ``aoi``
Assigns attributes to ``results``: ``times``, ``weather``,
``solar_position``, ``airmass``, ``total_irrad``, ``aoi``

See also
--------
Expand Down Expand Up @@ -1421,9 +1432,9 @@ def prepare_inputs(self, weather):
self.results.solar_position['azimuth'])

self.results.total_irrad = get_irradiance(
_tuple_from_dfs(self.weather, 'dni'),
_tuple_from_dfs(self.weather, 'ghi'),
_tuple_from_dfs(self.weather, 'dhi'),
_tuple_from_dfs(self.results.weather, 'dni'),
_tuple_from_dfs(self.results.weather, 'ghi'),
_tuple_from_dfs(self.results.weather, 'dhi'),
airmass=self.results.airmass['airmass_relative'],
model=self.transposition_model
)
Expand Down Expand Up @@ -1482,8 +1493,8 @@ def prepare_inputs_from_poa(self, data):

Notes
-----
Assigns attributes: ``weather``, ``total_irrad``, ``solar_position``,
``airmass``, ``aoi``.
Assigns attributes to ``results``: ``times``, ``weather``,
``total_irrad``, ``solar_position``, ``airmass``, ``aoi``.

See also
--------
Expand Down Expand Up @@ -1642,10 +1653,12 @@ def run_model(self, weather):

Notes
-----
Assigns attributes: ``solar_position``, ``airmass``, ``weather``,
``total_irrad``, ``aoi``, ``aoi_modifier``, ``spectral_modifier``,
and ``effective_irradiance``, ``cell_temperature``, ``dc``, ``ac``,
``losses``, ``diode_params`` (if dc_model is a single diode model).
Assigns attributes to ``results``: ``times``, ``weather``,
``solar_position``, ``airmass``, ``total_irrad``, ``aoi``,
``aoi_modifier``, ``spectral_modifier``, and
``effective_irradiance``, ``cell_temperature``, ``dc``, ``ac``,
``losses``, ``diode_params`` (if dc_model is a single diode
model).

See also
--------
Expand Down Expand Up @@ -1701,10 +1714,12 @@ def run_model_from_poa(self, data):

Notes
-----
Assigns attributes: ``solar_position``, ``airmass``, ``weather``,
``total_irrad``, ``aoi``, ``aoi_modifier``, ``spectral_modifier``,
and ``effective_irradiance``, ``cell_temperature``, ``dc``, ``ac``,
``losses``, ``diode_params`` (if dc_model is a single diode model).
Assigns attributes to results: ``times``, ``weather``,
``solar_position``, ``airmass``, ``total_irrad``, ``aoi``,
``aoi_modifier``, ``spectral_modifier``, and
``effective_irradiance``, ``cell_temperature``, ``dc``, ``ac``,
``losses``, ``diode_params`` (if dc_model is a single diode
model).

See also
--------
Expand Down Expand Up @@ -1798,7 +1813,7 @@ def run_model_from_effective_irradiance(self, data=None):
If optional column ``'poa_global'`` is present, these data are used.
If ``'poa_global'`` is not present, ``'effective_irradiance'`` is used.

Assigns attributes: ``weather``, ``total_irrad``,
Assigns attributes to results: ``times``, ``weather``, ``total_irrad``,
``effective_irradiance``, ``cell_temperature``, ``dc``, ``ac``,
``losses``, ``diode_params`` (if dc_model is a single diode model).

Expand Down
Loading