
%calcIO.m   This calcs the IO matrix
%Run this after running loadIO.m

s = size(M,1);
ratio = repmat(sum(M),s,1);
MS = M ./ ratio;

IO = MS * U;

% Now write to a wk1 file
wk1write('Aggregate IO table.wk1',IO,1,1);
   
   
   
   
   
