TELEMAC-MASCARET Forum: F2 in Subroutine SUSPENSION_FLUX

Welcome, Guest

TOPIC: F2 in Subroutine SUSPENSION_FLUX_MIXTE

F2 in Subroutine SUSPENSION_FLUX_MIXTE 11 years 2 weeks ago #10866

qilong
OFFLINE
Expert Boarder
Posts: 340
Thank you received: 33

Hello,

I'm doing the mixed sediment modelling and recently I've found a problem in the subroutine SUSPENSION_FLUX_MIXTE. I used the debugging mode and I found, in my case, the variable F2 (percentage of the mass of mud) could be NaN if some layers are completely eroded.

So I changed the expression

F2= MS_VASE(I, J)/(MS_VASE(I, J) + MS_SABLE(I, J))

to

F2= MS_VASE(I, J)/MAX(MS_VASE(I, J) + MS_SABLE(I, J),1.D-04)

Maybe it's safer to use sand fraction instead of mass percentage when the layer thickness is zero.

Kind regards,
Qilong

The administrator has disabled public write access.

The following user(s) said Thank You: Pablo

F2 in Subroutine SUSPENSION_FLUX_MIXTE 11 years 2 weeks ago #10872

cath OFFLINE Moderator Posts: 52 Thank you received: 24	Thanks ... We will correct this in the next release. Cath
	The administrator has disabled public write access. The following user(s) said Thank You: qilong