.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "auto_examples/plot_solar_radiation.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code .. rst-class:: sphx-glr-example-title .. _sphx_glr_auto_examples_plot_solar_radiation.py: Exploration of heating from solar radiation ------------------------------------------- In this example, we see how the solar radiation varies with the solar altitude angle, :math:`H_s`, the difference between the solar azimuth and the span azimuth (or bearing), :math:`\left|\gamma_c - \gamma_s\right|`, and the albedo, :math:`F`. .. GENERATED FROM PYTHON SOURCE LINES 11-13 Imports and utilities ^^^^^^^^^^^^^^^^^^^^^ .. GENERATED FROM PYTHON SOURCE LINES 13-20 .. code-block:: Python import matplotlib.cm as cm import matplotlib.pyplot as plt import numpy as np from linerate.equations import cigre601 .. GENERATED FROM PYTHON SOURCE LINES 21-23 Simulation parameters ^^^^^^^^^^^^^^^^^^^^^ .. GENERATED FROM PYTHON SOURCE LINES 23-29 .. code-block:: Python solar_altitude = np.linspace(0, np.pi / 2)[:, np.newaxis] azimuth_difference = np.linspace(0, np.pi / 2, 10)[np.newaxis, :] albedo = 0.2 height_above_sea_level = 0 .. GENERATED FROM PYTHON SOURCE LINES 30-32 Intermediate calculations ^^^^^^^^^^^^^^^^^^^^^^^^^ .. GENERATED FROM PYTHON SOURCE LINES 32-37 .. code-block:: Python sin_H_s = np.sin(solar_altitude) cos_eta = np.cos(solar_altitude) * np.cos(azimuth_difference) sin_eta = np.sqrt(1 - cos_eta**2) .. GENERATED FROM PYTHON SOURCE LINES 38-40 Compute the solar radiation components ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ .. GENERATED FROM PYTHON SOURCE LINES 40-52 .. code-block:: Python I_B = cigre601.solar_heating.compute_direct_solar_radiation( sin_H_s, clearness_ratio=1, height_above_sea_level=height_above_sea_level ) I_d = cigre601.solar_heating.compute_diffuse_sky_radiation(I_B, sin_H_s) I_T_F0 = cigre601.solar_heating.compute_global_radiation_intensity( I_B, I_d, albedo=0.0, sin_angle_of_sun_on_line=sin_eta, sin_solar_altitude=sin_H_s ) I_T_F = cigre601.solar_heating.compute_global_radiation_intensity( I_B, I_d, albedo=albedo, sin_angle_of_sun_on_line=sin_eta, sin_solar_altitude=sin_H_s ) .. GENERATED FROM PYTHON SOURCE LINES 53-57 Create visualisation ^^^^^^^^^^^^^^^^^^^^ This section can be skipped, since most of the code is just there to create pretty plots. .. GENERATED FROM PYTHON SOURCE LINES 57-112 .. code-block:: Python # Setup figure and axes fig = plt.figure(figsize=(11, 1.9)) axes = [fig.add_axes([0.07, 0.27, 0.19, 0.7])] axes += [ fig.add_axes([0.30, 0.27, 0.18, 0.7], sharey=axes[0]), fig.add_axes([0.53, 0.27, 0.18, 0.7], sharey=axes[0]), fig.add_axes([0.76, 0.27, 0.18, 0.7], sharey=axes[0]), ] cbar_ax = fig.add_axes([0.955, 0.27, 0.015, 0.7]) # Add plots axes[0].plot(np.degrees(solar_altitude), I_B, color="k") axes[1].plot(np.degrees(solar_altitude), I_d, color="k") for i, d_gamma in enumerate(azimuth_difference.ravel()): color = cm.cividis(d_gamma / azimuth_difference.max()) d_gamma = np.degrees(d_gamma) axes[2].plot(np.degrees(solar_altitude), I_T_F0[:, i], color=color) axes[3].plot(np.degrees(solar_altitude), I_T_F[:, i], color=color) for ax in axes: # Setup y-axes ax.set_xlabel(r"Solar altitude $[^\circ]$") ax.set_xlim(0, 90) ax.set_xticks([0, 30, 60, 90]) # Setup y-axes to be shared ax.set_ylim(0, I_T_F.max() * 1.05) ax.set_yticks([0, 500, 1000, 1360]) # Include tick for solar constant ax.set_yticklabels([0, 500, 1000, "$G_{SC}$"]) ax.axhline(1360, color="k", linestyle="--") # Add dashed line for solar constant # Removed shared ticks axes[1].tick_params(labelleft=False) axes[2].tick_params(labelleft=False) axes[3].tick_params(labelleft=False) # Setup labels axes[0].set_ylabel(r"$I_B~[\mathrm{W}~\mathrm{m}^{-1}]$", labelpad=-1) axes[1].set_ylabel(r"$I_d~[\mathrm{W}~\mathrm{m}^{-1}]$") axes[2].set_ylabel("$F=0$\n$I_T~[\mathrm{W}~\mathrm{m}^{-1}]$") # noqa axes[3].set_ylabel(f"$F={albedo}$\n$I_T~[\mathrm{{W}}~\mathrm{{m}}^{{-1}}]$") # noqa # Colorbar cbar_ax.imshow(azimuth_difference.T, aspect="auto", cmap="cividis") cbar_ax.yaxis.set_label_position("right") cbar_ax.set_ylabel(r"$\left|\gamma_c - \gamma_s\right|~[^\circ]$", labelpad=-10) cbar_ax.yaxis.tick_right() cbar_ax.set_yticks(cbar_ax.get_ylim()) cbar_ax.set_yticklabels([0, 90]) cbar_ax.set_xticks([]) plt.show() .. image-sg:: /auto_examples/images/sphx_glr_plot_solar_radiation_001.png :alt: plot solar radiation :srcset: /auto_examples/images/sphx_glr_plot_solar_radiation_001.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none /home/runner/work/linerate/linerate/examples/plot_solar_radiation.py:43: SyntaxWarning: invalid escape sequence '\m' ) /home/runner/work/linerate/linerate/examples/plot_solar_radiation.py:44: SyntaxWarning: invalid escape sequence '\m' I_d = cigre601.solar_heating.compute_diffuse_sky_radiation(I_B, sin_H_s) /home/runner/work/linerate/linerate/examples/plot_solar_radiation.py:44: SyntaxWarning: invalid escape sequence '\m' I_d = cigre601.solar_heating.compute_diffuse_sky_radiation(I_B, sin_H_s) .. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 0.163 seconds) .. _sphx_glr_download_auto_examples_plot_solar_radiation.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: plot_solar_radiation.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: plot_solar_radiation.py ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_