sediment-laden river plume: sediment point source

[jamesmao]

Dear ROMS users:

I am trying to simulate sediment-laden river plume. The model grid is similar to the test case ---- riverplume in ROMS. But, I don't know how to introduce sediment from the river. Is it similar to TS_psource? Does ROMS has option of sediment_psource? or I need to add it by myself?

Has anyone experience setting this up ? Thanks a lot.

Any help welcome!


Peng Cheng

[abever]

Here is what we changed to get riverine sediment.

added the river sediment to varinfo.dat
'river_mud1'
...
...
'river_sediment, scalar, series'
'river_time'
'idRtrc(idmud)'
'nulvar'
1.0d0

in get_data.F add
# ifdef RIVER_SEDIMENT
ELSE IF (SOURCES(ng)%Fsrc(is).eq.4.0_r8) THEN
SOURCES(ng)%Lsrc(is,itemp)=.true.
SOURCES(ng)%Lsrc(is,isalt)=.true.
SOURCES(ng)%Ltracer(itemp)=.true.
SOURCES(ng)%Ltracer(isalt)=.true.
DO i=1,NCS
SOURCES(ng)%Lsrc(is,idmud(i))=.true.
SOURCES(ng)%Ltracer(idmud(i))=.true.
END DO
# endif

in mod_ncparam.F add
CASE ('idRtrc(idmud)')
DO i=1,NCS
idRtrc(idmud(i))=varid-NCS+i
END DO

Make sure to set the flag to 4, as we have specified in get_data.F to include the sediment tracer. The ifdef RIVER... is simply so you can turn it on or off.

If this doesn't seem to work let me know.
Aaron

[jamesmao]

Aaron:

Thanks a lot for your help. I am using an analytical gird which is similar to the test case --- river plume. So, I added several lines in ana_psource and it works (I am using ROMS 2.1).

In ana_psource

DO is=1,Nsrc
Lsrc(is,itemp)=.true.
Lsrc(is,isalt)=.true.
DO ised=1,NST
Lsrc(is,idsed(ised))=.true.
END DO
END DO

...

DO k=1,N(ng)
DO is=1,Nsrc
Tsrc(is,k,itemp)=20.0_r8
Tsrc(is,k,isalt)=0.0_r8
DO ised=1,NST
Tsrc(is,k,idsed(ised))=50.0_r8
END DO
END DO
END DO

I also defined sediment bed layer properties in ana_sediment.

Thanks again, Aaron. I will try your method using netcdf input files.

Peng