function area = storage_integrate(q)

DT = diff(q.data(:,1));
DT = DT * ones(1,size(q.data,2)-1);

area = sum(0.5*(DT.*(q.data(1:(end-1),2:end)+q.data(2:end,2:end))));
