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 SUSPLOAD
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
It is supposed to read the river_mud_01. So what's wrong with the case?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