Hello,
Recently I looked at the code for computing erosion flux for mixed sediment. But I'm a little confused about the following piece of code in the subroutine "suspension_computation.f" from line 578 to line 591:
IF(.NOT.SEDCO) THEN
IF(DEBUG > 0) WRITE(LU,*) 'SUSPENSION_FLUX_MIXTE'
CALL SUSPENSION_FLUX_MIXTE(T4,HN,FDM,NPOIN,CHARR,XMVE,XMVS,
& VCE,GRAV,HMIN,XWC,ZERO,
& PARTHENIADES,FLUER,FLUER_VASE,
& ZREF,AC,CSTAEQ,QS_C,ICQ,DEBUG,
& AVAIL,NSICLA,ES,TOCE_VASE,
& TOCE_SABLE,NOMBLAY,
& DT,TOCE_MIXTE%R,MS_SABLE,
& MS_VASE)
IF (DEBUG > 0) WRITE(LU,*) 'END_SUSPENSION_FLUX_MOY'
ENDIF
IF(SEDCO) CALL OS('X=Y ',X=FLUER, Y=FLUER_VASE)
I think in the subroutine "suspension_flux_mixte.f", FLUER_VASE and FLUER_SABLE are calculated. But in the above code, only for cohesive sediment it sets FLUER=FLUER_VASE. What about non-cohesive sediment? Is there a line missing to set FLUER=FLUER_SABLE if (.NOT. SEDCO)?
Kind regards,
Qilong