Novel numerical approach for 1D variable density shallow flows over uneven rigid and erodible beds