Hi,
The way the model works is that the critical erosion stresses are constant per mud layer and specified by keywords.
Bed layer thicknesses may vary in depth and also horizontally depending on the consolidation state of the mud. The initial layer thicknesses are set in subroutine init_compo_coh.f.
In the default file, layer thickesses are uniform, but in order to specify some more rigid bed in some areas, the easiest way would be to empty the top layers in those areas. You could program this in init_compo_coh.f
In the 3D model we are programming a spatially varying bed shear strength which can be read from the geometry file. This will be for the next release 7.0 and could easily be programed in Sisyphe.
Best luck, Cath