Qdata = csvread('./missouri_data.csv',1,0);
discharge_data = 0.0283168466 .* Qdata(:,2); % convert cfs column to m3/s
Q_ws(1:run_time/timestep) = 0;
for i = 1:run_time/timestep;
    pull = round(rand() * numel(discharge_data)); % date to pull
    Q_ws(i) = discharge_data(pull);
end