function [rivers,jd]=getriv(start,stop,adjust);
%[rivers,jd]=getriv(start,stop,[adjust]);
% discharge from Merrimack, Saco, Androscoggin+Kennebec and Penobscot
% if [adjust] is present, correct for downstream area not gauged
[merr,jd]=getts('/usr/users/rps/massbay/runoff/merrimack.cdf','merrimack',start,stop);
[andro,jd]=getts('/usr/users/rps/massbay/runoff/androscoggin.cdf',...
     'androscoggin',start,stop);
[kenn,jd]=getts('/usr/users/rps/massbay/runoff/kennebec.cdf','kennebec',start,stop);
[penob,jd]=getts('/usr/users/rps/massbay/runoff/penobscot.cdf',...
      'penobscot',start,stop);
[saco,jd]=getts('/usr/users/rps/massbay/runoff/saco.cdf','saco',start,stop);
if (exist('adjust')),
  rivers=merr(:)*1.08;
  rivers(:,2)=saco(:)*1.31;
  rivers(:,3)=(andro(:)+kenn(:))*1.15;
  rivers(:,4)=penob(:)*1.10;
else 
  rivers=merr(:);
  rivers(:,2)=saco(:);
  rivers(:,3)=andro(:)+kenn(:);
  rivers(:,4)=penob(:);
end
