How to control sediment resuspension in ROMS(version783)?

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
zsurenj
Posts: 9
Joined: Fri Apr 24, 2015 4:06 pm
Location: Zhongshan university

How to control sediment resuspension in ROMS(version783)?

#1 Unread post by 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.
Attachments
yees.h
(1.93 KiB) Downloaded 252 times
sediment_yees.in
(40.36 KiB) Downloaded 223 times

jcwarner
Posts: 1172
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: How to control sediment resuspension in ROMS(version783)

#2 Unread post by 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
Posts: 9
Joined: Fri Apr 24, 2015 4:06 pm
Location: Zhongshan university

Re: How to control sediment resuspension in ROMS(version783)

#3 Unread post by 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.
Attachments
2Fig_Sca_bstr.pdf
(838.4 KiB) Downloaded 235 times

jcwarner
Posts: 1172
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: How to control sediment resuspension in ROMS(version783)

#4 Unread post by 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?

Post Reply