Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Wed Sep 19, 2018 1:26 am




Post new topic Reply to topic  [ 4 posts ] 

All times are UTC

Author Message
PostPosted: Fri Oct 14, 2016 1:11 pm 
Offline

Joined: Wed Apr 07, 2004 10:48 pm
Posts: 52
Location: 2nd Institute of Oceanography,SOA
Hi, ROMS users,

Has anyone succeeded in introduce the river sediments into the model? The effects of rivers' temperature and salinity can be seen clearly in my model but not the sediment.

I intended to add a mud to the present circulation model, which has already included the Yangtze River successfully (could be recognized from salinity). MUD_Ltsrc = T. and cpp file was defined like this:
#define SEDIMENT
#ifdef SEDIMENT
# define SUSPLOAD
#endif
#define ANA_SEDIMENT

The model can run for a year until I stopped it. I plotted the model output mud_01, it seems to have no signals of riverine mud_01 in the estuary. I also compared the case with MUD_Ltsrc = F, but nothing is different.

sediment_inp.h will read in a new variable related to river sediment: LtracerSrc
and I checked the LtracerSrc in the Build/get_data.f90, and found
Quote:
DO i=1,NT(ng)
IF (LtracerSrc(i,ng)) THEN
CALL get_ngfld (ng, iNLM, idRtrc(i), SSF(ng)%ncid, &
& 1, SSF(ng), update(1), &
& 1, Nsrc(ng), N(ng), 2, 1, Nsrc(ng), N(ng), &
& SOURCES(ng) % TsrcG(:,:,:,i))
IF (exit_flag.ne.NoError) RETURN
END IF
END DO


It is supposed to read the river_mud_01. So what's wrong with the case?
Any comments?

Feng


Last edited by FengZhou on Fri Oct 14, 2016 1:25 pm, edited 1 time in total.

Top
 Profile  
Reply with quote  
PostPosted: Fri Oct 14, 2016 1:22 pm 
Offline

Joined: Wed Dec 31, 2003 6:16 pm
Posts: 704
Location: USGS, USA
There could be a couple of issues, so lets try to figure it out.
- Did you set NCS (or NNS) to 1 at top of ocean.in?
- how are you adding sediment to the river> netcdf file? If you look at the stdout, does it show that it is trying to read the file and getting mud_01 values?
- you could just make the magnitude of the mud be very large (1000000) and run model for just a few time steps to see if the mud is coming in for testing.
- set Lsediment == T in the sediment.in file.
- set MUD_Ltsrc == T

-john


Top
 Profile  
Reply with quote  
PostPosted: Fri Oct 14, 2016 2:08 pm 
Offline

Joined: Wed Apr 07, 2004 10:48 pm
Posts: 52
Location: 2nd Institute of Oceanography,SOA
Hi, John,

Thank you very much! You are right. The river sediment concentration is one order lower than the local value caused by erosion and resuspension, that is why it cannot be seen.

regards,

Feng


Top
 Profile  
Reply with quote  
PostPosted: Fri Oct 14, 2016 2:15 pm 
Offline

Joined: Wed Dec 31, 2003 6:16 pm
Posts: 704
Location: USGS, USA
great. keep us posted if you any other issues.


Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 4 posts ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 1 guest


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group