Discussion :

[saida29]

Bonjour
Dans mon code lors de calcule de la vitesse de fluid j'utilise l'expression suivante:
"fluid_u(i,j)=sqrt(fluid_ux(i,j)**2.d0+fluid_uy(i,j)**2.d0)"
mon probleme est lorsque j'affiche fluid_u(i,j) je trouve 0 alors que l'affichage de sqrt(fluid_ux(i,j)**2.d0+fluid_uy(i,j)**2.d0) donne des valeurs non nuls
(print*, fluid_u(i,j),sqrt(fluid_ux(i,j)**2.d0+fluid_uy(i,j)**2.d0))
comment ça peut arriver et à quoi est dit?
aidez moi svp
Merci bien

[François L.]

Difficile de répondre à ça sans avoir un exemple minimal de code, mais je suggère de vérifier que fluid_u est bien déclaré comme un nombre à virgule flottante (REAL).

Par ailleurs, je déconseille l'usage du **2.d0, à remplacer par **2. Je ne sais pas si c'est le cas sur tous les compilateurs, mais x**2.d0 plante si x est négatif (ce qui est normal s'il passe par de l'exponentielle de log). Mais bon, ce n'était pas la question