In order of invocation in Hydro code:

sed_mod - parameters for sediment model

read_sed_input - subroutine to read the sediment model inputs written in sediment.in file; optionally
                 calculate settling vel and critical shear stress;

sed_init - subroutine to set initial conditions for sediment tracer fields and initial bed conditions

sed_roughness (inside sed_friction.F90) - calculate bottom roughness

sediment - routine to compute the sediment source and sinks. It includes the vertical settling of 
           the sediment, erosional and depositional flux, transport of multiple grain sizes, bedload 
           based on Meyer Peter and Mueller (not active yet) or van Rijn. Morphological model is found inside.
           Vertical settling is handled with semi-Lagrangian PPM and eventually cast into a body force for
           the transport solver.

[other routines:
sed_bedload - bedload routines (including call to solve_jcg() for morphol. model)
sed_friction - bottom formulations
sed_filter - morphological filters
sed_misc_subs - misc routines
]

----------------------------------------------------------------------------------------------------
Joseph's notes on doing SED simulations:
(1) itr_met=3 has issue when the susp. load becomes very large; always use upwind (itr_met=1 and tvd.prop=0) first
   and fine tune .prop later;
(2) using a very small dzb_min (i.e. unlimit the Cd) will lead to reduced near -bottom vel and stress -> little erosion!
(3) From erosion & deposition tests, the suspended load is sensitive to vertical diffusion. Try itur=3: k-kl and diffmin=1.e-4.
