normvi
C NORMVI SOURCE CHAT 05/01/13 02:00:47 5004 *************************************************************************** *** SP 'NORMVI' : calcule de maniere explicite la vitesse, sa norme et *** 'UXY' (cf. ci-desssous) au pt de depart. Module valable dans element *** reference ou element reel. *** *** APPELES 1 = aucun *** APPELES 2 = aucun *** *** E = 'NDIM' dimension de l'espace *** 'UE' vitesse dans element considere *** *** S = XNORM' norme de la vitesse *** 'UXY' -> racine carre de UX*UX+UY*UY (CAS 3D) *** *** Auteur Cyril Nou **************************************************************************** IMPLICIT INTEGER(I-N) IMPLICIT REAL*8 (A-H,O-Z) DIMENSION UE(3) UXY=0.D0 XNORM=0.D0 *** cas 2D DO 10 I=1,NDIM XNORM=XNORM+UE(I)**2 10 CONTINUE XNORM=SQRT(XNORM) *** cas 3D IF (NDIM.EQ.3) THEN DO 20 I=1,NDIM-1 UXY=UXY+UE(I)**2 20 CONTINUE UXY=SQRT(UXY) ENDIF RETURN END
© Cast3M 2003 - Tous droits réservés.
Mentions légales