This program calculates the average flux (3D distribution) for each energy bin across all times and viewing directions. The resulting energy spectrum is then plotted on a double-logarithmic graph.