open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi all,
I'm still dealing with suspension in 3D. I'm working with the V6P2 (it will change soon).

I can input sediment but as soon as i specify the settling velocity "VITESSE DE CHUTE CONSTANT" with " INFLUENCE DE LA TURBULENCE SUR LA VITESSE DE CHUTE = NON " i've got NaN.In the telemacv6p2.dico the default values for the settling velocity is set to 0.

The settling velocity is computed in VITCHU.f in V6P2. (ALL OS( 'X=C ' , X=W_SED , C=WCHU_CONST ))
WC_TURB.f computed the values for 'INFLUENCE DE LA TURBULENCE SUR LA VITESSE DE CHUTE = OUI', not my case.

My .cas files :
/

---

/ SEDIMENTS
/

---

NOMBRE DE TRACEURS : 1
NOM DU TRACEUR= 'SEDIMENT G/L'
SEDIMENT : OUI
MASSE VOLUMIQUE DU SEDIMENT : 1800
VALEURS INITIALES DES TRACEURS : 0.0
DIAMETRE MOYEN DES GRAINS : 0.0001 /100um
INFLUENCE DE LA TURBULENCE SUR LA VITESSE DE CHUTE = NON
VITESSE DE CHUTE CONSTANT : 0.0077 /VPD 0
SOLVEUR POUR LA DIFFUSION DU SEDIMENT : 7
OPTION DU SOLVEUR POUR LA DIFFUSION DU SEDIMENT : 5
SCHEMA POUR LA CONVECTION DES TRACEURS : 1
PROFILS DES TRACEURS SUR LA VERTICALE : 1;1;1;1;1;1
/1=profil constant; 2=profil de Rouse
VALEURS IMPOSEES DES TRACEURS : 999;999;999;999;2;9999 /gl

Anyone know what happens when I fixed this settling velocity (with a Soulsby law) ?
(fyi the maximun value at the bottom is around 0.1 m/s, so thanks to the Rouse number particules with 100µm are transported by suspension only)

Hi,
It's ok for now.
No formulae are defined si V6P2 and the settling velocity has to be negative.
Sorry for this quick question.
Best regards
P.N