modified power_stats.py so that it can compute a maximum
power by hour 24 element timeseries for all sites. Then made function
that compute load diversity factor