about the erosion flux

[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]

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

[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.