question about transport equation of sediment

[zhaoqian]

Hi all,
When dealing with transport equation of sediment, (1) Should I define SUSPLOAD and TS_MPDATA as well as SEDIMENT? (2)What's the function of subroutine sediment.F and step3d_t.F?
Thanks in advance !

[jcwarner]

these are different but related questions.

TS_MPDATA is a choice for an advection scheme for the tracers. I typically use this advection scheme when using sediment since mpdata is positive definite.

for the sedimnet, the main required cpp is
#define SEDIMENT
with that defined, you can add other choices like
#defined SUSP_LOAD
and / or
#define BEDLOAD_MPM or BEDLOAD_SOULSBY.
the choices are listed in ROMS/Include/cppdefs.h
-j

[zhaoqian]

Thank you for solving my first question.
I am confused with the lower boundary of sediment transport. Where is it defined , in set_tide.F ( which calls ana_btflux or the forcing given ) or somewhere in sediment.F ?

[jcwarner]

i am not quite clear about what u mean 'lower boundary of sediment'
may i suggest that you look in sediment.F and see how that is laid out. it calls different routines depending on what was #defined. There is also some info on the wiki and we have a basic paper in Comp and Geosci (Warner et al 2008) that explains some of these issues.
-j

[zhaoqian]

i am not quite clear about what u mean 'lower boundary of sediment'

Sorry. I means lower boundary condition of sediment. I don't know where the code is, in set_tide.F ( which calls ana_btflux or the forcing given ) or somewhere in sediment.F ?.