Example draft

Author: echedey-ls

docs/examples/spectrum/plot_martin_ruiz_spectral_modifier_figs4to6.py

@@ -0,0 +1,52 @@
+"""
+Modeling N. Martin and J. Ruiz Spectral Modifiers
+=================================================
+
+Mimic Figures 4, 5 & 6 from paper [1]_.
+Note raw data is unavailable, so we are only plotting the line given from the
+model.
+
+References
+----------
+.. [1] Martín, N. and Ruiz, J.M. (1999), A new method for the spectral
+   characterisation of PV modules. Prog. Photovolt: Res. Appl., 7: 299-310.
+   :doi:10.1002/(SICI)1099-159X(199907/08)7:4<299::AID-PIP260>3.0.CO;2-0
+"""
+
+from pvlib.spectrum.mismatch import martin_ruiz_spectral_modifier
+from pvlib.location import Location
+from pvlib.iotools import get_pvgis_tmy
+from pvlib.irradiance import get_extra_radiation, clearness_index
+from datetime import datetime, timedelta
+import numpy as np
+import pandas as pd
+import matplotlib.pyplot as plt
+
+# clearness_index = np.linspace(0.56, 0.82, 10)
+
+site = Location(40.4534, -3.7270, altitude=664, name='IES-UPM, Madrid')
+
+time = pd.date_range(start=datetime(2020, 1, 1), end=datetime(2020, 12, 31),
+                     freq=timedelta(hours=1))
+
+tmy_data, _, _, _ = get_pvgis_tmy(site.latitude, site.longitude,
+                                  map_variables=True)
+tmy_data = tmy_data.shift(freq=timedelta(minutes=-30))
+
+solar_pos = site.get_solarposition(tmy_data.index)
+
+extra_rad = get_extra_radiation(tmy_data.index)
+
+clearness = clearness_index(ghi=tmy_data['ghi'],
+                            solar_zenith=solar_pos['zenith'],
+                            extra_radiation=extra_rad)
+
+print('clearness')
+print(clearness)
+clearness.plot()
+
+plt.show()
+
+print('pressure')
+pressure = site.get_airmass(model='kasten1966')
+print(pressure)