subroutine river_read_input(subbasin_input_file_id)
use input, only : read_real_column
integer, intent(in) :: subbasin_input_file_id
call read_real_column(subbasin_input_file_id, "chl", chl(1, :), range=(/0., 1e4/), closed="n")
chl(2, :) = chl(1, :)
call read_real_column(subbasin_input_file_id, "chs", chs, range=(/0., 3./), closed="u")
chss = chs
call read_real_column(subbasin_input_file_id, "chw", chw(1, :), range=(/0., 1e4/), closed="n")
chw(2, :) = chw(1, :)
chw(2, :) = chw(2, :) * chwc0
call read_real_column(subbasin_input_file_id, "chd", chd)
call read_real_column(subbasin_input_file_id, "chk", chk(1, :), 0.37_dp)
call read_real_column(subbasin_input_file_id, "chk", chk(2, :), 0.37_dp)
call read_real_column(subbasin_input_file_id, "chn", chn, 0.075_dp)
call read_real_column(subbasin_input_file_id, "ovn", ovn, 0.15_dp)
call read_real_column(subbasin_input_file_id, "rt", rt, 0.0_dp)
call read_real_column(subbasin_input_file_id, "sdtsav", sdtsav, 0.0_dp)
end subroutine river_read_input