about the erosion flux

Sediment modeling collaborators: issues, applications, information exchange

Moderators: arango, robertson, rsignell

Post Reply
Message
Author
duckweed

about the erosion flux

#1 Unread post by duckweed »

Hi,
i found in my his.nc,the erosion flux is eual to 0 all the time.
However, i changed in the sed_fluxes.F. I write the erosion flux to log file.

Code: Select all

!  Compute erosion, ero_flux (kg/m2).
!
            cff1=(1.0_r8-bed(i,j,1,iporo))*bed_frac(i,j,1,ised)
            cff2=dt(ng)*Erate(ised,ng)*cff1
            cff3=Srho(ised,ng)*cff1
            cff4=bed_mass(i,j,1,bnew,ised)
            ero_flux(i,j,ised)=                                         &
     &               MIN(MAX(0.0_r8,cff2*(cff*tau_w(i,j)-1.0_r8)),      &
     &                   MIN(cff3*bottom(i,j,iactv),cff4)+              &
     &                   settling_flux(i,j,ised))
!my adjust			
			write(*,*) ero_flux(60,40,1) 
In the logfile the ero_flux is not 0.It's a small value such as 0.001.However,in his.nc,the ero_flux is 0.
Why?

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

Re: about the erosion flux

#2 Unread post by jcwarner »

really not sure. is the entire field 0? does any sediment get resupseded? are the mud_01 or sand01 greater than zero?

duckweed

Re: about the erosion flux

#3 Unread post by duckweed »

Thank you !
I checked the his.nc. The entire erosion is 0. And the dep_net is 0.
The sand_01 is greater than 0, is about 0.008.
I have a little confused about resuspension. which case can cause resuspension? 
I know, if ero_flux<settling_flux,deposit.

Code: Select all

            IF ((ero_flux(i,j,ised)-settling_flux(i,j,ised)).lt.        &
     &           0.0_r8) THEN
So, i think, if ero_flux>settling_flux,resuspension.However, i can't find it in the code.

Post Reply