How to control sediment resuspension in ROMS(version783)?

[zsurenj]

Dear all:
Does anyone have similar problem like this:when the sediment modular activated defined as following in app.h:
#define SEDIMENT
#define SUSPLOAD

while sediment parameters defined as following in sediment.in:

MUD_Ltsrc == T
! Surface erosion rate (kg/m2/s).

MUD_ERATE == 5.0d-2

! Critical shear for erosion and deposition (N/m2).

MUD_TAU_CE == 0.1d-5
MUD_TAU_CD == 0.1d-5
That is say that critical shear for erosion is very small, and surface erosion rate is big.
then the model didn't get any sediment resuspension,So what's the problem?
while I change MUD_TAU_CE from large to small value, the model still can't get sediment resuspended, but sediment(mud and sand) from river can appear.

[jcwarner]

I think you should add a bunch of variables to your history file and see what is happening. activate write out of bottom stresses, ssc, bed thickness, bed mass, etc etc. then you can look at computed bottom stress as compared to tau_crit erosion and bedmass etc to get a better picture of what is happening.
-j

[zsurenj]

Jcwarner,Thanks your advice.In fact,I already checked history file, you can find the bottom stress is right in the attached picture,which is greater than tau_crit erosion in the most model region, but there are not sediment resuspended, so I am confused about it. By the way, I can get right results using ROMS module from COAWST.

[jcwarner]

still could be a few things. I saw in your .h file

#define SEDIMENT
#if def SEDIMENT
# define SUSPLOAD
# define SED_MORPH
# undef SED_DENS
#endif

that second line looks funny. make it be
#ifdef SEDIMENT

also, can you post your stdout?
what it is the bed thickness?