function q = ri1(data)
% 2nd Riemann invariant for acoustics

q = data(:,1) + 2*data(:,2);
end
