Hello,
Two points:
(1) please first check that you problem is within the range of instabilities in which bars are formed (the beta coefficient beta=h/L, with L the channel width see papers from Zolezzi, Seminara, etc.).
(2) As you are starting a problem from scratch, you'll need to (i) have a 'clean' hydrodynamics that will serve you as hotstart for your (ii) morphodynamics/bar problem.
Please proceed step by step, moving towards the more complex aspects of your application. You can also take a look to our article appeared last year in JHE (Mendoza et al., 2016).
Hope it helps,
Pablo