function avg = storage_average(q)

area = storage_integrate(q);
avg = area / (q.data(end,1)-q.data(1,1));
